Methylome Dynamics of Bovine Gametes and in vivo Early Embryos

A review of the role of transcription factors in regulating adipogenesis and lipogenesis in beef cattle

In the past few decades, genomic selection and other refined strategies have been used to increase the growth rate and lean meat production of beef cattle. Nevertheless, the fast growth rates of cattle breeds are often accompanied by a reduction in intramuscular fat (IMF) deposition, impairing meat quality. Transcription factors play vital roles in regulating adipogenesis and lipogenesis in beef cattle. Meanwhile, understanding the role of transcription factors in regulating adipogenesis and lipogenesis in beef cattle has gained significant attention to increase IMF deposition and meat quality. Therefore, the aim of this paper was to provide a comprehensive summary and valuable insight into the complex role of transcription factors in adipogenesis and lipogenesis in beef cattle. This review summarizes the contemporary studies in transcription factors in adipogenesis and lipogenesis, genome-wide analysis of transcription factors, epigenetic regulation of transcription factors, nutritional regulation of transcription factors, metabolic signalling pathways, functional genomics methods, transcriptomic profiling of adipose tissues, transcription factors and meat quality and comparative genomics with other livestock species. In conclusion, transcription factors play a crucial role in promoting adipocyte development and fatty acid biosynthesis in beef cattle. They control adipose tissue formation and metabolism, thereby improving meat quality and maintaining metabolic balance. Understanding the processes by which these transcription factors regulate adipose tissue deposition and lipid metabolism will simplify the development of marbling or IMF composition in beef cattle.

In Vitro Culture Alters Cell Lineage Composition and Cellular Metabolism of Bovine Blastocyst

Profiling bovine blastocyst transcriptome at the single-cell level has enabled us to reveal the first cell lineage segregation, during which the inner cell mass (ICM), trophectoderm (TE), and an undefined population of transitional cells were identified. By comparing the transcriptome of blastocysts derived in vivo (IVV), in vitro from a conventional culture medium (IVC), and in vitro from an optimized reduced nutrient culture medium (IVR), we found a delay of the cell fate commitment to ICM in the IVC and IVR embryos. Developmental potential differences between IVV, IVC, and IVR embryos were mainly contributed by ICM and transitional cells. Pathway analysis of these non-TE cells between groups revealed highly active metabolic and biosynthetic processes, reduced cellular signaling, and reduced transmembrane transport activities in IVC embryos that may lead to reduced developmental potential. IVR embryos had lower activities in metabolic and biosynthetic processes but increased cellular signaling and transmembrane transport, suggesting these cellular mechanisms may contribute to improved blastocyst development compared to IVC embryos. However, the IVR embryos had compromised development compared to IVV embryos with notably over-active transmembrane transport activities that impaired ion homeostasis.

Pleomorphic adenoma gene1 in reproduction and implication for embryonic survival in cattle: a review

The pleomorphic adenoma gene1 (PLAG1) encodes a DNA-binding, C2H2 zinc-finger protein which acts as a transcription factor that regulates the expression of diverse genes across different organs and tissues; hence, the name pleomorphic. Rearrangements of the PLAG1 gene, and/or overexpression, are associated with benign tumors and cancers in a variety of tissues. This is best described for pleomorphic adenoma of the salivary glands in humans. The most notable expression of PLAG1 occurs during embryonic and fetal development, with lesser expression after birth. Evidence has accumulated of a role for PLAG1 protein in normal early embryonic development and placentation in mammals. PLAG1 protein influences the expression of the ike growth factor 2 (IGF2) gene and production of IGF2 protein. IGF2 is an important mitogen in ovarian follicles/oocytes, embryos, and fetuses. The PLAG1-IGF2 axis, therefore, provides one pathway whereby PLAG1 protein can influence embryonic survival and pregnancy. PLAG1 also influences over 1,000 other genes in embryos including those associated with ribosomal assembly and proteins. Brahman (Bos indicus) heifers homozygous for the PLAG1 variant, rs109815800 (G > T), show greater fertility than contemporary heifers with either one, or no copy, of the variant. Greater fertility in heifers homozygous for rs109815800 could be the result of early puberty and/or greater embryonic survival. The present review first looks at the broader roles of the PLAG1 gene and PLAG1 protein and then focuses on the emerging role of PLAG1/PLAG1 in embryonic development and pregnancy. A deeper understanding of factors which influence embryonic development is required for the next transformational increase in embryonic survival and successful pregnancy for both in vivo and in vitro derived embryos in cattle.

Molecular and cellular programs underlying the development of bovine pre-implantation embryos

Early embryonic mortality is a major cause of infertility in cattle, yet the underlying molecular causes remain a mystery. Over the past half century, assisted reproductive technologies such as in vitro fertilisation and somatic cell nuclear transfer have been used to improve cattle reproductive efficiency; however, reduced embryo developmental potential is seen compared to their in vivo counterparts. Recent years have seen exciting progress across bovine embryo research, including genomic profiling of embryogenesis, new methods for improving embryo competence, and experimenting on building bovine embryos from stem cell cultures. These advances are beginning to define bovine embryo molecular and cellular programs and could potentially lead to improved embryo health. Here, I highlight the current status of molecular determinants and cellular programs of bovine embryo development and new opportunities to improve the bovine embryo health.

Maternal restricted- and over- feeding during gestation perturb offspring sperm epigenome in sheep

In brief

Inadequate maternal nutrition during gestation can have immediate and lifelong effects on offspring. This study shows that maternal restricted - and over- nutrition during gestation do not affect semen characteristics in F1 male offspring but alters offspring sperm sncRNA profiles and DNA methylome in sheep.

Abstract

