Comparative genome-wide methylation analysis of longissimus dorsi muscles between Japanese black (Wagyu) and Chinese Red Steppes cattle

Comprehensive Analysis of Methylome and Transcriptome to Identify Potential Genes Regulating Porcine Testis Development

DNA methylation plays a critical role in regulating gene expression during testicular development. However, few studies report on candidate genes related to the DNA methylation regulation of porcine testicular development. This study examined the differentially expressed genes (DEGs) and their methylation levels in testicular tissues from pigs at 60 days of age (60 d) and 180 days of age (180 d) using RNA-Seq and whole genome bisulfite sequencing (WGBS). It was determined that DNA methylation primarily occurs in the cytosine-guanine (CG) context, and the analysis identified 106,282 differentially methylated regions (DMRs) corresponding to 12,385 differentially methylated genes (DMGs). Further integrated analysis of RNA-Seq and WGBS data revealed 1083 DMGs negatively correlated with the expression of DEGs. GO analysis showed that these genes were significantly enriched in spermatogenesis, germ cell development, and spermatid differentiation. The screening of enriched genes revealed that hyper-methylation repressed ADAM30, ADAM3A, DPY19L2, H2BC1, MAK, RPL10L, SPATA16, and YBX2, while hypo-methylation elevated CACNA1I, CADM1, CTNNB1, JAM2, and PAFAH1B3 expression. Additionally, the methylation status of the key genes ADAM3A, ADAM30, YBX2, JAM2, PAFAH1B3, and CTNNB1 was detected by bisulfite sequencing PCR (BSP). This study offers insights into the epigenetic regulation mechanisms underlying porcine testicular development.

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.

DNA methylation dynamics during yak adipocyte differentiation

In mammals, epigenetic modifications involving DNA methylation are necessary for the completion of the cell differentiation process. However, the global DNA methylation landscape and its dynamics during yak adipocyte differentiation remain unexplored. Here, we performed whole-genome bisulfite sequencing (WGBS) to asses DNA methylation in yak preadipocytes and adipocytes, combining these results with those of our previous studies on changes in chromatin accessibility and gene expression during yak adipogenesis. The results showed that CG methylation levels were lower in promoter than in exon and intron, and gradually decreasing from the distal regions to transcription start site (TSS). There was a significant negative correlation between CG methylation levels located in promoter and gene expression levels. The 46 genes shared by CG differentially methylated regions (DMRs) and differential chromatin accessibility were significantly enriched in Hedgehog and PI3K-Akt signaling pathways. ATAC-seq peaks with high chromatin accessibility located in both promoter (≤ 2 kb from TSS) and distal (> 2 kb from TSS) regions corresponded to low methylation levels. Taken together, these findings demonstrated that DNA methylation and its interactions with chromatin accessibility regulate gene expression during yak adipocyte differentiation, contributing to the understanding of mechanisms of various epigenetic modifications and their interactions in adipogenesis.

Molecular and Cellular Mechanisms of Intramuscular Fat Development and Growth in Cattle

Intramuscular fat, also referred to as marbling fat, is the white fat deposited within skeletal muscle tissue. The content of intramuscular fat in the skeletal muscle, particularly the longissimus dorsi muscle, of cattle is a critical determinant of beef quality and value. In this review, we summarize the process of intramuscular fat development and growth, the factors that affect this process, and the molecular and epigenetic mechanisms that mediate this process in cattle. Compared to other species, cattle have a remarkable ability to accumulate intramuscular fat, partly attributed to the abundance of sources of fatty acids for synthesizing triglycerides. Compared to other adipose depots such as subcutaneous fat, intramuscular fat develops later and grows more slowly. The commitment and differentiation of adipose precursor cells into adipocytes as well as the maturation of adipocytes are crucial steps in intramuscular fat development and growth in cattle. Each of these steps is controlled by various factors, underscoring the complexity of the regulatory network governing adipogenesis in the skeletal muscle. These factors include genetics, epigenetics, nutrition (including maternal nutrition), rumen microbiome, vitamins, hormones, weaning age, slaughter age, slaughter weight, and stress. Many of these factors seem to affect intramuscular fat deposition through the transcriptional or epigenetic regulation of genes directly involved in the development and growth of intramuscular fat. A better understanding of the molecular and cellular mechanisms by which intramuscular fat develops and grows in cattle will help us develop more effective strategies to optimize intramuscular fat deposition in cattle, thereby maximizing the quality and value of beef meat.

