Genome-wide analysis of DNA methylation patterns in horse

Different expression patterns of DNA methyltransferases during horse testis development

DNA methyltransferases (DNMTs) are important epigenetic modification during spermatogenesis. To further evaluate the pattern of DNMTs in horse testes during development, we investigated the expression and localization of DNMT1, DNMT3a and DNMT3b at different time points. The qRT-PCR results showed that DNMT1 expression was maintained in testes tissue from 6-month-old (0.5y) to 2-year-old (2y) of age and decreased after 3-year-old (3y) (P < 0.01). The expression levels of DNMT3a and DNMT3b peaked in testes tissue at 3y (P < 0.01). At 4-year-old (4y), the expression of DNMT3a and DNMT3b was decreased and became similar to that at 0.5y. Immunofluorescence of DNMT1, DNMT3a and DNMT3b on testis samples confirmed the differential expression and localization of these three DNA methylation transferases during horse development. Further molecular biological studies are needed to understand the implications of the expression patterns of these DNMTs in horse testes.

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

Introduction

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

Methods

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

Results

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

Conclusion

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

Dynamics of the Equine Placental DNA Methylome and Transcriptome from Mid- to Late Gestation

The placenta is a temporary organ that is essential for the survival of the fetus, with a lifelong effect on the health of both the offspring and the dam. The functions of the placenta are controlled by its dynamic gene expression during gestation. In this study, we aimed to investigate the equine placental DNA methylome as one of the fundamental mechanisms that controls the gene expression dynamic. Chorioallantois samples from four (4M), six (6M), and ten (10M) months of gestation were used to map the methylation pattern of the placenta. Globally, methylation levels increased toward the end of gestation. We identified 921 differentially methylated regions (DMRs) between 4M and 6M, 1225 DMRs between 4M and 10M, and 1026 DMRs between 6M and 10M. A total of 817 genes carried DMRs comparing 4M and 6M, 978 comparing 4M and 10M, and 804 comparing 6M and 10M. We compared the transcriptomes between the samples and found 1381 differentially expressed genes (DEGs) when comparing 4M and 6M, 1428 DEGs between 4M and 10M, and 741 DEGs between 6M and 10M. Finally, we overlapped the DEGs and genes carrying DMRs (DMRs-DEGs). Genes exhibiting (a) higher expression, low methylation and (b) low expression, high methylation at different time points were identified. The majority of these DMRs-DEGs were located in introns (48.4%), promoters (25.8%), and exons (17.7%) and were involved in changes in the extracellular matrix; regulation of epithelial cell migration; vascularization; and regulation of minerals, glucose, and metabolites, among other factors. Overall, this is the first report highlighting the dynamics in the equine placenta methylome during normal pregnancy. The findings presented serve as a foundation for future studies on the impact of abnormal methylation on the outcomes of equine pregnancies.

Integrating genome-wide methylation and transcriptome-wide analyses to reveal the genetic mechanism of milk traits in Kazakh horses

Horse Milk has important quantitative characteristics and high economic value. However, the DNA methylation regulators involved in horse milk traits have not been clarified. To explore the important role of genome-wide DNA methylation in regulating equine milk yield, this study systematically investigated the genome-wide DNA methylation profiles of Kazakh horse blood by comparing a high-production group (HP, average daily milk yield of 7.5 kg) and low-production group (LP, average daily milk yield of 3.2 kg) using deep whole-genome bisulfite sequencing. First, both groups showed similar proportions of methylation at CpG sites. Subsequently, we identified 26,677 differential methylated regions (DMRs) of CG, 15 DMRs of CHG, 480 DMRs of CHH and 8268 DMR-related genes (DMGs). GO and KEGG analyses revealed that some DMGs were involved in regulating milk and milk component formation. By combining the WGBS-seq and the previous RNA-seq data, a total of 94 overlapping genes were obtained. Finally, we found that 9 DMGs are likely involved in milk production by Kazakh horses.

The Innovative Informatics Approaches of High-Throughput Technologies in Livestock: Spearheading the Sustainability and Resiliency of Agrigenomics Research

