Genome-wide DNA methylation profiles and their relationships with mRNA and the microRNA transcriptome in bovine muscle tissue (Bos taurine)

MEF2A promoter methylation negatively regulates mRNA transcription and affects myoblast physiological function in cattle

This study investigates the regulatory effects of methylation in the promoter region of the bovine MEF2A gene on its transcription levels and the impact on bovine myoblasts. Transcription levels and promoter methylation status of MEF2A in the same tissues of calves and adult cattle were assessed using qRT-PCR and BSP methods. The results indicated that MEF2A expression levels in calves were significantly lower than those in adult cattle (P < 0.05), while the methylation rate of MEF2A was significantly higher in calves (P < 0.05), suggesting a correlation between high methylation levels and reduced gene expression. Subsequently, MEF2A overexpression and interference vectors were transfected into bovine myoblasts to examine the effects of altered MEF2A expression on its promoter methylation status. The findings revealed that MEF2A overexpression significantly reduced the methylation rate (P < 0.01), whereas MEF2A interference increased the methylation rate (P < 0.01), aligning with the expression trends of DNMT1. Furthermore, bovine myoblasts were treated with varying concentrations of the methylation inhibitor 5-Aza-dC to evaluate changes in MEF2A promoter methylation and mRNA levels. The effects on cell cycle progression, apoptosis, and other growth parameters were assessed using flow cytometry, ELISA, and qRT-PCR. Results showed that a concentration of 1 μM 5-Aza-dC effectively reduced MEF2A promoter methylation and significantly upregulated MEF2A expression, leading to accelerated cell cycle progression and increased secretion levels of GH and INS, all differences being statistically significant (P < 0.01). Additionally, 1 μM of 5-Aza-dC promoted apoptosis, with qRT-PCR results for relevant genes supporting this finding. In conclusion, methylation of the MEF2A promoter negatively regulates its mRNA transcription levels, thereby impacting the growth and development of Guanling cattle myoblasts. These results provide valuable insights for the genetic improvement of cattle through marker-assisted selection.

Intergenerational effects of maternal rate of body weight gain on the multi-omics hepatic profiles of bovine fetuses

Maternal periconceptual nutrition affects the growth trajectory of developing fetuses by modulating gene expression. The regulatory mechanisms and their role in fetal development remain underexplored in livestock models. Herein, we investigated the effects of maternal rate of body weight (BW) gain during early gestation on the DNA methylation, microRNA profiles, and their interaction with the hepatic gene expression in female fetuses. At breeding, 36 crossbred beef heifers (∼13 months of age) were randomly assigned to a nutritional plane to gain Low (0.28 kg/day; LG, n = 18) or Moderate (0.79 kg/day; MG, n = 18) BW through the first 83 days of gestation. A subset of pregnant heifers (n = 17) was selected, and fetal liver samples were collected on day 83 of gestation for DNA methylation and miRNA-Sequencing. After data quality control, miRDeep2 and Bismark tools were used to analyze miRNA and methylation data, respectively. The bta-miR-206 was the only differentially expressed miRNA (FDR = 0.02). Eight differentially methylated genes were identified (DMGs, FDR < 0.1). The over-represented pathways and biological processes (adj. p < 0.05) for bta-miR-206 targeted genes were associated with embryonic development, energy metabolism, and mineral transport, whereas the DMGs regulated anatomical structural development and transcriptional regulation. Our results show that key genes involved with liver metabolism, tissue structure, and function were regulated by DNA methylation and the miR-206. However, further investigation is warranted to determine physiological responses and long-term consequences on animal performance.

In search of epigenetic hallmarks of different tissues: an integrative omics study of horse liver, lung, and heart

DNA methylation and microRNA (miRNA) expression are epigenetic mechanisms essential for regulating tissue-specific gene expression and metabolic processes. However, high-resolution transcriptome, methylome, or miRNAome data is only available for a few model organisms and selected tissues. Up to date, only a few studies have reported on gene expression, DNA methylation, or miRNA expression in adult equine tissues at the genome-wide level. In the present study, we used RNA-Seq, miRNA-seq, and reduced representation bisulfite sequencing (RRBS) data from the heart, lung, and liver tissues of healthy cold-blooded horses to identify differentially expressed genes (DEGs), differentially expressed miRNA (DE miRNA) and differentially methylated sites (DMSs) between three types of horse tissues. Additionally, based on integrative omics analysis, we described the observed interactions of epigenetic mechanisms with tissue-specific gene expression alterations. The obtained data allowed identification from 4067 to 6143 DMSs, 9733 to 11,263 mRNAs, and 155 to 185 microRNAs, differentially expressed between various tissues. We pointed out specific genes whose expression level displayed a negative correlation with the level of CpG methylation and miRNA expression and revealed biological processes that they enrich. Furthermore, we confirmed and validated the accuracy of the Next-Generation Sequencing (NGS) results with bisulfite sequencing PCR (BSP) and quantitative PCR (qPCR). This comprehensive analysis forms a strong foundation for exploring the epigenetic mechanisms involved in tissue differentiation, especially the growth and development of the equine heart, lungs, and liver.

Long non-coding RNA (LncRNA) and epigenetic factors: their role in regulating the adipocytes in bovine

Investigating the involvement of long non-coding RNAs (lncRNAs) and epigenetic processes in bovine adipocytes can provide valuable new insights into controlling adipogenesis in livestock. Long non-coding RNAs have been associated with forming chromatin loops that facilitate enhancer-promoter interactions during adipogenesis, as well as regulating important adipogenic transcription factors like C/EBPα and PPARγ. They significantly influence gene expression regulation at the post-transcriptional level and are extensively researched for their diverse roles in cellular functions. Epigenetic modifications such as chromatin reorganization, histone alterations, and DNA methylation subsequently affect the activation of genes related to adipogenesis and the progression of adipocyte differentiation. By investigating how fat deposition is epigenetically regulated in beef cattle, scientists aim to unravel molecular mechanisms, identify key regulatory genes and pathways, and develop targeted strategies for modifying fat deposition to enhance desirable traits such as marbling and meat tenderness. This review paper delves into lncRNAs and epigenetic factors and their role in regulating bovine adipocytes while focusing on their potential as targets for genetic improvement to increase production efficiency. Recent genomics advancements, including molecular markers and genetic variations, can boost animal productivity, meeting global demands for high-quality meat products. This review establishes a foundation for future research on understanding regulatory networks linked to lncRNAs and epigenetic changes, contributing to both scholarly knowledge advancement and practical applications within animal agriculture.

An Appraisal for Providing Charge Transfer (CT) Through Synthetic Porous Frameworks for their Semiconductor Applications

In recent years, there has been considerable focus on the development of charge transfer (CT) complex formation as a means to modify the band gaps of organic materials. In particular, CT complexes alternate layers of aromatic molecules with donor (D) and acceptor (A) properties to provide inherent electrical conductivity. In particular, the synthetic porous frameworks as attractive D-A components have been extensively studied in recent years in comparison to existing D-A materials. Therefore, in this work, the synthetic porous frameworks are classified into conjugated microporous polymers (CMPs), covalent organic frameworks (COFs), and metal-organic frameworks (MOFs) and compare high-quality materials for CT in semiconductors. This work updates the overview of the above porous frameworks for CT, starting with their early history regarding their semiconductor applications, and lists CT concepts and selected key developments in their CT complexes and CT composites. In addition, the network formation methods and their functionalization are discussed to provide access to a variety of potential applications. Furthermore, several theoretical investigations, efficiency improvement techniques, and a discussion of the electrical conductivity of the porous frameworks are also highlighted. Finally, a perspective of synthetic porous framework studies on CT performance is provided along with some comparisons.

Identification of Susceptibility Genes Underlying Bovine Respiratory Disease in Xinjiang Brown Cattle Based on DNA Methylation