Genome-wide epigenetic dynamics during postnatal skeletal muscle growth in Hu sheep

Hypertrophy and fiber transformation are two prominent features of postnatal skeletal muscle development. However, the role of epigenetic modifications is less understood. ATAC-seq, whole genome bisulfite sequencing, and RNA-seq were applied to investigate the epigenetic dynamics of muscle in Hu sheep at 3 days, 3 months, 6 months, and 12 months after birth. All 6865 differentially expressed genes were assigned into three distinct tendencies, highlighting the balanced protein synthesis, accumulated immune activities, and restrained cell division in postnatal development. We identified 3742 differentially accessible regions and 11799 differentially methylated regions that were associated with muscle-development-related pathways in certain stages, like D3-M6. Transcription factor network analysis, based on genomic loci with high chromatin accessibility and low methylation, showed that ARID5B, MYOG, and ENO1 were associated with muscle hypertrophy, while NR1D1, FADS1, ZFP36L2, and SLC25A1 were associated with muscle fiber transformation. Taken together, these results suggest that DNA methylation and chromatin accessibility contributed toward regulating the growth and fiber transformation of postnatal skeletal muscle in Hu sheep.

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.

Assessing Genomic Diversity and Signatures of Selection in Chinese Red Steppe Cattle Using High-Density SNP Array

Chinese Red Steppe Cattle (CRS), a composite cattle breed, is well known for its milk production, high slaughter rate, carcass traits, and meat quality. Nowadays, it is widely bred in Jilin and Hebei Province and the Inner Mongolia Autonomous region. However, the population structure and the genetic basis of prominent characteristics of CRS are still unknown. In this study, we systematically describe their population structure, genetic diversity, and selection signature based on genotyping data from 61 CRS individuals with GGP Bovine 100 K chip. The results showed that CRS cattle had low inbreeding levels and had formed a unique genetic structure feature. Using two complementary methods (including comprehensive haplotype score and complex likelihood ratio), we identified 1291 and 1285 potentially selected genes, respectively. There were 141 genes annotated in common 106 overlapping genomic regions covered 5.62 Mb, including PLAG1, PRKG2, DGAT1, PARP10, TONSL, ADCK5, and BMP3, most of which were enriched in pathways related to muscle growth and differentiation, milk production, and lipid metabolism. This study will contribute to understanding the genetic mechanism behind artificial selection and give an extensive reference for subsequent breeding.

Profiling the immune epigenome across global cattle breeds

BACKGROUND

Understanding the variation between well and poorly adapted cattle breeds to local environments and pathogens is essential for breeding cattle with improved climate and disease-resistant phenotypes. Although considerable progress has been made towards identifying genetic differences between breeds, variation at the epigenetic and chromatin levels remains poorly characterized. Here, we generate, sequence and analyse over 150 libraries at base-pair resolution to explore the dynamics of DNA methylation and chromatin accessibility of the bovine immune system across three distinct cattle lineages.

RESULTS

We find extensive epigenetic divergence between the taurine and indicine cattle breeds across immune cell types, which is linked to the levels of local DNA sequence divergence between the two cattle sub-species. The unique cell type profiles enable the deconvolution of complex cellular mixtures using digital cytometry approaches. Finally, we show distinct sub-categories of CpG islands based on their chromatin and methylation profiles that discriminate between classes of distal and gene proximal islands linked to discrete transcriptional states.

CONCLUSIONS

Our study provides a comprehensive resource of DNA methylation, chromatin accessibility and RNA expression profiles of three diverse cattle populations. The findings have important implications, from understanding how genetic editing across breeds, and consequently regulatory backgrounds, may have distinct impacts to designing effective cattle epigenome-wide association studies in non-European breeds.