For more than a decade, next-generation sequencing (NGS) has been emerging as the mainstay of agrigenomics research. High-throughput technologies have made it feasible to facilitate research at the scale and cost required for using this data in livestock research. Scale frameworks of sequencing for agricultural and livestock improvement, management, and conservation are partly attributable to innovative informatics methodologies and advancements in sequencing practices. Genome-wide sequence-based investigations are often conducted worldwide, and several databases have been created to discover the connections between worldwide scientific accomplishments. Such studies are beginning to provide revolutionary insights into a new era of genomic prediction and selection capabilities of various domesticated livestock species. In this concise review, we provide selected examples of the current state of sequencing methods, many of which are already being used in animal genomic studies, and summarize the state of the positive attributes of genome-based research for cattle (Bos taurus), sheep (Ovis aries), pigs (Sus scrofa domesticus), horses (Equus caballus), chickens (Gallus gallus domesticus), and ducks (Anas platyrhyncos). This review also emphasizes the advantageous features of sequencing technologies in monitoring and detecting infectious zoonotic diseases. In the coming years, the continued advancement of sequencing technologies in livestock agrigenomics will significantly influence the sustained momentum toward regulatory approaches that encourage innovation to ensure continued access to a safe, abundant, and affordable food supplies for future generations.

Dynamic changes of genomic methylation profiles at different growth stages in Chinese Tan sheep

Background

Tan sheep, an important local sheep breed in China, is famous for their fur quality. One-month-old Tan sheep have white, curly hair with beautiful flower spikes, commonly known as “nine bends”, which has high economic value. However, the “nine bends” characteristic gradually disappears with age; consequently, the economic value of the Tan sheep decreases. Age-related changes in DNA methylation have been reported and may be responsible for age-induced changes in gene expression. Until now, no genome-wide surveys have been conducted to identify potential DNA methylation sites involved in different sheep growth stages. In this study we investigated the dynamic changes of genome-wide DNA methylation profiles in Tan sheep using DNA from skin and deep whole-genome bisulfite sequencing, and compared the DNA methylation levels at three different growth stages: 1, 24, and 48 months old (mon1, mon24, and mon48, respectively).

Results

In this study, 11 skin samples from three growth stages (four for mon1, four for mon24, and three for mon48) were used for DNA methylation analysis and gene expression profiling. There were 52, 288 and 236 differentially methylated genes (DMGs) identified between mon1 and mon24, mon1 and mon48, and mon24 and mon48, respectively. Of the differentially methylated regions, 1.11%, 7.61%, and 7.65% were in the promoter in mon1 vs. mon24, mon24 vs. mon48, and mon1 vs. mon48, respectively. DMGs were enriched in the MAPK and WNT signaling pathways, which are related to age growth and hair follicle morphogenesis processes. There were 51 DMGs associated with age growth and curly fleece formation. Four DMGs between mon1 and mon48 (KRT71, CD44, ROR2 and ZDHHC13) were further validated by bisulfite sequencing.

Conclusions

This study revealed dynamic changes in the genomic methylation profiles of mon1, mon24, and mon48 sheep, and the percentages of methylated cytosines were 3.38%, 2.85% and 4.17%, respectively. Of the DMGs, KRT71 and CD44 were highly methylated in mon1, and ROR2 and ZDHHC13 were highly methylated in mon48. These findings provide foundational information that may be used to develop strategies for potentially retaining the lamb fur and thus improving the economic value of Tan sheep.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-021-00632-9.

Differential Methylation and Transcriptome Integration Analysis Identified Differential Methylation Annotation Genes and Functional Research Related to Hair Follicle Development in Sheep