DNA methylation is a form of epigenetic regulation, having pivotal parts in controlling cellular expansion and expression levels within genes. Although blood DNA methylation has been studied in humans and other species, its prominence in cattle is largely unknown. This study aimed to methodically probe the genomic methylation map of Xinjiang brown (XJB) cattle suffering from bovine respiratory disease (BRD), consequently widening cattle blood methylome ranges. Genome-wide DNA methylation profiling of the XJB blood was investigated through whole-genome bisulfite sequencing (WGBS). Many differentially methylated regions (DMRs) obtained by comparing the cases and controls groups were found within the CG, CHG, and CHH (where H is A, T, or C) sequences (16,765, 7502, and 2656, respectively), encompassing 4334 differentially methylated genes (DMGs). Furthermore, GO/KEGG analyses showed that some DMGs were involved within immune response pathways. Combining WGBS-Seq data and existing RNA-Seq data, we identified 71 significantly differentially methylated (DMGs) and expressed (DEGs) genes (p < 0.05). Next, complementary analyses identified nine DMGs (LTA, STAT3, IKBKG, IRAK1, NOD2, TLR2, TNFRSF1A, and IKBKB) that might be involved in the immune response of XJB cattle infected with respiratory diseases. Although further investigations are needed to confirm their exact implication in the involved immune processes, these genes could potentially be used for a marker-assisted selection of animals resistant to BRD. This study also provides new knowledge regarding epigenetic control for the bovine respiratory immune process.

Diverse WGBS profiles of longissimus dorsi muscle in Hainan black goats and hybrid goats

BACKGROUND

Goat products have played a crucial role in meeting the dietary demands of people since the Neolithic era, giving rise to a multitude of goat breeds globally with varying characteristics and meat qualities. The primary objective of this study is to pinpoint the pivotal genes and their functions responsible for regulating muscle fiber growth in the longissimus dorsi muscle (LDM) through DNA methylation modifications in Hainan black goats and hybrid goats.

METHODS

Whole-genome bisulfite sequencing (WGBS) was employed to scrutinize the impact of methylation on LDM growth. This was accomplished by comparing methylation differences, gene expression, and their associations with growth-related traits.

RESULTS

In this study, we identified a total of 3,269 genes from differentially methylated regions (DMR), and detected 189 differentially expressed genes (DEGs) through RNA-seq analysis. Hypo DMR genes were primarily enriched in KEGG terms associated with muscle development, such as MAPK and PI3K-Akt signaling pathways. We selected 11 hub genes from the network that intersected the gene sets within DMR and DEGs, and nine genes exhibited significant correlation with one or more of the three LDM growth traits, namely area, height, and weight of loin eye muscle. Particularly, PRKG1 demonstrated a negative correlation with all three traits. The top five most crucial genes played vital roles in muscle fiber growth: FOXO3 safeguarded the myofiber's immune environment, FOXO6 was involved in myotube development and differentiation, and PRKG1 facilitated vasodilatation to release more glucose. This, in turn, accelerated the transfer of glucose from blood vessels to myofibers, regulated by ADCY5 and AKT2, ultimately ensuring glycogen storage and energy provision in muscle fibers.

CONCLUSION

This study delved into the diverse methylation modifications affecting critical genes, which collectively contribute to the maintenance of glycogen storage around myofibers, ultimately supporting muscle fiber growth.

Effects of DNA Methylation on Gene Expression and Phenotypic Traits in Cattle: A Review

Gene expression in cells is determined by the epigenetic state of chromatin. Therefore, the study of epigenetic changes is very important to understand the regulatory mechanism of genes at the molecular, cellular, tissue and organ levels. DNA methylation is one of the most studied epigenetic modifications, which plays an important role in maintaining genome stability and ensuring normal growth and development. Studies have shown that methylation levels in bovine primordial germ cells, the rearrangement of methylation during embryonic development and abnormal methylation during placental development are all closely related to their reproductive processes. In addition, the application of bovine male sterility and assisted reproductive technology is also related to DNA methylation. This review introduces the principle, development of detection methods and application conditions of DNA methylation, with emphasis on the relationship between DNA methylation dynamics and bovine spermatogenesis, embryonic development, disease resistance and muscle and fat development, in order to provide theoretical basis for the application of DNA methylation in cattle breeding in the future.

Whole-genome DNA methylation profiling reveals epigenetic signatures in developing muscle in Tan and Hu sheep and their offspring

Introduction

The Tan sheep is a popular local breed in China because of its tenderness and flavor. The Hu sheep breed is also famous for its high litter size, and its muscle growth rate is faster than that of Tan sheep. However, the epigenetic mechanism behind these muscle-related phenotypes is unknown.

Methods

In this study, the longissimus dorsi tissue from 18 6 month-old Tan sheep, Hu sheep, and Tan-Hu F2 generation (6 sheep per population) were collected. After genomic DNA extraction, whole-genome bisulfite sequencing (WGBS) and bioinformatics analysis were performed to construct genome-wide DNA methylome maps for the Tan sheep, Hu sheep and their Tan-Hu F2 generation.

Results

Distinct genome-wide DNA methylation patterns were observed between Tan sheep and Hu sheep. Moreover, DNA methylated regions were significantly increased in the skeletal muscle from Tan sheep vs. the F2 generation compared to the Hu sheep vs. F2 generation and the Tan sheep vs. Hu sheep. Compared with Hu sheep, the methylation levels of actin alpha 1 (ACTA1), myosin heavy chain 11 (MYH11), Wiskott-Aldrich syndrome protein (WAS), vav guanine nucleotide exchange factor 1 (VAV1), fibronectin 1 (FN1) and Rho-associated protein kinase 2 (ROCK2) genes were markedly distinct in the Tan sheep. Furthermore, Gene Ontology analysis indicated that these genes were involved in myotube differentiation, myotube cell development, smooth muscle cell differentiation and striated muscle cell differentiation.

Conclusion

The findings from this study, in addition to data from previous research, demonstrated that the ACTA1, MYH11, WAS, VAV1, FN1, and ROCK2 genes may exert regulatory effects on muscle development.

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.

Global genomic methylation related to the degree of parasitism in cattle

The objective of the present study was to characterize a herd of 72 ½ Angus × ½ Nellore heifers, identify the resistant, resilient and susceptible animals to parasites, relate the overall DNA methylation of these animals with the degree of parasitism, evaluated by the egg count per gram of feces (EPG), Haematobia irritans count (horn fly) and Rhipicephalus microplus count (bovine tick). The experiment was carried out in a completely randomized design, containing 72 treatments, with each animal considered a treatment, and 11 repetitions, with each collection within a year considered a repetition. The data obtained from the counts of the evaluated parasites were subjected to statistical analysis using the SISVAR program, to classify heifers according to the degree of parasitism in low (resistant), intermediary (resilient) and high (susceptible) parasite load for infection by nematodes, infestation by ticks and flies. Addition the animals in these three groups, by hierarchical grouping using the GENES program, heifers were classified as to the degree of parasitism by the three parasites along with the DNA methylation content of the animals in each group. A negative relationship was observed between resistance and methylated DNA content in both classifications, with the resistant, resilient, and susceptible animals showing the highest, intermediate, and lowest methylated DNA quantifications, respectively. Thus, the methodologies used herein enabled the classification of 72 heifers according to the degree of collective infection by gastrointestinal nematodes and infestation by ticks and horn flies, thereby establishing a link between the degree of parasitic resistance in cattle and the global methylated DNA quantification.

Nanodiamond-Quantum Sensors Reveal Temperature Variation Associated to Hippocampal Neurons Firing

Temperature is one of the most relevant parameters for the regulation of intracellular processes. Measuring localized subcellular temperature gradients is fundamental for a deeper understanding of cell function, such as the genesis of action potentials, and cell metabolism. Notwithstanding several proposed techniques, at the moment detection of temperature fluctuations at the subcellular level still represents an ongoing challenge. Here, for the first time, temperature variations (1 °C) associated with potentiation and inhibition of neuronal firing is detected, by exploiting a nanoscale thermometer based on optically detected magnetic resonance in nanodiamonds. The results demonstrate that nitrogen-vacancy centers in nanodiamonds provide a tool for assessing various levels of neuronal spiking activity, since they are suitable for monitoring different temperature variations, respectively, associated with the spontaneous firing of hippocampal neurons, the disinhibition of GABAergic transmission and the silencing of the network. Conjugated with the high sensitivity of this technique (in perspective sensitive to < 0.1 °C variations), nanodiamonds pave the way to a systematic study of the generation of localized temperature gradients under physiological and pathological conditions. Furthermore, they prompt further studies explaining in detail the physiological mechanism originating this effect.