Comparative genome-wide methylation analysis reveals epigenetic regulation of muscle development in grass carp (Ctenopharyngodon idellus) fed with whole faba bean

Crisp grass carp (CGC), the most representative improved varieties of grass carp (GC), features higher muscle hardness after feeding faba bean (Vicia faba L.) for 90-120 days. DNA methylation, a most widely studied epigenetic modification, plays an essential role in muscle development. Previous studies have identified numerous differentially expressed genes (DEGs) between CGC and GC. However, it remains unknown if the expression levels of these DEGs are influenced by DNA methylation. In the present study, we performed a comprehensive analysis of DNA methylation profiles between CGC and GC, and identified important candidate genes related to muscle development coupled with the transcriptome sequencing data. A total of 9,318 differentially methylated genes (DMGs) corresponding to 155,760 differentially methylated regions (DMRs) were identified between the two groups under the CG context in promoter regions. Combined with the transcriptome sequencing data, 14 key genes related to muscle development were identified, eight of which (gsk3b, wnt8a, wnt11, axin2, stat1, stat2, jak2, hsp90) were involved in muscle fiber hyperplasia, six of which (tgf-β1, col1a1, col1a2, col1a3, col4a1, col18a1) were associated with collagen synthesis in crisp grass carp. The difference of methylation levels in the key genes might lead to the expression difference, further resulting in the increase of muscle hardness in crisp grass carp. Overall, this study can help further understand how faba bean modulates muscle development by the epigenetic modifications, providing novel insights into the texture quality improvement in other aquaculture fish species by nutritional programming.

Genome-Wide DNA Methylation Patterns of Muscle and Tail-Fat in DairyMeade Sheep and Mongolian Sheep

Simple Summary

DNA methylation is an important epigenetic modification and plays an important role in the regulation of gene expression. The study of DNA methylation will help to explore the effects of epigenetic modifications, other than DNA sequence variation, on biological phenotypes and physiological functions, promoting the revolution of livestock selection and breeding practices. DairyMeade sheep (thin-tailed) and Mongolian sheep (fat-tailed) have large differences in their dairy and meat production performances, as well as their tail phenotype, thus providing us with good materials for genomic DNA methylation studies. The below results provided a genome-wide DNA methylation landscape of muscle and tail-fat tissues between DairyMeade sheep and Mongolian sheep and a series of differentiated methylation regions (DMRs) in which CAMK2D may play a crucial role in fat metabolism and meat quality traits. These results may help us to understand how DNA methylation affects the economic traits of domestic animals like sheep.

Abstract

This study aimed to explore the genome-wide DNA methylation differences between muscle and tail-fat tissues of DairyMeade sheep (thin-tailed, lean carcass) and Mongolian sheep (fat-tailed, fat-deposited carcass). Whole-genome bisulfite sequencing (WGBS) was conducted and the global DNA methylation dynamics were mapped. Generally, CGs had a higher DNA methylation level than CHHs and CHGs, and tail-fat tissues had higher CG methylation levels than muscle tissues. For DNA repeat elements, SINE had the highest methylation level, while Simple had the lowest. When dividing the gene promoter region into small bins (200 bp per bin), the bins near the transcription start site (±200 bp) had the highest CG count per bin but the lowest DNA methylation levels. A series of DMRs were identified in muscle and tail-fat tissues between the two breeds. Among them, the introns of gene CAMK2D (calcium/calmodulin-dependent protein kinase II δ) demonstrated significant DNA methylation level differences between the two breeds in both muscle and tail-fat tissues, and it may play a crucial role in fat metabolism and meat quality traits. This study may provide basic datasets and references for further epigenetic modification studies during sheep genetic improvement.

A comparative methylome analysis reveals conservation and divergence of DNA methylation patterns and functions in vertebrates

Background

