Genome-wide DNA methylation profiles changes associated with constant heat stress in pigs as measured by bisulfite sequencing

Functional transcriptome analysis revealed upregulation of MAPK-SMAD signalling pathways in chronic heat stress in crossbred cattle

Animal geneticists and breeders have the impending challenge of enhancing the resilience of Indian livestock to heat stress through better selection strategies. Climate change's impact on livestock is more intense in tropical countries like India where dairy cattle crossbreeds are more sensitive to heat stress. The main reason for this study was to find the missing relative changes in transcript levels in thermo-neutral and heat stress conditions in crossbred cattle through whole-transcriptome analysis of RNA-Seq data. Differentially expressed genes (DEGs) identified based on the minimum log twofold change value and false discovery rate 0.05 revealed 468 up-regulated genes and 2273 down-regulated significant genes. Functional annotation and pathway analysis of these significant DEGs were compared based on Gene Ontology (Biological process), Kyoto Encyclopedia of Genes and Genome (KEGG), and Reactome pathways using g: Profiler, ShinyGO v0.76, and iDEP.951 web tools. On finding network visualization, the most over-represented and correlated pathways were neuronal and sensory organ development, calcium signalling pathway, Mitogen-activated protein kinase (MAPK) and Smad signalling pathway, Ras-proximate-1, or Ras-related protein 1 (Rap 1) signalling pathway, apoptosis, and oxidative stress. Similarly, down-regulated genes were most expressed in mRNA processing, immune system, B-cell receptor signalling pathway, Nucleotide oligomerization domain (NOD)-like receptors (NLRs) signalling pathway and nonsense-mediated decay (NMD) pathway. The heat stress-responsive genes identified in this study will facilitate our understanding of the molecular basis for climate resilience and heat tolerance in Indian dairy crossbreeds.

Elucidation of population stratifying markers and selective sweeps in crossbred Landlly pig population using genome-wide SNP data

The present study was aimed at the identification of population stratifying markers from the commercial porcine SNP 60K array and elucidate the genome-wide selective sweeps in the crossbred Landlly pig population. Original genotyping data, generated on Landlly pigs, was merged in various combinations with global suid breeds that were grouped as exotic (global pig breeds excluding Indian and Chinese), Chinese (Chinese pig breeds only), and outgroup pig populations. Post quality control, the genome-wide SNPs were ranked for their stratifying power within each dataset in TRES (using three different criteria) and FIFS programs and top-ranked SNPs (0.5K, 1K, 2K, 3K, and 4K densities) were selected. PCA plots were used to assess the stratification power of low-density panels. Selective sweeps were elucidated in the Landlly population using intra- and inter-population haplotype statistics. Additionally, Tajima's D-statistics were calculated to determine the status of balancing selection in the Landlly population. PCA plots showed 0.5K marker density to effectively stratify Landlly from other pig populations. The A-score in DAPC program revealed the Delta statistic of marker selection to outperform other methods (informativeness and FST methods) and that 3000-marker density was suitable for stratification of Landlly animals from exotic pig populations. The results from selective sweep analysis revealed the Landlly population to be under selection for mammary (NAV2), reproductive efficiency (JMY, SERGEF, and MAP3K20), body conformation (FHIT, WNT2, ASRB, DMGDH, and BHMT), feed efficiency (CSRNP1 and ADRA1A), and immunity (U6, MYO3B, RBMS3, and FAM78B) traits. More than two methods suggested sweeps for immunity and feed efficiency traits, thus giving a strong indication for selection in this direction. The study is the first of its kind in Indian pig breeds with a comparison against global breeds. In conclusion, 500 markers were able to effectively stratify the breeds. Different traits under selective sweeps (natural or artificial selection) can be exploited for further improvement.

Joint Analysis of Genome-wide DNA Methylation and Transcription Sequencing Identifies the Role of BAX Gene in Heat Stress-Induced-Sertoli Cells Apoptosis

The problem of male infertility is a global health crisis and poses a serious threat to the well-being of families. Under heat stress (HS), the reduction of Sertoli cells (SCs) inhibits energy transport and nutrient supply to germ cells, leading to spermatogenesis failure. DNA methylation of genes is a central epigenetic regulatory mechanism in mammalian reproduction. However, it remains unclear how DNA methylation regulates gene expression in heat-stressed SCs. In this study, we investigated whether the decrease in SC levels during HS could be related to epigenetic DNA modifications. The cells exposed to HS showed changes in differential methylation cytosines and regions (DMCs/DMRs) and differential expression genes (DEGs), but not in global DNA methylations. One of the most important biological processes affected by HS is cell apoptosis induced by the intrinsic apoptotic signaling pathway (GO: 2,001,244, P < 0.05) by enrichment in the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). The joint analysis showed that several gene expressions in RNA-seq and WGBS overlapped and the shortlisted genes BAX, HSPH1, HSF1B, and BAG were strongly correlated with stress response and apoptosis. Methylation-specific PCR (MSP) and flow cytometry (FCM) analyzes showed that reduced promoter methylation and enhanced gene expression of BAX with a consequence of apoptosis. The activity of BAX, as well as an increase in its expression, is likely to result in a reduction of SCs population which could further impair ATP supply and adversely affect membrane integrity. These findings provide novel insights into the molecular mechanisms through which stressors cause male reproductive dysfunction and a new molecular etiology of male infertility.

Genome-Wide DNA Methylation Analysis and Functional Validation of Litter Size Traits in Jining Grey Goats

DNA methylation (DNAm) is associated with the reproductive system. However, the genetic mechanism through which DNAm regulates gene expression and thus affects litter size in goats is unclear. Therefore, in the present work, genome-wide DNAm profiles of HP and LP Jining Grey goat ovary tissues were comprehensively analyzed via WGBS, and RNA-Seq data were combined to identify candidate genes associated with litter size traits in the Jining Grey goat. Finally, BSP and RT-qPCR were used to verify the sequencing results of the key genes. Notably, the DNMT genes were downregulated at the expression level in the HP group. Both groups exhibited comparable levels of methylation. A total of 976 differentially methylated regions (DMRs) (973 DMRs for CG and 3 DMRs for CHG) and 310 differentially methylated genes (DMGs) were identified in this study. Through integration of WGBS and RNA-Seq data, we identified 59 differentially methylated and differentially expressed genes (DEGs) and ultimately screened 8 key DMGs (9 DMRS) associated with litter size traits in Jining Grey goats (SERPINB2: chr24_62258801_62259000, NDRG4: chr18_27599201_27599400, CFAP43: chr26_27046601_27046800, LRP1B. chr2_79720201_79720400, EPHA6: chr1_40088601_40088800, TTC29: chr17_59385801_59386000, PDE11A: chr2_117418601_117418800 and PGF: chr10_ 16913801_16914000 and chr10_16916401_16916600). In summary, our research comprehensively analyzed the genome-wide DNAm profiles of HP and LP Jining Grey goat ovary tissues. The data findings suggest that DNAm in goat ovaries may play an important role in determining litter size.

Influence of e-waste exposure on DNA damage and DNA methylation in people living near recycling sites

The association between long-term exposure to e-waste and poor health is well established, but how e-waste exposure affects DNA methylation is understudied. In this study, we measured the DNA damage levels and the alternation of DNA methylation in peripheral blood mononuclear cells (PBMCs) collected from a population exposed to e-waste. The concentration of 28 PCB congeners in the blood samples of e-waste recycling workers was elevated than those of the reference group. DNA damage levels were significantly higher than that of samples from the reference group by detecting the SCGE, CA, and CBMN assays. Eventually, we found that the methylation level of 1233 gene loci was changed in the exposure group. Bioinformatic analysis of differential genes revealed that the hypermethylated genes were enriched in cell component movement and regulation of cell function, and hypomethylated genes were involved in the cellular metabolic process. Among the 30 genes we tested, 14 genes showed a negative correlation between methylation level and expression level. Therefore, e-waste exposure potentially increased the levels of DNA damage and alters DNA methylation, which would likely impact human health.

Integrated Analysis of DNA Methylome and Transcriptome Reveals Epigenetic Regulation of Cold Tolerance in Litopenaeus vannamei