Deep Small RNA Sequencing Reveals Important miRNAs Related to Muscle Development and Intramuscular Fat Deposition in Longissimus dorsi Muscle From Different Goat Breeds

MicroRNAs (miRNAs) are a class of small non-coding RNAs that have been shown to play important post-transcriptional regulatory roles in the growth and development of skeletal muscle tissues. However, limited research into the effect of miRNAs on muscle development in goats has been reported. In this study, Liaoning cashmere (LC) goats and Ziwuling black (ZB) goats with significant phenotype difference in meat production performance were selected and the difference in Longissimus dorsi muscle tissue expression profile of miRNAs between the two goat breeds was then compared using small RNA sequencing. A total of 1,623 miRNAs were identified in Longissimus dorsi muscle tissues of the two goat breeds, including 410 known caprine miRNAs, 928 known species-conserved miRNAs and 285 novel miRNAs. Of these, 1,142 were co-expressed in both breeds, while 230 and 251 miRNAs were only expressed in LC and ZB goats, respectively. Compared with ZB goats, 24 up-regulated miRNAs and 135 miRNAs down-regulated were screened in LC goats. A miRNA-mRNA interaction network showed that the differentially expressed miRNAs would target important functional genes associated with muscle development and intramuscular fat deposition. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that the target genes of differentially expressed miRNAs were significantly enriched in Ras, Rap 1, FoxO, and Hippo signaling pathways. This study suggested that these differentially expressed miRNAs may be responsible for the phenotype differences in meat production performance between the two goat breeds, thereby providing an improved understanding of the roles of miRNAs in muscle tissue of goats.

Feedomics provides bidirectional omics strategies between genetics and nutrition for improved production in cattle

Increasing the efficiency and sustainability of cattle production is an effective way to produce valuable animal proteins for a growing human population. Genetics and nutrition are the 2 major research topics in selecting cattle with beneficial phenotypes and developing genetic potentials for improved performance. There is an inextricable link between genetics and nutrition, which urgently requires researchers to uncover the underlying molecular mechanisms to optimize cattle production. Feedomics integrates a range of omic techniques to reveal the mechanisms at different molecular levels related to animal production and health, which can provide novel insights into the relationships of genes and nutrition/nutrients. In this review, we summarized the applications of feedomics techniques to reveal the effect of genetic elements on the response to nutrition and investigate how nutrients affect the functional genome of cattle from the perspective of both nutrigenetics and nutrigenomics. We highlighted the roles of rumen microbiome in the interactions between host genes and nutrition. Herein, we discuss the importance of feedomics in cattle nutrition research, with a view to ensure that cattle exhibit the best production traits for human consumption from both genetic and nutritional aspects.

Characterization of tRNA expression profiles in large offspring syndrome

Background

Assisted Reproductive Technologies (ART) use can increase the risk of congenital overgrowth syndromes, such as large offspring syndrome (LOS) in ruminants. Epigenetic variations are known to influence gene expression and differentially methylated regions (DMRs) were previously determined to be associated with LOS in cattle. We observed DMRs overlapping tRNA clusters which could affect tRNA abundance and be associated with tissue specificity or overgrowth. Variations in tRNA expression have been identified in several disease pathways suggesting an important role in the regulation of biological processes. Understanding the role of tRNA expression in cattle offers an opportunity to reveal mechanisms of regulation at the translational level. We analyzed tRNA expression in the skeletal muscle and liver tissues of day 105 artificial insemination-conceived, ART-conceived with a normal body weight, and ART-conceived bovine fetuses with a body weight above the 97^th percentile compared to Control-AI.

Results

Despite the centrality of tRNAs to translation, in silico predictions have revealed dramatic differences in the number of tRNA genes between humans and cattle (597 vs 1,659). Consistent with reports in human, only a fraction of predicted tRNA genes are expressed. We detected the expression of 474 and 487 bovine tRNA genes in the muscle and liver with the remainder being unexpressed. 193 and 198 unique tRNA sequences were expressed in all treatment groups within muscle and liver respectively. In addition, an average of 193 tRNA sequences were expressed within the same treatment group in different tissues. Some tRNA isodecoders were differentially expressed between treatment groups. In the skeletal muscle and liver, we categorized 11 tRNA isoacceptors with undetected expression as well as an isodecoder that was unexpressed in the liver (Ser^GGA). Our results identified variation in the proportion of tRNA gene copies expressed between tissues and differences in the highest contributing tRNA anticodon within an amino acid family due to treatment and tissue type. Out of all amino acid families, roughly half of the most highly expressed tRNA isoacceptors correlated to their most frequent codon in the bovine genome.

Conclusion

Although the number of bovine tRNA genes is nearly triple of that of the tRNA genes in human, there is a shared occurrence of transcriptionally inactive tRNA genes in both species. We detected differential expression of tRNA genes as well as tissue- and treatment- specific tRNA transcripts with unique sequence variations that could modulate translation during protein homeostasis or cellular stress, and give rise to regulatory products targeting genes related to overgrowth in the skeletal muscle and/or tumor development in the liver of LOS individuals. While the absence of certain isodecoders may be relieved by wobble base pairing, missing tRNA species could increase the likelihood of mistranslation or mRNA degradation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08496-7.

Global DNA methylation profiles of buffalo (Bubalus bubalis) preimplantation embryos produced by handmade cloning and in vitro fertilization

Somatic cell nuclear transfer technique (SCNT) has proved to be an outstanding method of multiplication of elite animals but accompanied with low efficiency and live birth rate of cloned animals. Epigenetic alterations of DNA has been one of the culprits behind this issue. Cloned embryos are found to deviate slightly from regular pattern of demethylation and re-methylation at the time of nuclear reprogramming and embryonic development when compared with embryos produced by in vitro fertilization (IVF). Thus, the present study was aimed at evaluating global DNA methylation profiles of cloned embryos at 2-cell, 8-cell and blastocyst stages and compare it with corresponding stages of embryos produced by IVF by using MeDIP-Sequencing on Illumina-based platform. We found out that cloned embryos exhibited significantly different DNA methylation pattern as compared to IVF embryos with respect to distribution of differentially methylated regions in different components of genome, CpG islands distribution and methylation status, gene ontological profiles and pathways affected throughout the developmental stages. The data generated from MeDIP-Seq was validated at blastocyst stage cloned and IVF embryos by bisulfite-sequencing PCR on five randomly selected gene regions.

Genome-Wide DNA Methylome and Transcriptome Analysis of Porcine Testicular Cells Infected With Transmissible Gastroenteritis Virus

Transmissible gastroenteritis virus (TGEV) is a porcine pathogen causing highly communicable gastrointestinal infection that are lethal for suckling piglets. In an attempt to delineate the pathogenic mechanism of TGEV-infected porcine testicular cells (ST cells), we conducted a whole genome analysis of DNA methylation and expression in ST cells through reduced bisulfate-seq and RNA-seq. We examined alterations in the methylation patterns and recognized 1764 distinct methylation sites. 385 differentially expressed genes (DEGs) were enriched in the viral defense and ribosome biogenesis pathways. Integrative analysis identified two crucial genes (EMILIN2, RIPOR3), these two genes expression were negatively correlated to promoter methylation. In conclusion, alterations in DNA methylation and differential expression of genes reveal that their potential functional interactions in TGEV infection. Our data highlights the epigenetic and transcriptomic landscapes in TGEV-infected ST cells and provides a reliable dataset for screening TGEV resistance genes and genetic markers.

Genes Related to Fat Metabolism in Pigs and Intramuscular Fat Content of Pork: A Focus on Nutrigenetics and Nutrigenomics

Simple Summary