Cytosine DNA methylation is a heritable epigenetic mark present in most eukaryotic groups. While the patterns and functions of DNA methylation have been extensively studied in mouse and human, their conservation in other vertebrates remains poorly explored. In this study, we interrogated the distribution and function of DNA methylation in primary fibroblasts of seven vertebrate species including bio-medical models and livestock species (human, mouse, rabbit, dog, cow, pig, and chicken).

Results

Our data highlight both divergence and conservation of DNA methylation patterns and functions. We show that the chicken genome is hypomethylated compared to other vertebrates. Furthermore, compared to mouse, other species show a higher frequency of methylation of CpG-rich DNA. We reveal the conservation of large unmethylated valleys and patterns of DNA methylation associated with X-chromosome inactivation through vertebrate evolution and make predictions of conserved sets of imprinted genes across mammals. Finally, using chemical inhibition of DNA methylation, we show that the silencing of germline genes and endogenous retroviruses (ERVs) are conserved functions of DNA methylation in vertebrates.

Conclusions

Our data highlight conserved properties of DNA methylation in vertebrate genomes but at the same time point to differences between mouse and other vertebrate species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01270-x.

Genome-Wide DNA Methylation and Its Effect on Gene Expression During Subclinical Mastitis in Water Buffalo

Subclinical mastitis (SCM) in buffalo is one of the most challenging paradoxes for the dairy sector with very significant milk production losses and poses an imminent danger to milch animal’s milk-producing ability. We present here the genome-wide methylation specific to SCM in water buffalo and its consequential effect on the gene expression landscape for the first time. Whole-genome DNA methylation profiles from peripheral blood lymphocytes and gene expression profiles from milk somatic cells of healthy and SCM cases were catalogued from the MeDIP-Seq and RNA-Seq data. The average methylation in healthy buffaloes was found to be higher than that in the SCM-infected buffaloes. DNA methylation was abundant in the intergenic region followed by the intronic region in both healthy control and SCM groups. A total of 3,950 differentially methylated regions (DMRs) were identified and annotated to 370 differentially methylated genes (DMGs), most of which were enriched in the promoter region. Several important pathways were activated due to hypomethylation and belonged to the Staphylococcus aureus infection, Th17 cell differentiation, and antigen processing and presentation pathways along with others of defense responses. DNA methylome was compared with transcriptome to understand the regulatory role of DNA methylation on gene expression specific to SCM in buffaloes. A total of 4,778 significant differentially expressed genes (DEGs) were extracted in response to SCM, out of which 67 DMGs were also found to be differentially expressed, suggesting that during SCM, DNA methylation could be one of the epigenetic regulatory mechanisms of gene expression. Genes like CSF2RB, LOC102408349, C3 and PZP like, and CPAMD8 were found to be downregulated in our study, which are known to be involved in the immune response to SCM. Association of DNA methylation with transposable elements, miRNAs, and lncRNAs was also studied. The present study reports a buffalo SCM web resource (BSCM2TDb) available at http://webtom.cabgrid.res.in/BSCM2TDb that catalogues all the mastitis-related information of the analyses results of this study in a single place. This will be of immense use to buffalo researchers to understand the host–pathogen interaction involving SCM, which is required in endeavors of mastitis control and management.

Integrative Genome-Wide DNA Methylome and Transcriptome Analysis of Ovaries from Hu Sheep with High and Low Prolific

DNA methylation plays an important role in biological processes by affecting gene expression. However, how DNA methylation regulates phenotypic variation in Hu sheep remains unclear. Therefore, we generated genome-wide DNA methylation and transcriptomic profiles in the ovaries of Hu sheep with different prolificacies and genotypes (FecBB and FecB+). Results showed that ovary DNA methylome and transcriptome were significantly different between high prolificacy and low prolificacy Hu sheep. Comparative methylome analyses identified 10,644, 9,594, and 12,214 differentially methylated regions and 87, 1,121, and 2,375 genes, respectively, showing differential expression levels in three different comparison groups. Female reproduction-associated differentially methylated regions-related genes and differentially expressed genes were enriched, thereby the respective interaction networks were constructed. Furthermore, systematical integrative analyses revealed a negative correlation between DNA methylation around the transcriptional start site and gene expression levels, which was confirmed by testing the expression of integrin β2 subunit (ITGB2) and lysosome-associated protein transmembrane-4 beta (LAPTM4B) in vivo and in vitro. These findings demonstrated that DNA methylation influences the propensity for prolificacy by affecting gene expression in the ovaries, which may contribute to a greater understanding of the epigenome and transcriptome that will be useful for animal breeding.