There is a growing body of evidence that inadequate maternal nutrition during gestation can have immediate and lifelong effects on offspring. However, little is known about the effects of maternal nutrition during gestation on male offspring reproduction. Here, using a sheep model of maternal restricted - and over - nutrition (60 or 140% of the National Research Council requirements) during gestation, we found that maternal restricted - and over - nutrition do not affect semen characteristics (i.e. volume, sperm concentration, pH, sperm motility, sperm morphology) or scrotal circumference in male F1 offspring. However, using small RNA sequencing analysis, we demonstrated that both restricted - and over - nutrition during gestation induced marked changes in composition and expression of sperm small noncoding RNAs (sncRNAs) subpopulations including in male F1 offspring. Whole-genome bisulfite sequencing analysis further identified specific genomic loci where poor maternal nutrition resulted in alterations in DNA methylation. These findings indicate that maternal restricted - and over - nutrition during gestation induce epigenetic modifications in sperm of F1 offspring sperm in sheep, which may contribute to environmentally influenced phenotypes in ruminants.

Third-Generation Sequencing Reveals the Adaptive Role of the Epigenome in Three Deep-Sea Polychaetes

The roles of DNA methylation in invertebrates are poorly characterized, and critical data are missing for the phylum Annelida. We fill this knowledge gap by conducting the first genome-wide survey of DNA methylation in the deep-sea polychaetes dominant in deep-sea vents and seeps: Paraescarpia echinospica, Ridgeia piscesae, and Paralvinella palmiformis. DNA methylation calls were inferred from Oxford Nanopore sequencing after assembling high-quality genomes of these animals. The genomes of these worms encode all the key enzymes of the DNA methylation metabolism and possess a mosaic methylome similar to that of other invertebrates. Transcriptomic data of these polychaetes support the hypotheses that gene body methylation strengthens the expression of housekeeping genes and that promoter methylation acts as a silencing mechanism but not the hypothesis that DNA methylation suppresses the activity of transposable elements. The conserved epigenetic profiles of genes responsible for maintaining homeostasis under extreme hydrostatic pressure suggest DNA methylation plays an important adaptive role in these worms.

Pre-Implantation Bovine Embryo Evaluation-From Optics to Omics and Beyond

Approximately 80% of the ~1.5 million bovine embryos transferred in 2021 were in vitro produced. However, only ~27% of the transferred IVP embryos will result in live births. The ~73% pregnancy failures are partly due to transferring poor-quality embryos, a result of erroneous stereomicroscopy-based morphological evaluation, the current method of choice for pre-transfer embryo evaluation. Numerous microscopic (e.g., differential interference contrast, electron, fluorescent, time-lapse, and artificial-intelligence-based microscopy) and non-microscopic (e.g., genomics, transcriptomics, epigenomics, proteomics, metabolomics, and nuclear magnetic resonance) methodologies have been tested to find an embryo evaluation technique that is superior to morphologic evaluation. Many of these research tools can accurately determine embryo quality/viability; however, most are invasive, expensive, laborious, technically sophisticated, and/or time-consuming, making them futile in the context of in-field embryo evaluation. However accurate they may be, using complex methods, such as RNA sequencing, SNP chips, mass spectrometry, and multiphoton microscopy, at thousands of embryo production/collection facilities is impractical. Therefore, future research is warranted to innovate field-friendly, simple benchtop tests using findings already available, particularly from omics-based research methodologies. Time-lapse monitoring and artificial-intelligence-based automated image analysis also have the potential for accurate embryo evaluation; however, further research is warranted to innovate economically feasible options for in-field applications.

Establishment of bovine trophoblast stem cells

Here, we report that a chemical cocktail (LCDM: leukemia inhibitory factor [LIF], CHIR99021, dimethinedene maleate [DiM], minocycline hydrochloride), previously developed for extended pluripotent stem cells (EPSCs) in mice and humans, enables de novo derivation and long-term culture of bovine trophoblast stem cells (TSCs). Bovine TSCs retain developmental potency to differentiate into mature trophoblast cells and exhibit transcriptomic and epigenetic (chromatin accessibility and DNA methylome) features characteristic of trophectoderm cells from early bovine embryos. The bovine TSCs established in this study will provide a model to study bovine placentation and early pregnancy failure.

Intergenic transcription in in vivo developed bovine oocytes and pre-implantation embryos

Background

Intergenic transcription, either failure to terminate at the transcription end site (TES), or transcription initiation at other intergenic regions, is present in cultured cells and enhanced in the presence of stressors such as viral infection. Transcription termination failure has not been characterized in natural biological samples such as pre-implantation embryos which express more than 10,000 genes and undergo drastic changes in DNA methylation.

Results

Using Automatic Readthrough Transcription Detection (ARTDeco) and data of in vivo developed bovine oocytes and embryos, we found abundant intergenic transcripts that we termed as read-outs (transcribed from 5 to 15 kb after TES) and read-ins (transcribed 1 kb up-stream of reference genes, extending up to 15 kb up-stream). Read-throughs (continued transcription from TES of expressed reference genes, 4-15 kb in length), however, were much fewer. For example, the numbers of read-outs and read-ins ranged from 3,084 to 6,565 or 33.36-66.67% of expressed reference genes at different stages of embryo development. The less copious read-throughs were at an average of 10% and significantly correlated with reference gene expression (P < 0.05). Interestingly, intergenic transcription did not seem to be random because many intergenic transcripts (1,504 read-outs, 1,045 read-ins, and 1,021 read-throughs) were associated with common reference genes across all stages of pre-implantation development. Their expression also seemed to be regulated by developmental stages because many were differentially expressed (log2 fold change ≥ 2, P < 0.05). Additionally, while gradual but un-patterned decreases in DNA methylation densities 10 kb both up- and down-stream of the intergenic transcribed regions were observed, the correlation between intergenic transcription and DNA methylation was insignificant. Finally, transcription factor binding motifs and polyadenylation signals were found in 27.2% and 12.15% of intergenic transcripts, respectively, suggesting considerable novel transcription initiation and RNA processing.

Conclusion

In summary, in vivo developed oocytes and pre-implantation embryos express large numbers of intergenic transcripts, which are not related to the overall DNA methylation profiles either up- or down-stream.

Differentially methylated regions identified in bovine embryos are not observed in adulthood

