A comparative analysis of chromatin accessibility in cattle, pig, and mouse tissues

A review of the role of epigenetic studies for intramuscular fat deposition in beef cattle

Intramuscular fat (IMF) deposition profoundly influences meat quality and economic value in beef cattle production. Meanwhile, contemporary developments in epigenetics have opened new outlooks for understanding the molecular basics of IMF regulation, and it has become a key area of research for world scholars. Therefore, the aim of this paper was to provide insight and synthesis into the intricate relationship between epigenetic mechanisms and IMF deposition in beef cattle. The methodology involves a thorough analysis of existing literature, including pertinent books, academic journals, and online resources, to provide a comprehensive overview of the role of epigenetic studies in IMF deposition in beef cattle. This review summarizes the contemporary studies in epigenetic mechanisms in IMF regulation, high-resolution epigenomic mapping, single-cell epigenomics, multi-omics integration, epigenome editing approaches, longitudinal studies in cattle growth, environmental epigenetics, machine learning in epigenetics, ethical and regulatory considerations, and translation to industry practices from perspectives of IMF deposition in beef cattle. Moreover, this paper highlights DNA methylation, histone modifications, acetylation, phosphorylation, ubiquitylation, non-coding RNAs, DNA hydroxymethylation, epigenetic readers, writers, and erasers, chromatin immunoprecipitation followed by sequencing, whole genome bisulfite sequencing, epigenome-wide association studies, and their profound impact on the expression of crucial genes governing adipogenesis and lipid metabolism. Nutrition and stress also have significant influences on epigenetic modifications and IMF deposition. The key findings underscore the pivotal role of epigenetic studies in understanding and enhancing IMF deposition in beef cattle, with implications for precision livestock farming and ethical livestock management. In conclusion, this review highlights the crucial significance of epigenetic pathways and environmental factors in affecting IMF deposition in beef cattle, providing insightful information for improving the economics and meat quality of cattle production.

Chromatin profiling and state predictions reveal insights into epigenetic regulation during early porcine development

BACKGROUND

Given their physiological similarities to humans, pigs are increasingly used as model organisms in human-oriented biomedical studies. Additionally, their value to animal agriculture across the globe has led to the development of numerous studies to investigate how to improve livestock welfare and production efficiency. As such, pigs are uniquely poised as compelling models that can yield findings with potential implications in both human and animal contexts. Despite this, many gaps remain in our knowledge about the foundational mechanisms that govern gene expression in swine across different developmental stages, particularly in early development. To address some of these gaps, we profiled the histone marks H3K4me3, H3K27ac, and H3K27me3 and the SWI/SNF central ATPase BRG1 in two porcine cell lines representing discrete early developmental time points and used the resulting information to construct predicted chromatin state maps for these cells. We combined this approach with analysis of publicly available RNA-seq data to examine the relationship between epigenetic status and gene expression in these cell types.

RESULTS

In porcine fetal fibroblast (PFF) and trophectoderm cells (PTr2), we saw expected patterns of enrichment for each of the profiled epigenetic features relative to specific genomic regions. H3K4me3 was primarily enriched at and around global gene promoters, H3K27ac was enriched in promoter and intergenic regions, H3K27me3 had broad stretches of enrichment across the genome and narrower enrichment patterns in and around the promoter regions of some genes, and BRG1 primarily had detectable enrichment at and around promoter regions and in intergenic stretches, with many instances of H3K27ac co-enrichment. We used this information to perform genome-wide chromatin state predictions for 10 different states using ChromHMM. Using the predicted chromatin state maps, we identified a subset of genomic regions marked by broad H3K4me3 enrichment, and annotation of these regions revealed that they were highly associated with essential developmental processes and consisted largely of expressed genes. We then compared the identities of the genes marked by these regions to genes identified as cell-type-specific using transcriptome data and saw that a subset of broad H3K4me3-marked genes was also specifically expressed in either PFF or PTr2 cells.

CONCLUSIONS

These findings enhance our understanding of the epigenetic landscape present in early swine development and provide insight into how variabilities in chromatin state are linked to cell identity. Furthermore, this data captures foundational epigenetic details in two valuable porcine cell lines and contributes to the growing body of knowledge surrounding the epigenetic landscape in this species.

Xenopus tropicalis osteoblast-specific open chromatin regions reveal promoters and enhancers involved in human skeletal phenotypes and shed light on early vertebrate evolution

While understanding the genetic underpinnings of osteogenesis has far-reaching implications for skeletal diseases and evolution, a comprehensive characterization of the osteoblastic regulatory landscape in non-mammalian vertebrates is still lacking. Here, we compared the ATAC-Seq profile of Xenopus tropicalis (Xt) osteoblasts to a variety of non mineralizing control tissues, and identified osteoblast-specific nucleosome free regions (NFRs) at 527 promoters and 6747 distal regions. Sequence analyses, Gene Ontology, RNA-Seq and ChIP-Seq against four key histone marks confirmed that the distal regions correspond to bona fide osteogenic transcriptional enhancers exhibiting a shared regulatory logic with mammals. We report 425 regulatory regions conserved with human and globally associated to skeletogenic genes. Of these, 35 regions have been shown to impact human skeletal phenotypes by GWAS, including one trps1 enhancer and the runx2 promoter, two genes which are respectively involved in trichorhinophalangeal syndrome type I and cleidocranial dysplasia. Intriguingly, 60 osteoblastic NFRs also align to the genome of the elephant shark, a species lacking osteoblasts and bone tissue. To tackle this paradox, we chose to focus on dlx5 because its conserved promoter, known to integrate regulatory inputs during mammalian osteogenesis, harbours an osteoblast-specific NFR in both frog and human. Hence, we show that dlx5 is expressed in Xt and elephant shark odontoblasts, supporting a common cellular and genetic origin of bone and dentine. Taken together, our work (i) unravels the Xt osteogenic regulatory landscape, (ii) illustrates how cross-species comparisons harvest data relevant to human biology and (iii) reveals that a set of genes including bnc2, dlx5, ebf3, mir199a, nfia, runx2 and zfhx4 drove the development of a primitive form of mineralized skeletal tissue deep in the vertebrate lineage.