Genetic Variations and Differential DNA Methylation to Face Contrasted Climates in Small Ruminants: An Analysis on Traditionally-Managed Sheep and Goats

The way in which living organisms mobilize a combination of long-term adaptive mechanisms and short-term phenotypic plasticity to face environmental variations is still largely unknown. In the context of climate change, understanding the genetic and epigenetic bases for adaptation and plasticity is a major stake for preserving genomic resources and the resilience capacity of livestock populations. We characterized both epigenetic and genetic variations by contrasting 22 sheep and 21 goats from both sides of a climate gradient, focusing on free-ranging populations from Morocco. We produced for each individual Whole-Genome Sequence at 12X coverage and MeDIP-Seq data, to identify regions under selection and those differentially methylated. For both species, the analysis of genetic differences (F_ST) along the genome between animals from localities with high vs. low temperature annual variations detected candidate genes under selection in relation to environmental perception (5 genes), immunity (4 genes), reproduction (8 genes) and production (11 genes). Moreover, we found for each species one differentially methylated gene, namely AGPTA4 in goat and SLIT3 in sheep, which were both related, among other functions, to milk production and muscle development. In both sheep and goats, the comparison between genomic regions impacted by genetic and epigenetic variations suggests that climatic variations impacted similar biological pathways but different genes.

Differential DNA methylation analysis reveals key genes in Chinese Qingyu and Landrace pigs

The Chinese Qingyu pig is a typical domestic fatty pig breed and an invaluable indigenous genetic resource in China. Compared with the Landrace pig, the Qingyu pig has unique meat characteristics, including muscle development, intramuscular fat, and other meat quality traits. At present, few studies have explored epigenetic differences due to DNA methylation between the Qingyu pig and the Landrace pig. In this study, 30 Qingyu pigs and 31 Landrace pigs were subjected to reduced representation bisulfite sequencing (RRBS). Genome-wide differential DNA methylation analysis was conducted. Six genomic regions, including regions on Sus scrofa chromosome (SSC) 1: 266.09-274.23 Mb, SSC5: 0.88-10.68 Mb, SSC8: 41.23-48.51 Mb, SSC12: 45.43-54.38 Mb, SSC13: 202.15-207.95 Mb, and SSC14: 126.43-139.85 Mb, were regarded as key regions that may be associated with phenotypic differences between the Qingyu pig and the Landrace pig. Furthermore, according to further analysis, five differentially methylated genes (ADCY1, FUBP3, GRIN2B, KIT, and PIK3R6) were identified as key candidate genes that might be associated with meat characteristics. Our findings provide new insights into the differences in DNA methylation between the Qingyu pig and the Landrace pig. These results enrich the epigenetic research of the Chinese Qingyu pig.

Analysis of DNA Methylation Profiles in Mandibular Condyle of Chicks With Crossed Beaks Using Whole-Genome Bisulfite Sequencing

Crossed beaks have been observed in at least 12 chicken strains around the world, which severely impairs their growth and welfare. To explore the intrinsic factor causing crossed beaks, this study measured the length of bilateral mandibular ramus of affected birds, and investigated the genome-wide DNA methylation profiles of normal and affected sides of mandibular condyle. Results showed that the trait was caused by impaired development of unilateral mandibular ramus, which is extended through calcification of mandibular condyle. The methylation levels in the CG contexts were higher than that of CHG and CHH, with the highest methylation level of gene body region, followed by transcription termination sites and downstream. Subsequently, we identified 1,568 differentially methylated regions and 1,317 differentially methylated genes in CG contexts. Functional annotation analysis of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes showed that these genes were involved in bone mineralization and bone morphogenesis. Furthermore, by combining the WGBS and previous RNA-Seq data, 11 overlapped genes were regulated by both long non-coding RNA and DNA methylation. Among them, FIGNL1 is an important gene in calcification of mandibular condyle. Generally, because the affected genes play key roles in maintaining mandibular calcification, these changes may be pivotal factors of crossed beaks.

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.