The establishment of epigenetic marks during the reprogramming window is susceptible to environmental influences, and stimuli during this critical stage can cause altered DNA methylation in offspring. In a previous study, we found that low levels of sulphur and cobalt (low S/Co) in the diet offered to oocyte donors altered the DNA methylome of bovine embryos. However, due to the extensive epigenetic reprogramming that occurs during embryogenesis, we hypothesized that the different methylation regions (DMRs) identified in the blastocysts may not maintain in adulthood. Here, we aimed to characterize DMRs previously identified in embryos, in the blood and sperm of adult progenies of two groups of heifers (low S/Co and control). We used six bulls and characterized the DNA methylation levels of KDM2A, KDM5A, KMT2D, and DOT1L genes. Our results showed that all DMRs analysed in both groups and tissues were hypermethylated unlike that noticed in the embryonic methylome profiles. These results suggest that embryo DMRs were reprogrammed during the final stages of de novo methylation during embryogenesis or later in development. Therefore, due to the highly dynamic epigenetic state during early embryonic development, we suggest that is essential to validate the DMRs found in embryos in adult individuals.

Current status of genome-wide epigenetic profiling of mammalian preimplantation embryos

Background

Genome-wide information on epigenetic modifications in mammalian preimplantation embryos was an unexplored sanctuary of valuable research insights protected by the difficulty of its analysis. However, that is no longer the case, and many epigenome maps are now available for sightseeing there.

Methods

This review overviews the current status of genome-wide epigenetic profiling in terms of DNA methylome and histone modifications in mammalian preimplantation embryos.

Main findings

As the sensitivity of methods for analyzing epigenetic modifications increased, pioneering work began to explore the genome-wide epigenetic landscape in the mid-2010s, first for DNA methylation and then for histone modifications. Since then, a huge amount of data has accumulated, revealing typical epigenetic profiles in preimplantation development and, more recently, changes in response to environmental interventions.

Conclusions

These accumulating data may be used to improve the quality of preimplantation embryos, both in terms of their short-term developmental competence and their subsequent long-term health implications.

One-carbon metabolite supplementation to heifers for the first 14 d of the estrous cycle alters the plasma and hepatic one-carbon metabolite pool and methionine-folate cycle enzyme transcript abundance in a dose-dependent manner

The objective of this study was to determine the dose of folate and vitamin B12 in beef heifers fed rumen protected methionine and choline required to maintain increased B12 levels and intermediates of the methionine-folate cycle in circulation. Angus heifers (n = 30; BW = 392.6 ± 12.6 kg) were individually fed and assigned to one of five treatments: 0XNEG: Total mixed ration (TMR) and saline injections at day 0 and 7 of the estrous cycle, 0XPOS: TMR, rumen protected methionine (MET) fed at 0.08% of the diet DM, rumen protected choline (CHOL) fed at 60 g/d, and saline injections at day 0 and 7, 0.5X: TMR, MET, CHOL, 5 mg B12, and 80 mg folate at day 0 and 7, 1X: TMR, MET CHOL, 10 mg vitamin B12, and 160 mg folate at day 0 and 7, and 2X: TMR, MET, CHOL, 20 mg B12, and 320 mg folate at day 0 and 7. All heifers were estrus synchronized but not bred, and blood was collected on day 0, 2, 5, 7, 9, 12, and 14 of a synchronized estrous cycle. Heifers were slaughtered on day 14 of the estrous cycle for liver collection. Serum B12 concentrations were greater in the 0.5X, 1X, and 2X, compared with 0XNEG and 0XPOS on all days after treatment initiation (P < 0.0001). Serum folate concentrations were greater for the 2X treatment at day 5, 7, and 9 of the cycle compared with all other treatments (P ≤ 0.05). There were no differences (P ≥ 0.19) in hepatic methionine-cycle or choline analyte concentrations by treatment. Concentrations of hepatic folate cycle intermediates were always greater (P ≤ 0.04) in the 2X treatment compared with the 0XNEG and 0XPOS heifers. Serum methionine was greater (P = 0.04) in the 0.5X and 2X heifers compared with 0XNEG, and S-adenosylhomocysteine (SAH) tended (P = 0.06) to be greater in the 0.5X heifers and the S-adenosylmethionine (SAM):SAH ratio was decreased (P = 0.05) in the 0.5X treatment compared with the 0XNEG, 0XPOS, and 2X heifers. The hepatic transcript abundance of MAT2A and MAT2B were decreased (P ≤ 0.02) in the 0.5X heifers compared with the 0XNEG, 0XPOS, and 2X heifers. These data support that beef heifers fed rumen protected methionine and choline require 20 mg B12 and 320 mg folate once weekly to maintain increased concentrations of B12 and folate in serum. Furthermore, these data demonstrate that not all supplementation levels are equal in providing positive responses, and that some levels, such as the 0.5X, may result in a stoichiometric imbalance in the one-carbon metabolism pathway that results in a decreased SAM:SAH ratio.

Metabolism-epigenetic interactions on in vitro produced embryos

Metabolism and epigenetics, which reciprocally regulate each other in different cell types, are fundamental aspects of cellular adaptation to the environment. Evidence in cancer and stem cells has shown that the metabolic status modifies the epigenome while epigenetic mechanisms regulate the expression of genes involved in metabolic processes, thereby altering the metabolome. This crosstalk occurs as many metabolites serve as substrates or cofactors of chromatin-modifying enzymes. If we consider the intense metabolic dynamic and the epigenetic remodelling of the embryo, the comprehension of these regulatory networks will be important not only for understanding early embryonic development, but also to determine in vitro culture conditions that support embryo development and may insert positive regulatory marks that may persist until adult life. In this review, we focus on how metabolism may affect epigenetic reprogramming of the early stages of development, in particular acetylation and methylation of histone and DNA. We also present other metabolic modifications in bovine embryos, such as lactylation, highlighting the promising epigenetic and metabolic targets to improve conditions for in vitro embryo development.

Metabolic enhancement of the one carbon metabolism (OCM) in bovine oocytes IVM increases the blastocyst rate: evidences for a OCM checkpoint