Cross-Species Comparative DNA Methylation Reveals Novel Insights into Complex Trait Genetics among Cattle, Sheep, and Goats

The cross-species characterization of evolutionary changes in the functional genome can facilitate the translation of genetic findings across species and the interpretation of the evolutionary basis underlying complex phenotypes. Yet, this has not been fully explored between cattle, sheep, goats, and other mammals. Here, we systematically characterized the evolutionary dynamics of DNA methylation and gene expression in 3 somatic tissues (i.e. brain, liver, and skeletal muscle) and sperm across 7 mammalian species, including 3 ruminant livestock species (cattle, sheep, and goats), humans, pigs, mice, and dogs, by generating and integrating 160 DNA methylation and transcriptomic data sets. We demonstrate dynamic changes of DNA hypomethylated regions and hypermethylated regions in tissue-type manner across cattle, sheep, and goats. Specifically, based on the phylo-epigenetic model of DNA methylome, we identified a total of 25,074 hypomethylated region extension events specific to cattle, which participated in rewiring tissue-specific regulatory network. Furthermore, by integrating genome-wide association studies of 50 cattle traits, we provided novel insights into the genetic and evolutionary basis of complex phenotypes in cattle. Overall, our study provides a valuable resource for exploring the evolutionary dynamics of the functional genome and highlights the importance of cross-species characterization of multiomics data sets for the evolutionary interpretation of complex phenotypes in cattle livestock.

A perspective from the EU: unintended genetic changes in plants caused by NGT-their relevance for a comprehensive molecular characterisation and risk assessment

Several regions in the world are currently holding discussions in regard to the regulation of new genomic techniques (NGTs) and their application in agriculture. The European Commission, for instance, is proposing the introduction of specific regulation for NGT plants. Various questions need to be answered including e.g., the extent to which NGT-induced intended and unintended genetic modifications must be subjected to a mandatory risk assessment as part of an approval procedure. This review mostly focuses on findings in regard to unintended genetic changes that can be caused by the application of NGTs. More specifically, the review deals with the application of the nuclease CRISPR/Cas, which is currently the most important tool for developing NGT plants, and its potential to introduce double strand breaks (DSBs) at a targeted DNA sequence. For this purpose, we identified the differences in comparison to non-targeted mutagenesis methods used in conventional breeding. The review concludes that unintended genetic changes caused by NGT processes are relevant to risk assessment. Due to the technical characteristics of NGTs, the sites of the unintended changes, their genomic context and their frequency (in regard to specific sites) mean that the resulting gene combinations (intended or unintended) may be unlikely to occur with conventional methods. This, in turn, implies that the biological effects (phenotypes) can also be different and may cause risks to health and the environment. Therefore, we conclude that the assessment of intended as well as unintended genetic changes should be part of a mandatory comprehensive molecular characterisation and risk assessment of NGT plants that are meant for environmental releases or for market authorisation.

Characterization of transcriptional enhancers in the chicken genome using CRISPR-mediated activation

DNA regulatory elements intricately control when, where, and how genes are activated. Therefore, understanding the function of these elements could unveil the complexity of the genetic regulation network. Genome-wide significant variants are predominantly found in non-coding regions of DNA, so comprehending the predicted functional regulatory elements is crucial for understanding the biological context of these genomic markers, which can be incorporated into breeding programs. The emergence of CRISPR technology has provided a powerful tool for studying non-coding regulatory elements in genomes. In this study, we leveraged epigenetic data from the Functional Annotation of Animal Genomes project to identify promoter and putative enhancer regions associated with three genes (HBBA, IRF7, and PPARG) in the chicken genome. To identify the enhancer regions, we designed guide RNAs targeting the promoter and candidate enhancer regions and utilized CRISPR activation (CRISPRa) with dCas9-p300 and dCas9-VPR as transcriptional activators in chicken DF-1 cells. By comparing the expression levels of target genes between the promoter activation and the co-activation of the promoter and putative enhancers, we were able to identify functional enhancers that exhibited augmented upregulation. In conclusion, our findings demonstrate the remarkable efficiency of CRISPRa in precisely manipulating the expression of endogenous genes by targeting regulatory elements in the chicken genome, highlighting its potential for functional validation of non-coding regions.

An organism-wide ATAC-seq peak catalog for the bovine and its use to identify regulatory variants