Hair follicle growth and development are a complex and long-term physiological process, which is regulated by a variety of physical factors and signal pathways. Increasing the understanding of the epigenetic regulation and function of candidate genes related to hair follicle development will help to better understand the molecular regulatory mechanisms of hair follicle development. In this study, the methylated DNA immunoprecipitation sequencing (MeDIP-seq) was used to obtain the genome-wide methylation map of the hair follicular development of Super Merino sheep in six stages (fetal skin tissue at 65d, 85d, 105d, 135d, 7d, and 30d after birth). Combined with the results of previous RNA-sequencing, 65 genes were screened out that were both differential methylation and differential expression, including EDN1, LAMC2, NR1D1, RORB, MyOZ3, and WNT2 gene. Differential methylation genes were enriched in Wnt, TNF, TGF-beta, and other signaling pathways related to hair follicle development. The bisulfite sequencing PCR results and MeDIP-seq were basically consistent, indicating that the sequencing results were accurate. As a key gene in the Wnt signaling pathway, both differential methylation and expression gene identified by MeDIP-seq and RNA-seq, further exploration of the function of WNT2 gene revealed that the DNA methylation of exon 5 (CpG11 site) promoted the expression of WNT2 gene. The overexpression vector of lentivirus pLEX-MCS-WNT2 was constructed, and WNT2 gene effectively promoted the proliferation of sheep skin fibroblasts. The results showed that WNT2 gene could promote the growth and development of skin and hair follicles. The results of this study will provide a theoretical basis for further research on sheep hair follicle development and gene regulation mechanisms.

DNA methylation and histone acetylation are involved in Wnt10b expression during the secondary hair follicle cycle in Angora rabbits

Secondary hair follicles (SHFs) in the Angora rabbit exhibit classic cyclic hair development, but the multiple molecular signals involved in hair cycling are yet to be explored in detail. In the present study, we investigated the expression pattern, methylation and histone H3 acetylation status of Wnt10b, as a molecular signal participating in hair cycling, during the SHF cycle in the Angora rabbit. Expression of Wnt10b at the anagen phase was significantly higher than that at both the telogen and catagen phases, suggesting that Wnt10b might serve as a critical activator during cyclic transition of SHFs. Methylation frequency of the fifth CpG site (CpG5-175 bp) in CpG islands at the anagen phase was lower than that at both the catagen and telogen phases. The methylation status of the CpG5 site was negatively correlated with Wnt10b expression. This indicated that the methylation of CpG5 might participate in Wnt10b transcriptional suppression in SHFs. Furthermore, histone H3 acetylation status in the regions-256~-11 bp and 98 ~ 361 bp were significantly lower at both the catagen and telogen phases than at the anagen phase. The histone H3 acetylation level was significantly positively correlated with Wnt10b expression. This confirmed that histone acetylation was likely involved in upregulating Wnt10b transcription in SHFs. Additionally, potential binding to the transcription factors ZF57 and HDBP was predicted within the CpG5 site. In conclusion, our findings reveal the epigenetic mechanism of Wnt10b transcription and provide a new insight into epigenetic regulation during the SHF cycle in the Angora rabbit.

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

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

Next-Generation Sequencing in Equine Genomics

The sequencing and assembly of a reference genome for the horse has been revolutionary for investigation of horse health and performance. Next-generation sequencing (NGS) methods represent a second revolution in equine genomics. Researchers can align and compare DNA and RNA sequencing data to the reference genome to explore variation that may contribute or be attributed to disease. NGS has also facilitated the translation of research discovery to clinically relevant applications. This article discusses the history and development of NGS, details some of the available sequencing platforms, and describes currently available applications in the context of both discovery and clinical settings.

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.

Genome-Wide DNA Methylation and RNA Analysis Reveal Potential Mechanism of Resistance to Streptococcus agalactiae in GIFT Strain of Nile Tilapia (Oreochromis niloticus )