Genome-wide DNA methylation analysis in Chinese Chenghua and Yorkshire pigs

Background

The Chinese Chenghua pig (CHP) is a typical Chinese domestic fatty pig breed with superior meat quality characteristics, while the Yorkshire pig (YP) has the characteristics of fast growth and a high rate of lean meat. Long term natural selection and artificial selection resulted in great phenotypic differences between the two breeds, including growth, development, production performance, meat quality, and coat color. However, genome-wide DNA methylation differences between CHP and YP remain unclear.

Results

DNA methylation data were generated for muscle tissues of CHP and YP using reduced representation bisulfite sequencing (RRBS). In this study, a total of 2,416,211 CpG sites were identified. Besides, the genome-wide DNA methylation analysis revealed 722 differentially methylated regions (DMRs) and 466 differentially methylated genes (DMGs) in pairwise CHP vs. YP comparison. Six key genomic regions (Sus scrofa chromosome (SSC)1:253.47–274.23 Mb, SSC6:148.71–169.49 Mb, SSC7:0.25–9.86 Mb, SSC12:43.06–61.49 Mb, SSC14:126.43–140.95 Mb, and SSC18:49.17–54.54 Mb) containing multiple DMRs were identified, and differences of methylation patterns in these regions may be related to phenotypic differences between CHP and YP. Based on the functional analysis of DMGs, 8 DMGs (ADCY1, AGBL4, EXOC2, FUBP3, PAPPA2, PIK3R1, MGMT and MYH8) were considered as important candidate genes associated with muscle development and meat quality traits in pigs.

Conclusions

This study explored the difference in meat quality between CHP and YP from the epigenetic point of view, which has important reference significance for the local pork industry and pork food processing.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-00977-0.

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.

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.

Putative Epigenetic Biomarkers of Stress in Red Blood Cells of Chickens Reared Across Different Biomes

Production animals are constantly subjected to early adverse environmental conditions that influence the adult phenotype and produce epigenetic effects. CpG dinucleotide methylation in red blood cells (RBC) could be a useful epigenetic biomarker to identify animals subjected to chronic stress in the production environment. Here we compared a reduced fraction of the RBC methylome of chickens exposed to social isolation to non-exposed. These experiments were performed in two different locations: Brazil and Sweden. The aim was to identify stress-associated DNA methylation profiles in RBC across these populations, in spite of the variable conditions to which birds are exposed in each facility and their different lineages. Birds were increasingly exposed to a social isolation treatment, combined with food and water deprivation, at random periods of the day from weeks 1-4 after hatching. We then collected the RBC DNA from individuals and compared a reduced fraction of their methylome between the experimental groups using two bioinformatic approaches to identify differentially methylated regions (DMRs): one using fixed-size windows and another that preselected differential peaks with MACS2. Three levels of significance were used (P ≤ 0.05, P ≤ 0.005, and P ≤ 0.0005) to identify DMRs between experimental groups, which were then used for different analyses. With both of the approaches more DMRs reached the defined significance thresholds in BR individuals compared to SW. However, more DMRs had higher fold change values in SW compared to BR individuals. Interestingly, ChrZ was enriched above expectancy for the presence of DMRs. Additionally, when analyzing the locations of these DMRs in relation to the transcription starting site (TSS), we found three peaks with high DMR presence: 10 kb upstream, the TSS itself, and 20-40 kb downstream. Interestingly, these peaks had DMRs with a high presence (>50%) of specific transcription factor binding sites. Three overlapping DMRs were found between the BR and SW population using the most relaxed p-value (P ≤ 0.05). With the most stringent p-value (P ≤ 0.0005), we found 7 and 4 DMRs between treatments in the BR and SW populations, respectively. This study is the first approximation to identify epigenetic biomarkers of long-term exposure to stress in different lineages of production animals.

