Machine learning to identify endometrial biomarkers predictive of pregnancy success following artificial insemination in dairy cows.

The objective was to identify a set of genes whose transcript abundance is predictive of a cow's ability to become pregnant following artificial insemination (AI). Endometrial epithelial cells from the uterine body were collected for RNA sequencing using the cytobrush method from 193 first-service Holstein cows at estrus prior to AI (day 0). A group of 253 first-service cows not used for cytobrush collection were controls. There was no effect of cytobrush collection on pregnancy outcomes at day 30 or 70 or on pregnancy loss between day 30 and 70. There were 2 upregulated and 214 downregulated genes (FDR < 0.05, absolute fold change >2-fold) for cows pregnant at day 30 versus those that were not pregnant. Functional terms overrepresented in the downregulated genes included those related to immune and inflammatory responses. Machine learning for fertility biomarkers with the R package BORUTA resulted in identification of 57 biomarkers that predicted pregnancy outcome at day 30 with an average accuracy of 77%. Thus, machine learning can identify predictive biomarkers of pregnancy in endometrium with high accuracy. Moreover, sampling of endometrial epithelium using the cytobrush can help understand functional characteristics of the endometrium at AI without compromising cow fertility. Functional characteristics of the genes comprising the set of biomarkers is indicative that a major determinant of cow fertility, at least for first insemination after calving, is immune status of the uterus, which, in turn, is likely to reflect the previous history of uterine disease.

Effect of Holstein genotype on ex-vivo interleukin-1ß response to lipopolysaccharide (LPS), lipoteichoic acid (LTA) and heat-killed Gram-negative and Gram-positive bacteria.

Effects of Holstein genotype on interleukin-1ß response were assessed by ex-vivo stimulation of whole blood with lipopolysaccharide (LPS), lipoteichoic acid (LTA), or sonicated, heat-killed Gram-negative or Gram-positive bacteria. Holstein genotypes were unselected Holsteins (UH, n = 14) not subjected to selection pressures since the mid-1960s and contemporary Holsteins (CH, n = 13). Milk yield of UH and CH cows differ by more than 4500 kg/lactation. Whole blood was mixed with 0.01 µg LPS, 10 µg LTA or 2.5 × 106 CFU of sonicated, heat-killed E. coli, K. pneumoniae, S. marcescens, S. aureus, S. dysgalactiae, or S. uberis per mL of blood and incubated (4 h, 37 °C). Plasma IL-1ß was quantified by ELISA and log10-transformed concentrations analyzed with a multivariate linear mixed effects model. Responses to bacteria were greater than responses to LPS or LTA. Responses to LPS, LTA and the Gram-negative stimulants were greater in UH than in CH cows while responses to the Gram-positive bacteria did not differ between Holstein genotypes. In both genotypes, strong correlations were detected among IL-1ß responses to the Gram-negative stimulants and to LTA. There were strong correlations among IL-1ß responses to the Gram-positive bacteria in CH cows but only between S. aureus and S. dysgalactiae in UH cows. The IL-1ß response to S. uberis was highly correlated with responses to all of the Gram-negative stimulants in CH cows but only with E. coli in the UH cows. The reduced immune response could make contemporary cows more susceptible to infection by Gram-negative bacteria. Results confirm selection practices since the mid-1960s have altered immune response in the Holstein, at least to Gram-negative bacteria, and validate the need for additional studies to further evaluate the impacts of these selection practices on immune function in contemporary Holsteins.

Truncation of IFT80 causes early embryonic loss in Holstein cattle associated with Holstein haplotype 2.