We report the generation of an organism-wide catalog of 976,813 cis-acting regulatory elements for the bovine detected by the assay for transposase accessible chromatin using sequencing (ATAC-seq). We regroup these regulatory elements in 16 components by nonnegative matrix factorization. Correlation between the genome-wide density of peaks and transcription start sites, correlation between peak accessibility and expression of neighboring genes, and enrichment in transcription factor binding motifs support their regulatory potential. Using a previously established catalog of 12,736,643 variants, we show that the proportion of single-nucleotide polymorphisms mapping to ATAC-seq peaks is higher than expected and that this is owing to an approximately 1.3-fold higher mutation rate within peaks. Their site frequency spectrum indicates that variants in ATAC-seq peaks are subject to purifying selection. We generate eQTL data sets for liver and blood and show that variants that drive eQTL fall into liver- and blood-specific ATAC-seq peaks more often than expected by chance. We combine ATAC-seq and eQTL data to estimate that the proportion of regulatory variants mapping to ATAC-seq peaks is approximately one in three and that the proportion of variants mapping to ATAC-seq peaks that are regulatory is approximately one in 25. We discuss the implication of these findings on the utility of ATAC-seq information to improve the accuracy of genomic selection.

Improving the annotation of the cattle genome by annotating transcription start sites in a diverse set of tissues and populations using Cap Analysis Gene Expression sequencing

Understanding the genomic control of tissue-specific gene expression and regulation can help to inform the application of genomic technologies in farm animal breeding programs. The fine mapping of promoters [transcription start sites (TSS)] and enhancers (divergent amplifying segments of the genome local to TSS) in different populations of cattle across a wide diversity of tissues provides information to locate and understand the genomic drivers of breed- and tissue-specific characteristics. To this aim, we used Cap Analysis Gene Expression (CAGE) sequencing, of 24 different tissues from 3 populations of cattle, to define TSS and their coexpressed short-range enhancers (<1 kb) in the ARS-UCD1.2_Btau5.0.1Y reference genome (1000bulls run9) and analyzed tissue and population specificity of expressed promoters. We identified 51,295 TSS and 2,328 TSS-Enhancer regions shared across the 3 populations (dairy, beef-dairy cross, and Canadian Kinsella composite cattle from 2 individuals, 1 of each sex, per population). Cross-species comparative analysis of CAGE data from 7 other species, including sheep, revealed a set of TSS and TSS-Enhancers that were specific to cattle. The CAGE data set will be combined with other transcriptomic information for the same tissues to create a new high-resolution map of transcript diversity across tissues and populations in cattle for the BovReg project. Here we provide the CAGE data set and annotation tracks for TSS and TSS-Enhancers in the cattle genome. This new annotation information will improve our understanding of the drivers of gene expression and regulation in cattle and help to inform the application of genomic technologies in breeding programs.

Profiling of Chromatin Accessibility in Pigs across Multiple Tissues and Developmental Stages

The study of chromatin accessibility across tissues and developmental stages is essential for elucidating the transcriptional regulation of various phenotypes and biological processes. However, the chromatin accessibility profiles of multiple tissues in newborn pigs and across porcine liver development remain poorly investigated. Here, we used ATAC-seq and rRNA-depleted RNA-seq to profile open chromatin maps and transcriptional features of heart, kidney, liver, lung, skeletal muscle, and spleen in newborn pigs and porcine liver tissue in the suckling and adult stages, respectively. Specifically, by analyzing a union set of protein-coding genes (PCGs) and two types of transcripts (lncRNAs and TUCPs), we obtained a comprehensive annotation of consensus ATAC-seq peaks for each tissue and developmental stage. As expected, the PCGs with tissue-specific accessible promoters had active transcription and were relevant to tissue-specific functions. In addition, other non-coding tissue-specific peaks were involved in both physical activity and the morphogenesis of neonatal tissues. We also characterized stage-specific peaks and observed a close association between dynamic chromatin accessibility and hepatic function transition during liver postnatal development. Overall, this study expands our current understanding of epigenetic regulation in mammalian tissues and organ development, which can benefit both economic trait improvement and improve the biomedical usage of pigs.

Micrococcal nuclease sequencing of porcine sperm suggests enriched co-location between retained histones and genomic regions related to semen quality and early embryo development

The mammalian spermatozoon has a unique chromatin structure in which the majority of histones are replaced by protamines during spermatogenesis and a small fraction of nucleosomes are retained at specific locations of the genome. The sperm's chromatin structure remains unresolved in most animal species, including the pig. However, mapping the genomic locations of retained nucleosomes in sperm could help understanding the molecular basis of both sperm development and function as well as embryo development. This information could then be useful to identify molecular markers for sperm quality and fertility traits. Here, micrococcal nuclease digestion coupled with high throughput sequencing was performed on pig sperm to map the genomic location of mono- and sub-nucleosomal chromatin fractions in relation to a set of diverse functional elements of the genome, some of which were related to semen quality and early embryogenesis. In particular, the investigated elements were promoters, the different sections of the gene body, coding and non-coding RNAs present in the pig sperm, potential transcription factor binding sites, genomic regions associated to semen quality traits and repeat elements. The analysis yielded 25,293 and 4,239 peaks in the mono- and sub-nucleosomal fractions, covering 0.3% and 0.02% of the porcine genome, respectively. A cross-species comparison revealed positional conservation of the nucleosome retention in sperm between the pig data and a human dataset that found nucleosome enrichment in genomic regions of importance in development. Both gene ontology analysis of the genes mapping nearby the mono-nucleosomal peaks and the identification of putative transcription factor binding motifs within the mono- and the sub- nucleosomal peaks showed enrichment for processes related to sperm function and embryo development. There was significant motif enrichment for Znf263, which in humans was suggested to be a key regulator of genes with paternal preferential expression during early embryogenesis. Moreover, enriched positional intersection was found in the genome between the mono-nucleosomal peaks and both the RNAs present in pig sperm and the RNAs related to sperm quality. There was no co-location between GWAS hits for semen quality in swine and the nucleosomal sites. Finally, the data evidenced depletion of mono-nucleosomes in long interspersed nuclear elements and enrichment of sub-nucleosomes in short interspersed repeat elements.These results suggest that retained nucleosomes in sperm could both mark regulatory elements or genes expressed during spermatogenesis linked to semen quality and fertility and act as transcriptional guides during early embryogenesis. The results of this study support the undertaking of ambitious research using a larger number of samples to robustly assess the positional relationship between histone retention in sperm and the reproductive ability of boars.