DNA methylation is an important epigenetic modification that has been shown to be associated with responses to non-biological stressors. However, there is currently no research on DNA methylation in response to environmental signals in shrimp. In this study, we conducted a comprehensive comparative analysis of DNA methylation profiles and differentially expressed genes between two strains of Litopenaeus vannamei with significantly different cold tolerance through whole genome bisulfite sequencing (WGBS) and transcriptome sequencing. Between Lv-C and Lv-T (constant temperature of 28 °C and low temperatures of 18 °C and 10 °C) under cytosine-guanine (CG) environments, 39,100 differentially methylated regions (DMRs) were identified, corresponding to 9302 DMR-related genes (DMRGs). The DMRs were mainly located in the gene body (exons and introns). Gene Ontology (GO) analysis showed that these DMRGs were significantly enriched in cell parts, catalytic activity, and metabolic processes. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed significant enrichment of these DMRGs in pathways such as proteasome (ko03050), oxidative phosphorylation (ko00190), mTOR signaling pathway (ko04150), fatty acid metabolism (ko01212), and fatty acid degradation (ko00071). The comprehensive results suggested that L. vannamei mainly regulates gene expression in response to low temperatures through hypermethylation or demethylation of some genes involved in thermogenesis, glycolysis, the autophagy pathway, the peroxisome, and drug metabolism pathways. These results provide important clues for studying DNA methylation patterns and identifying cold tolerance genes in shrimp.

Pig fetal skeletal muscle development is associated with genome-wide DNA hypomethylation and corresponding alterations in transcript and microRNA expression

Fetal myogenesis represents a critical period of porcine skeletal muscle development and requires coordinated expression of thousands of genes. Epigenetic mechanisms, including DNA methylation, drive transcriptional regulation during development; however, these processes are understudied in developing porcine tissues. We performed bisulfite sequencing to assess DNA methylation in pig longissimus dorsi muscle at 41- and 70-days gestation (dg), as well as RNA- and small RNA-sequencing to identify coordinated changes in methylation and expression between myogenic stages. We identified 45 739 differentially methylated regions (DMRs) between stages, and the majority (N = 34 232) were hypomethylated at 70 versus 41 dg. Integration of methylation and transcriptomic data revealed strong associations between differential gene methylation and expression. Differential miRNA methylation was significantly negatively correlated with abundance, and dynamic expression of assayed miRNAs persisted postnatally. Motif analysis revealed significant enrichment of myogenic regulatory factor motifs among hypomethylated regions, suggesting that DNA hypomethylation may function to increase accessibility of muscle-specific transcription factors. We show that developmental DMRs are enriched for GWAS SNPs for muscle- and meat-related traits, demonstrating the potential for epigenetic processes to influence phenotypic diversity. Our results enhance understanding of DNA methylation dynamics of porcine myogenesis and reveal putative cis-regulatory elements governed by epigenetic processes.

Transcriptome and DNA Methylation Responses in the Liver of Yellowfin Seabream Under Starvation Stress

Fish suffer from starvation due to environmental risks such as extreme weather in the wild and due to insufficient feedings in farms. Nutrient problems from short-term or long-term starvation conditions can result in stress-related health problems for fish. Yellowfin seabream (Acanthopagrus latus) is an important marine economic fish in China. Understanding the molecular responses to starvation stress is vital for propagation and culturing yellowfin seabream. In this study, the transcriptome and genome-wide DNA methylation levels in the livers of yellowfin seabream under 14-days starvation stress were analyzed. One hundred sixty differentially expressed genes (DEGs) by RNA-Seq analysis and 737 differentially methylated-related genes by whole genome bisulfite sequencing analysis were identified. GO and KEGG pathway enrichment analysis found that energy metabolism-related pathways such as glucose metabolism and lipid metabolism were in response to starvation. Using bisulfite sequencing PCR, we confirmed the presence of CpG methylation differences within the regulatory region of a DEG ppargc1a in response to 14-days starvation stress. This study revealed the molecular responses of livers in response to starvation stress at the transcriptomic and whole genome DNA methylation levels in yellowfin seabream.

RNA-Seq Analysis Identifies Differentially Expressed Genes in the Longissimus dorsi of Wagyu and Chinese Red Steppe Cattle

Meat quality has a close relationship with fat and connective tissue; therefore, screening and identifying functional genes related to lipid metabolism is essential for the production of high-grade beef. The transcriptomes of the Longissimus dorsi muscle in Wagyu and Chinese Red Steppe cattle, breeds with significant differences in meat quality and intramuscular fat deposition, were analyzed using RNA-seq to screen for candidate genes associated with beef quality traits. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the 388 differentially expressed genes (DEGs) were involved in biological processes such as short-chain fatty acid metabolism, regulation of fatty acid transport and the peroxisome proliferator-activated receptor (PPAR) signaling pathway. In addition, crystallin alpha B (CRYAB), ankyrin repeat domain 2 (ANKRD2), aldehyde dehydrogenase 9 family member A1 (ALDH9A1) and enoyl-CoA hydratase and 3-hydroxyacyl CoA dehydrogenase (EHHADH) were investigated for their effects on intracellular triglyceride and fatty acid content and their regulatory effects on genes in lipogenesis and fatty acid metabolism pathways. This study generated a dataset from transcriptome profiling of two cattle breeds, with differing capacities for fat-deposition in the muscle, and revealed molecular evidence that CRYAB, ANKRD2, ALDH9A1 and EHHADH are related to fat metabolism in bovine fetal fibroblasts (BFFs). The results provide potential functional genes for maker-assisted selection and molecular breeding to improve meat quality traits in beef cattle.

Thermoregulation mechanisms and perspectives for validating thermal windows in pigs with hypothermia and hyperthermia: An overview

Specific anatomical characteristics make the porcine species especially sensitive to extreme temperature changes, predisposing them to pathologies and even death due to thermal stress. Interest in improving animal welfare and porcine productivity has led to the development of various lines of research that seek to understand the effect of certain environmental conditions on productivity and the impact of implementing strategies designed to mitigate adverse effects. The non-invasive infrared thermography technique is one of the tools most widely used to carry out these studies, based on detecting changes in microcirculation. However, evaluations using this tool require reliable thermal windows; this can be challenging because several factors can affect the sensitivity and specificity of the regions selected. This review discusses the thermal windows used with domestic pigs and the association of thermal changes in these regions with the thermoregulatory capacity of piglets and hogs.

The effects of temperature and donor piglet age on the transcriptomic profile and energy metabolism of myoblasts

Rapid climate change is associated with frequent extreme heat events and the resulting thermal stress has consequences for the health, welfare, and growth of farm animals. The aim of this study was to characterize the transcriptional changes and the effects on energy metabolism in proliferating porcine myoblasts derived from piglets of different ages, representing differences in thermoregulatory abilities, and cultivated below (35°C) and above (39°C, 41°C) the standard cultivation temperature (37°C). Satellite cells originating from Musculus rhomboideus of piglets isolated on days 5 (P5, thermolabile) and 20 (P20, thermostable) of age were used. Our expression analyses highlighted differentially expressed genes in porcine myoblasts cultures under heat or cold induced stress. These gene sets showed enrichment for biological processes and pathways related to organelle fission, cell cycle, chromosome organization, and DNA replication. Culture at 35°C resulted in increased metabolic flux as well as a greater abundance of transcripts of the cold shock protein-encoding gene RBM3 and those of genes related to biological processes and signaling pathways, especially those involving the immune system (cytokine–cytokine receptor interaction, TNF and IL-17 signaling pathways). For cultivation at 39°C, differences in the expression of genes related to DNA replication and cell growth were identified. The highest glutathione index ratio was also found under 39°C. Meanwhile, cultivation at 41°C induced a heat stress response, including the upregulation of HSP70 expression and the downregulation of many biological processes and signaling pathways related to proliferative ability. Our analysis also identified differentially expressed genes between cells of donors with a not yet (P5) and already fully developed (P20) capacity for thermoregulation at different cultivation temperatures. When comparing P5 and P20, most of the changes in gene expression were detected at 37°C. At this optimal temperature, muscle cells can develop to their full capacity. Therefore, the most diverse molecular signaling pathways, including PI3K-Akt signaling, Wnt signaling, and EGFR tyrosine kinase inhibitor, were found and are more pronounced in muscle cells from 20-day-old piglets. These results contribute to a better understanding of the mechanisms underlying the adaptation of skeletal muscle cells to temperature stress in terms of their thermoregulatory ability.

Epigenomics of rats' liver and its cross-species functional annotation reveals key regulatory genes underlying short term heat-stress response