The one carbon metabolism (OCM) has a primary role in the process of oocyte maturation. In this study bovine oocytes were cultured for 24 h, up to MII stage, with standard medium supplemented or not with 8 metabolic enhancers of the OCM and the MII and blastocyst rate were compared. Additional analyses were performed on matured oocytes, cumulus cells, zygotes and blastocysts. The OCM supplementation increased the blastocyst rate derived from in vitro fertilization. The mitochondrial mass and DNMT3a protein expression were increased whereas DNA fragmentation decreased in matured oocytes. DNA methylation in female pronucleus of zygotes was increased. The supplementation did not directly affect the redox balance as ROS and GSH in matured oocytes and homocysteine in the spent medium were unchanged. The supplementation of the oocytes with metabolic enhancers of the OCM may increase the yield from the culture, likely due to improved DNA methylation and epigenetic programming. The lack of effects on MII rate with huge differences appearing at the blastocyst stage suggest the existence of a OCM metabolic check point that hampers oocytes progression to blastocyst post-fertilization, if they were not properly primed at the time of maturation.

High-resolution ribosome profiling reveals translational selectivity for transcripts in bovine preimplantation embryo development

High-resolution ribosome fractionation and low-input ribosome profiling of bovine oocytes and preimplantation embryos has enabled us to define the translational landscapes of early embryo development at an unprecedented level. We analyzed the transcriptome and the polysome- and non-polysome-bound RNA profiles of bovine oocytes (germinal vesicle and metaphase II stages) and early embryos at the two-cell, eight-cell, morula and blastocyst stages, and revealed four modes of translational selectivity: (1) selective translation of non-abundant mRNAs; (2) active, but modest translation of a selection of highly expressed mRNAs; (3) translationally suppressed abundant to moderately abundant mRNAs; and (4) mRNAs associated specifically with monosomes. A strong translational selection of low-abundance transcripts involved in metabolic pathways and lysosomes was found throughout bovine embryonic development. Notably, genes involved in mitochondrial function were prioritized for translation. We found that translation largely reflected transcription in oocytes and two-cell embryos, but observed a marked shift in the translational control in eight-cell embryos that was associated with the main phase of embryonic genome activation. Subsequently, transcription and translation become more synchronized in morulae and blastocysts. Taken together, these data reveal a unique spatiotemporal translational regulation that accompanies bovine preimplantation development.

Completing Single-Cell DNA Methylome Profiles via Transfer Learning Together With KL-Divergence

The high level of sparsity in methylome profiles obtained using whole-genome bisulfite sequencing in the case of low biological material amount limits its value in the study of systems in which large samples are difficult to assemble, such as mammalian preimplantation embryonic development. The recently developed computational methods for addressing the sparsity by imputing missing have their limits when the required minimum data coverage or profiles of the same tissue in other modalities are not available. In this study, we explored the use of transfer learning together with Kullback-Leibler (KL) divergence to train predictive models for completing methylome profiles with very low coverage (below 2%). Transfer learning was used to leverage less sparse profiles that are typically available for different tissues for the same species, while KL divergence was employed to maximize the usage of information carried in the input data. A deep neural network was adopted to extract both DNA sequence and local methylation patterns for imputation. Our study of training models for completing methylome profiles of bovine oocytes and early embryos demonstrates the effectiveness of transfer learning and KL divergence, with individual increase of 29.98 and 29.43%, respectively, in prediction performance and 38.70% increase when the two were used together. The drastically increased data coverage (43.80-73.6%) after imputation powers downstream analyses involving methylomes that cannot be effectively done using the very low coverage profiles (0.06-1.47%) before imputation.

An Interplay between Epigenetics and Translation in Oocyte Maturation and Embryo Development: Assisted Reproduction Perspective

Germ cell quality is a key prerequisite for successful fertilization and early embryo development. The quality is determined by the fine regulation of transcriptomic and proteomic profiles, which are prone to alteration by assisted reproduction technology (ART)-introduced in vitro methods. Gaining evidence shows the ART can influence preset epigenetic modifications within cultured oocytes or early embryos and affect their developmental competency. The aim of this review is to describe ART-determined epigenetic changes related to the oogenesis, early embryogenesis, and further in utero development. We confront the latest epigenetic, related epitranscriptomic, and translational regulation findings with the processes of meiotic maturation, fertilization, and early embryogenesis that impact the developmental competency and embryo quality. Post-ART embryo transfer, in utero implantation, and development (placentation, fetal development) are influenced by environmental and lifestyle factors. The review is emphasizing their epigenetic and ART contribution to fetal development. An epigenetic parallel among mouse, porcine, and bovine animal models and human ART is drawn to illustrate possible future mechanisms of infertility management as well as increase the awareness of the underlying mechanisms governing oocyte and embryo developmental complexity under ART conditions.

Sperm DNA methylation patterns at discrete CpGs and genes involved in embryonic development are related to bull fertility

Background

Despite a multifactorial approach being taken for the evaluation of bull semen quality in many animal breeding centres worldwide, reliable prediction of bull fertility is still a challenge. Recently, attention has turned to molecular mechanisms, which could uncover potential biomarkers of fertility. One of these mechanisms is DNA methylation, which together with other epigenetic mechanisms is essential for the fertilising sperm to drive normal embryo development and establish a viable pregnancy. In this study, we hypothesised that bull sperm DNA methylation patterns are related to bull fertility. We therefore investigated DNA methylation patterns from bulls used in artificial insemination with contrasting fertility scores.

Results

The DNA methylation patterns were obtained by reduced representative bisulphite sequencing from 10 high-fertility bulls and 10 low-fertility bulls, having average fertility scores of − 6.6 and + 6.5%, respectively (mean of the population was zero). Hierarchical clustering analysis did not distinguish bulls based on fertility but did highlight individual differences. Despite this, using stringent criteria (DNA methylation difference ≥ 35% and a q-value < 0.001), we identified 661 differently methylated cytosines (DMCs). DMCs were preferentially located in intergenic regions, introns, gene downstream regions, repetitive elements, open sea, shores and shelves of CpG islands. We also identified 10 differently methylated regions, covered by 7 unique genes (SFRP1, STXBP4, BCR, PSMG4, ARSG, ATP11A, RXRA), which are involved in spermatogenesis and early embryonic development.

Conclusion

This study demonstrated that at specific CpG sites, sperm DNA methylation status is related to bull fertility, and identified seven differently methylated genes in sperm of subfertile bulls that may lead to altered gene expression and potentially influence embryo development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08614-5.

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.