Recessive alleles represent genetic risk in populations that have undergone bottleneck events. We present a comprehensive framework for identification and validation of these genetic defects, including haplotype-based detection, variant selection from sequence data, and validation using knockout embryos. Holstein haplotype 2 (HH2), which causes embryonic death, was used to demonstrate the approach. Holstein haplotype 2 was identified using a deficiency-of-homozygotes approach and confirmed to negatively affect conception rate and stillbirths. Five carriers were present in a group of 183 sequenced Holstein bulls selected to maximize the coverage of unique haplotypes. Three variants concordant with haplotype calls were found in HH2: a high-priority frameshift mutation resulting, and 2 low-priority variants (1 synonymous variant, 1 premature stop codon). The frameshift in intraflagellar 80 (IFT80) was confirmed in a separate group of Holsteins from the 1000 Bull Genomes Project that shared no animals with the discovery set. IFT80-null embryos were generated by truncating the IFT80 transcript at exon 2 or 11 using a CRISPR-Cas9 system. Abattoir-derived oocytes were fertilized in vitro, and zygotes were injected at the one-cell stage either with a guide RNA and CAS9 mRNA complex (n = 100) or Cas9 mRNA (control, n = 100) before return to culture, and replicated 3 times. IFT80 is activated at the 8-cell stage, and IFT80-null embryos arrested at this stage of development, which is consistent with data from mouse hypomorphs and HH2 carrier-to-carrier matings. This frameshift in IFT80 on chromosome 1 at 107,172,615 bp (p.Leu381fs) disrupts WNT and hedgehog signaling, and is responsible for the death of homozygous embryos.

Genomic characterization of autozygosity and recent inbreeding trends in all major breeds of US dairy cattle.

Maintaining a genetically diverse dairy cattle population is critical to preserving adaptability to future breeding goals and avoiding declines in fitness. This study characterized the genomic landscape of autozygosity and assessed trends in genetic diversity in 5 breeds of US dairy cattle. We analyzed a sizable genomic data set containing 4,173,679 pedigreed and genotyped animals of the Ayrshire, Brown Swiss, Guernsey, Holstein, and Jersey breeds. Runs of homozygosity (ROH) of 2 Mb or longer in length were identified in each animal. The within-breed means for number and the combined length of ROH were highest in Jerseys (62.66 ± 8.29 ROH and 426.24 ± 83.40 Mb, respectively; mean ± SD) and lowest in Ayrshires (37.24 ± 8.27 ROH and 265.05 ± 85.00 Mb, respectively). Short ROH were the most abundant, but moderate to large ROH made up the largest proportion of genome autozygosity in all breeds. In addition, we identified ROH islands in each breed. This revealed selection patterns for milk production, productive life, health, and reproduction in most breeds and evidence for parallel selective pressure for loci on chromosome 6 between Ayrshire and Brown Swiss and for loci on chromosome 20 between Holstein and Jersey. We calculated inbreeding coefficients using 3 different approaches, pedigree-based (FPED), marker-based using a genomic relationship matrix (FGRM), and segment-based using ROH (FROH). The average inbreeding coefficient ranged from 0.06 in Ayrshires and Brown Swiss to 0.08 in Jerseys and Holsteins using FPED, from 0.22 in Holsteins to 0.29 in Guernsey and Jerseys using FGRM, and from 0.11 in Ayrshires to 0.17 in Jerseys using FROH. In addition, the effective population size at past generations (5-100 generations ago), the yearly rate of inbreeding, and the effective population size in 3 recent periods (2000-2009, 2010-2014, and 2015-2018) were determined in each breed to ascertain current and historical trends of genetic diversity. We found a historical trend of decreasing effective population size in the last 100 generations in all breeds and breed differences in the effect of the recent implementation of genomic selection on inbreeding accumulation.

A multi-tissue atlas of regulatory variants in cattle.

Characterization of genetic regulatory variants acting on livestock gene expression is essential for interpreting the molecular mechanisms underlying traits of economic value and for increasing the rate of genetic gain through artificial selection. Here we build a Cattle Genotype-Tissue Expression atlas (CattleGTEx) as part of the pilot phase of the Farm animal GTEx (FarmGTEx) project for the research community based on 7,180 publicly available RNA-sequencing (RNA-seq) samples. We describe the transcriptomic landscape of more than 100 tissues/cell types and report hundreds of thousands of genetic associations with gene expression and alternative splicing for 23 distinct tissues. We evaluate the tissue-sharing patterns of these genetic regulatory effects, and functionally annotate them using multiomics data. Finally, we link gene expression in different tissues to 43 economically important traits using both transcriptome-wide association and colocalization analyses to decipher the molecular regulatory mechanisms underpinning such agronomic traits in cattle.

Comparative transcriptome in large-scale human and cattle populations.