Molecular responses to heat stress are multifaceted and under a complex cellular post-transcriptional control. This study explores the epigenetic and transcriptional alterations induced by heat stress (42 °C for 120 min) in the liver of rats, by integrating ATAC-seq, RNA-Seq, and WGBS information. Out of 2586 differential ATAC-seq peaks induced by heat stress, 36 up-regulated and 22 down-regulated transcript factors (TFs) are predicted, such as Cebpα, Foxa2, Foxo4, Nfya and Sp3. Furthermore, 150,189 differentially methylated regions represent 2571 differentially expressed genes (DEGs). By integrating all data, 45 DEGs are concluded as potential heat stress response markers in rats. To comprehensively annotate and narrow down predicted markers, they are integrated with GWAS results of heat stress parameters in cows, and PheWAS data in humans. Besides better understanding of heat stress responses in mammals, INSR, MAPK8, RHPN2 and BTBD7 are proposed as candidate markers for heat stress in mammals.

Effect of acute heat shock on stress gene expression and DNA methylation in zebu (Bos indicus) and crossbred (Bos indicus × Bos taurus) dairy cattle

Environmental temperature is one of the major factors to affect health and productivity of dairy cattle. Gene expression networks within the cells and tissues coordinate stress response, metabolism, and milk production in dairy cattle. Epigenetic DNA methylations were found to mediate the effect of environment by regulating gene expression patterns. In the present study, we compared three Indian native zebu cattle, Bos indicus (Sahiwal, Tharparkar, and Hariana) and one crossbred Bos indicus × Bos taurus (Vrindavani) for stress gene expression and differences in the DNA methylation patterns. The results indicated acute heat shock to cultured PBMC affected their proliferation, stress gene expression, and DNA methylation. Interestingly, expressions of HSP70, HSP90, and STIP1 were found more pronounced in zebu cattle than the crossbred cattle. However, no significant changes were observed in global DNA methylation due to acute heat shock, even though variations were observed in the expression patterns of DNA methyltransferases (DNMT1, DNMT3a) and demethylases (TET1, TET2, and TET3) genes. The treatment 5-AzaC (5-azacitidine) that inhibit DNA methylation in proliferating PBMC caused significant increase in heat shock-induced HSP70 and STIP1 expression indicating that hypomethylation facilitated stress gene expression. Further targeted analysis DNA methylation in the promoter regions revealed no significant differences for HSP70, HSP90, and STIP1. However, there was a significant hypomethylation for BDNF in both zebu and crossbred cattle. Similarly, NR3C1 promoter region showed hypomethylation alone in crossbred cattle. Overall, the results indicated that tropically adapted zebu cattle had comparatively higher expression of stress genes than the crossbred cattle. Furthermore, DNA methylation may play a role in regulating expression of certain genes involved in stress response pathways.

Gestational heat stress alters skeletal muscle gene expression profiles and vascularity in fetal pigs in a sexually dimorphic manner

BACKGROUND

There is evidence that sow heat stress (HS) during gestation affects fetal development with implications for impaired muscle growth. We have previously demonstrated that maternal HS during early to mid-gestation compromised muscle fibre hyperplasia in developing fetal pigs. Thus, we hypothesised these phenotypic changes are associated with a change in expression of genes regulating fetal skeletal muscle development and metabolism. To test this, at d 60 of gestation, RNA sequencing and immunohistochemistry were performed on fetal longissimus dorsi (LD) muscle biopsies collected from pregnant gilts that had experienced either thermoneutral control (CON, 20 °C, n = 7 gilts, 18 LD samples) or controlled HS (cyclic 28 to 33 °C, n = 8 gilts, 23 LD samples) conditions for 3 weeks.

RESULTS

A total of 282 genes were differentially expressed between the HS and CON groups in female LD muscles (false discovery rate (FDR) ≤ 0.05), whereas no differentially expressed genes were detected in male LD muscles between the two groups (FDR > 0.05). Gestational HS increased the expression of genes associated with transcription corepressor activity, adipogenesis cascades, negative regulation of angiogenesis and pro-inflammatory signalling in female LD muscles. Immunohistochemical analyses revealed a decreased muscle vascularity density in fetuses from HS group for both sexes compared to those from the CON group (P = 0.004).

CONCLUSIONS

These results reveal gilt HS during early to mid-gestation altered gene expression profiles in fetal LD muscles in a sexually dimorphic manner. The molecular responses, including transcription and angiogenesis repressions and enhanced adipogenesis cascades, were exclusively observed in females. However, the associated reductions in muscle vascularity were observed independently of sexes. Collectively this may indicate female fetal pigs are more adaptive to gestational HS in terms of gene expression changes, and/or there may be sexually dimorphic differences with respect to the timing of muscle molecular responses to gestational HS.

A draft genome of Drung cattle reveals clues to its chromosomal fusion and environmental adaptation

Drung cattle (Bos frontalis) have 58 chromosomes, differing from the Bos taurus 2n = 60 karyotype. To date, its origin and evolution history have not been proven conclusively, and the mechanisms of chromosome fusion and environmental adaptation have not been clearly elucidated. Here, we assembled a high integrity and good contiguity genome of Drung cattle with 13.7-fold contig N50 and 4.1-fold scaffold N50 improvements over the recently published Indian mithun assembly, respectively. Speciation time estimation and phylogenetic analysis showed that Drung cattle diverged from Bos taurus into an independent evolutionary clade. Sequence evidence of centromere regions provides clues to the breakpoints in BTA2 and BTA28 centromere satellites. We furthermore integrated a circulation and contraction-related biological process involving 43 evolutionary genes that participated in pathways associated with the evolution of the cardiovascular system. These findings may have important implications for understanding the molecular mechanisms of chromosome fusion, alpine valleys adaptability and cardiovascular function.

Exploring crucial molecular events in pearl oyster after pre-grafting conditioning by genome-wide bisulfite sequencing for DNA methylation analysis

Pre-grafting condition is an important method to promote recovery from transplant surgery during pearl production. In the present study, we constructed two DNA methylomes from pearl oysters with and without conditioning to investigate the molecular mechanism of the pearl oyster Pinctada fucata martensii underlying the pre-grafting condition. A total of 4,594,997 and 4,930,813 methyl CG in the control (Con) and pre-grafting group (PT) were detected, resulting in the whole genome methylation profile and methylation pattern in P. f. martensii. Results reveal that the promoter, especially the CpG island-rich region, was more infrequently methylated than the gene function elements in P. f. martensii. A total of 51,957 differently methylated regions (DMRs) between Con and PT were obtained, including 3789 DMR in the promoter and 16,021 in the gene body. Based on gene ontology and pathway enrichment analyses, these DMRs were mainly related to "cellular process", "metabolic process", "Epstein-Barr virus infection", and "Fanconi anemia pathway". The methylation site in the promoter region may be associated with the promoter activity and transcription factor binding. These results help our understanding of the mechanism of pre-grafting condition, thereby providing key information in guiding to improve the conditioning methods for enhanced pearl oyster survival rate after transplantation.

A cross-species comparative approach to assessing multi- and transgenerational effects of endocrine disrupting chemicals

A wide range of chemicals have been identified as endocrine disrupting chemicals (EDCs) in vertebrate species. Most studies of EDCs have focused on exposure of both male and female adults to these chemicals; however, there is clear evidence that EDCs have dramatic effects when mature or developing gametes are exposed, and consequently are associated with in multigenerational and transgenerational effects. Several publications have reviewed such actions of EDCs in subgroups of species, e.g., fish or rodents. In this review, we take a holistic approach synthesizing knowledge of the effects of EDCs across vertebrate species, including fish, anurans, birds, and mammals, and discuss the potential mechanism(s) mediating such multi- and transgenerational effects. We also propose a series of recommendations aimed at moving the field forward in a structured and coherent manner.

B1 repetitive sequence methylation enhances wound healing of second‑degree burns in rats

The accumulation of DNA damage in burn wounds delays wound healing. DNA methylation by short interspersed nuclear element (SINE) small interfering (si)RNA prevents DNA damage and promotes cell proliferation. Therefore, SINE siRNA may be able to promote burn wound healing. Here, a SINE B1 siRNA was used to treat burn wounds in rats. Second-degree burn wounds were introduced on the backs of rats. The rats were then divided into three groups: a B1 siRNA-treated, saline-treated control, and saline + calcium phosphate-nanoparticle-treated control group (n=15/group). The wounds were imaged on days 0, 7, 14, 21 and 28 post-injury. The tissue sections were processed for methylation, histological and immunohistochemical examination, and scored based on the overall expression of histone H2AX phosphorylated on serine 139 (γH2AX) and 8-hydroxy-2'-deoxyguanosine (8-OHdG). Burn wound closure improved in the B1 siRNA-treated group compared with that in the control group, especially from days 14-28 post-injury (P<0.001). The overall pathological score and degree of B1 methylation in the B1 siRNA-treated group improved significantly at days 14-28 post-injury, with the maximum improvement observed on day 14 (P<0.01) compared with the NSS and Ca-P nanoparticle groups. Immunohistochemical staining revealed lower expression of γH2AX and 8-OHdG in the B1 siRNA-treated group than in the control groups at days 14-28 post-injury; the maximum improvement was observed on days 14 and 21. These data imply that administering SINE siRNA is a promising therapeutic option for managing second-degree burns.