Streptococcus agalactiae is an important pathogenic bacterium causing great economic loss in Nile tilapia (Oreochromis niloticus) culture. Resistant and susceptible groups sharing the same genome showed significantly different resistance to S. agalactiae in the genetically improved farmed tilapia strain of Nile tilapia. The resistance mechanism is unclear. We determined genome-wide DNA methylation profiles in spleen of resistant and susceptible O. niloticus at 5 h postinfection with S. agalactiae using whole-genome bisulfite sequencing. The methylation status was higher in the spleen samples from resistant fish than in the susceptible group. A total of 10,177 differentially methylated regions were identified in the two groups, including 3725 differentially methylated genes (DMGs) (3129 hyper-DMGs and 596 hypo-DMGs). The RNA sequencing showed 2374 differentially expressed genes (DEGs), including 1483 upregulated and 891 downregulated. Integrated analysis showed 337 overlapping DEGs and DMGs and 82 overlapping DEGs and differentially methylated region promoters. By integrating promoter DNA methylation with gene expression, we revealed four immune-related genes (Arnt2, Nhr38, Pcdh10, and Ccdc158) as key factors in epigenetic mechanisms contributing to pathogen resistance. Our study provided systematic methylome maps to explore the epigenetic mechanism and reveal the methylation loci of pathogen resistance and identified methylation-regulated genes that are potentially involved in defense against pathogens.

piRNA-63076 contributes to pulmonary arterial smooth muscle cell proliferation through acyl-CoA dehydrogenase

Piwi-interacting RNAs (piRNAs) are thought to be germline-specific and to be involved in maintaining genome stability during development. Recently, piRNA expression has been identified in somatic cells in diverse organisms. However, the roles of piRNAs in pulmonary arterial smooth muscle cell (PASMC) proliferation and the molecular mechanism underlying the hypoxia-regulated pathological process of pulmonary hypertension are not well understood. Using hypoxic animal models, cell and molecular biology, we obtained the first evidence that the expression of piRNA-63076 was up-regulated in hypoxia and was positively correlated with cell proliferation. Subsequently, we showed that acyl-CoA dehydrogenase (Acadm), which is negatively regulated by piRNA-63076 and interacts with Piwi proteins, was involved in hypoxic PASMC proliferation. Finally, Acadm inhibition under hypoxia was partly attributed to DNA methylation of the Acadm promoter region mediated by piRNA-63076. Overall, these findings represent invaluable resources for better understanding the role of epigenetics in pulmonary hypertension associated with piRNAs.

Comprehensive genome and transcriptome analyses reveal genetic relationship, selection signature, and transcriptome landscape of small-sized Korean native Jeju horse

The Jeju horse, indigenous to the Jeju Island in Korea may have originated from Mongolian horses. Adaptations to the local harsh environment have conferred Jeju horse with unique traits such as small-sized body, stocky head, and shorter limbs. These characteristics have not been studied previously at the genomic level. Therefore, we sequenced and compared the genome of 41 horses belonging to 6 breeds. We identified numerous breed-specific non-synonymous SNPs and loss-of-function mutants. Demographic and admixture analyses showed that, though Jeju horse is genetically the closest to the Mongolian breeds, its genetic ancestry is independent of that of the Mongolian breeds. Genome wide selection signature analysis revealed that genes such as LCORL, MSTN, HMGA2, ZFAT, LASP1, PDK4, and ACTN2, were positively selected in the Jeju horse. RNAseq analysis showed that several of these genes were also differentially expressed in Jeju horse compared to Thoroughbred horse. Comparative muscle fiber analysis showed that, the type I muscle fibre content was substantially higher in Jeju horse compared to Thoroughbred horse. Our results provide insights about the selection of complex phenotypic traits in the small-sized Jeju horse and the novel SNPs identified will aid in designing high-density SNP chip for studying other native horse breeds.

Ten years of the horse reference genome: insights into equine biology, domestication and population dynamics in the post-genome era