The intramuscular fat (IMF) or marbling is an essential pork sensory quality that influences the preference of the consumers and premiums for pork. IMF is the streak of visible fat intermixed with the lean within a muscle fibre and determines sensorial qualities of pork such as flavour, tenderness and juiciness. Fat metabolism and IMF development are controlled by dietary nutrients, genes, and their metabolic pathways in the pig. Nutrigenetics explains how the genetic make-up of an individual pig influences the pig’s response to dietary nutrient intake. Differently, nutrigenomics is the analysis of how the entire genome of an individual pig is affected by dietary nutrient intake. The knowledge of nutrigenetics and nutrigenomics, when harmonized, is a powerful tool in estimating nutrient requirements for swine and programming dietary nutrient supply according to an individual pig’s genetic make-up. The current paper aimed to highlight the roles of nutrigenetics and nutrigenomics in elucidating the underlying mechanisms of fat metabolism and IMF deposition in pigs. This knowledge is essential in redefining nutritional intervention for swine production and the improvement of some economically important traits such as growth performance, backfat thickness, IMF accretion, disease resistance etc., in animals.

Abstract

Fat metabolism and intramuscular fat (IMF) are qualitative traits in pigs whose development are influenced by several genes and metabolic pathways. Nutrigenetics and nutrigenomics offer prospects in estimating nutrients required by a pig. Application of these emerging fields in nutritional science provides an opportunity for matching nutrients based on the genetic make-up of the pig for trait improvements. Today, integration of high throughput “omics” technologies into nutritional genomic research has revealed many quantitative trait loci (QTLs) and single nucleotide polymorphisms (SNPs) for the mutation(s) of key genes directly or indirectly involved in fat metabolism and IMF deposition in pigs. Nutrient–gene interaction and the underlying molecular mechanisms involved in fatty acid synthesis and marbling in pigs is difficult to unravel. While existing knowledge on QTLs and SNPs of genes related to fat metabolism and IMF development is yet to be harmonized, the scientific explanations behind the nature of the existing correlation between the nutrients, the genes and the environment remain unclear, being inconclusive or lacking precision. This paper aimed to: (1) discuss nutrigenetics, nutrigenomics and epigenetic mechanisms controlling fat metabolism and IMF accretion in pigs; (2) highlight the potentials of these concepts in pig nutritional programming and research.

MeDIP-sequencing for profiling global DNA methylation in buffalo embryos produced by in vitro fertilization

Assisted reproductive technique like in vitro fertilization has contributed immensely in producing genetically improved livestock. Production of embryos under in vitro conditions can affect global DNA methylation pattern during the course of embryonic development. The present study is aimed at the generation and comparison of global DNA methylome of embryos at 2-cell, 8-cell and blastocyst stage of buffalo embryos produced by in vitro fertilization using MeDIP-Sequencing. It is observed that there is a profound difference in the global DNA methylation profile of IVF embryos at different developmental stages. These differences are manifested throughout the course of embryonic development. Pathways like Wnt signaling pathway, gonadotropin-releasing hormone receptor pathway and integrin signaling were found to be majorly affected by hypermethylation of DNA in IVF embryos throughout the development.

A 104-bp Structural Variation of the ADPRHL1 Gene Is Associated With Growth Traits in Chickens

Analyzing marker-assisted breeding is an important method utilized in modern molecular breeding. Recent studies have determined that a large number of molecular markers appear to explain the impact of "lost heritability" on human height. Therefore, it is necessary to locate molecular marker sites in poultry and investigate the possible molecular mechanisms governing their effects. In this study, we found a 104-bp insertion/deletion polymorphism in the 5'UTR of the ADPRHL1 gene through resequencing. In cross-designed F₂ resource groups, the indel was significantly associated with weight at 0, 2, 4, 6, and 10 weeks and a number of other traits [carcass weight (CW), semi-evisceration weight (SEW), evisceration weight (EW), claw weight (CLW), wings weight (DWW), gizzard weight (GW), pancreas weight (PW), chest muscle weight (CMW), leg weight (LW), leg muscle weight (LMW), shedding Weight (SW), liver rate (LR), and leg muscle rate (LMR)] (P < 0.05). In brief, the insertion-insertion (II) genotype was significantly associated with the greatest growth traits and meat quality traits, whereas the values associated with the insertion-deletion (ID) genotype were the lowest in the F₂ reciprocal cross chickens. The mutation sites were genotyped in 4,526 individuals from 12 different chicken breeds and cross-designed F₂ resource groups. The II genotype is the most important genotype in commercial broilers, and the I allele frequency observed in these breeds is relatively high. Deletion mutations tend to be fixed in commercial broilers. However, there is still considerable great potential for breeding in dual-purpose chickens and commercial laying hens. A luciferase reporter assay showed that the II genotype of the ADPRHL1 gene possessed 2.49-fold higher promoter activity than the DD genotype (P < 0.05). We hypothesized that this indel might affect the transcriptional activity of ADPRHL1, thereby affecting the growth traits of chickens. These findings may help to elucidate the function of the ADPRHL1 gene and facilitate enhanced reproduction in the chicken industry.

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 DNA methylation and transcription analysis in tongue and biceps femoris muscles of cloned pigs with macroglossia

Cloned animals are prone to abnormal phenotypes such as enlarged tongue, fetal oversize, and progeria. In the present study, whole-genome bisulfite sequencing and mRNA sequencing were performed on tongue and biceps femoris muscles of cloned piglets with and without macroglossia, in an attempt to elucidate the epigenetic causes of the macroglossia phenotype. We identified 14 958 and 18 752 differentially methylated regions in the tongue and biceps femoris muscles, respectively, of macroglossia piglets and these correspond to 4574 and 4772 differentially methylated genes compared with the control group (piglets without macroglossia). Larger methylation difference was found in tongue muscle than in biceps femoris muscle. In total, 114 genes in tongue and 72 genes in biceps femoris muscles were found to be differentially expressed between the two groups. Of these differentially expressed genes in tongue muscle, 31 were also differentially methylated genes, among which DIO3 and ZIC1 were imprinting or predicted imprinting genes. These two and another six overlapping genes (ALDH1A2, MKX, MAB21L2, CA3, RANBP3L, and MYL10) are crucial factors involved in embryonic development or tissue and organ development. GO enrichment analysis suggested possible alteration of these processes. Our study provides novel molecular insights into the formation of macroglossia in cloned pigs.

Multi-Omics Analysis of Key microRNA-mRNA Metabolic Regulatory Networks in Skeletal Muscle of Obese Rabbits

microRNAs (miRNAs), small non-coding RNA with a length of about 22 nucleotides, are involved in the energy metabolism of skeletal muscle cells. However, their molecular mechanism of metabolism in rabbit skeletal muscle is still unclear. In this study, 16 rabbits, 8 in the control group (CON-G) and 8 in the experimental group (HFD-G), were chosen to construct an obese model induced by a high-fat diet fed from 35 to 70 days of age. Subsequently, 54 differentially expressed miRNAs, 248 differentially expressed mRNAs, and 108 differentially expressed proteins related to the metabolism of skeletal muscle were detected and analyzed with three sequencing techniques (small RNA sequencing, transcriptome sequencing, and tandem mass tab (TMT) protein technology). It was found that 12 miRNAs and 12 core genes (e.g., CRYL1, VDAC3 and APIP) were significantly different in skeletal muscle from rabbits in the two groups. The network analysis showed that seven miRNA-mRNA pairs were involved in metabolism. Importantly, two miRNAs (miR-92a-3p and miR-30a/c/d-5p) regulated three transcription factors (MYBL2, STAT1 and IKZF1) that may be essential for lipid metabolism. These results enhance our understanding of molecular mechanisms associated with rabbit skeletal muscle metabolism and provide a basis for future studies in the metabolic diseases of human obesity.

DNA methylation status of SERPINA3 gene involved in mRNA expression in three cattle breeds

DNA methylation could take part in the gene expression and acts an important role in muscle development. In this study, DNA methylation and expression in adipose and muscle tissues were examined at the same time to evaluate the extent of epigenetic modifications and gene expression on the differentially methylated region (DMR) in SERPINA3. Chain reaction of bisulfite sequencing polymerase (BSP) was used to compared difference among DNA methylation patterns. The result of quantitative real-time PCR (qPCR) analysis showed that there was an extensive expression of SERPINA3 gene in tissue and there was a significant difference existing in muscle and adipose between Jiaxian cattle and individual of other breeds with increasing hybridization (p < 0.05). The statistic analyses indicated that DNA methylation patterns had a significant influence to the level of mRNA in tissue of fat and muscle. This study may be an important reference for investigating development of muscle tissue in cattle, and may promote the process of cattle molecular breeding.

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.