BACKGROUND: Cross-species comparison of transcriptomes is important for elucidating evolutionary molecular mechanisms underpinning phenotypic variation between and within species, yet to date it has been essentially limited to model organisms with relatively small sample sizes. RESULTS: Here, we systematically analyze and compare 10,830 and 4866 publicly available RNA-seq samples in humans and cattle, respectively, representing 20 common tissues. Focusing on 17,315 orthologous genes, we demonstrate that mean/median gene expression, inter-individual variation of expression, expression quantitative trait loci, and gene co-expression networks are generally conserved between humans and cattle. By examining large-scale genome-wide association studies for 46 human traits (average n = 327,973) and 45 cattle traits (average n = 24,635), we reveal that the heritability of complex traits in both species is significantly more enriched in transcriptionally conserved than diverged genes across tissues. CONCLUSIONS: In summary, our study provides a comprehensive comparison of transcriptomes between humans and cattle, which might help decipher the genetic and evolutionary basis of complex traits in both species.

Effect of Holstein genotype on ex-vivo cytokine response to lipopolysaccharide (LPS) and lipoteichoic acid (LTA) during the periparturient period.

Effects of Holstein genotype on innate immune response were assessed with ex-vivo lipopolysaccharide (LPS) and lipoteichoic acid (LTA) stimulation of whole blood from unselected (UH, n = 10) and contemporary (CH, n = 11) Holsteins that differ in production by more than 4,500 kg/lactation. Blood was collected at -14, 7, 28, and 49 days in milk (DIM), mixed with a pathogen-associated molecular pattern (PAMP) molecule (0.01 or 1.0 µg LPS or 10 or 100 µg LTA per mL blood) and incubated (4 h, 37 °C). Plasma cytokines were quantified by ELISA, log10-transformed and analyzed by repeated measures with DIM as the repeated effect. Cytokine responses increased with PAMP dose and decreased as DIM increased. There was a genotype by LPS dose interaction for IL-1ß as response to the low dose was greater in UH but did not differ between genotypes for the high dose. The IL-1ß response was greater while the IL-6 response to LTA tended to be greater in UH than in CH cows. The more negative energy balance of CH cows did not impact genotype difference in cytokine responses. Results indicate selection since the mid-1960s has decreased ex-vivo, whole blood cytokine response of CH cows to LPS and to LTA.

Genome-wide recombination map construction from single sperm sequencing in cattle.

BACKGROUND: Meiotic recombination is one of the important phenomena contributing to gamete genome diversity. However, except for human and a few model organisms, it is not well studied in livestock, including cattle. RESULTS: To investigate their distributions in the cattle sperm genome, we sequenced 143 single sperms from two Holstein bulls. We mapped meiotic recombination events at high resolution based on phased heterozygous single nucleotide polymorphism (SNP). In the absence of evolutionary selection pressure in fertilization and survival, recombination events in sperm are enriched near distal chromosomal ends, revealing that such a pattern is intrinsic to the molecular mechanism of meiosis. Furthermore, we further validated these findings in single sperms with results derived from sequencing its family trio of diploid genomes and our previous studies of recombination in cattle. CONCLUSIONS: To our knowledge, this is the first large-scale single sperm whole-genome sequencing effort in livestock, which provided useful information for future studies of recombination, genome instability, and male infertility.

Functional annotation of regulatory elements in cattle genome reveals the roles of extracellular interaction and dynamic change of chromatin states in rumen development during weaning.

We profiled landscapes of bovine regulatory elements and explored dynamic changes of chromatin states in rumen development during weaning. The regulatory elements (15 chromatin states) and their coordinated activities in cattle were defined through genome-wide profiling of four histone modifications, CTCF-binding, DNA accessibility, DNA methylation, and transcriptome in rumen epithelial tissues. Each chromatin state presented specific enrichment for sequence ontology, methylation, trait-associated variants, transcription, gene expression-associated variants, selection signatures, and evolutionarily conserved elements. During weaning, weak enhancers and flanking active transcriptional start sites (TSS) were the most dynamic chromatin states and occurred in tandem with significant variations in gene expression and DNA methylation, significantly associated with stature, production, and reproduction economic traits. By comparing with in vitro cultured epithelial cells and in vivo rumen tissues, we showed the commonness and uniqueness of these results, especially the roles of cell interactions and mitochondrial activities in tissue development.