The horse reference genome from the Thoroughbred mare Twilight has been available for a decade and, together with advances in genomics technologies, has led to unparalleled developments in equine genomics. At the core of this progress is the continuing improvement of the quality, contiguity and completeness of the reference genome, and its functional annotation. Recent achievements include the release of the next version of the reference genome (EquCab3.0) and generation of a reference sequence for the Y chromosome. Horse satellite‐free centromeres provide unique models for mammalian centromere research. Despite extremely low genetic diversity of the Y chromosome, it has been possible to trace patrilines of breeds and pedigrees and show that Y variation was lost in the past approximately 2300 years owing to selective breeding. The high‐quality reference genome has led to the development of three different SNP arrays and WGSs of almost 2000 modern individual horses. The collection of WGS of hundreds of ancient horses is unique and not available for any other domestic species. These tools and resources have led to global population studies dissecting the natural history of the species and genetic makeup and ancestry of modern breeds. Most importantly, the available tools and resources, together with the discovery of functional elements, are dissecting molecular causes of a growing number of Mendelian and complex traits. The improved understanding of molecular underpinnings of various traits continues to benefit the health and performance of the horse whereas also serving as a model for complex disease across species.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

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

Background

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

Results

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

Conclusions

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

Functional Annotation of Animal Genomes (FAANG): Current Achievements and Roadmap

Functional annotation of genomes is a prerequisite for contemporary basic and applied genomic research, yet farmed animal genomics is deficient in such annotation. To address this, the FAANG (Functional Annotation of Animal Genomes) Consortium is producing genome-wide data sets on RNA expression, DNA methylation, and chromatin modification, as well as chromatin accessibility and interactions. In addition to informing our understanding of genome function, including comparative approaches to elucidate constrained sequence or epigenetic elements, these annotation maps will improve the precision and sensitivity of genomic selection for animal improvement. A scientific community-driven effort has already created a coordinated data collection and analysis enterprise crucial for the success of this global effort. Although it is early in this continuing process, functional data have already been produced and application to genetic improvement reported. The functional annotation delivered by the FAANG initiative will add value and utility to the greatly improved genome sequences being established for domesticated animal species.

Whole-genome bisulfite sequencing of goat skins identifies signatures associated with hair cycling

Background

Hair follicles (HFs), upon development, undergo repetitive cycles of growth (anagen), regression (catagen), and rest (telogen). The transition between the stages is determined by multiple molecular signals, including DNA methylation, which plays important roles in mammalian cellular identity and is essential for the development of HFs. Secondary hair follicles (SHFs) in cashmere goat exhibit classic cyclic hair development, and little has been done on a genome-wide scale to examine potentially methylated genes involved in the hair cyclic transition.

Results

Genome-wide DNA methylation profiles between skin tissues sampled during the anagen and telogen stages in cashmere goats were investigated using whole-genome bisulfite sequencing (WGBS). The methylation status was observed to be higher in the skin samples with HFs in the telogen than those in the anagen stage. A total of 1311 differentially methylated regions (DMRs) were identified between the two groups, which contained 493 fully annotated DMR-related genes (DMGs) (269 Hyper- DMGs and 224 Hypo-DMGs). Furthermore, a significant over-representation of the functional categories for DMGs related to immune response and intercellular crosstalk during hair cycling was observed. By integrating DNA methylation and mRNA expression data, we revealed that four genes (FMN1, PCOLCE, SPTLC3, and COL5A1) are crucial factors for elucidating epigenetic mechanisms contributing to the telogen-to-anagen transition.

Conclusion

Our study provided systematic methylome maps pertaining to the hair cycling stages (anagen vs telogen) at a single-base resolution, and revealed stage-specific methylation loci during cashmere growth or quiescence. Furthermore, we identified epigenetically regulated genes that are potentially involved in HF development and growth in cashmere goats, and likely in other mammal species.

Electronic supplementary material

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

Transcriptome analysis reveals long intergenic non-coding RNAs involved in skeletal muscle growth and development in pig