DNA methylation studies in cattle

Investigation of the role of epigenetics in cattle breeding is gaining importance. DNA methylation represents an epigenetic modification which is essential for genomic stability and maintenance of development. Recently, DNA methylation research in cattle has intensified. The studies focus on the definition of methylomes in various organs and tissues in relation to the expression of genes underlying economically important traits, and explore methylome changes under developmental, environmental, disease, and diet influences. The investigations further characterize the methylation patterns of gametes in connection with their quality, and study methylome alterations in the developing naturally or assisted produced zygotes, embryos, and fetuses, considering their viability. A wide array of technologies developed for accurate and precise analysis of DNA methylation patterns is employed for both single-gene and genome-wide studies. Overall, the research is directed towards the identification of single methylation markers or their combinations which may be useful in the selection and breeding of animals to ensure cattle improvement.

Comprehensive analysis of lncRNA expression profiles in cytopathic biotype BVDV-infected MDBK cells provides an insight into biological contexts of host-BVDV interactions

Bovine viral diarrhea virus (BVDV) is the causative agent of bovine viral diarrhea-mucosal disease, which significantly affects the production performance of cattle, causing serious economic losses to the cattle industries worldwide. Up to now, some mechanisms involved in host–BVDV interaction are still not fully understood. The discovery of long non-coding RNAs (lncRNAs) has provided a new perspective on gene regulation in diverse biological contexts, particularly in viral infection and host immune responses. However, little is known about the profiles and functions of lncRNAs in host cells in response to BVDV infection. Here, we utilized Illumina sequencing to explore lncRNAs profiles in cytopathic (CP) biotype BVDV-infected MDBK cells to further reveal the potential roles of lncRNAs in BVDV infection and host–BVDV interaction with integrated analysis of lncRNAs and mRNA expression profiles. A total of 1747 significantly differentially expressed genes, DEGs (156 lncRNAs and 1591 mRNAs) were obtained via RNA-seq in BVDV-infected MDBK cells compared to mock-infected cells. Next, these DE lncRNAs and mRNAs were subjected to construct lncRNAs-mRNAs co-expression network followed by the prediction of potential functions of the DE lncRNAs. Co-expression network analysis elucidated that DE lncRNAs were significant enrichment in NOD-like receptor, TNF, NF-ĸB, ErbB, Ras, apoptosis, and fatty acid biosynthesis pathways, indicating that DE lncRNAs play important roles in host–BVDV interactions. Our data give an overview of changes in transcriptome and potential roles of lncRNAs, providing molecular biology basis for further exploring the mechanisms of host–BVDV interaction.

Comparative genome-wide methylation analysis of longissimus dorsi muscles in Yorkshire and Wannanhua pigs

DNA methylation was one of the earliest discovered epigenetic modifications in vertebrates, and is an important epigenetic mechanism involved in the expression of genes in many biological processes, including muscle growth and development. Its effects on economically important traits are evidenced in reported differences in meat quality traits between Chinese indigenous pig breeds (Wannanhua pig) and Western commercial pig breeds (Yorkshire pig), and this presents a unique model for analyzing the effects of DNA methylation on these traits. In the present study, a whole genome DNA methylation analysis was performed on the two breeds using methylated DNA immunoprecipitation. GO functional enrichment and pathway enrichment analyses identified differentially methylated genes primarily associated with fatty acid metabolism, biological processes of muscle development and signaling pathways related to muscle development and pork quality. Differentially methylated genes were verified by sodium pyrosequencing, and the results were consistent with the sequencing results. The results of the integrative analysis between DNA methylation and gene expression revealed that the DNA methylation levels showed a significantly negative correlation with gene expression levels around the transcription start site of genes. In total, 41 genes were both differentially expressed and methylated; these genes were related to fat metabolism, lipid metabolism and skeletal muscle development. This study could help further explore the molecular mechanisms and phenotypic differences in pig growth and development among different breeds.

The Profiling of DNA Methylation and Its Regulation on Divergent Tenderness in Angus Beef Cattle

Beef is an essential food source in the world. Beef quality, especially tenderness, has a significant impact on consumer satisfaction and industry profit. Many types of research to date have focused on the exploration of physiological and developmental mechanisms of beef tenderness. Still, the role and impact of DNA methylation status on beef tenderness have yet to be elucidated. In this study, we exhaustively analyzed the DNA methylation status in divergent tenderness observed in Angus beef. We characterized the methylation profiles related to beef tenderness and explored methylation distributions on the whole genome. As a result, differentially methylated regions (DMRs) associated with tenderness and toughness of beef were identified. Importantly, we annotated these DMRs on the bovine genome and explored bio-pathways of underlying genes and methylation biomarkers in beef quality. Specifically, we observed that the ATP binding cassette subfamily and myosin-related genes were highly methylated gene sets, and generation of neurons, regulation of GTPase activity, ion transport and anion transport, etc., were the significant pathways related with beef tenderness. Moreover, we explored the relationship between DNA methylation and gene expression in DMRs. Some methylated genes were identified as candidate biomarkers for beef tenderness. These results provide not only novel epigenetic information associated with beef quality but offer more significant insights into meat science, which will further help us explore the mechanism of muscle biology.

Comparative whole genome DNA methylation profiling across cattle tissues reveals global and tissue-specific methylation patterns

BACKGROUND

Efforts to improve animal health, and understand genetic bases for production, may benefit from a comprehensive analysis of animal genomes and epigenomes. Although DNA methylation has been well studied in humans and other model species, its distribution patterns and regulatory impacts in cattle are still largely unknown. Here, we present the largest collection of cattle DNA methylation epigenomic data to date.

RESULTS

Using Holstein cattle, we generated 29 whole genome bisulfite sequencing (WGBS) datasets for 16 tissues, 47 corresponding RNA-seq datasets, and 2 whole genome sequencing datasets. We did read mapping and DNA methylation calling based on two different cattle assemblies, demonstrating the high quality of the long-read-based assembly markedly improved DNA methylation results. We observed large differences across cattle tissues in the methylation patterns of global CpG sites, partially methylated domains (PMDs), hypomethylated regions (HMRs), CG islands (CGIs), and common repeats. We detected that each tissue had a distinct set of PMDs, which showed tissue-specific patterns. Similar to human PMD, cattle PMDs were often linked to a general decrease of gene expression and a decrease in active histone marks and related to long-range chromatin organizations, like topologically associated domains (TADs). We tested a classification of the HMRs based on their distributions relative to transcription start sites (TSSs) and detected tissue-specific TSS-HMRs and genes that showed strong tissue effects. When performing cross-species comparisons of paired genes (two opposite strand genes with their TSS located in the same HMR), we found out they were more consistently co-expressed among human, mouse, sheep, goat, yak, pig, and chicken, but showed lower consistent ratios in more divergent species. We further used these WGBS data to detect 50,023 experimentally supported CGIs across bovine tissues and found that they might function as a guard against C-to-T mutations for TSS-HMRs. Although common repeats were often heavily methylated, some young Bov-A2 repeats were hypomethylated in sperm and could affect the promoter structures by exposing potential transcription factor binding sites.

CONCLUSIONS

This study provides a comprehensive resource for bovine epigenomic research and enables new discoveries about DNA methylation and its role in complex traits.

Combined Analysis of DNA Methylome and Transcriptome Reveal Novel Candidate Genes Related to Porcine Escherichia coli F4ab/ac-Induced Diarrhea