Characterization of Chromatin Accessibility in Fetal Bovine Chondrocytes

Despite significant advances of the bovine epigenome investigation, new evidence for the epigenetic basis of fetal cartilage development remains lacking. In this study, the chondrocytes were isolated from long bone tissues of bovine fetuses at 90 days. The Assay for Transposase-Accessible Chromatin with high throughput sequencing (ATAC-seq) and transcriptome sequencing (RNA-seq) were used to characterize gene expression and chromatin accessibility profile in bovine chondrocytes. A total of 9686 open chromatin regions in bovine fetal chondrocytes were identified and 45% of the peaks were enriched in the promoter regions. Then, all peaks were annotated to the nearest gene for Gene Ontology (GO) and Kyoto Encylopaedia of Genes and Genomes (KEGG) analysis. Growth and development-related processes such as amide biosynthesis process (GO: 0043604) and translation regulation (GO: 006417) were enriched in the GO analysis. The KEGG analysis enriched endoplasmic reticulum protein processing signal pathway, TGF-β signaling pathway and cell cycle pathway, which are closely related to protein synthesis and processing during cell proliferation. Active transcription factors (TFs) were enriched by ATAC-seq, and were fully verified with gene expression levels obtained by RNA-seq. Among the top50 TFs from footprint analysis, known or potential cartilage development-related transcription factors FOS, FOSL2 and NFY were found. Overall, our data provide a theoretical basis for further determining the regulatory mechanism of cartilage development in bovine.

Profiling the immune epigenome across global cattle breeds

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Relating enhancer genetic variation across mammals to complex phenotypes using machine learning

Protein-coding differences between species often fail to explain phenotypic diversity, suggesting the involvement of genomic elements that regulate gene expression such as enhancers. Identifying associations between enhancers and phenotypes is challenging because enhancer activity can be tissue-dependent and functionally conserved despite low sequence conservation. We developed the Tissue-Aware Conservation Inference Toolkit (TACIT) to associate candidate enhancers with species' phenotypes using predictions from machine learning models trained on specific tissues. Applying TACIT to associate motor cortex and parvalbumin-positive interneuron enhancers with neurological phenotypes revealed dozens of enhancer-phenotype associations, including brain size-associated enhancers that interact with genes implicated in microcephaly or macrocephaly. TACIT provides a foundation for identifying enhancers associated with the evolution of any convergently evolved phenotype in any large group of species with aligned genomes.

Functional annotation of the animal genomes: An integrated annotation resource for the horse

The genomic sequence of the horse has been available since 2009, providing critical resources for discovering important genomic variants regarding both animal health and population structures. However, to fully understand the functional implications of these variants, detailed annotation of the horse genome is required. Due to the limited availability of functional data for the equine genome, as well as the technical limitations of short-read RNA-seq, existing annotation of the equine genome contains limited information about important aspects of gene regulation, such as alternate isoforms and regulatory elements, which are either not transcribed or transcribed at a very low level. To solve above problems, the Functional Annotation of the Animal Genomes (FAANG) project proposed a systemic approach to tissue collection, phenotyping, and data generation, adopting the blueprint laid out by the Encyclopedia of DNA Elements (ENCODE) project. Here we detail the first comprehensive overview of gene expression and regulation in the horse, presenting 39,625 novel transcripts, 84,613 candidate cis-regulatory elements (CRE) and their target genes, 332,115 open chromatin regions genome wide across a diverse set of tissues. We showed substantial concordance between chromatin accessibility, chromatin states in different genic features and gene expression. This comprehensive and expanded set of genomics resources will provide the equine research community ample opportunities for studies of complex traits in the horse.

Chicken chromatin accessibility atlas accelerates epigenetic annotation of birds and gene fine-mapping associated with growth traits