Long intergenic non-coding RNAs (lincRNAs) play essential roles in numerous biological processes and are widely studied. The skeletal muscle is an important tissue that plays an essential role in individual movement ability. However, lincRNAs in pig skeletal muscles are largely undiscovered and their biological functions remain elusive. In this study, we assembled transcriptomes using RNA-seq data published in previous studies of our laboratory group and identified 323 lincRNAs in porcine leg muscle. We found that these lincRNAs have shorter transcript length, fewer exons and lower expression level than protein-coding genes. Gene ontology and pathway analyses indicated that many potential target genes (PTGs) of lincRNAs were involved in skeletal-muscle-related processes, such as muscle contraction and muscle system process. Combined our previous studies, we found a potential regulatory mechanism in which the promoter methylation of lincRNAs can negatively regulate lincRNA expression and then positively regulate PTG expression, which can finally result in abnormal phenotypes of cloned piglets through a certain unknown pathway. This work detailed a number of lincRNAs and their target genes involved in skeletal muscle growth and development and can facilitate future studies on their roles in skeletal muscle growth and development.

Genome-wide analysis of DNA methylation in hypothalamus and ovary of Capra hircus

BACKGROUND

DNA methylation is a frequently studied epigenetic modification due to its role in regulating gene expression and hence in biological processes and in determining phenotypic plasticity in organisms. Rudimentary DNA methylation patterns for some livestock species are publically available: among these, goat methylome deserves to be further explored.

RESULTS

Genome-wide DNA methylation maps of the hypothalamus and ovary from Saanen goats were generated using Methyl-CpG binding domain protein sequencing (MBD-seq). Analysis of DNA methylation patterns indicate that the majority of methylation peaks found within genes are located gene body regions, for both organs. Analysis of the distribution of methylated sites per chromosome showed that chromosome X had the lowest number of methylation peaks. The X chromosome has one of the highest percentages of methylated CpG islands in both organs, and approximately 50% of the CpG islands in the goat epigenome are methylated in hypothalamus and ovary. Organ-specific Differentially Methylated Genes (DMGs) were correlated with the expression levels.

CONCLUSIONS

The comparison between transcriptome and methylome in hypothalamus and ovary showed that a higher level of methylation is not accompanied by a higher gene suppression. The genome-wide DNA methylation map for two goat organs produced here is a valuable starting point for studying the involvement of epigenetic modifications in regulating goat reproduction performance.

Comparative analysis of DNA methylome and transcriptome of skeletal muscle in lean-, obese-, and mini-type pigs

DNA methylation plays a pivotal role in biological processes by affecting gene expression. However, how DNA methylation mediates phenotype difference of skeletal muscle between lean-, obese-, and mini-type pigs remains unclear. We systematically carried out comparative analysis of skeletal muscle by integrating analysis of genome-wide DNA methylation, mRNA, lncRNA and miRNA profiles in three different pig breeds (obese-type Tongcheng, lean-type Landrace, and mini-type Wuzhishan pigs). We found that the differentially methylated genes (DMGs) were significantly associated with lipid metabolism, oxidative stress and muscle development. Among the identified DMGs, 253 genes were related to body-size and obesity. A set of lncRNAs and mRNAs including UCP3, FHL1, ANK1, HDAC4, and HDAC5 exhibited inversely changed DNA methylation and expression level; these genes were associated with oxidation reduction, fatty acid metabolism and cell proliferation. Gene regulatory networks involved in phenotypic variation of skeletal muscle were related to lipid metabolism, cellular movement, skeletal muscle development, and the p38 MAPK signaling pathway. DNA methylation potentially influences the propensity for obesity and body size by affecting gene expression in skeletal muscle. Our findings provide an abundant information of epigenome and transcriptome that will be useful for animal breeding and biomedical research.

Reduced representation bisulphite sequencing of ten bovine somatic tissues reveals DNA methylation patterns and their impacts on gene expression

Background

As a major epigenetic component, DNA methylation plays important functions in individual development and various diseases. DNA methylation has been well studied in human and model organisms, but only limited data exist in economically important animals like cattle.

Results

Using reduced representation bisulphite sequencing (RRBS), we obtained single-base-resolution maps of bovine DNA methylation from ten somatic tissues. In total, we evaluated 1,868,049 cytosines in CG-enriched regions. While we found slightly low methylation levels (29.87 to 38.06 %) in cattle, the methylation contexts (CGs and non-CGs) of cattle showed similar methylation patterns to other species. Non-CG methylation was detected but methylation levels in somatic tissues were significantly lower than in pluripotent cells. To study the potential function of the methylation, we detected 10,794 differentially methylated cytosines (DMCs) and 836 differentially methylated CG islands (DMIs). Further analyses in the same tissues revealed many DMCs (including non-CGs) and DMIs, which were highly correlated with the expression of genes involved in tissue development.