Enterotoxigenic Escherichia coli (ETEC) that express F4 (K88) fimbriae are the principal microorganisms responsible for bacterial diarrhea in neonatal and pre-weaning piglets. To better understand the molecular effects of ETEC F4ab/ac infection, we performed a genome-wide comparison of the changes in DNA methylation and gene expression in ETEC F4ab/ac infected porcine intestinal epithelial cells. We characterized the pattern of changes in methylation and found 3297 and 1593 differentially methylated regions in cells infected with F4ab and F4ac, respectively. Moreover, 606 and 780 differentially expressed genes (DEGs) in ETEC F4ab and F4ac infected cells were detected and these genes were highly enriched in immune/defense response related pathways. Integrative analysis identified 27 and 10 genes showing inverse correlations between promoter methylation and expression with ETEC F4ab/ac infection. Altered DNA methylation and expression of various genes suggested their roles and potential functional interactions upon ETEC F4ab/ac infection. Further functional analyses revealed that three DEGs (S100A9, SGO1, and ESPL1) in F4ab infected cells and three DEGs (MAP3K21, PAK6, and MPZL1) in F4ac infected cells are likely involved in the host cells response to ETEC infection. Our data provides further insight into the epigenetic and transcriptomic alterations of ETEC F4ab/ac infected porcine intestinal epithelial cells, and may advance the identification of biomarkers and drug targets for predicting susceptibility to and controlling ETEC F4ab/ac induced diarrhea.

Analysis of DNA methylation profiles during sheep skeletal muscle development using whole-genome bisulfite sequencing

BACKGROUND

DNA methylation is an epigenetic regulatory form that plays an important role in regulating the gene expression and the tissues development.. However, DNA methylation regulators involved in sheep muscle development remain unclear. To explore the functional importance of genome-scale DNA methylation during sheep muscle growth, this study systematically investigated the genome-wide DNA methylation profiles at key stages of Hu sheep developmental (fetus and adult) using deep whole-genome bisulfite sequencing (WGBS).

RESULTS

Our study found that the expression levels of DNA methyltransferase (DNMT)-related genes were lower in fetal muscle than in the muscle of adults. The methylation levels in the CG context were higher than those in the CHG and CHH contexts, and methylation levels were highest in introns, followed by exons and downstream regions. Subsequently, we identified 48,491, 17, and 135 differentially methylated regions (DMRs) in the CG, CHG, and CHH sequence contexts and 11,522 differentially methylated genes (DMGs). The results of bisulfite sequencing PCR (BSP) correlated well with the WGBS-Seq data. Moreover, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional annotation analysis revealed that some DMGs were involved in regulating skeletal muscle development and fatty acid metabolism. By combining the WGBS-Seq and previous RNA-Seq data, a total of 159 overlap genes were obtained between differentially expressed genes (DEGs) and DMGs (FPKM > 10 and fold change > 4). Finally, we found that 9 DMGs were likely to be involved in muscle growth and metabolism of Hu sheep.

CONCLUSIONS

We systemically studied the global DNA methylation patterns of fetal and adult muscle development in Hu sheep, which provided new insights into a better understanding of the epigenetic regulation of sheep muscle development.

Expression profiles of the MXD3 gene and association of sequence variants with growth traits in Xianan and Qinchuan cattle

Max dimerization protein 3 (MXD3) belongs to the MYC superfamily of basic helix‐loop‐helix leucine zipper transcription factors, and MXD3‐MAX heterodimers can bind to promoters of target genes to modulate their expression. The aim of this study was to determine the MXD3 mRNA expression levels in various cattle tissues comprising heart, liver, spleen, lung, kidney, Longissimus dorsi muscle and subcutaneous fat in Chinese Qinchuan and Xianan cattle breeds. The RT‐qPCR data showed that the MXD3 gene was variably expressed between all tissues and at levels that were significantly different between two breeds (p < .05). We used the polymerase chain reaction‐restriction fragment length polymorphism (PCR‐RFLP) method to investigate the possible association between single‐nucleotide polymorphisms (SNP) within the MXD3 gene and five different growth traits in cattle. We found two intronic SNPs (g.2694 C>T and g.3801 T>C) and one SNP in 3′untranslated region (3′UTR) (g.6263 G>A) of MXD3 gene. Association analysis revealed strong associations between pairwise and triple SNP combinations and the growth traits. Based on these results, we suggest that MXD3 polymorphisms could be useful as molecular markers in the Chinese beef cattle breeding program.

Transcriptome and DNA Methylation Analyses of the Molecular Mechanisms Underlying with Longissimus dorsi Muscles at Different Stages of Development in the Polled Yak

DNA methylation modifications are implicated in many biological processes. As the most common epigenetic mechanism DNA methylation also affects muscle growth and development. The majority of previous studies have focused on different varieties of yak, but little is known about the epigenetic regulation mechanisms in different age groups of animals. The development of muscles in the different stages of yak growth remains unclear. In this study, we selected the longissimus dorsi muscle tissue at three different growth stages of the yak, namely, 90-day-old fetuses (group E), six months old (group M), and three years old (group A). Using RNA-Seq transcriptome sequencing and methyl-RAD whole-genome methylation sequencing technology, changes in gene expression levels and DNA methylation status throughout the genome were investigated during the stages of yak development. Each group was represented by three biological replicates. The intersections of expression patterns of 7694 differentially expressed genes (DEGs) were identified (padj < 0.01, |log2FC| > 1.2) at each of the three developmental periods. Time-series expression profile clustering analysis indicated that the DEGs were significantly arranged into eight clusters which could be divided into two classes (padj < 0.05), class I profiles that were downregulated and class II profiles that were upregulated. Based on this cluster analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that DEGs from class I profiles were significantly (padj < 0.05) enriched in 21 pathways, the most enriched pathway being the Axon guidance signaling pathway. DEGs from the class II profile were significantly enriched in 58 pathways, the pathway most strongly enriched being Metabolic pathway. After establishing the methylation profiles of the whole genomes, and using two groups of comparisons, the three combinations of groups (M-vs.-E, M-vs.-A, A-vs.-E) were found to have 1344, 822, and 420 genes, respectively, that were differentially methylated at CCGG sites and 2282, 3056, and 537 genes, respectively, at CCWGG sites. The two sets of data were integrated and the negative correlations between DEGs and differentially methylated promoters (DMPs) analyzed, which confirmed that TMEM8C, IGF2, CACNA1S and MUSTN1 were methylated in the promoter region and that expression of the modified genes was negatively correlated. Interestingly, these four genes, from what was mentioned above, perform vital roles in yak muscle growth and represent a reference for future genomic and epigenomic studies in muscle development, in addition to enabling marker-assisted selection of growth traits.

Genome-Wide DNA Methylome and Transcriptome Analysis of Porcine Intestinal Epithelial Cells upon Deoxynivalenol Exposure

Deoxynivalenol (DON) is a type of mycotoxin that is disruptive to intestinal and immune systems. To better understand the molecular effects of DON exposure, we performed genome-wide comparisons of DNA methylation and gene expression from porcine intestinal epithelial cell IPEC-J2 upon DON exposure using reduced representation bisulfite sequencing and RNA-seq technologies. We characterized the methylation pattern changes and found 3030 differentially methylated regions. Moreover, 3226 genes showing differential expression were enriched in pathways of protein and nucleic acid synthesis and ribosome biogenesis. Integrative analysis identified 29 genes showing inverse correlations between promoter methylation and expression. Altered DNA methylation and expression of various genes suggested their roles and potential functional interactions upon DON exposure. Our data provided new insights into epigenetic and transcriptomic alterations of intestinal epithelial cells upon DON exposure and may advance the identification of biomarkers and drug targets for predicting and controlling the toxic effects of this common mycotoxin.

The landscape of DNA methylation associated with the transcriptomic network in layers and broilers generates insight into embryonic muscle development in chicken

Scope: As DNA methylation is one of the key epigenetic mechanisms involved in embryonic muscle development, elucidating its relationship with non-coding RNAs and genes is essential for understanding early muscle development. The methylome profiles of pre-hatching chicken across multiple developmental stages remain incomplete although several related studies have been reported. Methods: In this study, we performed single-base-resolution bisulfite sequencing together with RNA-seq of broilers and layers in different embryonic development points (E10, E13, E16 and E19) to explore the genetic basis of embryonic muscle development in chicken. The differential methylated regions and novel lncRNAs were identified for association analyses. Through genomic position and correlation analysis between DMRs and lncRNAs, the target lncRNAs were detected to participate in the embryonic muscle formation and the results were then verified in vitro experiments. Results: Comparison of methylome profiles between two chicken lines revealed that lower methylation in broilers might contribute to muscle development in embryonic period. Differential methylated region analysis showed that the majority of differential methylated regions were hypo-DMRs for broilers. Differential methylated genes were significantly enriched in muscle development-related terms at E13 and E19. Furthermore, we identified a long non-coding RNA MyH1-AS that potentially regulated embryonic muscle development, proved by the regulatory network construction and further in vitro experiments. Conclusion: Our study revealed an integrative landscape of middle- to late-stage of embryonic myogenesis in chicken, gave rise to a comprehensive understanding of epigenetic and transcriptional regulation in muscle development. Moreover, we provided a reliable data resource for further embryonic muscle development studies.