Exposure of pigs to glyphosate affects gene-specific DNA methylation and gene expression

Glyphosate (N-(phosphonomethyl)glycine) is a broad-spectrum systemic herbicide and crop desiccant. Glyphosate has long been suspected of leading to the development of cancer and of compromising fertility. Herbicides have been increasingly recognized as epigenetic modifiers, and the impact of glyphosate on human and animal health might be mediated by epigenetic modifications. This article presents the results from an animal study where pigs were exposed to glyphosate while feeding. The experimental setup included a control group with no glyphosate added to the feed and two groups of pigs with 20 ppm and 200 ppm of glyphosate added to the feed, respectively. After exposure, the pigs were dissected, and tissues of the small intestine, liver, and kidney were used for DNA methylation and gene expression analyses. No significant change in global DNA methylation was found in the small intestine, kidney, or liver. Methylation status was determined for selected genes involved in various functions such as DNA repair and immune defense. In a CpG island of the promoter for IL18, we observed significantly reduced DNA methylation for certain individual CpG positions. However, this change in DNA methylation had no influence on IL18 mRNA expression. The expression of the DNA methylation enzymes DNMT1, DNMT3A, and DNMT3B was measured in the small intestine, kidney, and liver of pigs exposed to glyphosate. No significant changes in relative gene expression were found for these enzymes following dietary exposure to 20 and 200 ppm glyphosate. In contrast, a significant increase in expression of the enzyme TET3, responsible for demethylation, was observed in kidneys exposed to 200 ppm glyphosate.

•

A large animal study with exposure of pigs to glyphosate is presented here.

•

Pigs were exposed to 20 ppm and 200 ppm in the diet.

•

No significant changes in global DNA methylation was observed.

•

A significantly reduced DNA methylation was found in the porcine IL18 promoter.

•

A significant increase in TET3 expression was seen in porcine kidneys exposed to 200 ppm glyphosate.

How Epigenetics Can Enhance Pig Welfare?

Epigenetics works as an interface between the individual and its environment to provide phenotypic plasticity to increase individual adaptation capabilities. Recently, a wide variety of epi-genetic findings have indicated evidence for its application in the development of putative epi-biomarkers of stress in farm animals. The purpose of this study was to evaluate previously reported stress epi-biomarkers in swine and encourage researchers to investigate potential paths for the development of a robust molecular tool for animal welfare certification. In this literature review, we report on the scientific concerns in the swine production chain, the management carried out on the farms, and the potential implications of these practices for the animals' welfare and their epigenome. To assess reported epi-biomarkers, we identified, from previous studies, potentially stress-related genes surrounding epi-biomarkers. With those genes, we carried out a functional enrichment analysis of differentially methylated regions (DMRs) of the DNA of swine subjected to different stress-related conditions (e.g., heat stress, intrauterine insult, and sanitary challenges). We identified potential epi-biomarkers for target analysis, which could be added to the current guidelines and certification schemes to guarantee and certify animal welfare on farms. We believe that this technology may have the power to increase consumers' trust in animal welfare.

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.

Whole Genome DNA Methylation Variations in Mammary Gland Tissues from Holstein Cattle Producing Milk with Various Fat and Protein Contents

Milk fat and protein contents are among key elements of milk quality, and they are attracting more attention in response to consumers' demand for high-quality dairy products. To investigate the potential regulatory roles of DNA methylation underlying milk component yield, whole genome bisulfite sequencing was employed to profile the global DNA methylation patterns of mammary gland tissues from 17 Canada Holstein cows with various milk fat and protein contents. A total of 706, 2420 and 1645 differentially methylated CpG sites (DMCs) were found between high vs. low milk fat (HMF vs. LMF), high vs. low milk protein (HMP vs. LMP), and high vs. low milk fat and protein (HMFP vs. LMFP) groups, respectively (q value < 0.1). Twenty-seven, 56 and 67 genes harboring DMCs in gene regions (denoted DMC genes) were identified for HMF vs. LMF, HMP vs. LMP and HMFP vs. LMFP, respectively. DMC genes from HMP vs. LMP and HMFP vs. LMFP comparisons were significantly overrepresented in GO terms related to aerobic electron transport chain and/or mitochondrial ATP (adenosine triphosphate) synthesis coupled electron transport. A total of 83 (HMF vs. LMF), 708 (HMP vs. LMP) and 408 (HMFP vs. LMFP) DMCs were co-located with 87, 147 and 158 quantitative trait loci (QTL) for milk component and yield traits, respectively. In conclusion, the identified methylation changes are potentially involved in the regulation of milk fat and protein yields, as well as the variation in reported co-located QTLs.

Population genomic, olfactory, dietary, and gut microbiota analyses demonstrate the unique evolutionary trajectory of feral pigs

Domestication is an intriguing evolutionary process. Many domestic populations are subjected to strong human-mediated selection, and when some individuals return to the wild, they are again subjected to selective forces associated with new environments. Generally, these feral populations evolve into something different from their wild predecessors and their members typically possess a combination of both wild and human selected traits. Feralisation can manifest in different forms on a spectrum from a wild to a domestic phenotype. This depends on how the rewilded domesticated populations can readapt to natural environments based on how much potential and flexibility the ancestral genome retains after its domestication signature. Whether feralisation leads to the evolution of new traits that do not exist in the wild or to convergence with wild forms, however, remains unclear. To address this question, we performed population genomic, olfactory, dietary, and gut microbiota analyses on different populations of Sus scrofa (wild boar, hybrid, feral and several domestic pig breeds). Porcine single nucleotide polymorphisms (SNPs) analysis shows that the feral population represents a cluster distinctly separate from all others. Its members display signatures of past artificial selection, as demonstrated by values of F_ST in specific regions of the genome and bottleneck signature, such as the number and length of runs of homozygosity. Generalised F_ST values, reacquired olfactory abilities, diet, and gut microbiota variation show current responses to natural selection. Our results suggest that feral pigs are an independent evolutionary unit which can persist so long as levels of human intervention remain unchanged.

Genome-wide DNA methylation profiles provide insight into epigenetic regulation of red and white muscle development in Chinese perch Siniperca chuatsi

Fish skeletal muscles are composed of spatially well-separated fiber types, namely, red and white muscles with different physiological functions and metabolism. To compare the DNA methylation profiles of the two types of muscle tissues and identify potential candidate genes for the muscle growth and development under epigenetic regulation, genome-wide DNA methylation of the red and white muscle in Chinese perch Siniperca chuatsi were comparatively analyzed using bisulfate sequencing methods. An average of 0.9 billion 150-bp paired-end reads were obtained, of which 86% were uniquely mapped to the genome. Methylation mostly occurred at CG sites at a ratio of 94.43% in the red muscle and 93.16% in the white muscle. The mean methylation levels at C-sites were 5.95% in red muscle and 5.83% in white muscle, whereas the mean methylation levels of CG, CHG, and CHH were 73.23%, 0.62%, and 0.67% in red muscle, and 71.01%, 0.62%, and 0.67% in white muscle, respectively. A total of 4192 differentially methylated genes (DMGs) were identified significantly enriched in cell signaling pathways related to skeletal muscle differentiation and growth. Various muscle-related genes, including myosin gene isoforms and regulatory factors, are differentially methylated in the promoter region between the red and white muscles. Further analysis of the transcriptional expression of these genes showed that the muscle regulatory factors (myf5, myog, pax3, pax7, and twitst2) and myosin genes (myh10, myh16, myo18a, myo7a, myo9a, and myl3) were differentially expressed between the two kinds of muscles, consistent with the DNA methylation analysis results. ELISA assays confirmed that the level of 5mC in red muscle was significantly higher than in white muscle (P < 0.05). The RT-qPCR assays revealed that the expression levels of the three DNA methylation transferase (dnmt) subtypes, dnmt1, dnmt3ab, and dnmt3bb1, were significantly higher in red muscle than in white muscle. The higher DNA methylation levels in the red muscle may result from higher DNA methylation transferase expression in the red muscles. Thus, this study might provide a theoretical foundation to better understand epigenetic regulation in the growth and development of red and white muscles in animals, at least in Chinese perch fish.