Mining RNAseq data reveals dynamic metaboloepigenetic profiles in human, mouse and bovine pre-implantation embryos

Metaboloepigenetic regulation has been reported in stem cells, germ cells, and tumor cells. Embryonic metaboloepigenetics, however, have just begun to be described. Here we analyzed RNAseq data to characterize the metaboloepigenetic profiles of human, mouse, and bovine pre-implantation embryos. In embryos, metaboloepigenetic reprogramming was species-specific, varied with the developmental stage and was disrupted with in vitro culture. Metabolic pathways and gene expressions were strongly correlated with early embryo DNA methylation and were changed with in vitro culture. Although the idea that the in vitro environment may influence development is not new, there has been little progress on improving pregnancy rates after decades using in vitro fertilization. Hence, the present data will contribute to understanding how the in vitro manipulation affects the metaboloepigenetic status of early embryos, which can be used to establish culture strategies aimed at improving the in vitro environment and, consequently, pregnancy rates and offspring health.

•

Embryonic metaboloepigenetic reprogramming is stage- and species-specific

•

In vitro culture disrupts the in vivo embryonic metaboloepigenetic reprogramming

•

Metabolic genes and pathways are highly correlated with embryo methylome

Low guanine content and biased nucleotide distribution in vertebrate mtDNA can cause overestimation of non-CpG methylation

Mitochondrial DNA (mtDNA) methylation in vertebrates has been hotly debated for over 40 years. Most contrasting results have been reported following bisulfite sequencing (BS-seq) analyses. We addressed whether BS-seq experimental and analysis conditions influenced the estimation of the levels of methylation in specific mtDNA sequences. We found false positive non-CpG methylation in the CHH context (fpCHH) using unmethylated Sus scrofa mtDNA BS-seq data. fpCHH methylation was detected on the top/plus strand of mtDNA within low guanine content regions. These top/plus strand sequences of fpCHH regions would become extremely AT-rich sequences after BS-conversion, whilst bottom/minus strand sequences remained almost unchanged. These unique sequences caused BS-seq aligners to falsely assign the origin of each strand in fpCHH regions, resulting in false methylation calls. fpCHH methylation detection was enhanced by short sequence reads, short library inserts, skewed top/bottom read ratios and non-directional read mapping modes. We confirmed no detectable CHH methylation in fpCHH regions by BS-amplicon sequencing. The fpCHH peaks were located in the D-loop, ATP6, ND2, ND4L, ND5 and ND6 regions and identified in our S. scrofa ovary and oocyte data and human BS-seq data sets. We conclude that non-CpG methylation could potentially be overestimated in specific sequence regions by BS-seq analysis.

Whole-Genome DNA Methylation Dynamics of Sheep Preimplantation Embryo Investigated by Single-Cell DNA Methylome Sequencing

The early stages of mammalian embryonic development involve the participation and cooperation of numerous complex processes, including nutritional, genetic, and epigenetic mechanisms. However, in embryos cultured in vitro, a developmental block occurs that affects embryo development and the efficiency of culture. Although the block period is reported to involve the transcriptional repression of maternal genes and transcriptional activation of zygotic genes, how epigenetic factors regulate developmental block is still unclear. In this study, we systematically analyzed whole-genome methylation levels during five stages of sheep oocyte and preimplantation embryo development using single-cell level whole genome bisulphite sequencing (SC-WGBS) technology. Then, we examined several million CpG sites in individual cells at each evaluated developmental stage to identify the methylation changes that take place during the development of sheep preimplantation embryos. Our results showed that two strong waves of methylation changes occurred, namely, demethylation at the 8-cell to 16-cell stage and methylation at the 16-cell to 32-cell stage. Analysis of DNA methylation patterns in different functional regions revealed a stable hypermethylation status in 3'UTRs and gene bodies; however, significant differences were observed in intergenic and promoter regions at different developmental stages. Changes in methylation at different stages of preimplantation embryo development were also compared to investigate the molecular mechanisms involved in sheep embryo development at the methylation level. In conclusion, we report a detailed analysis of the DNA methylation dynamics during the development of sheep preimplantation embryos. Our results provide an explanation for the complex regulatory mechanisms underlying the embryo developmental block based on changes in DNA methylation levels.

Whole Genome Methylation Analysis Reveals Role of DNA Methylation in Cow's Ileal and Ileal Lymph Node Responses to Mycobacterium avium subsp. paratuberculosis Infection

Johne's Disease (JD), caused by Mycobacterium avium subsp paratuberculosis (MAP), is an incurable disease of ruminants and other animal species and is characterized by an imbalance of gut immunity. The role of MAP infection on the epigenetic modeling of gut immunity during the progression of JD is still unknown. This study investigated the DNA methylation patterns in ileal (IL) and ileal lymph node (ILLN) tissues from cows diagnosed with persistent subclinical MAP infection over a one to 4 years period. DNA samples from IL and ILLN tissues from cows negative (MAPneg) (n = 3) or positive for MAP infection (MAPinf) (n = 4) were subjected to whole genome bisulfite sequencing. A total of 11,263 and 62,459 differentially methylated cytosines (DMCs), and 1259 and 8086 differentially methylated regions (DMRs) (FDR<0.1) were found between MAPinf and MAPneg IL and ILLN tissues, respectively. The DMRs were found on 394 genes (denoted DMR genes) in the IL and on 1305 genes in the ILLN. DMR genes with hypermethylated promoters/5'UTR [3 (IL) and 88 (ILLN)] or hypomethylated promoters/5'UTR [10 (IL) and 25 (ILLN)] and having multiple functions including response to stimulus/immune response (BLK, BTC, CCL21, AVPR1A, CHRNG, GABRA4, TDGF1), cellular processes (H2AC20, TEX101, GLA, NCKAP5L, RBM27, SLC18A1, H2AC20BARHL2, NLGN3, SUV39H1, GABRA4, PPA1, UBE2D2) and metabolic processes (GSTO2, H2AC20, SUV39H1, PPA1, UBE2D2) are potential DNA methylation candidate genes of MAP infection. The ILLN DMR genes were enriched for more biological process (BP) gene ontology (GO) terms (n = 374), most of which were related to cellular processes (27.6%), biological regulation (16.6%), metabolic processes (15.4%) and response to stimulus/immune response (8.2%) compared to 75 BP GO terms (related to cellular processes, metabolic processes and transport, and system development) enriched for IL DMR genes. ILLN DMR genes were enriched for more pathways (n = 47) including 13 disease pathways compared with 36 enriched pathways, including 7 disease/immune pathways for IL DMR genes. In conclusion, the results show tissue specific responses to MAP infection with more epigenetic changes (DMCs and DMRs) in the ILLN than in the IL tissue, suggesting that the ILLN and immune processes were more responsive to regulation by methylation of DNA relative to IL tissue. Our data is the first to demonstrate a potential role for DNA methylation in the pathogenesis of MAP infection in dairy cattle.