Effect of Temperature and Maternal Age on Recombination Rate in Cattle.

Meiotic recombination is a fundamental biological process that facilitates meiotic division and promotes genetic diversity. Recombination is phenotypically plastic and affected by both intrinsic and extrinsic factors. The effect of maternal age on recombination rates has been characterized in a wide range of species, but the effect's direction remains inconclusive. Additionally, the characterization of temperature effects on recombination has been limited to model organisms. Here we seek to comprehensively determine the impact of genetic and environmental factors on recombination rate in dairy cattle. Using a large cattle pedigree, we identified maternal recombination events within 305,545 three-generation families. By comparing recombination rate between parents of different ages, we found a quadratic trend between maternal age and recombination rate in cattle. In contrast to either an increasing or decreasing trend in humans, cattle recombination rate decreased with maternal age until 65 months and then increased afterward. Combining recombination data with temperature information from public databases, we found a positive correlation between environmental temperature during fetal development of offspring and recombination rate in female parents. Finally, we fitted a full recombination rate model on all related factors, including genetics, maternal age, and environmental temperatures. Based on the final model, we confirmed the effect of maternal age and environmental temperature during fetal development of offspring on recombination rate with an estimated heritability of 10% (SE = 0.03) in cattle. Collectively, we characterized the maternal age and temperature effects on recombination rate and suggested the adaptation of meiotic recombination to environmental stimuli in cattle. Our results provided first-hand information regarding the plastic nature of meiotic recombination in a mammalian species.

Trends in genetic diversity and the effect of inbreeding in American Angus cattle under genomic selection.

BACKGROUND: While the adoption of genomic evaluations in livestock has increased genetic gain rates, its effects on genetic diversity and accumulation of inbreeding have raised concerns in cattle populations. Increased inbreeding may affect fitness and decrease the mean performance for economically important traits, such as fertility and growth in beef cattle, with the age of inbreeding having a possible effect on the magnitude of inbreeding depression. The purpose of this study was to determine changes in genetic diversity as a result of the implementation of genomic selection in Angus cattle and quantify potential inbreeding depression effects of total pedigree and genomic inbreeding, and also to investigate the impact of recent and ancient inbreeding. RESULTS: We found that the yearly rate of inbreeding accumulation remained similar in sires and decreased significantly in dams since the implementation of genomic selection. Other measures such as effective population size and the effective number of chromosome segments show little evidence of a detrimental effect of using genomic selection strategies on the genetic diversity of beef cattle. We also quantified pedigree and genomic inbreeding depression for fertility and growth. While inbreeding did not affect fertility, an increase in pedigree or genomic inbreeding was associated with decreased birth weight, weaning weight, and post-weaning gain in both sexes. We also measured the impact of the age of inbreeding and found that recent inbreeding had a larger depressive effect on growth than ancient inbreeding. CONCLUSIONS: In this study, we sought to quantify and understand the possible consequences of genomic selection on the genetic diversity of American Angus cattle. In both sires and dams, we found that, generally, genomic selection resulted in decreased rates of pedigree and genomic inbreeding accumulation and increased or sustained effective population sizes and number of independently segregating chromosome segments. We also found significant depressive effects of inbreeding accumulation on economically important growth traits, particularly with genomic and recent inbreeding.

Single-cell transcriptomic analyses of dairy cattle ruminal epithelial cells during weaning.

Using the 10× Genomics Chromium Controller, we obtained scRNA-seq data of 5064 and 1372 individual cells from two Holstein calf ruminal epithelial tissues before and after weaning, respectively. We detected six distinct cell clusters, designated their cell types, and reported their marker genes. We then examined these clusters' underlining cell types and relationships by performing cell cycle, pseudotime trajectory, regulatory network, weighted gene co-expression network and gene ontology analyses. By integrating these cell marker genes with Holstein GWAS signals, we found they were enriched for animal production and body conformation traits. Finally, we confirmed their cell identities by comparing them with human and mouse stomach epithelial cells. This study presents an initial effort to implement single-cell transcriptomic analysis in cattle, and demonstrates ruminal tissue epithelial cell types and their developments during weaning, opening the door for new discoveries about tissue/cell type roles in complex traits at single-cell resolution.