Genome-wide profile in DNA methylation in goat ovaries of two different litter size populations

Although some studies have investigated the DNA methylation modification in goat ovaries, it is not understood DNA methylation related to goat litter size. This investigation was designed to explore the DNA methylation status in the ovaries of high litter size and low litter size groups using whole-genome bisulfite sequencing (WGBS). We found that there was global difference on DNA methylation in high litter size and low litter size goat ovaries. Many differentially methylated region-related genes (DMGs) were found in the ovaries of these two different goat populations. Moreover, enrichment analysis discovered that many DMGs were involved in gamete development, reproductive system development, wingless-type MMTV integration site family (WNT) signalling pathways and mitogen-activated protein kinase 1 (MAPK) signalling pathways. The data indicated that DNA methylation in goat ovaries may play important roles in the folliculogenesis, the oocyte ovulation rate and finally the litter size. This study provides a comprehensive analysis of genome-wide DNA methylation patterns in ovaries of high and low litter size goat which helps the understanding of ovarian DNA methylation in relation to goat fertility capability.

Comparative CpG methylation kinetic patterns of cis-regulatory regions of heat stress-related genes in Sahiwal and Frieswal cattle upon persistent heat stress

The kinetic patterns of CpG methylation of the cis-regulatory region of heat stress-related genes on exposed to heat stress (at 42 °C) between the Sahiwal and Frieswal cattle was compared in the present study. Using an in vitro whole blood culture model, cells were continuously exposed to heat stress (at 42 °C) for 6 h. Methylation levels of five genes, viz., GPX1, HSP70, HSP90, c-FOS, and JUN were estimated by SyberGreen-based quantitative methylation-specific PCR (qMSP) assay. CpG methylation kinetics at different time points of heat stress (0.5, 1, 2, 4, 6 h) were analyzed using mixed ANOVA. The initial methylation level, estimated at 37 °C, of HSP70 was significantly high in the Sahiwal breed. A significant (p<0.001) time-dependent hypomethylation of an antioxidant gene (GPX1) CpG islands was detected at the acute phase of the stress. Heat shock protein gene (HSP70) showed a similar CpG methylation kinetics where the hypomethylation was prominent from 1 h and persisted up to 4 h. The heat stress responses of both Sahiwal and Frieswal cattle were identical as there was no distinctiveness in the methylation kinetics of CpG islands of studied genes. The acclimatization of Frieswal cattle-a breed developed in India over the years to the tropical climatic conditions, maybe one of the reasons for this similarity. Thus, the present study results could pave a path to understand the molecular mechanism of heat stress and adaptation of indigenous and crossbred cattle populations to the changing scenario in tropical climate conditions.

Effects of temperature on proliferation of myoblasts from donor piglets with different thermoregulatory maturities

Background

Climate change and the associated risk for the occurrence of extreme temperature events or permanent changes in ambient temperature are important in the husbandry of farm animals. The aim of our study was to investigate the effects of permanent cultivation temperatures below (35 °C) and above (39 °C, 41 °C) the standard cultivation temperature (37 °C) on porcine muscle development. Therefore, we used our porcine primary muscle cell culture derived from satellite cells as an in vitro model. Neonatal piglets have limited thermoregulatory stability, and several days after birth are required to maintain their body temperature. To consider this developmental step, we used myoblasts originating from thermolabile (five days of age) and thermostable piglets (twenty days of age).

Results

The efficiency of myoblast proliferation using real-time monitoring via electrical impedance was comparable at all temperatures with no difference in the cell index, slope or doubling time. Both temperatures of 37 °C and 39 °C led to similar biochemical growth properties and cell viability. Only differences in the mRNA expression of myogenesis-associated genes were found at 39 °C compared to 37 °C with less MYF5, MYOD and MSTN and more MYH3 mRNA. Myoblasts grown at 35 °C are smaller, exhibit higher DNA synthesis and express higher amounts of the satellite cell marker PAX7, muscle growth inhibitor MSTN and metabolic coactivator PPARGC1A. Only permanent cultivation at 41 °C resulted in higher HSP expression at the mRNA and protein levels. Interactions between the temperature and donor age showed that MYOD, MYOG, MYH3 and SMPX mRNAs were temperature-dependently expressed in myoblasts of thermolabile but not thermostable piglets.

Conclusions

We conclude that 37 °C to 39 °C is the best physiological temperature range for adequate porcine myoblast development. Corresponding to the body temperatures of piglets, it is therefore possible to culture primary muscle cells at 39 °C. Only the highest temperature of 41 °C acts as a thermal stressor for myoblasts with increased HSP expression, but it also accelerates myogenic development. Cultivation at 35 °C, however, leads to less differentiated myoblasts with distinct thermogenetic activity. The adaptive behavior of derived primary muscle cells to different cultivation temperatures seems to be determined by the thermoregulatory stability of the donor piglets.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-021-00376-4.

Genome-wide DNA methylation profiling of high-fat emulsion-induced hyperlipidemia mice intervened by a polysaccharide from Cyclocarya paliurus (Batal) Iljinskaja

Genome-wide DNA methylation was used to study the lipid-lowering effect of Cyclocarya paliurus (Batal) Iljinskaja polysaccharide (CPP). The objective of this study was to investigate the hypolipidemic effects and the potential underlying mechanisms of action of CPP-2 in high-fat emulsion (HFE)-induced mice. The results showed that CPP-2 reduced the level of genome-wide DNA methylation in the liver of HFE-induced mice, which had a lipid-lowering effect by regulating the AMP-activated protein kinase (AMPK) signaling-, fatty acid metabolism-, fatty acid biosynthesis- and adipocytokine signaling pathways. A series of lipid metabolism genes were screened out by conjoint analysis of the Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Hereafter, fatty acid synthase (FAS) and peroxisome proliferators-activated receptor α (PPARα) as target genes were selected to validate the accuracy of the results. The findings demonstrated that CPP-2 might be effective in lowering the lipid content, thereby protecting against HFE-induced hyperlipidemia.

Transcriptomics analysis reveals the effect of Broussonetia papyrifera L. fermented feed on meat quality traits in fattening lamb

To date, utilization of feed grains is increasing, which competes for human food. It is imperative to develop and utilize unconventional feed materials. Broussonetia papyrifera L. (B. papyrifera) is a good feeding material with high crude protein, crude fat, and low crude fiber, which is widely distributed in China. In this study, 12 Dorper ♂×Hu ♀  crossbred weaned male lambs were seleted into four groups based on the feed that ratio of the B. papyrifera fermented feed in the total mixed diet (0%, 6%, 18%, and 100%), to character the lambs' longissimus dorsi (LD) fatty acids, morphology and transcriptome. Results showed that the muscle fiber's diameter and area were the smallest in the 100% group. The highest content of beneficial fatty acids and the lowest content of harmful fatty acids in group 18%. RNA-seq identified 443 differentially expressed genes (DEGs) in the LD of lambs from 4 groups. Among these genes, 169 (38.1%) were up-regulated and 274 (61.9%) were down-regulated. The DEGs were mostly enriched in in fatty acid metabolism, arginine and proline metabolism, and PPAR signaling pathways. Our results provide knowledge to understand effect of different ratios of B. papyrifera fermented feed on sheep meat quality traits, also a basis for understanding of the molecular regulation mechanism of B. papyrifera fermented feed affecting on sheep meat quality.

Analysis of long non-coding RNAs in skeletal muscle of Bama Xiang pigs in response to heat stress

The aim of this study was to identify differentially expressed long non-coding RNAs (lncRNAs) molecules and predict their target genes related to muscle development and lipid metabolism in longissimus dorsi (LD) muscles of Bama Xiang pigs under constant heat stress. Ten male Bama Xiang pigs with an average initial body weight of 14 kg were randomly divided into control group (22°C) and heat stress (35 °C) group. The experiment lasted for 28 days. All the pigs were slaughtered at the end of the experiment, and LD muscles were collected for muscle quality analysis and transcriptome sequencing. Heat stress reduced meat quality of Bama Xiang pigs. lncRNAs in LD were identified systematically by deep RNA sequencing between the two groups. The results showed that 365 lncRNAs from the LD were identified, including 128 intergenic lncRNAs, 82 intronic lncRNAs, and 155 anti-sense lncRNAs. The differences lie in transcript of length, number of exons and wider size distribution, and expression level per KB fragment in three subtypes of lncRNAs. The three types of transposable elements coverage, including Line/L1, SINE/tRNA, and LTR/ERVL-MaLR, are the highest in mRNA and the three subtypes of lncRNAs in pigs. lncRNAs and mRNAs were different in comparison of features. The results predicted the target genes of the significant differentially expressed lncRNAs related to muscle development and lipid metabolism. This is the first study to expand the knowledge about muscle-related lncRNAs biology in Bama Xiang pigs under heat stress and will contribute to the development of alleviating the adverse effects of heat stress on pork quality targeting lncRNAs.