Conclusions

In summary, our study provides a baseline dataset and essential information for DNA methylation profiles of cattle.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3116-1) contains supplementary material, which is available to authorized users.

Genome-wide analysis of DNA methylation in obese, lean, and miniature pig breeds

DNA methylation is a crucial epigenetic modification involved in diverse biological processes. There is significant phenotypic variance between Chinese indigenous and western pig breeds. Here, we surveyed the genome-wide DNA methylation profiles of blood leukocytes from three pig breeds (Tongcheng, Landrace, and Wuzhishan) by methylated DNA immunoprecipitation sequencing. The results showed that DNA methylation was enriched in gene body regions and repetitive sequences. LINE/L1 and SINE/tRNA-Glu were the predominant methylated repeats in pigs. The methylation level in the gene body regions was higher than in the 5' and 3' flanking regions of genes. About 15% of CpG islands were methylated in the pig genomes. Additionally, 2,807, 2,969, and 5,547 differentially methylated genes (DMGs) were identified in the Tongcheng vs. Landrace, Tongcheng vs. Wuzhishan, and Landrace vs. Wuzhishan comparisons, respectively. A total of 868 DMGs were shared by the three contrasts. The DMGs were significantly enriched in development- and metabolism-related biological processes and pathways. Finally, we identified 32 candidate DMGs associated with phenotype variance in pigs. Our research provides a DNA methylome resource for pigs and furthers understanding of epigenetically regulated phenotype variance in mammals.

Profiling the genome-wide DNA methylation pattern of porcine ovaries using reduced representation bisulfite sequencing

Substantial evidence has shown that DNA methylation regulates the initiation of ovarian and sexual maturation. Here, we investigated the genome-wide profile of DNA methylation in porcine ovaries at single-base resolution using reduced representation bisulfite sequencing. The biological variation was minimal among the three ovarian replicates. We found hypermethylation frequently occurred in regions with low gene abundance, while hypomethylation in regions with high gene abundance. The DNA methylation around transcriptional start sites was negatively correlated with their own CpG content. Additionally, the methylation level in the bodies of genes was higher than that in their 5′ and 3′ flanking regions. The DNA methylation pattern of the low CpG content promoter genes differed obviously from that of the high CpG content promoter genes. The DNA methylation level of the porcine ovary was higher than that of the porcine intestine. Analyses of the genome-wide DNA methylation in porcine ovaries would advance the knowledge and understanding of the porcine ovarian methylome.

Genome-wide DNA methylation profile of developing deciduous tooth germ in miniature pigs

Background

DNA methylation is an important epigenetic modification critical to the regulation of gene expression during development. To date, little is known about the role of DNA methylation in tooth development in large animal models. Thus, we carried out a comparative genomic analysis of genome-wide DNA methylation profiles in E50 and E60 tooth germ from miniature pigs using methylated DNA immunoprecipitation-sequencing (MeDIP-seq).

Results

We observed different DNA methylation patterns during the different developmental stages of pig tooth germ. A total of 2469 differentially methylated genes were identified. Functional analysis identified several signaling pathways and 104 genes that may be potential key regulators of pig tooth development from E50 to E60.

Conclusions

The present study provided a comprehensive analysis of the global DNA methylation pattern of tooth germ in miniature pigs and identified candidate genes that potentially regulate tooth development from E50 to E60.

Electronic supplementary material

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

Single Base-Resolution Methylome of the Dizygotic Sheep