Evolutionary epigenomic analyses in mammalian early embryos reveal species-specific innovations and conserved principles of imprinting

While mouse remains the most popular model, the conservation of parental-to-embryonic epigenetic transition across mammals is poorly defined. Through analysis of oocytes and early embryos in human, bovine, porcine, rat, and mouse, we revealed remarkable species-specific innovations as no single animal model fully recapitulates the human epigenetic transition. In rodent oocytes, transcription-dependent DNA methylation allows methylation of maternal imprints but not intergenic paternal imprints. Unexpectedly, prevalent DNA hypermethylation, paralleled by H3K36me2/3, also occurs in nontranscribed regions in porcine and bovine oocytes, except for megabase-long “CpG continents (CGCs)” where imprinting control regions preferentially reside. Broad H3K4me3 and H3K27me3 domains exist in nonhuman oocytes, yet only rodent H3K27me3 survives beyond genome activation. Coincidently, regulatory elements preferentially evade H3K27me3 in rodent oocytes, and failure to do so causes aberrant embryonic gene repression. Hence, the diverse mammalian innovations of parental-to-embryonic transition center on a delicate “to-methylate-or-not” balance in establishing imprints while protecting other regulatory regions.

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.

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.

Dysregulated Gene Expression of Imprinted and X-Linked Genes: A Link to Poor Development of Bovine Haploid Androgenetic Embryos

Mammalian uniparental embryos are efficient models for genome imprinting research and allow studies on the contribution of the paternal and maternal genomes to early embryonic development. In this study, we analyzed different methods for production of bovine haploid androgenetic embryos (hAE) to elucidate the causes behind their poor developmental potential. Results indicate that hAE can be efficiently generated by using intracytoplasmic sperm injection and oocyte enucleation at telophase II. Although androgenetic haploidy does not disturb early development up to around the 8-cell stage, androgenetic development is disturbed after the time of zygote genome activation and hAE that reach the morula stage are less capable to reach the blastocyst stage of development. Karyotypic comparisons to parthenogenetic- and ICSI-derived embryos excluded chromosomal segregation errors as causes of the developmental constraints of hAE. However, analysis of gene expression indicated abnormal levels of transcripts for key long non-coding RNAs involved in X chromosome inactivation and genomic imprinting of the KCNQ1 locus, suggesting an association with X chromosome and some imprinted loci. Moreover, transcript levels of methyltransferase 3B were significantly downregulated, suggesting potential anomalies in hAE establishing de novo methylation. Finally, the methylation status of imprinted control regions for XIST and KCNQ1OT1 genes remained hypomethylated in hAE at the morula and blastocyst stages, confirming their origin from spermatozoa. Thus, our results exclude micromanipulation and chromosomal abnormalities as major factors disturbing the normal development of bovine haploid androgenotes. In addition, although the cause of the arrest remains unclear, we have shown that the inefficient development of haploid androgenetic bovine embryos to develop to the blastocyst stage is associated with abnormal expression of key factors involved in X chromosome activity and genomic imprinting.

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.

The Epigenetics of Gametes and Early Embryos and Potential Long-Range Consequences in Livestock Species-Filling in the Picture With Epigenomic Analyses

The epigenome is dynamic and forged by epigenetic mechanisms, such as DNA methylation, histone modifications, chromatin remodeling, and non-coding RNA species. Increasing lines of evidence support the concept that certain acquired traits are derived from environmental exposure during early embryonic and fetal development, i.e., fetal programming, and can even be "memorized" in the germline as epigenetic information and transmitted to future generations. Advances in technology are now driving the global profiling and precise editing of germline and embryonic epigenomes, thereby improving our understanding of epigenetic regulation and inheritance. These achievements open new avenues for the development of technologies or potential management interventions to counteract adverse conditions or improve performance in livestock species. In this article, we review the epigenetic analyses (DNA methylation, histone modification, chromatin remodeling, and non-coding RNAs) of germ cells and embryos in mammalian livestock species (cattle, sheep, goats, and pigs) and the epigenetic determinants of gamete and embryo viability. We also discuss the effects of parental environmental exposures on the epigenetics of gametes and the early embryo, and evidence for transgenerational inheritance in livestock.

Dynamics of paternal contributions to early embryo development in large animals

This review focuses on current knowledge of paternal contributions to preimplantation embryonic development with particular emphasis on large animals. Specifically, the included content aims to summarize genomic and epigenomic contributions of paternally expressed genes, their regulation, and chromatin structure that are indispensable for early embryo development. The accumulation of current knowledge will summarize conserved allelic function among species to include functional molecular and genomic studies across large domestic animals in context with reference to founding experimental models.

Tricarboxylic Acid Cycle Metabolites as Mediators of DNA Methylation Reprogramming in Bovine Preimplantation Embryos