The development of epigenetic maps, such as the ENCODE project in humans, provides resources for gene regulation studies and a reference for research of disease-related regulatory elements. However, epigenetic information, such as a bird-specific chromatin accessibility atlas, is currently lacking for the thousands of bird species currently described. The major genomic difference between birds and mammals is their shorter introns and intergenic distances, which seriously hinders the use of humans and mice as a reference for studying the function of important regulatory regions in birds. In this study, using chicken as a model bird species, we systematically compiled a chicken chromatin accessibility atlas using 53 Assay of Transposase Accessible Chromatin sequencing (ATAC-seq) samples across 11 tissues. An average of 50 796 open chromatin regions were identified per sample, cumulatively accounting for 20.36% of the chicken genome. Tissue specificity was largely reflected by differences in intergenic and intronic peaks, with specific functional regulation achieved by two mechanisms: recruitment of several sequence-specific transcription factors and direct regulation of adjacent functional genes. By integrating data from genome-wide association studies, our results suggest that chicken body weight is driven by different regulatory variants active in growth-relevant tissues. We propose CAB39L (active in the duodenum), RCBTB1 (muscle and liver), and novel long non-coding RNA ENSGALG00000053256 (bone) as candidate genes regulating chicken body weight. Overall, this study demonstrates the value of epigenetic data in fine-mapping functional variants and provides a compendium of resources for further research on the epigenetics and evolution of birds and mammals.

Butyrate Induces Modifications of the CTCF-Binding Landscape in Cattle Cells

Butyrate is produced in the rumen from microbial fermentation and is related to several functions, including cell differentiation and proliferation. Butyrate supplementation in calves can accelerate rumen development. DNA-protein interactions, such as the CCCTC-binding factor (CTCF), play essential roles in chromatin organization and gene expression regulation. Although CTCF-binding sites have been identified recently in cattle, a deeper characterization, including differentially CTCF-binding sites (DCBS), is vital for a better understanding of butyrate’s role in the chromatin landscape. This study aimed to identify CTCF-binding regions and DCBS under a butyrate-induced condition using ChIP-seq in bovine cells; 61,915 CTCF peaks were identified in the butyrate and 51,347 in the control. From these regions, 2265 DCBS were obtained for the butyrate vs. control comparison, comprising ~90% of induced sites. Most of the butyrate DCBS were in distal intergenic regions, showing a potential role as insulators. Gene ontology enrichment showed crucial terms for the induced DCBS, mainly related to cellular proliferation, cell adhesion, and growth regulation. Interestingly, the ECM-receptor interaction pathway was observed for the induced DCBS. Motif enrichment analysis further identified transcription factors, including CTCF, BORIS, TGIF2, and ZIC3. When DCBS was integrated with RNA-seq data, putative genes were identified for the repressed DCBS, including GATA4. Our study revealed promising candidate genes in bovine cells by a butyrate-induced condition that might be related to the regulation of rumen development, such as integrins, keratins, and collagens. These results provide a better understanding of the function of butyrate in cattle rumen development and chromatin landscape regulation.

Differentially CTCF-Binding Sites in Cattle Rumen Tissue during Weaning

The weaning transition in calves is characterized by major structural changes such as an increase in the rumen capacity and surface area due to diet changes. Studies evaluating rumen development in calves are vital to identify genetic mechanisms affected by weaning. This study aimed to provide a genome-wide characterization of CTCF-binding sites and differentially CTCF-binding sites (DCBS) in rumen tissue during the weaning transition of four Holstein calves to uncover regulatory elements in rumen epithelial tissue using ChIP-seq. Our study generated 67,280 CTCF peaks for the before weaning (BW) and 39,891 for after weaning (AW). Then, 7401 DCBS were identified for the AW vs. BW comparison representing 0.15% of the cattle genome, comprising ~54% of induced DCBS and ~46% of repressed DCBS. Most of the induced and repressed DCBS were in distal intergenic regions, showing a potential role as insulators. Gene ontology enrichment revealed many shared GO terms for the induced and the repressed DCBS, mainly related to cellular migration, proliferation, growth, differentiation, cellular adhesion, digestive tract morphogenesis, and response to TGFβ. In addition, shared KEGG pathways were obtained for adherens junction and focal adhesion. Interestingly, other relevant KEGG pathways were observed for the induced DCBS like gastric acid secretion, salivary secretion, bacterial invasion of epithelial cells, apelin signaling, and mucin-type O-glycan biosynthesis. IPA analysis further revealed pathways with potential roles in rumen development during weaning, including TGFβ, Integrin-linked kinase, and Integrin signaling. When DCBS were further integrated with RNA-seq data, 36 putative target genes were identified for the repressed DCBS, including KRT84, COL9A2, MATN3, TSPAN1, and AJM1. This study successfully identified DCBS in cattle rumen tissue after weaning on a genome-wide scale and revealed several candidate target genes that may have a role in rumen development, such as TGFβ, integrins, keratins, and SMADs. The information generated in this preliminary study provides new insights into bovine genome regulation and chromatin landscape.

Inferring mammalian tissue-specific regulatory conservation by predicting tissue-specific differences in open chromatin

BACKGROUND

Evolutionary conservation is an invaluable tool for inferring functional significance in the genome, including regions that are crucial across many species and those that have undergone convergent evolution. Computational methods to test for sequence conservation are dominated by algorithms that examine the ability of one or more nucleotides to align across large evolutionary distances. While these nucleotide alignment-based approaches have proven powerful for protein-coding genes and some non-coding elements, they fail to capture conservation of many enhancers, distal regulatory elements that control spatial and temporal patterns of gene expression. The function of enhancers is governed by a complex, often tissue- and cell type-specific code that links combinations of transcription factor binding sites and other regulation-related sequence patterns to regulatory activity. Thus, function of orthologous enhancer regions can be conserved across large evolutionary distances, even when nucleotide turnover is high.

RESULTS