Earlier demethylation of myogenic genes contributes to embryonic precocious terminal differentiation of myoblasts in miniature pigs

Here, we performed whole-genome bisulfite sequencing of longissimus dorsi muscle from Landrace and Wuzhishan (WZS) miniature pigs during 18, 21, and 28 d postcoitum. It was uncovered that in regulatory regions only around transcription start sites (TSSs), gene expression and methylation showed negative correlation, whereas in gene bodies, positive correlation occurred. Furthermore, earlier myogenic gene demethylation around TSSs and earlier hypermethylation of myogenic genes in gene bodies were considered to trigger their earlier expression in miniature pigs. Furthermore, by analyzing the methylation pattern of the myogenic differentiation 1(MyoD) promoter and distal enhancer, we found that earlier demethylation of the MyoD distal enhancer in WZSs contributes to its earlier expression. Moreover, DNA demethylase Tet1 was found to be involved in the demethylation of the myogenin promoter and promoted immortalized mouse myoblast cell line (C2C12) and porcine embryonic myogenic cell differentiation. This study reveals that earlier demethylation of myogenic genes contributes to precocious terminal differentiation of myoblasts in miniature pigs.-Zhang, X., Nie, Y., Cai, S., Ding, S., Fu, B., Wei, H., Chen, L., Liu, X., Liu, M., Yuan, R., Qiu, B., He, Z., Cong, P., Chen, Y., Mo, D. Earlier demethylation of myogenic genes contributes to embryonic precocious terminal differentiation of myoblasts in miniature pigs.

Survey of allele specific expression in bovine muscle

Allelic imbalance is a common phenomenon in mammals that plays an important role in gene regulation. An Allele Specific Expression (ASE) approach can be used to detect variants with a cis-regulatory effect on gene expression. In cattle, this type of study has only been done once in Holstein. In our study we performed a genome-wide analysis of ASE in 19 Limousine muscle samples. We identified 5,658 ASE SNPs (Single Nucleotide Polymorphisms showing allele specific expression) in 13% of genes with detectable expression in the Longissimus thoraci muscle. Interestingly we found allelic imbalance in AOX1, PALLD and CAST genes. We also found 2,107 ASE SNPs located within genomic regions associated with meat or carcass traits. In order to identify causative cis-regulatory variants explaining ASE we searched for SNPs altering binding sites of transcription factors or microRNAs. We identified one SNP in the 3’UTR region of PRNP that could be a causal regulatory variant modifying binding sites of several miRNAs. We showed that ASE is frequent within our muscle samples. Our data could be used to elucidate the molecular mechanisms underlying gene expression imbalance.

Role of Maternal Dietary Protein and Amino Acids on Fetal Programming, Early Neonatal Development, and Lactation in Swine

Simple Summary

Dietary protein is an important nutrient source for sows, necessary for not only growth and production, but also other physiological functions. Protein limitations in maternal diets have the potential to impair fetal myogenesis, while excess maternal dietary protein appears to only have minor effects on early fetal muscle formation. Effects of maternal protein deficiency on increased fat deposition in porcine neonates is inconsistent with gene expressions in the neonates. Sufficient maternal dietary protein can enhance porcine milk protein and fat concentration. Understanding the function of protein and amino acids in sows and the effects on their offspring can provide rational approaches for the regulation of piglet growth and further improvements in meat quality in the future.

Abstract

Maternal nutrition plays a vital role in fetal development, early development of neonates, and lactation and regulates the lifetime productivity of offspring. During pregnancy, maternal nutrition alters expression of the fetal genome and the development of tissues and organs via fetal programming. After parturition, maternal nutrition continues to regulate growth and development of piglets through maternal milk, which contains carbohydrates, lipids, proteins and oligosaccharides. Thus, deficiencies in maternal nutrition are detrimental to development of piglets, which can lead to inefficient growth and decreased carcass merit. Protein is an important nutritional component for sows, which not only functions in muscle development, but also plays a vital role in embryonic and neonatal development and lactation. Although effects of maternal undernutrition on neonatal development have been widely studied in sows, the function of different maternal dietary protein levels on fetal development, neonatal growth and lactation performance of sows is largely unknown. Determination of the effects and underlying mechanisms of maternal dietary protein levels on development of piglets is vital to the pork industry. Therefore, we summarized recent reports regarding mechanisms of effects of maternal protein levels on regulation of conceptus growth and early postnatal development though uterine fetal programming and lactation in swine.

Epigenetic regulation of gene expression in Chinese Hamster Ovary cells in response to the changing environment of a batch culture

The existence of dynamic cellular phenotypes in changing environmental conditions is of major interest for cell biologists who aim to understand the mechanism and sequence of regulation of gene expression. In the context of therapeutic protein production by Chinese Hamster Ovary (CHO) cells, a detailed temporal understanding of cell‐line behavior and control is necessary to achieve a more predictable and reliable process performance. Of particular interest are data on dynamic, temporally resolved transcriptional regulation of genes in response to altered substrate availability and culture conditions. In this study, the gene transcription dynamics throughout a 9‐day batch culture of CHO cells was examined by analyzing histone modifications and gene expression profiles in regular 12‐ and 24‐hr intervals, respectively. Three levels of regulation were observed: (a) the presence or absence of DNA methylation in the promoter region provides an ON/OFF switch; (b) a temporally resolved correlation is observed between the presence of active transcription‐ and promoter‐specific histone marks and the expression level of the respective genes; and (c) a major mechanism of gene regulation is identified by interaction of coding genes with long non‐coding RNA (lncRNA), as observed in the regulation of the expression level of both neighboring coding/lnc gene pairs and of gene pairs where the lncRNA is able to form RNA–DNA–DNA triplexes. Such triplex‐forming regions were predominantly found in the promoter or enhancer region of the targeted coding gene. Significantly, the coding genes with the highest degree of variation in expression during the batch culture are characterized by a larger number of possible triplex‐forming interactions with differentially expressed lncRNAs. This indicates a specific role of lncRNA‐triplexes in enabling rapid and large changes in transcription. A more comprehensive understanding of these regulatory mechanisms will provide an opportunity for new tools to control cellular behavior and to engineer enhanced phenotypes.

Contrasting effects of acute and chronic stress on the transcriptome, epigenome, and immune response of Atlantic salmon

Stress experienced during early life may have lasting effects on the immune system, with impacts on health and disease dependent on the nature and duration of the stressor. The epigenome is especially sensitive to environmental stimuli during early life and represents a potential mechanism through which stress may cause long-lasting health effects. However, the extent to which the epigenome responds differently to chronic vs acute stressors is unclear, especially for non-mammalian species. We examined the effects of acute stress (cold-shock during embryogenesis) and chronic stress (absence of tank enrichment during larval-stage) on global gene expression (using RNA-seq) and DNA methylation (using RRBS) in the gills of Atlantic salmon (Salmo salar) four months after hatching. Chronic stress induced pronounced transcriptional differences, while acute stress caused few lasting transcriptional effects. However, both acute and chronic stress caused lasting and contrasting changes in the methylome. Crucially, we found that acute stress enhanced transcriptional immune response to a pathogenic challenge (bacterial lipopolysaccharide, LPS), while chronic stress suppressed it. We identified stress-induced changes in promoter and gene-body methylation that were associated with altered expression for a small proportion of immune-related genes, and evidence of wider epigenetic regulation within signalling pathways involved in immune response. Our results suggest that stress can affect immuno-competence through epigenetic mechanisms, and highlight the markedly different effects of chronic larval and acute embryonic stress. This knowledge could be used to harness the stimulatory effects of acute stress on immunity, paving the way for improved stress and disease management through epigenetic conditioning.

Genome-wide differences in DNA methylation changes in caprine ovaries between oestrous and dioestrous phases