Effects of Incorporating Dry Matter Intake and Residual Feed Intake into a Selection Index for Dairy Cattle Using Deterministic Modeling.

The inclusion of feed efficiency in the breeding goal for dairy cattle has been discussed for many years. The effects of incorporating feed efficiency into a selection index were assessed by indirect selection (dry matter intake) and direct selection (residual feed intake) using deterministic modeling. Both traits were investigated in three ways: (1) restricting the trait genetic gain to zero, (2) applying negative selection pressure, and (3) applying positive selection pressure. Changes in response to selection from economic and genetic gain perspectives were used to evaluate the impact of including feed efficiency with direct or indirect selection in an index. Improving feed efficiency through direct selection on residual feed intake was the best scenario analyzed, with the highest overall economic response including favorable responses to selection for production and feed efficiency. Over time, the response to selection is cumulative, with the potential for animals to reduce consumption by 0.16 kg to 2.7 kg of dry matter per day while maintaining production. As the selection pressure increased on residual feed intake, the response to selection for production, health, and fertility traits and body condition score became increasingly less favorable. This work provides insight into the potential long-term effects of selecting for feed efficiency as residual feed intake.

Comparison of Gene Editing Versus Conventional Breeding to Introgress the POLLED Allele Into the Tropically Adapted Australian Beef Cattle Population.

Dehorning is the process of physically removing horns to protect animals and humans from injury, but the process is costly, unpleasant, and faces increasing public scrutiny. Genetic selection for polled (hornless), which is genetically dominant to horned, is a long-term solution to eliminate the need for dehorning. However, due to the limited number of polled Australian Brahman bulls, the northern Australian beef cattle population remains predominantly horned. The potential to use gene editing to produce high-genetic-merit polled cattle was recently demonstrated. To further explore the concept, this study simulated introgression of the POLLED allele into a tropically adapted Australian beef cattle population via conventional breeding or gene editing (top 1% or 10% of seedstock bulls/year) for 3 polled mating schemes and compared results to baseline selection on genetic merit (Japan Ox selection index, $JapOx) alone, over the course of 20 years. The baseline scenario did not significantly decrease the 20-year HORNED allele frequency (80%), but resulted in one of the fastest rates of genetic gain ($8.00/year). Compared to the baseline, the conventional breeding scenarios where polled bulls were preferentially used for breeding, regardless of their genetic merit, significantly decreased the 20-year HORNED allele frequency (30%), but resulted in a significantly slower rate of genetic gain ($6.70/year, P ≤ 0.05). The mating scheme that required the exclusive use of homozygous polled bulls, resulted in the lowest 20-year HORNED allele frequency (8%), but this conventional breeding scenario resulted in the slowest rate of genetic gain ($5.50/year). The addition of gene editing the top 1% or 10% of seedstock bull calves/year to each conventional breeding scenario resulted in significantly faster rates of genetic gain (up to $8.10/year, P ≤ 0.05). Overall, our study demonstrates that, due to the limited number of polled Australian Brahman bulls, strong selection pressure on polled will be necessary to meaningfully increase the number of polled animals in this population. Moreover, these scenarios illustrate how gene editing could be a tool for accelerating the development of high-genetic-merit homozygous polled sires to mitigate the current trade-off of slower genetic gain associated with decreasing HORNED allele frequency in the Australian Brahman population.

Invited review: The future of selection decisions and breeding programs: What are we breeding for, and who decides?