Integrative analysis of genome-wide DNA methylation and gene expression profiles reveals important epigenetic genes related to milk production traits in dairy cattle

Epigenetic modification plays a critical role in establishing and maintaining cell differentiation, embryo development, tumorigenesis and many complex diseases. However, little is known about the epigenetic regulatory mechanisms for milk production in dairy cattle. Here, we conducted an epigenome-wide study, together with gene expression profiles to identify important epigenetic candidate genes related to the milk production traits in dairy cattle. Whole-genome bisulphite sequencing and RNA sequencing were employed to detect differentially methylated genes (DMG) and differentially expressed genes (DEG) in blood samples in dry period and lactation period between two groups of cows with extremely high and low milk production performance. A total of 10,877 and 6,617 differentially methylated regions were identified between the two groups in the two periods, which corresponded to 3,601 and 2,802 DMGs, respectively. Furthermore, 156 DEGs overlap with DMGs in comparison of the two groups, and 131 DEGs overlap with DMGs in comparison of the two periods. By integrating methylome, transcriptome and GWAS data, some potential candidate genes for milk production traits in dairy cattle were suggested, such as DOCK1, PTK2 and PIK3R1. Our studies may contribute to a better understanding of epigenetic modification on milk production traits of dairy cattle.

Impacts of Epigenetic Processes on the Health and Productivity of Livestock

The dynamic changes in the epigenome resulting from the intricate interactions of genetic and environmental factors play crucial roles in individual growth and development. Numerous studies in plants, rodents, and humans have provided evidence of the regulatory roles of epigenetic processes in health and disease. There is increasing pressure to increase livestock production in light of increasing food needs of an expanding human population and environment challenges, but there is limited related epigenetic data on livestock to complement genomic information and support advances in improvement breeding and health management. This review examines the recent discoveries on epigenetic processes due to DNA methylation, histone modification, and chromatin remodeling and their impacts on health and production traits in farm animals, including bovine, swine, sheep, goat, and poultry species. Most of the reports focused on epigenome profiling at the genome-wide or specific genic regions in response to developmental processes, environmental stressors, nutrition, and disease pathogens. The bulk of available data mainly characterized the epigenetic markers in tissues/organs or in relation to traits and detection of epigenetic regulatory mechanisms underlying livestock phenotype diversity. However, available data is inadequate to support gainful exploitation of epigenetic processes for improved animal health and productivity management. Increased research effort, which is vital to elucidate how epigenetic mechanisms affect the health and productivity of livestock, is currently limited due to several factors including lack of adequate analytical tools. In this review, we (1) summarize available evidence of the impacts of epigenetic processes on livestock production and health traits, (2) discuss the application of epigenetics data in livestock production, and (3) present gaps in livestock epigenetics research. Knowledge of the epigenetic factors influencing livestock health and productivity is vital for the management and improvement of livestock productivity.

Weaning Induces Stress-Dependent DNA Methylation and Transcriptional Changes in Piglet PBMCs

Changes to the epigenome, including those to DNA methylation, have been proposed as mechanisms by which stress can induce long-term physiological changes in livestock species. Pig weaning is associated with dietary and social stress, both of which elicit an immune response and changes to the hypothalamic–pituitary–adrenal (HPA) axis. While differential methylation following stress has been assessed in model organisms, it remains poorly understood how the pig methylome is altered by stressors in production settings. We quantified changes in CpG methylation and transcript abundance in piglet peripheral blood mononuclear cells (PBMCs) following weaning and also assessed differential patterns in pigs exhibiting high and low stress response as measured by cortisol concentration and lesion scores. Blood was collected from nine gilt piglets 24 h before and after weaning, and whole-genome bisulfite sequencing (WGBS) and RNA-sequencing were performed on six and nine animals, respectively, at both time points. We identified 2,674 differentially methylated regions (DMRs) that were enriched within promoters of genes associated with lymphocyte stimulation and transcriptional regulation. Stress groups displayed unique differential methylation and expression patterns associated with activation and suppression of T cell immunity in low and high stress animals, respectively. Differential methylation was strongly associated with differential expression; specifically, upregulated genes were enriched among hypomethylated genes. We observed post-weaning hypermethylation of the glucocorticoid receptor (NR3C1) promoter and a significant decrease in NR3C1 expression (n = 9, p = 6.1 × 10^–3). Our results indicate that weaning-associated stress elicits genome-wide methylation changes associated with differential gene expression, reduced T cell activation, and an altered HPA axis response.

Analysis of Genome DNA Methylation at Inherited Coat Color Dilutions of Rex Rabbits

Background: The dilution of color in rabbits is associated with many different genetic mechanisms that form different color groups. A number of previous studies have revealed potential regulatory mechanisms by which epigenetics regulate pigmentation. However, the genome-wide DNA methylation involved in animal coat color dilution remains unknown. Results: We compared genome-wide DNA methylation profiles in Rex rabbit hair follicles in a Chinchilla group (Ch) and a diluted Chinchilla group (DCh) through whole-genome bisulfite sequencing (WGBS). Approximately 3.5% of the cytosine sites were methylated in both groups, of which the CG methylation type was in greatest abundance. In total, we identified 126,405 differentially methylated regions (DMRs) between the two groups, corresponding to 11,459 DMR-associated genes (DMGs). Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed that these DMGs were principally involved in developmental pigmentation and Wnt signaling pathways. In addition, two DMRs were randomly selected to verify that the WGBS data were reliable using bisulfite sequencing PCR, and seven DMGs were analyzed to establish the relationship between the level of DNA methylation and mRNA expression using qRT-PCR. Due to the limitation of small sample size, replication of the results with a larger sample size would be important in future studies. Conclusion: These findings provide evidence that there is an association between inherited color dilution and DNA methylation alterations in hair follicles, greatly contributing to our understanding of the epigenetic regulation of rabbit pigmentation.

Comparative genome-wide DNA methylation analysis reveals epigenomic differences in response to heat-humidity stress in Bombyx mori

DNA methylation is an important epigenetic modification and has been shown to be involved in the response to abiotic stress. However, there are few studies on DNA methylation in insect response to environmental signals. In this study, we conducted a comprehensive comparative analysis of DNA methylation profiles between two silkworm strains with significantly different resistance to heat and humidity by whole-genome bisulfite sequencing (WGBS). We identified, in total, 2934 differentially methylated regions (DMRs) between RT_48h (resistant strain with high-temperature/humidity treatment for 48 h) and ST_48h (sensitive strain with high-temperature/humidity treatment for 48 h) under cytosine context (CG), which corresponded to 1230 DMR-related genes (DMGs), and the DMRs were primarily located in the gene body (exon and intron) region. Gene ontology (GO) and KEGG analysis showed that these DMGs were most significantly enriched in binding, cellular metabolic process, and RNA transport pathways. Moreover, 10 DMGs have been revealed to be involved in the heat-humidity stress response in the silkworm. The results of this study indicated that DNA methylation plays crucial roles in silkworm response to environmental stressors and provides important clues to identify key resistance genes in silkworm under high-temperature/humidity stress response.

Genome-Wide Assessment of DNA Methylation in Chicken Cardiac Tissue Exposed to Different Incubation Temperatures and CO2 Levels