Sheep is an important livestock in the world for meat, dairy and wool production. The third version of sheep reference genome has been recently assembled, but sheep DNA methylome has not been profiled yet. In this study, we report the comprehensive sheep methylome with 94.38% cytosine coverage at single base resolution by sequencing DNA samples from Longissimus dorsi of dizygotic Sunit sheep, which were bred in different habitats. We also compared methylomes between the twin sheep. DNA methylation status at genome-scale differentially methylated regions (DMRs), functional genomic regions and 248 DMR-containing genes were identified between the twin sheep. Gene ontology (GO) and KEGG annotations of these genes were performed to discover computationally predicted function. Lipid metabolism, sexual maturity and tumor-associated categories were observed to significantly enrich DMR-containing genes. These findings could be used to illustrate the relationship between phenotypic variations and gene methylation patterns.

Genome-wide analysis of DNA methylation in pigs using reduced representation bisulfite sequencing

DNA methylation plays a major role in the epigenetic regulation of gene expression. Although a few DNA methylation profiling studies of porcine genome which is one of the important biomedical models for human diseases have been reported, the available data are still limited. We tried to study methylation patterns of diverse pig tissues as a study of the International Swine Methylome Consortium to generate the swine reference methylome map to extensively evaluate the methylation profile of the pig genome at a single base resolution. We generated and analysed the DNA methylome profiles of five different tissues and a cell line originated from pig. On average, 39.85 and 62.1% of cytosine and guanine dinucleotides (CpGs) of CpG islands and 2 kb upstream of transcription start sites were covered, respectively. We detected a low rate (an average of 1.67%) of non-CpG methylation in the six samples except for the neocortex (2.3%). The observed global CpG methylation patterns of pigs indicated high similarity to other mammals including humans. The percentage of CpG methylation associated with gene features was similar among the tissues but not for a 3D4/2 cell line. Our results provide essential information for future studies of the porcine epigenome.

HEpD: a database describing epigenetic differences between Thoroughbred and Jeju horses

With the advent of next-generation sequencing technology, genome-wide maps of DNA methylation are now available. The Thoroughbred horse is bred for racing, while the Jeju horse is a traditional Korean horse bred for racing or food. The methylation profiles of equine organs may provide genomic clues underlying their athletic traits. We have developed a database to elucidate genome-wide DNA methylation patterns of the cerebrum, lung, heart, and skeletal muscle from Thoroughbred and Jeju horses. Using MeDIP-Seq, our database provides information regarding significantly enriched methylated regions beyond a threshold, methylation density of a specific region, and differentially methylated regions (DMRs) for tissues from two equine breeds. It provided methylation patterns at 784 gene regions in the equine genome. This database can potentially help researchers identify DMRs in the tissues of these horse species and investigate the differences between the Thoroughbred and Jeju horse breeds.

Utilization of farm animal genetic resources in a changing agro-ecological environment in the Nordic countries

Livestock production is the most important component of northern European agriculture and contributes to and will be affected by climate change. Nevertheless, the role of farm animal genetic resources in the adaptation to new agro-ecological conditions and mitigation of animal production’s effects on climate change has been inadequately discussed despite there being several important associations between animal genetic resources and climate change issues. The sustainability of animal production systems and future food security require access to a wide diversity of animal genetic resources. There are several genetic questions that should be considered in strategies promoting adaptation to climate change and mitigation of environmental effects of livestock production. For example, it may become important to choose among breeds and even among farm animal species according to their suitability to a future with altered production systems. Some animals with useful phenotypes and genotypes may be more useful than others in the changing environment. Robust animal breeds with the potential to adapt to new agro-ecological conditions and tolerate new diseases will be needed. The key issue in mitigation of harmful greenhouse gas effects induced by livestock production is the reduction of methane (CH₄) emissions from ruminants. There are differences in CH₄ emissions among breeds and among individual animals within breeds that suggest a potential for improvement in the trait through genetic selection. Characterization of breeds and individuals with modern genomic tools should be applied to identify breeds that have genetically adapted to marginal conditions and to get critical information for breeding and conservation programs for farm animal genetic resources. We conclude that phenotyping and genomic technologies and adoption of new breeding approaches, such as genomic selection introgression, will promote breeding for useful characters in livestock species.