In many cell types, epigenetic changes are partially regulated by the availability of metabolites involved in the activity of chromatin-modifying enzymes. Even so, the association between metabolism and the typical epigenetic reprogramming that occurs during preimplantation embryo development remains poorly understood. In this work, we explore the link between energy metabolism, more specifically the tricarboxylic acid cycle (TCA), and epigenetic regulation in bovine preimplantation embryos. Using a morphokinetics model of embryonic development (fast- and slow-developing embryos), we show that DNA methylation (5mC) and hydroxymethylation (5hmC) are dynamically regulated and altered by the speed of the first cleavages. More specifically, slow-developing embryos fail to perform the typical reprogramming that is necessary to ensure the generation of blastocysts with higher ability to establish specific cell lineages. Transcriptome analysis revealed that such differences were mainly associated with enzymes involved in the TCA cycle rather than specific writers/erasers of DNA methylation marks. This relationship was later confirmed by disturbing the embryonic metabolism through changes in α-ketoglutarate or succinate availability in culture media. This was sufficient to interfere with the DNA methylation dynamics despite the fact that blastocyst rates and total cell number were not quite affected. These results provide the first evidence of a relationship between epigenetic reprogramming and energy metabolism in bovine embryos. Likewise, levels of metabolites in culture media may be crucial for precise epigenetic reprogramming, with possible further consequences in the molecular control and differentiation of cells.

Locus-specific analysis of DNA methylation patterns in cloned and in vitro fertilized porcine embryos

Porcine somatic cell nuclear transfer (SCNT) is currently inefficient, as 1–3.95% of reconstructed embryos survive to term; inadequate or erroneous epigenetic reprogramming of the specialized donor somatic nucleus could be a primary reason. Therefore, a locus-specific analysis of DNA methylation dynamics in embryogenesis and the DNA methylation status of gametes and donor cells used for SCNT were conducted in the following developmentally important gene loci: POU5F1, NANOG, SOX2, H19, IGF2, IGF2R, XIST; and the retrotransposon LINE-1. There were significant epigenetic differences between the gametes and the somatic donor cells. Three gamete-specific differentially methylated regions (DMRs) in POU5F1, XIST, and LINE-1 were identified. A delayed demethylation process at POU5F1 and LINE-1 loci occurred after three successive cleavages, compared to the in vitro fertilized (IVF) embryos. Although cloned embryos could undergo de-methylation and re-methylation dynamics at the DMRs of imprinted genes (H19,IGF2R, and XIST), the re-methylation process was compromised, unlike in fertilized embryos. LINE-1 loci are widely dispersed across the whole genome, and LINE-1 DMR might be a potential porcine nuclear reprogramming epi-marker. Data from observations in our present and previous studies, and two published articles were pooled to produce a schematic diagram of locus-specific, DNA methylation dynamics of cloned and IVF embryos during porcine early embryogenesis. This also indicated aberrant DNA methylation reprogramming events, including inadequate DNA demethylation and insufficient re-methylation in cloned embryos. Further research should focus on mechanisms underlying demethylation during the early cleavage of embryos and de novo DNA methylation at the blastocyst stage.

Sperm DNA Hypomethylation Proximal to Reproduction Pathway Genes in Maturing Elite Norwegian Red Bulls

Genomic selection in modern farming demands sufficient semen production in young bulls. Factors affecting semen quality and production capacity in young bulls are not well understood; DNA methylation, a complicated phenomenon in sperm cells, is one such factors. In this study, fresh and frozen-thawed semen samples from the same Norwegian Red (NR) bulls at both 14 and 17 months of age were examined for sperm chromatin integrity parameters, ATP content, viability, and motility. Furthermore, reduced representation bisulfite libraries constructed according to two protocols, the Ovation^® RRBS Methyl-Seq System (Ovation method) and a previously optimized gel-free method and were sequenced to study the sperm DNA methylome in frozen-thawed semen samples. Sperm quality analyses indicated that sperm concentration, total motility and progressivity in fresh semen from 17 months old NR bulls were significantly higher compared to individuals at 14 months of age. The percentage of DNA fragmented sperm cells significantly decreased in both fresh and frozen-thawed semen samples in bulls with increasing age. Libraries from the Ovation method exhibited a greater percentage of read loss and shorter read size following trimming. Downstream analyses for reads obtained from the gel-free method revealed similar global sperm DNA methylation but differentially methylated regions (DMRs) between 14- and 17 months old NR bulls. The majority of identified DMRs were hypomethylated in 14 months old bulls. Most of the identified DMRs (69%) exhibited a less than 10% methylation difference while only 1.5% of DMRs exceeded a 25% methylation difference. Pathway analysis showed that genes annotated with DMRs having low methylation differences (less than 10%) and DMRs having between 10 and 25% methylation differences, could be associated with important hormonal signaling and sperm function relevant pathways, respectively. The current research shows that RRBS in parallel with routine sperm quality analyses could be informative in reproductive capacity of young NR bulls. Although global sperm DNA methylation levels in 14 and 17 months old NR bulls were similar, regions with low and varying levels of DNA methylation differences can be identified and linked with important sperm function and hormonal pathways.

Erasing gametes to write blastocysts: metabolism as the new player in epigenetic reprogramming

Understanding preimplantation embryonic development is crucial for the improvement of assisted reproductive technologies and animal production. To achieve this goal, it is important to consider that gametes and embryos are highly susceptible to environmental changes. Beyond the metabolic adaptation, the dynamic status imposed during follicular growth and early embryogenesis may create marks that will guide the molecular regulation during prenatal development, and consequently impact the offspring phenotype. In this context, metaboloepigenetics has gained attention, as it investigates the crosstalk between metabolism and molecular control, i.e., how substrates generated by metabolic pathways may also act as players of epigenetic modifications. In this review, we present the main metabolic and epigenetic events of pre-implantation development, and how these systems connect to open possibilities for targeted manipulation of reproductive technologies and animal production systems.

The landscape of accessible chromatin in bovine oocytes and early embryos