Temperature and CO₂ concentration during incubation have profound effects on broiler chick development, and numerous studies have identified significant effects on hatch heart weight (HW) as a result of differences in these parameters. Early life environment has also been shown to affect broiler performance later in life; it has thus been suggested that epigenetic mechanisms may mediate long-term physiological changes induced by environmental stimuli. DNA methylation is an epigenetic modification that can confer heritable changes in gene expression. Using reduced-representation bisulfite sequencing (RRBS), we assessed DNA methylation patterns in cardiac tissue of 84 broiler hatchlings incubated at two egg shell temperatures (EST; 37.8°C and 38.9°C) and three CO₂ concentrations (0.1%, 0.4%, and 0.8%) from day 8 of incubation onward. We assessed differential methylation between EST treatments and identified 2,175 differentially methylated (DM) CpGs (1,121 hypermethylated, 1,054 hypomethylated at 38.9° vs. 37.8°) in 269 gene promoters and 949 intragenic regions. DM genes (DMGs) were associated with heart developmental processes, including cardiomyocyte proliferation and differentiation. We identified enriched binding motifs among DM loci, including those for transcription factors associated with cell proliferation and heart development among hypomethylated CpGs that suggest increased binding ability at higher EST. We identified 9,823 DM CpGs between at least two CO₂ treatments, with the greatest difference observed between 0.8 and 0.1% CO₂ that disproportionately impacted genes involved in cardiac muscle development and response to low oxygen levels. Using HW measurements from the same chicks, we performed an epigenome-wide association study (EWAS) for HW, and identified 23 significantly associated CpGs, nine of which were also DM between ESTs. We found corresponding differences in transcript abundance between ESTs in three DMGs (ABLIM2, PITX2, and THRSP). Hypomethylation of an exonic CpG in PITX2 at 38.9°C was associated with increased expression, and suggests increased cell proliferation in broiler hatchlings incubated at higher temperatures. Overall, these results identified numerous epigenetic associations between chick incubation factors and heart development that may manifest in long-term differences in animal performance.

Genome-Wide DNA Methylation Signatures of Sea Cucumber Apostichopus japonicus during Environmental Induced Aestivation

Organisms respond to severe environmental changes by entering into hypometabolic states, minimizing their metabolic rates, suspending development and reproduction, and surviving critical ecological changes. They come back to an active lifestyle once the environmental conditions are conducive. Marine invertebrates live in the aquatic environment and adapt to environmental changes in their whole life. Sea cucumbers and sponges are only two recently known types of marine organisms that aestivate in response to temperature change. Sea cucumber has become an excellent model organism for studies of environmentally-induced aestivation by marine invertebrates. DNA methylation, the most widely considered epigenetic marks, has been reported to contribute to phenotypic plasticity in response to environmental stress in aquatic organisms. Most of methylation-related enzymes, including DNA methyltransferases, Methyl-CpG binding domain proteins, and DNA demethylases, were up-regulated during aestivation. We conducted high-resolution whole-genome bisulfite sequencing of the intestine from sea cucumber at non-aestivation and deep-aestivation stages. Further DNA methylation profile analysis was also conducted across the distinct genomic features and entire transcriptional units. A different elevation in methylation level at internal exons was observed with clear demarcation of intron/exon boundaries during transcriptional unit scanning. The lowest methylation level occurs in the first exons, followed by the last exons and the internal exons. A significant increase in non-CpG methylation (CHG and CHH) was observed within the intron and mRNA regions in aestivation groups. A total of 1393 genes were annotated within hypermethylated DMRs (differentially methylated regions), and 749 genes were annotated within hypomethylated DMRs. Differentially methylated genes were enriched in the mRNA surveillance pathway, metabolic pathway, and RNA transport. Then, 24 hypermethylated genes and 15 hypomethylated genes were Retrovirus-related Pol polyprotein from transposon (RPPT) genes. This study provides further understanding of epigenetic control on environmental induced hypometabolism in aquatic organisms.

Use of Different Cooling Methods in Pig Facilities to Alleviate the Effects of Heat Stress-A Review

An increase in the frequency of hot periods, which has been observed over the past decades, determines the novel approach to livestock facilities improvement. The effects of heat stress are revealed in disorders in physiological processes, impaired immunity, changes in behaviour and decreases in animal production, thus implementation of cooling technologies is a key factor for alleviating these negative consequences. In pig facilities, various cooling methods have been implemented. Air temperature may be decreased by using adiabatic cooling technology such as a high-pressure fogging system or evaporative pads. In modern-type buildings large-surface evaporative pads may support a tunnel ventilation system. Currently a lot of attention has also been paid to developing energy- and water-saving cooling methods, using for example an earth-air or earth-to-water heat exchanger. The pigs' skin surface may be cooled by using sprinkling nozzles, high-velocity air stream or conductive cooling pads. The effectiveness of these technologies is discussed in this article, taking into consideration the indicators of animal welfare such as respiratory rate, skin surface and body core temperature, performance parameters and behavioural changes.

Time-lagged and acute impact of heat stress on production and fertility traits in the local dual-purpose cattle breed "Rotes Höhenvieh" under pasture-based conditions

Climate change causes rising temperatures and extreme weather events worldwide, with possible detrimental time-lagged and acute impact on production and functional traits of cattle kept in outdoor production systems. The aim of the present study was to infer the influence of mean daily temperature humidity index (mTHI) and number of heat stress days (nHS) from different recording periods on birth weight (BWT), 200 d- and 365 d-weight gain (200 dg, 365 dg) of calves, and on the probability of stillbirth (SB), and calving interval (CINT) of their dams. Data recording included 4,362 observations for BWT, 3,136 observations for 200 dg, 2,502 observations for 365 dg, 9,293 observations for the birth status, and 2,811 observations for CINT of the local dual-purpose cattle breed "Rotes Höhenvieh" (RHV). Trait responses on mTHI and nHS were studied via generalized linear mixed model applications with identity link functions for Gaussian traits (BWT, 200 dg, 365 dg, CINT) and logit link functions for binary SB. High mTHI and high nHS before autumn births had strongest detrimental impact on BWT across all antepartum- (a.p.) periods (34.4 ± 0.79 kg maximum). Prolonged CINT was observed when cows suffered heat stress (HS) before or after spring calvings, with maximum length of 391.6 ± 3.82 d (56 d a.p.-period). High mTHI and high nHS during the 42 d- and 56 d a.p.-period implied increased probabilities for SB. We found a significant (P < 0.05) seasonal effect on SB in model 3 across all a.p.-periods, with the highest probability in autumn (maximum of 5.4 ± 0.82% in the 7 d a.p.-period). Weight gains of calves (200 dg and 365 dg) showed strongest HS response for mTHI and nHS measurements from the long-term postnatal periods (42 d- and 56 d-periods), with minimum 200 dg of 194.2 ± 4.15 kg (nHS of 31 to 42 d in the 42 d-period) or minimum 365 dg of 323.8 ± 3.82 kg (mTHI ≥ 60 in the 42 d-period). Calves born in summer, combined with high mTHI or high nHS pre- or postnatal, had lower weight gains, compared with calves born in other calving seasons or under cooler conditions. Highest BWT, weight gains, and shortest CINT mostly were detected under cool to moderate climate conditions for mTHI, and small to moderate nHS. Results indicate acute and time-lagged HS effects and address possible HS-induced epigenetic modifications of the bovine genome across generations and limited acclimatization processes to heat, especially when heat occurs during the cooler spring and autumn months.

Genome-Wide DNA Methylation and Gene Expression Profiles in Cows Subjected to Different Stress Level as Assessed by Cortisol in Milk

Dairy cattle health, wellbeing and productivity are deeply affected by stress. Its influence on metabolism and immune response is well known, but the underlying epigenetic mechanisms require further investigation. In this study, we compared DNA methylation and gene expression signatures between two dairy cattle populations falling in the high- and low-variant tails of the distribution of milk cortisol concentration (MC), a neuroendocrine marker of stress in dairy cows. Reduced Representation Bisulfite Sequencing was used to obtain a methylation map from blood samples of these animals. The high and low groups exhibited similar amounts of methylated CpGs, while we found differences among non-CpG sites. Significant methylation changes were detected in 248 genes. We also identified significant fold differences in the expression of 324 genes. KEGG and Gene Ontology (GO) analysis showed that genes of both groups act together in several pathways, such as nervous system activity, immune regulatory functions and glucocorticoid metabolism. These preliminary results suggest that, in livestock, cortisol secretion could act as a trigger for epigenetic regulation and that peripheral changes in methylation can provide an insight into central nervous system functions.

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.

Genome-Wide Transcriptome and Metabolome Analyses Provide Novel Insights and Suggest a Sex-Specific Response to Heat Stress in Pigs

Heat stress (HS) negatively impacts pig production and swine health. Therefore, to understand the genetic and metabolic responses of pigs to HS, we used RNA-Seq and high resolution magic angle spinning (HR-MAS) NMR analyses to compare the transcriptomes and metabolomes of Duroc pigs (n = 6, 3 barrows and 3 gilts) exposed to heat stress (33 °C and 60% RH) with a control group (25 °C and 60% RH). HS resulted in the differential expression of 552 (236 up, 316 down) and 879 (540 up, 339 down) genes and significant enrichment of 30 and 31 plasma metabolites in female and male pigs, respectively. Apoptosis, response to heat, Toll-like receptor signaling and oxidative stress were enriched among the up-regulated genes, while negative regulation of the immune response, ATP synthesis and the ribosomal pathway were enriched among down-regulated genes. Twelve and ten metabolic pathways were found to be enriched (among them, four metabolic pathways, including arginine and proline metabolism, and three metabolic pathways, including pantothenate and CoA biosynthesis), overlapping between the transcriptome and metabolome analyses in the female and male group respectively. The limited overlap between pathways enriched with differentially expressed genes and enriched plasma metabolites between the sexes suggests a sex-specific response to HS in pigs.