Genome-Wide DNA Methylation Analysis of Mammary Gland Tissues From Chinese Holstein Cows With Staphylococcus aureus Induced Mastitis

Staphylococcus aureus intramammary infection is one of the most common causes of chronic mastitis in dairy cows, whose development may be associated with epigenetic changes in the expression of important host defense genes. This study aimed to construct a genome-wide DNA methylation profile of the mammary gland of Chinese Holstein cows (n = 3) following experimentally induced S. aureus mastitis, and to explore the potential gene regulatory mechanisms affected by DNA methylation during S. aureus mastitis. DNA was extracted from S. aureus-positive (n = 3) and S. aureus-negative (n = 3) mammary gland quarters and subjected to methylation-dependent restriction-site associated DNA sequencing (Methyl-RAD Seq). Results showed that CmCGG/CmCWGG DNA methylation sites were unevenly distributed and concentrated on chromosomes 5, 11, and 19, and within intergenic regions and intron regions of genes. Compared with healthy control quarters, 9,181 significantly differentially methylated (DM) CmCGG sites and 1,790 DM CmCWGG sites were found in the S. aureus-positive quarters (P < 0.05, |log2FC| > 1). Furthermore, 363 CmCGG differently methylated genes (DMGs) and 301 CmCWGG DMGs (adjusted P < 0.05, |log2FC| > 1) were identified. Gene ontology and KEGG enrichment analysis indicated that CmCGG DMGs are involved in immune response pathways, while the CmCWGG DMGs were mainly enriched in gene ontology terms related to metabolism. The mRNAs of 526 differentially methylated CmCGG genes and 124 differentially methylated CmCWGG genes were also significantly differentially expressed (RNA-Seq data) in the same samples, herein denoted differentially methylated and expressed genes (DMEGs) (P < 0.05). Functional enrichment analysis of DMEGs revealed roles related to biological processes, especially the regulation of immune response to diseases. CmCGG DMEGs like IL6R, TNF, BTK, IL1R2, and TNFSF8 enriched in several immune-related GO terms and pathways indicated their important roles in host immune response and their potential as candidate genes for S. aureus mastitis. These results suggest potential regulatory roles for DNA methylation in bovine mammary gland processes during S. aureus mastitis and serves as a reference for future epigenetic regulation and mechanistic studies.

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.

Biomarkers of health and welfare: A One Health perspective from the laboratory side

A biomarker is any measurement reflecting an interaction between a biological system and a potential hazard, which may be chemical, physical, or biological. The One World, One Health concept established that human and animal health and the environmental state are highly interconnected, sharing common aspects that can be applied globally in these three components. In this paper, we review how the concept of One Health can be applied to biomarkers of health and welfare, with a special focus on five points that can be applied to any biomarker when it is expected to be used to evaluate the human, animal or environmental health. Three of these points are: (1) the different biomarkers that can be used, (2) the different sample types where the biomarkers can be analysed, and (3) the main methods that can be used for their measurement. In addition, we will evaluate two key points needed for adequate use of a biomarker in any situation: (4) a proper analytical validation in the sample that it is going to be used, and (5) a correct selection of the biomarker. It is expected that this knowledge will help to have a broader idea about the use of biomarkers of health and welfare and also will contribute to a better and more accurate use of these biomarkers having in mind their One Health perspective.

Comparative Analysis of the Liver and Spleen Transcriptomes between Holstein and Yunnan Humped Cattle

Simple Summary

Cattle are important agricultural animals that provide essential sources of meat, milk, fertilizer for crops, clothing, and animal traction for human use, and the demand for these products has increased in recent years. There are existing differences in disease resistance between different cattle breeds. However, the genetic basis underlying disease resistance differences is poorly understood and requires further investigation. In this study, many immune- and disease-relevant genes and pathways were identified between Holstein and Yunnan humped cattle using RNA-sequencing. The novel findings regarding the genetic basis underlying disease resistance differences between zebu cattle and taurine cattle will provide a scientific basis and key technical support for disease-resistant breeding of domestic cattle, and thus have important social and economic significance.