We present a new machine learning-based approach for evaluating enhancer conservation that leverages the combinatorial sequence code of enhancer activity rather than relying on the alignment of individual nucleotides. We first train a convolutional neural network model that can predict tissue-specific open chromatin, a proxy for enhancer activity, across mammals. Next, we apply that model to distinguish instances where the genome sequence would predict conserved function versus a loss of regulatory activity in that tissue. We present criteria for systematically evaluating model performance for this task and use them to demonstrate that our models accurately predict tissue-specific conservation and divergence in open chromatin between primate and rodent species, vastly out-performing leading nucleotide alignment-based approaches. We then apply our models to predict open chromatin at orthologs of brain and liver open chromatin regions across hundreds of mammals and find that brain enhancers associated with neuron activity have a stronger tendency than the general population to have predicted lineage-specific open chromatin.

CONCLUSION

The framework presented here provides a mechanism to annotate tissue-specific regulatory function across hundreds of genomes and to study enhancer evolution using predicted regulatory differences rather than nucleotide-level conservation measurements.

Characterization of Accessible Chromatin Regions in Cattle Rumen Epithelial Tissue during Weaning

Weaning in ruminants is characterized by the transition from a milk-based diet to a solid diet, which drives a critical gastrointestinal tract transformation. Understanding the regulatory control of this transformation during weaning can help to identify strategies to improve rumen health. This study aimed to identify regions of accessible chromatin in rumen epithelial tissue in pre- and post-weaning calves and investigate differentially accessible regions (DARs) to uncover regulatory elements in cattle rumen development using the ATAC-seq approach. A total of 126,071 peaks were identified, covering 1.15% of the cattle genome. From these accessible regions, 2766 DARs were discovered. Gene ontology enrichment resulted in GO terms related to the cell adhesion, anchoring junction, growth, cell migration, motility, and morphogenesis. In addition, putative regulatory canonical pathways were identified (TGFβ, integrin-linked kinase, integrin signaling, and regulation of the epithelial–mesenchymal transition). Canonical pathways integrated with co-expression results showed that TGFβ and ILK signaling pathways play essential roles in rumen development through the regulation of cellular adhesions. In this study, DARs during weaning were identified, revealing enhancers, transcription factors, and candidate target genes that represent potential biomarkers for the bovine rumen development, which will serve as a molecular tool for rumen development studies.

Profiling of open chromatin in developing pig (Sus scrofa) muscle to identify regulatory regions

There is very little information about how the genome is regulated in domestic pigs (Sus scrofa). This lack of knowledge hinders efforts to define and predict the effects of genetic variants in pig breeding programs. To address this knowledge gap, we need to identify regulatory sequences in the pig genome starting with regions of open chromatin. We used the "Improved Protocol for the Assay for Transposase-Accessible Chromatin (Omni-ATAC-Seq)" to identify putative regulatory regions in flash-frozen semitendinosus muscle from 24 male piglets. We collected samples from the smallest-, average-, and largest-sized male piglets from each litter through five developmental time points. Of the 4661 ATAC-Seq peaks identified that represent regions of open chromatin, >50% were within 1 kb of known transcription start sites. Differential read count analysis revealed 377 ATAC-Seq defined genomic regions where chromatin accessibility differed significantly across developmental time points. We found regions of open chromatin associated with downregulation of genes involved in muscle development that were present in small-sized fetal piglets but absent in large-sized fetal piglets at day 90 of gestation. The dataset that we have generated provides a resource for studies of genome regulation in pigs and contributes valuable functional annotation information to filter genetic variants for use in genomic selection in pig breeding programs.

Decoding the Equine Genome: Lessons from ENCODE

The horse reference genome assemblies, EquCab2.0 and EquCab3.0, have enabled great advancements in the equine genomics field, from tools to novel discoveries. However, significant gaps of knowledge regarding genome function remain, hindering the study of complex traits in horses. In an effort to address these gaps and with inspiration from the Encyclopedia of DNA Elements (ENCODE) project, the equine Functional Annotation of Animal Genome (FAANG) initiative was proposed to bridge the gap between genome and gene expression, providing further insights into functional regulation within the horse genome. Three years after launching the initiative, the equine FAANG group has generated data from more than 400 experiments using over 50 tissues, targeting a variety of regulatory features of the equine genome. In this review, we examine how valuable lessons learned from the ENCODE project informed our decisions in the equine FAANG project. We report the current state of the equine FAANG project and discuss how FAANG can serve as a template for future expansion of functional annotation in the equine genome and be used as a reference for studies of complex traits in horse. A well-annotated reference functional atlas will also help advance equine genetics in the pan-genome and precision medicine era.

Pig genome functional annotation enhances the biological interpretation of complex traits and human disease

The functional annotation of livestock genomes is crucial for understanding the molecular mechanisms that underpin complex traits of economic importance, adaptive evolution and comparative genomics. Here, we provide the most comprehensive catalogue to date of regulatory elements in the pig (Sus scrofa) by integrating 223 epigenomic and transcriptomic data sets, representing 14 biologically important tissues. We systematically describe the dynamic epigenetic landscape across tissues by functionally annotating 15 different chromatin states and defining their tissue-specific regulatory activities. We demonstrate that genomic variants associated with complex traits and adaptive evolution in pig are significantly enriched in active promoters and enhancers. Furthermore, we reveal distinct tissue-specific regulatory selection between Asian and European pig domestication processes. Compared with human and mouse epigenomes, we show that porcine regulatory elements are more conserved in DNA sequence, under both rapid and slow evolution, than those under neutral evolution across pig, mouse, and human. Finally, we provide biological insights on tissue-specific regulatory conservation, and by integrating 47 human genome-wide association studies, we demonstrate that, depending on the traits, mouse or pig might be more appropriate biomedical models for different complex traits and diseases.