Cross-talk between energy metabolism and epigenetics during temperature stress response in C2C12 myoblasts

Objective: Environmental stress induces disturbances in cell energy metabolism and may cause epigenetic modifications. This study aimed to understand the possible impact of temperature stress (35 °C, 39 °C and 41 °C, compared to control 37 °C) on energy metabolism and epigenetic modifications, such as DNA methylation and histone H4 acetylation, as well as its effects on the expression of genes responsible for epigenetic changes, in mouse skeletal myoblasts (C2C12 cells). Methods: The results showed significantly reduced maximal respiration and spare respiratory capacity under heat stress (39 °C and 41 °C), suggesting that mitochondrial functions were compromised under these conditions. The glycolytic capacity and glycolysis markedly increased following low-temperature stress (35 °C). The results suggested that, under cold stress, cells prefer glycolysis as a rapid compensatory mechanism to meet energy requirements for adaptive thermogenic response. Results: Epigenetic changes (histone H4 acetylation and global DNA methylation) were observed under both heat and cold stress. Among the genes coding for DNA methyltransferases, the Dnmt3a was significantly increased under high-temperature conditions (39 °C and 41 °C), while Dnmt1 expression was significantly increased at low temperature (35 °C), indicating that under these conditions the cells preferred maintenance of methylation to de novo methylation activity. An expression pattern similar to Dnmt3a was observed for Gcn5, encoding for a histone acetyltransferase. The study revealed that temperature stress induced changes in the metabolic profiles, as well as epigenetic modifications, including the dynamics of the key enzymes. Conclusion: The results indicated the existence of crosstalk mechanisms between energy metabolism and epigenetics during cell stress response.

Differences in the fast muscle methylome provide insight into sex-specific epigenetic regulation of growth in Nile tilapia during early stages of domestication

Growth is a complex trait whose variability within a population cannot be explained solely by genetic variation. Epigenetic regulation is often suggested as an important factor shaping the phenotype, but its association with growth can be highly context- and species-dependent. Nevertheless, the mechanisms involved in epigenetic regulation of growth in fish are poorly understood. We have used reduced representation bisulphite sequencing to determine the genome-wide CpG methylation patterns in male and female Nile tilapia of different sizes but at the same early stage of domestication. The average CpG methylation level in the reduced genome representation was 63% across groups but many sites displayed group-specific methylation patterns. The number of differentially methylated (DM) CpGs was much higher when the interaction between sex and weight was included rather than when these factors were considered separately. There were 1128 DM CpGs between large and small females and 970 DM CpGs between large and small males. We have found many growth-related genes associated with DM CpGs, namely map3k5 and akt3 in females and gadd45g and ppargc1a in males. Only 5% of CpG locations associated with growth were common to both sexes. In particular, the autophagy-related gene atg14 displayed a high association of methylation with growth exclusively in males. The sexually dimorphic association between atg14 methylation and growth may uncover novel metabolic mechanisms at play during mouth brooding in Nile tilapia females. Taken together, our data suggest that epigenetic regulation of growth in Nile tilapia involves different gene networks in males and females.

Alteration in bile acids profile in Large White pigs during chronic heat exposure

Bile acids (BAs) are critical for cholesterol homeostasis and new roles in metabolism and endocrinology have been demonstrated recently. It remains unknown whether BA metabolism can be affected by heat stress (HS). The objective of this study was to describe the shifts in serum, hepatic and intestinal BA profiles induced by chronic HS. Twenty-seven Large White pigs weighing 40.8 ± 2.7 kg were assigned to one of the three treatments: a control group (CON, 23 °C), a HS group (33 °C), or a pair-fed group (PF, 23 °C and fed the same amount as HS group) for 21 d. The concentrations of taurine-conjugated BAs (TUDCA and THDCA in serum and TCDCA, TUDCA, THDCA and THCA in liver) were decreased in HS and PF pigs. However, in HS pigs, a reduction in taurine-conjugated BAs (TCBA) correlated with decreased liver genes expression of BA synthesis, conjugation and uptake transport. BA regulated-genes (FXR, TGR5 and FGFR4) in HS pigs and TGR5, FGFR4 and KLβ in PF pigs were down-regulated in liver. In ileum, total BAs and glycoursodeoxycholic acid concentrations were higher in HS pigs than other groups and PF group, respectively (P < 0.05). TCBA (P = 0.01) and tauroursodeoxycholic acid (P < 0.01) were decreased in PF group. BA transporters (OSTα and MRP3) were up-regulated in HS pigs compared with CON and PF pigs, respectively (P < 0.01). In cecum, ursodeoxycholic acid was higher in HS (P = 0.02) group than CON group. The expression of apical sodium-coupled bile acid transporter (P = 0.04) was lower in HS pigs than CON pigs, while OSTβ (P < 0.01) was greater in HS group than PF group. These results suggest that chronic HS suppressed liver activity of synthesis and uptake of TCBA, at least in part, which was independent of reduced feed intake.

The role of N6-methyladenosine RNA methylation in the heat stress response of sheep (Ovis aries)

With the intensive development of the sheep industry and increasing global temperatures, heat stress in sheep has become an increasingly severe and important issue in recent years. The level of N6-methyladenosine (m6A) RNA methylation changes in response to stress plays important roles in stress responses. However, the role of m6A in the heat stress response of sheep remains unclear. To explore this issue, we measured heat stress protein (HSP) expression, liver function indexes, m6A on RNA, m6A-related enzyme expression, and tissue damage in sheep that had been subjected to heat stress. At the transcriptome level, our results showed significant increases in m6A on RNA and increased mRNA levels of HSPs (HSP70, HSP90, and HSP110) and m6A-related enzymes [METTL3 (methyltransferase-like 3), METTL14 (methyltransferase-like 14), WTAP (wilms tumor 1-associated protein), FTO (fat mass and obesity-associated protein), ALKBH5 (alkB homologue 5), YTHDF1-3 (YTH domain family proteins), and YTHDC1-2 (YTH domain-containing proteins)] following heat stress. At the protein level, the expression of METTL3, YTHDF1-2, and YTHDC2 showed no significant differences following heat stress. However, in contrast to its mRNA level after heat stress, the protein expression of YTHDF3 was reduced, while the expression of HSPs (HSP70, HSP90, and HSP110), METTL14, WTAP, FTO, ALKBH5, YTHDF3, and YTHDC1 increased in line with their measured mRNA levels. Histological experiments revealed that heat stress caused varying degrees of damage to sheep liver tissue. Moreover, immunohistochemical staining indicated that the m6A-related enzymes were expressed in sheep hepatocytes, and differences in expression patterns were observed between the control and heat stress groups. In summary, differences in the level of m6A and the expression of m6A-related enzymes in the liver of sheep were observed after heat stress. This indicates that m6A is involved in the regulation of heat stress in sheep. Our findings provide a new avenue for studying the responses to heat stress in sheep.

Hypomethylation in the promoter region of ZPBP as a potential litter size indicator in Berkshire pigs

In pigs, litter size is typically defined as the total number of piglets born (TNB) or the number of piglets born alive (NBA). Increasing pig litter size is of great economic interest as a means to increase productivity. The capacity of the uterus is a critical component of litter size and may play a central role in prolificacy. In this study, we investigated litter-size-related epigenetic markers in uterine tissue from Berkshire pigs with smaller litter size groups (SLGs) and larger litter size groups (LLGs) using genome-wide bisulfite sequencing (GWBS). A total of 3269 differentially methylated regions (DMRs) were identified: 1566 were hypermethylated and 1703 hypomethylated in LLG compared to SLG. The zona pellucida binding protein (ZPBP) gene was significantly hypomethylated in the LLG promoter region, and its expression was significantly upregulated in uterine tissue. Thus, the methylation status of ZPBP gene was identified as a potential indicator of litter size. Furthermore, we verified its negative correlation with litter size traits (TNB and NBA) in whole blood samples from 172 Berkshire sows as a blood-based biomarker by a porcine methylation-specific restriction enzyme polymerase chain reaction (PMP) assay. The results suggest that the methylation status of the ZPBP gene can serve as a valuable epigenetic biomarker for hyperprolific sows.