Genetic selection has been a very successful tool for the long-term improvement of livestock populations, and the rapid adoption of genomic selection over the last decade has doubled the rate of gain in some populations. Breeding programs seek to identify genetically superior parents of the next generation, typically as a function of an index that combines information about many economically important traits into a single number. In the United States, the data that drive this system are collected through the national dairy herd improvement program that began more than a century ago. The resulting information about animal performance, pedigree, and genotype is used to compute genomic evaluations for comparing and ranking animals for selection. However, the full expression of genetic potential requires that animals are placed in environments that can support such performance. The Agricultural Research Service of the US Department of Agriculture and the Council on Dairy Cattle Breeding collaborate to deliver state-of-the-art genomic evaluations to the dairy industry. Today, most breeding stock are selected and marketed using the net merit dollars (NM$) selection index, which evolved from 2 traits in 1926 (milk and fat yield) to a combination of 36 individual traits following the last NM$ update in 2018. Updates to NM$ require the estimation of many different values, and it can be difficult to achieve consensus from stakeholders on what should be added to, or removed from, the index at each review, and how those traits should be weighted. Over time, the majority of the emphasis in the index has shifted from yield traits to fertility, health, and fitness traits. Phenotypes for some of these new traits are difficult or expensive to measure, or require changes to on-farm habits that have not been widely adopted. This is driving interest in sensor-based systems that provide continuous measurements of the farm environment, individual animal performance, and detailed milk composition. There is also a need to capture more detailed data about the environment in which animals perform, including information about feeding, housing, milking systems, and infectious and parasitic load. However, many challenges accompany these new technologies, including a lack of standardization or validation, need for high-speed internet connections, increased computational requirements, and interpretations that are often not backed by direct observations of biological phenomena. This work will describe how US selection objectives are developed, as well as discuss opportunities and challenges associated with new technologies for measuring and recording animal performance.

An alternative interpretation of residual feed intake by phenotypic recursive relationships in dairy cattle.

There has been increasing interest in residual feed intake (RFI) as a measure of net feed efficiency in dairy cattle. Residual feed intake phenotypes are obtained as residuals from linear regression encompassing relevant factors (i.e., energy sinks) to account for body tissue mobilization. By rearranging the single-trait linear regression, we showed a causal RFI interpretation underlying the linear regression for RFI. It postulates recursive effects in energy allocation from energy sinks on dry matter intake, but the feedback or simultaneous effects are nonexistent. A Bayesian recursive structural equation model was proposed for directly predicting RFI and energy sinks and estimating relevant genetic parameters simultaneously. A simplified Markov chain Monte Carlo algorithm was described. The recursive model is asymptotically equivalent to one-step linear regression for RFI, yet extends the analytical capacity to multiple-trait analysis.

Epigenomics and genotype-phenotype association analyses reveal conserved genetic architecture of complex traits in cattle and human.

BACKGROUND: Lack of comprehensive functional annotations across a wide range of tissues and cell types severely hinders the biological interpretations of phenotypic variation, adaptive evolution, and domestication in livestock. Here we used a combination of comparative epigenomics, genome-wide association study (GWAS), and selection signature analysis, to shed light on potential adaptive evolution in cattle. RESULTS: We cross-mapped 8 histone marks of 1300 samples from human to cattle, covering 178 unique tissues/cell types. By uniformly analyzing 723 RNA-seq and 40 whole genome bisulfite sequencing (WGBS) datasets in cattle, we validated that cross-mapped histone marks captured tissue-specific expression and methylation, reflecting tissue-relevant biology. Through integrating cross-mapped tissue-specific histone marks with large-scale GWAS and selection signature results, we for the first time detected relevant tissues and cell types for 45 economically important traits and artificial selection in cattle. For instance, immune tissues are significantly associated with health and reproduction traits, multiple tissues for milk production and body conformation traits (reflecting their highly polygenic architecture), and thyroid for the different selection between beef and dairy cattle. Similarly, we detected relevant tissues for 58 complex traits and diseases in humans and observed that immune and fertility traits in humans significantly correlated with those in cattle in terms of relevant tissues, which facilitated the identification of causal genes for such traits. For instance, PIK3CG, a gene highly specifically expressed in mononuclear cells, was significantly associated with both age-at-menopause in human and daughter-still-birth in cattle. ICAM, a T cell-specific gene, was significantly associated with both allergic diseases in human and metritis in cattle. CONCLUSION: Collectively, our results highlighted that comparative epigenomics in conjunction with GWAS and selection signature analyses could provide biological insights into the phenotypic variation and adaptive evolution. Cattle may serve as a model for human complex traits, by providing additional information beyond laboratory model organisms, particularly when more novel phenotypes become available in the near future.

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.