Background

DNA methylation plays a vital role in reproduction. Entire genome DNA methylation changes during the oestrous phase (ES) and dioestrous phase (DS) in the ovaries of Guanzhong dairy goats were investigated using bisulphite sequencing to understand the molecular biological mechanisms of these goats’ oestrous cycle.

Results

We discovered distinct genome-wide DNA methylation patterns in ES and DS ovaries. A total of 26,910 differentially methylated regions were upregulated and 21,453 differentially methylated regions were downregulated in the ES samples compared with the DS samples (P-values ≤0.05 and fold change of methylation ratios ≥2). Differentially methylated region analysis showed hypomethylation in the gene body regions and hypermethylation in the joining region between upstream regions and gene bodies. The methylation ratios of the STAR, FGF2, FGF12, BMP5 and SMAD6 genes in the ES samples were lower than those of the DS samples (P-values ≤0.05 and fold change of methylation ratios ≥2). Conversely, the methylation ratios of the EGFR, TGFBR2, IGF2BP1 and MMD2 genes increased in the ES samples compared with the DS samples. In addition, 223 differentially methylated genes were found in the GnRH signalling pathway (KO04912), ovarian steroidogenesis pathway (KO04913), oestrogen signalling pathway (KO04915), oxytocin signalling pathway (KO04921), insulin secretion pathway (KO04911) and MAPK signalling pathway (KO04010).

Conclusions

This study is the first large-scale comparison of the high-resolution DNA methylation landscapes of oestrous and dioestrous ovaries from dairy goats. Previous studies and our investigations have shown that the NR5A2, STAR, FGF2 and BMP5 genes might have potential application value in regulating caprine oestrus.

Electronic supplementary material

The online version of this article (10.1186/s40104-018-0301-x) contains supplementary material, which is available to authorized users.

Comparative whole genome DNA methylation profiling of cattle sperm and somatic tissues reveals striking hypomethylated patterns in sperm

Background

Although sperm DNA methylation has been studied in humans and other species, its status in cattle is largely unknown.

Results

Using whole-genome bisulfite sequencing (WGBS), we profiled the DNA methylome of cattle sperm through comparison with three somatic tissues (mammary gland, brain, and blood). Large differences between cattle sperm and somatic cells were observed in the methylation patterns of global CpGs, pericentromeric satellites, partially methylated domains (PMDs), hypomethylated regions (HMRs), and common repeats. As expected, we observed low methylation in the promoter regions and high methylation in the bodies of active genes. We detected selective hypomethylation of megabase domains of centromeric satellite clusters, which may be related to chromosome segregation during meiosis and their rapid transcriptional activation upon fertilization. We found more PMDs in sperm cells than in somatic cells and identified meiosis-related genes such asKIF2B and REPIN1, which are hypomethylated in sperm but hypermethylated in somatic cells. In addition to the common HMRs around gene promoters, which showed substantial differences between sperm and somatic cells, the sperm-specific HMRs also targeted to distinct spermatogenesis-related genes, including BOLL, MAEL, ASZ1, SYCP3, CTCFL, MND1, SPATA22, PLD6, DDX4, RBBP8, FKBP6, and SYCE1. Although common repeats were heavily methylated in both sperm and somatic cells, some young Bov-A2 repeats, which belong to the SINE family, were hypomethylated in sperm and could affect the promoter structures by introducing new regulatory elements.

Conclusions

Our study provides a comprehensive resource for bovine sperm epigenomic research and enables new discoveries about DNA methylation and its role in male fertility.

MTCH2 promotes adipogenesis in intramuscular preadipocytes via an m6A-YTHDF1-dependent mechanism

Intramuscular fat is considered a potential factor that is associated with meat quality in animal production and insulin resistance in humans. N⁶-methyladenosine (m⁶A) modification of mRNA plays an important role in regulating adipogenesis. However, the effects of m⁶A on the adipogenesis of intramuscular preadipocytes and associated mechanisms remain unknown. Here, we performed m⁶A sequencing to compare m⁶A methylome of the longissimus dorsi muscles (LDMs) between Landrace pigs (lean-type breed) and Jinhua pigs (obese-type breed with higher levels of intramuscular fat). Transcriptome-wide m⁶A profiling of porcine LDMs was highly conserved with humans and mice. Furthermore, we identified a unique methylated gene in Jinhua pigs named mitochondrial carrier homology 2 ( MTCH2). The m⁶A levels of MTCH2 mRNA were reduced by introducing a synonymous mutation, and adipogenesis test results showed that the MTCH2 mutant was inferior with regard to adipogenesis compared with the MTCH2 wild-type. We then found that MTCH2 protein expression was positively associated with m⁶A levels, and an YTH domain family protein 1-RNA immunoprecipitation-quantitative PCR assay indicated that MTCH2 mRNA was a target of the YTH domain family protein 1. This study provides comprehensive m⁶A profiles of LDM transcriptomes in pigs and suggests an essential role for m⁶A modification of MTCH2 in intramuscular fat regulation.-Jiang, Q., Sun, B., Liu, Q., Cai, M., Wu, R., Wang, F., Yao, Y., Wang, Y., Wang, X. MTCH2 promotes adipogenesis in intramuscular preadipocytes via an m⁶A-YTHDF1-dependent mechanism.

DNA methylation patterns in peripheral blood mononuclear cells from Holstein cattle with variable milk yield

Background

Milk yield for Holstein cows has doubled over five decades due to genetic selection and changes to management, but the molecular mechanisms that facilitated this increase are mostly unknown. Epigenetic modifications to the cattle genome are a plausible molecular mechanism to cause variation in milk yield and our objective was to describe genome-wide DNA methylation patterns in peripheral blood mononuclear cells (PBMC) from mature Holstein dairy cows with variable milk yield.

Results

Whole genome MeDIP-seq was performed following DNA extraction from PBMC of 6 lactating dairy cows from 4 different herds that varied in milk yield from 13,556 kg to 23,105 kg per 305 day lactation. We describe methylation across the genome and for 13,677 protein coding genes. Repetitive element reads were primarily mapped to satellite (36.4%), SINE (29.1%), and LINE (23.7%) regions and the majority (78.4%) of CpG sites were sequenced at least once. DNA methylation was generally low upstream of genes with the nadir occurring 95 bp prior to the transcription start site (TSS). Methylation was lower in the first exon than in later exons, was highest for introns near the intron-exon junctions, and declined downstream as the distance from the gene increased. We identified 72 differentially methylated regions (DMR) between high milk yield cows and their control, and 252 DMR across herd environments.

Conclusions

This reference methylome for cattle with extreme variation in milk yield phenotype provides a resource to more fully evaluate relationships between DNA methylation and phenotype in populations subject to selection. The detection of DMR in cows of varying milk yield suggests potential to exploit epigenetic variation in cattle improvement programs.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5124-9) contains supplementary material, which is available to authorized users.

Transposon-associated epigenetic silencing during Pleurotus ostreatus life cycle

Transposable elements constitute an important fraction of eukaryotic genomes. Given their mutagenic potential, host-genomes have evolved epigenetic defense mechanisms to limit their expansion. In fungi, epigenetic modifications have been widely studied in ascomycetes, although we lack a global picture of the epigenetic landscape in basidiomycetes. In this study, we analysed the genome-wide epigenetic and transcriptional patterns of the white-rot basidiomycete Pleurotus ostreatus throughout its life cycle. Our results performed by using high-throughput sequencing analyses revealed that strain-specific DNA methylation profiles are primarily involved in the repression of transposon activity and suggest that 21 nt small RNAs play a key role in transposon silencing. Furthermore, we provide evidence that transposon-associated DNA methylation, but not sRNA production, is directly involved in the silencing of genes surrounded by transposons. Remarkably, we found that nucleus-specific methylation levels varied in dikaryotic strains sharing identical genetic complement but different subculture conditions. Finally, we identified key genes activated in the fruiting process through the comparative analysis of transcriptomes. This study provides an integrated picture of epigenetic defense mechanisms leading to the transcriptional silencing of transposons and surrounding genes in basidiomycetes. Moreover, our findings suggest that transcriptional but not methylation reprogramming triggers fruitbody development in P. ostreatus.