Chromatin reorganization governs the regulation of gene expression during preimplantation development. However, the landscape of chromatin dynamics in this period has not been explored in bovine. In this study, we constructed a genome-wide map of accessible chromatin in bovine oocytes and early embryos using an improved assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) which revealed unique features of the accessible chromatin during bovine early embryo development. We found that chromatin accessibility is low in oocytes and 2-/4-cell embryos, followed by a significant increase in embryos during major embryonic genome activation (EGA), and peaked in elongating day 14 embryos. Genome-wide characteristics of open chromatin showed that ATAC-seq signals in both transcription start sites (TSS) and transcription end sites (TES) were strong. Additionally, the distal ATAC-seq peaks were enriched in repeat elements in a type-specific and stage-specific manner. We further unveiled a series of transcription factor (TF) motifs with distinct variation of enrichment from distal ATAC-seq peaks. By integrated analysis of chromatin accessibility with transcriptomes and DNA methylomes in bovine early embryos, we showed that promoter accessibility was positively correlated with gene expression, especially during major EGA, and was strongly correlated to DNA methylation and CpG density. Finally, we identified the critical chromatin signatures and TFs that differ between in vivo and in vitro derived blastocysts, which provides insights to the potential mechanisms leading to low quality of embryos produced in vitro. Together, this comprehensive analysis revealed critical features of chromatin landscape and epigenetic reprogramming during bovine preimplantation embryo development.

DNA methylation dynamics during zygotic genome activation in goat

DNA methylation is a crucial element in the epigenetic regulation of mammalian embryonic development. However, the subtle changes in DNA methylation differ in species, and, little information is known regarding the dynamics of DNA methylation at the single-base resolution in goat. In the present study, we studied the DNA methylation dynamics during goat zygotic genome activation (ZGA) at global and single-base resolution using immunostaining and reduced representation bisulfite sequencing, respectively. We showed that DNA methylation was decreased both at global and single-base resolution, and the expression of TET1 was increased while DNMT1 was decreased during ZGA in goat. We identified 51058 tiles of differential methylation regions (DMRs), which were enriched in the developmental process, the regulation of developmental process, AMPK signaling pathway, mTOR signaling pathway, autophagy, and lysosome, as revealed by GO and KEGG enrichment analysis. Furthermore, we found an association between the methylation level and the expression of imprinted genes (IGF2R, PEG3, and ZFP64), maternal genes (TRIM28, SETD1A, SIN3A, and NPM2), and zygotic genes (DUXA, IGF2BP1, WT1, and ZIM3), suggesting that DNA methylation is in the tight control of ZGA in goat by regulating the expression of the critical genes. Our data will help to understand the stochastic ZGA events to achieve better development of goat embryos in vitro and provide an excellent source for further ZGA studies.

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.

DNA methylation in the vertebrate germline: balancing memory and erasure

Cytosine methylation is a DNA modification that is critical for vertebrate development and provides a plastic yet stable information module in addition to the DNA code. DNA methylation memory establishment, maintenance and erasure is carefully balanced by molecular machinery highly conserved among vertebrates. In mammals, extensive erasure of epigenetic marks, including 5-methylcytosine (5mC), is a hallmark of early embryo and germline development. Conversely, global cytosine methylation patterns are preserved in at least some non-mammalian vertebrates over comparable developmental windows. The evolutionary mechanisms which drove this divergence are unknown, nevertheless a direct consequence of retaining epigenetic memory in the form of 5mC is the enhanced potential for transgenerational epigenetic inheritance (TEI). Given that DNA methylation dynamics remains underexplored in most vertebrate lineages, the extent of information transferred to offspring by epigenetic modification might be underestimated.

DNA methylation changes during preimplantation development reveal inter-species differences and reprogramming events at imprinted genes

Preimplantation embryos experience profound resetting of epigenetic information inherited from the gametes. Genome-wide analysis at single-base resolution has shown similarities but also species differences between human and mouse preimplantation embryos in DNA methylation patterns and reprogramming. Here, we have extended such analysis to two key livestock species, the pig and the cow. We generated genome-wide DNA methylation and whole-transcriptome datasets from gametes to blastocysts in both species. In oocytes from both species, a distinctive bimodal methylation landscape is present, with hypermethylated domains prevalent over hypomethylated domains, similar to human, while in the mouse the proportions are reversed.

An oocyte-like pattern of methylation persists in the cleavage stages, albeit with some reduction in methylation level, persisting to blastocysts in cow, while pig blastocysts have a highly hypomethylated landscape. In the pig, there was evidence of transient de novo methylation at the 8–16 cell stages of domains unmethylated in oocytes, revealing a complex dynamic of methylation reprogramming. The methylation datasets were used to identify germline differentially methylated regions (gDMRs) of known imprinted genes and for the basis of detection of novel imprinted loci. Strikingly in the pig, we detected a consistent reduction in gDMR methylation at the 8–16 cell stages, followed by recovery to the blastocyst stage, suggesting an active period of imprint stabilization in preimplantation embryos. Transcriptome analysis revealed absence of expression in oocytes of both species of ZFP57, a key factor in the mouse for gDMR methylation maintenance, but presence of the alternative imprint regulator ZNF445. In conclusion, our study reveals species differences in DNA methylation reprogramming and suggests that porcine or bovine models may be closer to human in key aspects than in the mouse model.

Metabolism-associated genome-wide epigenetic changes in bovine oocytes during early lactation

Dietary intake in early lactating cows is outmatched by milk production. These cows experience a negative energy balance, resulting in a distinct blood metabolism and poor reproductive function due to impaired ovulation and increased embryo loss. We hypothesize that oocytes from lactating cows undergoing transient metabolic stress exhibit a different epigenetic profile crucial for developmental competence. To investigate this, we collected oocytes from metabolically-profiled cows at early- and mid-postpartum stages and characterized their epigenetic landscape compared with control heifers using whole-genome bisulfite sequencing. Early-postpartum cows were metabolically deficient with a significantly lower energy balance and significantly higher concentrations of non-esterified fatty acids and beta-hydroxybutyrate than mid-postpartum animals and control heifers. Accordingly, 32,990 early-postpartum-specific differentially methylated regions (DMRs) were found in genes involved in metabolic pathways, carbon metabolism, and fatty acid metabolism, likely descriptive of the epigenetic regulation of metabolism in early-postpartum oocytes. DMRs found overlapping CpG islands and exons of imprinted genes such as MEST and GNAS in early-postpartum oocytes suggest that early lactation metabolic stress may affect imprint acquisition, which could explain the embryo loss. This whole-genome approach introduces potential candidate genes governing the link between metabolic stress and the reproductive outcome of oocytes.