Abstract

Previous studies have shown that Yunnan humped cattle have higher disease resistance than pure taurine cattle, such as Holsteins. However, there exists limited information about the molecular genetic basis underlying disease resistance differences between them. The objective of this study was to compare differentially expressed genes (DEGs) in the liver and spleen tissues of Holstein and Yunnan humped cattle through comparative transcriptome analysis, using RNA-sequencing. In total, 1564 (647 up- and 917 down-regulated genes) and 1530 (716 up- and 814 down-regulated genes) DEGs were obtained in the liver and spleen tissues of Holstein and Yunnan humped cattle comparison groups, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the DEGs were mainly associated with the RIG-I signaling pathway, immune responses, major histocompatibility complex (MHC) class I protein complex and complement activation, human T-cell lymphotropic virus type-I (HTLV-I) infection. Some genes related to immune function, such as C1QB, CD55, MASP2, C4BPA, MAVS, NOD2, and CD46, were up-regulated in Yunnan humped cattle, while C2, SERPING1, SERPINE1, TIRAP, TLR2, and TLR6 were down-regulated. The expression levels of 11 selected DEGs, analyzed by quantitative reverse-transcription polymerase chain reaction (RT-qPCR), were consistent with the deep sequencing results by RNA-sequencing. Our results will provide a scientific basis and key technical support for disease-resistant breeding of domestic cattle.

Detection of InDel variations within seven candidate genes and their associations with phenotypic traits in three cattle breeds

The Xinjiang brown cattle, Red steppe cattle, and Yunling cattle are indigenous cultivated cattle breeds in Chinese frontier provinces, and they produce high-grade beef and milk products, however, their genetic diversity in many important genes related to excellent meat and milk production is still unknown. Our previous studies have found that several candidate genes (e.g., SREBP1c and PAX7) were associated with bovine economically important phenotypic traits, but none has been reported in the above-mentioned three cattle breeds. Since the InDel (insertion/deletion) marker becomes a useful tool applied in the animal molecular breeding, herein, we firstly found that the InDel variations of seven candidate genes in these cattle. Results showed that the genotypic and allelic distributions of these seven genes were remarkably different among these three cattle (p < 0.05 or p < 0.01). Furthermore, the InDel variations of SREBP1c and PAX7 genes were significantly associated with eight phenotypic traits in Xinjiang brown cattle (p < 0.05 or p < 0.01), respectively, suggesting that they can become the useful DNA markers.

DNA methylation pattern of the goat PITX1 gene and its effects on milk performance

Paired-like homeodomain transcription factor 1 (PITX1) is a pivotal gene in the hypothalamic–pituitary–adrenal axis, which is a well-known pathway affecting lactation performance. The aim of this study was to analyze the DNA methylation profile of the PITX1 gene and its relevance to milk performance in Xinong Saanen dairy goats; thus, potential epigenetic markers of lactation performance were identified. A total of 267 goat blood samples were divided into “low” and “high” groups according to two milk traits: the average milk yield (AMY) and the average milk density (AMD). One CpG island in the 3′-flanking region of the PITX1 gene was identified as being related to milk performance. Fisher's exact test demonstrated that the methylation rates of the overall CpG island and the 3rd and 12th CpG-dinucleotide loci in the blood were significantly associated with the AMY, and the overall methylation rate of the high AMY group was relative hypomethylation compared with the low AMY group. The overall methylation rates of this CpG island in mammary gland tissue from dry and lactation periods again exhibited a significant difference: the lactation period showed relative hypomethylation compared with the dry period. Bioinformatic transcription factor binding site prediction identified some lactation performance related transcription factors in this CpG island, such as CTCF, STAT, SMAD, CDEF, SP1, and KLFS. Briefly, overall methylation changes of the CpG island in the PITX1 gene are relevant to lactation performance, which will be valuable for future studies and epigenetic marker-assisted selection (eMAS) in the breeding of goats with respect to lactation performance.