Chromatin accessibility and regulatory vocabulary across indicine cattle tissues

BACKGROUND

Spatiotemporal changes in the chromatin accessibility landscape are essential to cell differentiation, development, health, and disease. The quest of identifying regulatory elements in open chromatin regions across different tissues and developmental stages is led by large international collaborative efforts mostly focusing on model organisms, such as ENCODE. Recently, the Functional Annotation of Animal Genomes (FAANG) has been established to unravel the regulatory elements in non-model organisms, including cattle. Now, we can transition from prediction to validation by experimentally identifying the regulatory elements in tropical indicine cattle. The identification of regulatory elements, their annotation and comparison with the taurine counterpart, holds high promise to link regulatory regions to adaptability traits and improve animal productivity and welfare.

RESULTS

We generate open chromatin profiles for liver, muscle, and hypothalamus of indicine cattle through ATAC-seq. Using robust methods for motif discovery, motif enrichment and transcription factor binding sites, we identify potential master regulators of the epigenomic profile in these three tissues, namely HNF4, MEF2, and SOX factors, respectively. Integration with transcriptomic data allows us to confirm some of their target genes. Finally, by comparing our results with Bos taurus data we identify potential indicine-specific open chromatin regions and overlaps with indicine selective sweeps.

CONCLUSIONS

Our findings provide insights into the identification and analysis of regulatory elements in non-model organisms, the evolution of regulatory elements within two cattle subspecies as well as having an immediate impact on the animal genetics community in particular for a relevant productive species such as tropical cattle.

Opportunity to improve livestock traits using 3D genomics

The advent of high-throughput chromosome conformation capture and sequencing (Hi-C) has enabled researchers to probe the 3D architecture of the mammalian genome in a genome-wide manner. Simultaneously, advances in epigenomic assays, such as chromatin immunoprecipitation and sequencing (ChIP-seq) and DNase-seq, have enabled researchers to study cis-regulatory interactions and chromatin accessibility across the same genome-wide scale. The use of these data has revealed many unique insights into gene regulation and disease pathomechanisms in several model organisms. With the advent of these high-throughput sequencing technologies, there has been an ever-increasing number of datasets available for study; however, this is often limited to model organisms. Livestock species play critical roles in the economies of developing and developed nations alike. Despite this, they are greatly underrepresented in the 3D genomics space; Hi-C and related technologies have the potential to revolutionise livestock breeding by enabling a more comprehensive understanding of how production traits are controlled. The growth in human and model organism Hi-C data has seen a surge in the availability of computational tools for use in 3D genomics, with some tools using machine learning techniques to predict features and improve dataset quality. In this review, we provide an overview of the 3D genome and discuss the status of 3D genomics in livestock before delving into advancing the field by drawing inspiration from research in human and mouse. We end by offering future directions for livestock research in the field of 3D genomics.

Evidence for and localization of proposed causative variants in cattle and pig genomes

BACKGROUND

This paper reviews the localization of published potential causative variants in contemporary pig and cattle reference genomes, and the evidence for their causality. In spite of the difficulties inherent to the identification of causative variants from genetic mapping and genome-wide association studies, researchers in animal genetics have proposed putative causative variants for several traits relevant to livestock breeding.

RESULTS

For this review, we read the literature that supports potential causative variants in 13 genes (ABCG2, DGAT1, GHR, IGF2, MC4R, MSTN, NR6A1, PHGK1, PRKAG3, PLRL, RYR1, SYNGR2 and VRTN) in cattle and pigs, and localized them in contemporary reference genomes. We review the evidence for their causality, by aiming to separate the evidence for the locus, the proposed causative gene and the proposed causative variant, and report the bioinformatic searches and tactics needed to localize the sequence variants in the cattle or pig genome.

CONCLUSIONS

Taken together, there is usually good evidence for the association at the locus level, some evidence for a specific causative gene at eight of the loci, and some experimental evidence for a specific causative variant at six of the loci. We recommend that researchers who report new potential causative variants use referenced coordinate systems, show local sequence context, and submit variants to repositories.

Tissue Resources for the Functional Annotation of Animal Genomes

In order to generate an atlas of the functional elements driving genome expression in domestic animals, the Functional Annotation of Animal Genome (FAANG) strategy was to sample many tissues from a few animals of different species, sexes, ages, and production stages. This article presents the collection of tissue samples for four species produced by two pilot projects, at INRAE (National Research Institute for Agriculture, Food and Environment) and the University of California, Davis. There were three mammals (cattle, goat, and pig) and one bird (chicken). It describes the metadata characterizing these reference sets (1) for animals with origin and selection history, physiological status, and environmental conditions; (2) for samples with collection site and tissue/cell processing; (3) for quality control; and (4) for storage and further distribution. Three sets are identified: set 1 comprises tissues for which collection can be standardized and for which representative aliquots can be easily distributed (liver, spleen, lung, heart, fat depot, skin, muscle, and peripheral blood mononuclear cells); set 2 comprises tissues requiring special protocols because of their cellular heterogeneity (brain, digestive tract, secretory organs, gonads and gametes, reproductive tract, immune tissues, cartilage); set 3 comprises specific cell preparations (immune cells, tracheal epithelial cells). Dedicated sampling protocols were established and uploaded in https://data.faang.org/protocol/samples. Specificities between mammals and chicken are described when relevant. A total of 73 different tissues or tissue sections were collected, and 21 are common to the four species. Having a common set of tissues will facilitate the transfer of knowledge within and between species and will contribute to decrease animal experimentation. Combining data on the same samples will facilitate data integration. Quality control was performed on some tissues with RNA extraction and RNA quality control. More than 5,000 samples have been stored with unique identifiers, and more than 4,000 were uploaded onto the Biosamples database, provided that standard ontologies were available to describe the sample. Many tissues have already been used to implement FAANG assays, with published results. All samples are available without restriction for further assays. The requesting procedure is described. Members of FAANG are encouraged to apply a range of molecular assays to characterize the functional status of collected samples and share their results, in line with the FAIR (Findable, Accessible, Interoperable, and Reusable) data principles.

Successful ATAC-Seq From Snap-Frozen Equine Tissues

An assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) has become an increasingly popular method to assess genome-wide chromatin accessibility in isolated nuclei from fresh tissues. However, many biobanks contain only snap-frozen tissue samples. While ATAC-seq has been applied to frozen brain tissues in human, its applicability in a wide variety of tissues in horse remains unclear. The Functional Annotation of Animal Genome (FAANG) project is an international collaboration aimed to provide high quality functional annotation of animal genomes. The equine FAANG initiative has generated a biobank of over 80 tissues from two reference female animals and experiments to begin to characterize tissue specificity of genome function for prioritized tissues have been performed. Due to the logistics of tissue collection and storage, extracting nuclei from a large number of tissues for ATAC-seq at the time of collection is not always practical. To assess the feasibility of using stored frozen tissues for ATAC-seq and to provide a guideline for the equine FAANG project, we compared ATAC-seq results from nuclei isolated from frozen tissue to cryopreserved nuclei (CN) isolated at the time of tissue harvest in liver, a highly cellular homogenous tissue, and lamina, a relatively acellular tissue unique to the horse. We identified 20,000-33,000 accessible chromatin regions in lamina and 22-61,000 in liver, with consistently more peaks identified using CN isolated at time of tissue collection. Our results suggest that frozen tissues are an acceptable substitute when CN are not available. For more challenging tissues such as lamina, nuclei extraction at the time of tissue collection is still preferred for optimal results. Therefore, tissue type and accessibility to intact nuclei should be considered when designing ATAC-seq experiments.

The landscape of chromatin accessibility in skeletal muscle during embryonic development in pigs

BACKGROUND

The development of skeletal muscle in pigs during the embryonic stage is precisely regulated by transcriptional mechanisms, which depend on chromatin accessibility. However, how chromatin accessibility plays a regulatory role during embryonic skeletal muscle development in pigs has not been reported. To gain insight into the landscape of chromatin accessibility and the associated genome-wide transcriptome during embryonic muscle development, we performed ATAC-seq and RNA-seq analyses of skeletal muscle from pig embryos at 45, 70 and 100 days post coitus (dpc).

RESULTS

In total, 21,638, 35,447 and 60,181 unique regions (or peaks) were found across the embryos at 45 dpc (LW45), 70 dpc (LW70) and 100 dpc (LW100), respectively. More than 91% of the peaks were annotated within - 1 kb to 100 bp of transcription start sites (TSSs). First, widespread increases in specific accessible chromatin regions (ACRs) from embryos at 45 to 100 dpc suggested that the regulatory mechanisms became increasingly complicated during embryonic development. Second, the findings from integrated ATAC-seq and RNA-seq analyses showed that not only the numbers but also the intensities of ACRs could control the expression of associated genes. Moreover, the motif screening of stage-specific ACRs revealed some transcription factors that regulate muscle development-related genes, such as MyoG, Mef2c, and Mef2d. Several potential transcriptional repressors, including E2F6, OTX2 and CTCF, were identified among the genes that exhibited different regulation trends between the ATAC-seq and RNA-seq data.

CONCLUSIONS

This work indicates that chromatin accessibility plays an important regulatory role in the embryonic muscle development of pigs and regulates the temporal and spatial expression patterns of key genes in muscle development by influencing the binding of transcription factors. Our results contribute to a better understanding of the regulatory dynamics of genes involved in pig embryonic skeletal muscle development.

Functional annotations of three domestic animal genomes provide vital resources for comparative and agricultural research

Gene regulatory elements are central drivers of phenotypic variation and thus of critical importance towards understanding the genetics of complex traits. The Functional Annotation of Animal Genomes consortium was formed to collaboratively annotate the functional elements in animal genomes, starting with domesticated animals. Here we present an expansive collection of datasets from eight diverse tissues in three important agricultural species: chicken (Gallus gallus), pig (Sus scrofa), and cattle (Bos taurus). Comparative analysis of these datasets and those from the human and mouse Encyclopedia of DNA Elements projects reveal that a core set of regulatory elements are functionally conserved independent of divergence between species, and that tissue-specific transcription factor occupancy at regulatory elements and their predicted target genes are also conserved. These datasets represent a unique opportunity for the emerging field of comparative epigenomics, as well as the agricultural research community, including species that are globally important food resources.

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.