Genome-wide analysis of DNA Methylation profiles on sheep ovaries associated with prolificacy using whole-genome Bisulfite sequencing

Methods in DNA methylation array dataset analysis: A review

Understanding the intricate relationships between gene expression levels and epigenetic modifications in a genome is crucial to comprehending the pathogenic mechanisms of many diseases. With the advancement of DNA Methylome Profiling techniques, the emphasis on identifying Differentially Methylated Regions (DMRs/DMGs) has become crucial for biomarker discovery, offering new insights into the etiology of illnesses. This review surveys the current state of computational tools/algorithms for the analysis of microarray-based DNA methylation profiling datasets, focusing on key concepts underlying the diagnostic/prognostic CpG site extraction. It addresses methodological frameworks, algorithms, and pipelines employed by various authors, serving as a roadmap to address challenges and understand changing trends in the methodologies for analyzing array-based DNA methylation profiling datasets derived from diseased genomes. Additionally, it highlights the importance of integrating gene expression and methylation datasets for accurate biomarker identification, explores prognostic prediction models, and discusses molecular subtyping for disease classification. The review also emphasizes the contributions of machine learning, neural networks, and data mining to enhance diagnostic workflow development, thereby improving accuracy, precision, and robustness.

DNA Methylation of Postnatal Liver Development in Pigs

DNA methylation plays an important role in the development and tissue differentiation of eukaryotes. In this study, bisulfite sequencing (BS-seq) technology was used to analyze the DNA methylation profiles of liver tissues taken from Rongchang pigs at three postnatal feeding stages, including newborn, suckling, and adult. The DNA methylation pattern across the genomes or genic region showed little difference between the three stages. We observed 419 differentially methylated regions (DMRs) in promoters, corresponding to 323 genes between newborn and suckling stages, in addition to 288 DMRs, corresponding to 134 genes, between suckling and adult stages and 351 DMRs, corresponding to 293 genes, between newborn and adult stages. These genes with DMRs were mainly enriched in metabolic, immune-related functional processes. Correlation analysis showed that the methylation level of gene promoters was significantly negatively correlated with gene expression. Further, we found that genes related to nutritional metabolism, e.g., carbohydrate metabolism (FAHD1 and GUSB) or fatty acid metabolism (LPIN1 and ACOX2), lost DNA methylation in their promoter, with mRNA expression increased in newborn pigs compared with those in the suckling stage. A few fatty acid metabolism-related genes (SLC27A5, ACOX2) were hypomethylated and highly expressed in the newborn stage, which might satisfy the nutritional requirements of Rongchang pigs with high neonatal birth rates. In the adult stage, HMGCS2-which is related to fatty acid β-oxidation-was hypomethylated and highly expressed, which explains that the characteristics of high energy utilization in adult Rongchang pigs and their immune-related genes (CD68, STAT2) may be related to the establishment of liver immunity. This study provides a comprehensive analysis of genome-wide DNA methylation patterns in pig liver postnatal development and growth. Our findings will serve as a valuable resource in hepatic metabolic studies and the agricultural food industry.

Comparative transcriptome and methylome of polar bears, giant and red pandas reveal diet-driven adaptive evolution

Epigenetic regulation plays an important role in the evolution of species adaptations, yet little information is available on the epigenetic mechanisms underlying the adaptive evolution of bamboo-eating in both giant pandas (Ailuropoda melanoleuca) and red pandas (Ailurus fulgens). To investigate the potential contribution of epigenetic to the adaptive evolution of bamboo-eating in giant and red pandas, we performed hepatic comparative transcriptome and methylome analyses between bamboo-eating pandas and carnivorous polar bears (Ursus maritimus). We found that genes involved in carbohydrate, lipid, amino acid, and protein metabolism showed significant differences in methylation and expression levels between the two panda species and polar bears. Clustering analysis of gene expression revealed that giant pandas did not form a sister group with the more closely related polar bears, suggesting that the expression pattern of genes in livers of giant pandas and red pandas have evolved convergently driven by their similar diets. Compared to polar bears, some key genes involved in carbohydrate metabolism and biological oxidation and cholesterol synthesis showed hypomethylation and higher expression in giant and red pandas, while genes involved in fat digestion and absorption, fatty acid metabolism, lysine degradation, resistance to lipid peroxidation and detoxification showed hypermethylation and low expression. Our study elucidates the special nutrient utilization mechanism of giant pandas and red pandas and provides some insights into the molecular mechanism of their adaptive evolution of bamboo feeding. This has important implications for the breeding and conservation of giant pandas and red pandas.

Analysis of 206 whole-genome resequencing reveals selection signatures associated with breed-specific traits in Hu sheep

As an invaluable Chinese sheep germplasm resource, Hu sheep are renowned for their high fertility and beautiful wavy lambskins. Their distinctive characteristics have evolved over time through a combination of artificial and natural selection. Identifying selection signatures in Hu sheep can provide a straightforward insight into the mechanism of selection and further uncover the candidate genes associated with breed-specific traits subject to selection. Here, we conducted whole-genome resequencing on 206 Hu sheep individuals, each with an approximate 6-fold depth of coverage. And then we employed three complementary approaches, including composite likelihood ratio, integrated haplotype homozygosity score and the detection of runs of homozygosity, to detect selection signatures. In total, 10 candidate genomic regions displaying selection signatures were simultaneously identified by multiple methods, spanning 88.54 Mb. After annotating, these genomic regions harbored collectively 92 unique genes. Interestingly, 32 candidate genes associated with reproduction were distributed in nine genomic regions detected. Out of them, two stood out as star candidates: BMPR1B and GNRH2, both of which have documented associations with fertility, and a HOXA gene cluster (HOXA1-5, HOXA9, HOXA10, HOXA11 and HOXA13) had also been linked to fertility. Additionally, we identified other genes that are related to hair follicle development (LAMTOR3, EEF1A2), ear size (HOXA1, KCNQ2), fat tail formation (HOXA10, HOXA11), growth and development (FAF1, CCNDBP1, GJB2, GJA3), fat deposition (ACOXL, JAZF1, HOXA3, HOXA4, HOXA5, EBF4), immune (UBR1, FASTKD5) and feed intake (DAPP1, RNF17, NPBWR2). Our results offer novel insights into the genetic mechanisms underlying the selection of breed-specific traits in Hu sheep and provide a reference for sheep genetic improvement programs.

Identification of Susceptibility Genes Underlying Bovine Respiratory Disease in Xinjiang Brown Cattle Based on DNA Methylation

DNA methylation is a form of epigenetic regulation, having pivotal parts in controlling cellular expansion and expression levels within genes. Although blood DNA methylation has been studied in humans and other species, its prominence in cattle is largely unknown. This study aimed to methodically probe the genomic methylation map of Xinjiang brown (XJB) cattle suffering from bovine respiratory disease (BRD), consequently widening cattle blood methylome ranges. Genome-wide DNA methylation profiling of the XJB blood was investigated through whole-genome bisulfite sequencing (WGBS). Many differentially methylated regions (DMRs) obtained by comparing the cases and controls groups were found within the CG, CHG, and CHH (where H is A, T, or C) sequences (16,765, 7502, and 2656, respectively), encompassing 4334 differentially methylated genes (DMGs). Furthermore, GO/KEGG analyses showed that some DMGs were involved within immune response pathways. Combining WGBS-Seq data and existing RNA-Seq data, we identified 71 significantly differentially methylated (DMGs) and expressed (DEGs) genes (p < 0.05). Next, complementary analyses identified nine DMGs (LTA, STAT3, IKBKG, IRAK1, NOD2, TLR2, TNFRSF1A, and IKBKB) that might be involved in the immune response of XJB cattle infected with respiratory diseases. Although further investigations are needed to confirm their exact implication in the involved immune processes, these genes could potentially be used for a marker-assisted selection of animals resistant to BRD. This study also provides new knowledge regarding epigenetic control for the bovine respiratory immune process.

Genome-Wide DNA Methylation Analysis and Functional Validation of Litter Size Traits in Jining Grey Goats

DNA methylation (DNAm) is associated with the reproductive system. However, the genetic mechanism through which DNAm regulates gene expression and thus affects litter size in goats is unclear. Therefore, in the present work, genome-wide DNAm profiles of HP and LP Jining Grey goat ovary tissues were comprehensively analyzed via WGBS, and RNA-Seq data were combined to identify candidate genes associated with litter size traits in the Jining Grey goat. Finally, BSP and RT-qPCR were used to verify the sequencing results of the key genes. Notably, the DNMT genes were downregulated at the expression level in the HP group. Both groups exhibited comparable levels of methylation. A total of 976 differentially methylated regions (DMRs) (973 DMRs for CG and 3 DMRs for CHG) and 310 differentially methylated genes (DMGs) were identified in this study. Through integration of WGBS and RNA-Seq data, we identified 59 differentially methylated and differentially expressed genes (DEGs) and ultimately screened 8 key DMGs (9 DMRS) associated with litter size traits in Jining Grey goats (SERPINB2: chr24_62258801_62259000, NDRG4: chr18_27599201_27599400, CFAP43: chr26_27046601_27046800, LRP1B. chr2_79720201_79720400, EPHA6: chr1_40088601_40088800, TTC29: chr17_59385801_59386000, PDE11A: chr2_117418601_117418800 and PGF: chr10_ 16913801_16914000 and chr10_16916401_16916600). In summary, our research comprehensively analyzed the genome-wide DNAm profiles of HP and LP Jining Grey goat ovary tissues. The data findings suggest that DNAm in goat ovaries may play an important role in determining litter size.

Pan-Omics in Sheep: Unveiling Genetic Landscapes

Multi-omics-integrated analysis, known as panomics, represents an advanced methodology that harnesses various high-throughput technologies encompassing genomics, epigenomics, transcriptomics, proteomics, and metabolomics. Sheep, playing a pivotal role in agricultural sectors due to their substantial economic importance, have witnessed remarkable advancements in genetic breeding through the amalgamation of multiomics analyses, particularly with the evolution of high-throughput technologies. This integrative approach has established a robust theoretical foundation, enabling a deeper understanding of sheep genetics and fostering improvements in breeding strategies. The comprehensive insights obtained through this approach shed light on diverse facets of sheep development, including growth, reproduction, disease resistance, and the quality of livestock products. This review primarily focuses on the application of principal omics analysis technologies in sheep, emphasizing correlation studies between multiomics data and specific traits such as meat quality, wool characteristics, and reproductive features. Additionally, this paper anticipates forthcoming trends and potential developments in this field.

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.

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.

Comparison of genome-wide DNA methylation patterns between antler precartilage and cartilage

Deer antlers are the only mammalian organs that can fully regenerate after being lost and provide a valuable model for cartilage development. As one of the best-studied epigenetic mechanisms, DNA methylation is known to engage in organ and tissue development. This study aimed to investigate the role of DNA methylation in antler chondrogenesis by comparing whole-genome DNA methylation between precartilage and cartilage. Quantitative reverse transcription PCR (RT-qPCR) showed significant differences in the expression levels of DNA methyltransferase genes (DNMT1, DNMT3A, and DNMT3B) between precartilage and cartilage. Subsequently, we obtained DNA methylation profiles of antler precartilage and cartilage tissues by whole-genome bisulfite sequencing. Although sequencing data indicated that overall methylation levels at CpG and non-CpG sites were similar between precartilage and cartilage, 140,784 differentially methylated regions (DMRs, P < 0.05) and 3,941 DMR-related genes were identified. Gene ontology (GO) analysis of DMR-related genes demonstrated some significantly enriched GO terms (P < 0.05) related to chondrogenesis, including insulin receptor binding, collage trimer, integrin binding, and extracellular matrix structural constituent. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of DMR-related genes uncovered that the PI3K/AKT, cortisol synthesis and secretion, glycosaminoglycan biosynthesis-keratan sulfate, Hippo, and NF-κB signaling pathways might play a pivotal role in the transition of precartilage to cartilage. Moreover, we found that 25 DMR-related genes, including CD44, IGF1, ITGAV, ITGB1, RUNX1, COL2A1, COMP, and TAGLN, were most likely involved in antler chondrogenesis. In conclusion, this study revealed the genome-wide DNA methylation patterns of antler precartilage and cartilage, which may contribute to understanding the epigenetic regulation of antler chondrogenesis.

Transcriptome and DNA Methylation Responses in the Liver of Yellowfin Seabream Under Starvation Stress

Fish suffer from starvation due to environmental risks such as extreme weather in the wild and due to insufficient feedings in farms. Nutrient problems from short-term or long-term starvation conditions can result in stress-related health problems for fish. Yellowfin seabream (Acanthopagrus latus) is an important marine economic fish in China. Understanding the molecular responses to starvation stress is vital for propagation and culturing yellowfin seabream. In this study, the transcriptome and genome-wide DNA methylation levels in the livers of yellowfin seabream under 14-days starvation stress were analyzed. One hundred sixty differentially expressed genes (DEGs) by RNA-Seq analysis and 737 differentially methylated-related genes by whole genome bisulfite sequencing analysis were identified. GO and KEGG pathway enrichment analysis found that energy metabolism-related pathways such as glucose metabolism and lipid metabolism were in response to starvation. Using bisulfite sequencing PCR, we confirmed the presence of CpG methylation differences within the regulatory region of a DEG ppargc1a in response to 14-days starvation stress. This study revealed the molecular responses of livers in response to starvation stress at the transcriptomic and whole genome DNA methylation levels in yellowfin seabream.

Short variation of the sheep PDGFD gene is correlated with litter size

Platelet-derived growth factor (PDGF) family, exert plays a key role in embryonic development, cell proliferation, cell migration, angiogenesis and reproduction. Related studies about GWAS analyses have found that PDGFD significantly affected deposition of tail fat in sheep, but there are no studies on reproduction in animals. In this study, three breed of sheep were used to find insertion/deletion (indel) fragment polymorphism of PDGFD which including Australian white (AUW) sheep (Meat type, n = 932), Guiqian semi-fine wool (GSFW) sheep (wool type, n = 60) and East Friensian milk (EFM) sheep (dairy type, n = 60). Only a 18-bp variation was polymorphic in the study AUW sheep population and the genotypes of different sheep breed are also specific. Moreover, the association analysis indicated that this variant was associated with litter size of AUW sheep in the first parity (p < 0.05). The litter size of II genotype was significantly lower than other genotypes in the first parity (p < 0.05). We also revealed that the PDGFD gene was relatively conservative in eight species, PDGFD mRNA expression in 832 sheep samples implying this gene was related to reproduction traits. Hence, these finding demonstrated the one-cause multipotency of PDGFD gene. Collectively, these results suggest that this indel can be used as an effective marker for sheep breeding.

Comparative Transcriptomics and Methylomics Reveal Adaptive Responses of Digestive and Metabolic Genes to Dietary Shift in Giant and Red Pandas

Both the giant panda (Ailuropoda melanoleuca) and red panda (Ailurus fulgens) belong to the order Carnivora, but have changed their dietary habits to eating bamboo exclusively. The convergent evolution characteristics of their morphology, genome and gut flora have been found in the two pandas. However, the research on the convergent adaptation of their digestion and metabolism to the bamboo diet, mediated by the dietary shift of the two pandas at the gene-expression and epigenetic regulation levels, is still lacking. We therefore used RNA sequencing among five species (two pandas and three non-herbivore mammals) and bisulfite sequencing among three species (two pandas and a carnivore ferret) to sequence key digestion and metabolism tissues (stomach and small intestine). Our results provide evidence that the convergent differentially expressed genes (related to carbohydrate utilization, bile secretion, Lys and Arg metabolism, vitamin B12 utilization and cyanide detoxification) of the two pandas are adaptive responses to the bamboo diet containing low lipids, low Lys and Arg, low vitamin B12 and high cyanide. We also profiled the genome-wide methylome maps of giant panda, red panda and ferret, and the results indicated that the promoter methylation of the two pandas may regulate digestive and metabolic genes to adapt to sudden environmental changes, and then, transmit genetic information to future generations to evolve into bamboo eaters. Taken together, our study provides new insights into the molecular mechanisms of the dietary shift and the adaptation to a strict bamboo diet in both pandas using comparative transcriptomics and methylomics.

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

Simple Summary

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

Abstract

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

DNA methylation plays an important role in iron-overloaded Tibetans

The prevalence of iron overload in Tibetans in Tibet is higher than that in Han. DNA methylation (DNAm) is closely related to iron metabolism and iron level. Nevertheless, the epigenetic status of Tibetans with iron overload is unknown, and we therefore aimed to explore whether the phenomenon observed in the Tibetan population is regulated by epigenetics. The results showed that 2.26% of cytosine was methylated in the whole genome, and that the rate of CG cytosine methylation was higher in individuals in the iron overload (TH) group than in those in the iron normal (TL) group. We analyzed differentially methylated genes (DMGs) in whole-genome bisulfite sequencing data from the TH and TL groups of high-altitude Tibetans. Protein-protein interaction and pathway analyses of candidate DMGs related to iron uptake and transport showed that epigenetic changes in 10 candidate genes (ACO1, CYBRD1, FLVCR1, HFE, HMOX2, IREB2, NEDD8, SLC11A2, SLC40A1 and TFRC) are likely to relate to iron overload. This work reveals, for the first time, changes of DNAm in Tibetan people with iron overload, which suggest that DNAm is a mechanism underlying differences in iron content between individuals in the high-altitude Tibetan population. Our findings should contribute to the study of iron metabolism and the overall health status of Tibetans.

Epigenomic profiling indicates a role for DNA methylation in the postnatal liver and pancreas development of giant pandas

Giant pandas are unique within Carnivora with a strict bamboo diet. Here, the epigenomic profiles of giant panda liver and pancreas tissues collected from three important feeding stages were investigated using BS-seq. Few differences in DNA methylation profiles were exhibited between no feeding and suckling groups in both tissues. However, we observed a tendency toward a global loss of DNA methylation in the gene-body and promoter region of metabolism-related genes from newborn to adult. Correlation analysis revealed a significant negative correlation between the changes in methylation levels within gene promoters and gene expression. The majority of genes related to nutrition metabolism had lost DNA methylation with increased mRNA expression in adult giant pandas. The few galactose metabolism and unsaturated fatty acid metabolism related genes that were hypomethylated and highly-expressed at early stages of giant panda development may meet the nutritional requirement of this species' highly altricial neonates.

Variability in Global DNA Methylation Rate Across Tissues and Over Time in Sheep

Recent studies highlighted the influence of epigenetic marks in the variability of many complex traits, both in plants and animals. These studied focused only on specific sites of the genome having differentially methylated profiles among individuals and/or tissues. In contrast, we recently used the methylation rate of the entire genome as a unique measure considered as a novel quantitative phenotype in sheep. This phenotype named global DNA methylation rate (GDMR), measured by luminometric assay, integrates the methylation level of each CpG dinucleotide within the 6 million of CCGG sites along the ovine genome. GDMR measured in blood previously showed moderate heritability of 0.20 and provided evidence for a genetic determinism. The main objective of the present study was to better characterize the GDMR phenotype in various tissues and investigate its variability in several breeds of sheep reared in the same environment. GDMR was measured on blood samples collected monthly from 59 growing male and female lambs (24 Romane, 23 Blackbelly and 12 Charollais), between birth and 4 months of age. Blood GDMR was on average around 80% and was influenced by the sampling date (p &lt; 0.001), the breed (p = 0.002) and the sex (p = 0.002). In addition, GDMR was determined in 12 somatic (frontal lobe, pituitary gland, heart, lung, sub cutaneous and perirenal adipose tissue, skeletal muscle, liver, spleen, adrenal gland, medulla and cortical kidney) and 6 reproductive tissues (ovary, oviduct, uterus, testis, epididymis and seminal vesicle). GDMR was on average 70% in somatic tissues but marked variation was observed depending on the tissue. The GDMR measured in blood was higher than that measured in other somatic tissues, and is not a good proxy of less accessible tissues. Female reproductive tissues had a 10% higher GDMR than male reproductive tissues. We demonstrated a significant influence of the breed on blood GDMR, certainly reflecting the influence of different genetic backgrounds. The effect of the breed on GDMR may be related to their specific abilities to adapt to and live in different conditions.

Genome-Wide Methylation Profiling in the Thalamus of Scrapie Sheep

Scrapie is a neurodegenerative disorder belonging to the group of transmissible spongiform encephalopathy (TSE). Scrapie occurs in sheep and goats, which are considered good natural animal models of these TSE. Changes in DNA methylation occur in the central nervous system (CNS) of patients suffering from prion-like neurodegenerative diseases, such as Alzheimer's disease. Nevertheless, potential DNA methylation alterations have not yet been investigated in the CNS of any prion disease model or naturally infected cases, neither in humans nor in animals. Genome-wide DNA methylation patterns were studied in the thalamus obtained from sheep naturally infected with scrapie at a clinical stage (n = 4) and from controls (n = 4) by performing a whole-genome bisulfite sequencing (WGBS) analysis. Ewes carried the scrapie-susceptible ARQ/ARQ PRNP genotype and were sacrificed at a similar age (4–6 years). Although the average genomic methylation levels were similar between the control and the scrapie animals, we identified 8,907 significant differentially methylated regions (DMRs) and 39 promoters (DMPs). Gene Ontology analysis revealed that hypomethylated DMRs were enriched in genes involved in transmembrane transport and cell adhesion, whereas hypermethylated DMRs were related to intracellular signal transduction genes. Moreover, genes highly expressed in specific types of CNS cells and those previously described to be differentially expressed in scrapie brains contained DMRs. Finally, a quantitative PCR (qPCR) validation indicated differences in the expression of five genes (PCDH19, SNCG, WDR45B, PEX1, and CABIN1) that matched the methylation changes observed in the genomic study. Altogether, these results suggest a potential regulatory role of DNA methylation in prion neuropathology.

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

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

Whole Genome DNA Methylation Variations in Mammary Gland Tissues from Holstein Cattle Producing Milk with Various Fat and Protein Contents

Milk fat and protein contents are among key elements of milk quality, and they are attracting more attention in response to consumers' demand for high-quality dairy products. To investigate the potential regulatory roles of DNA methylation underlying milk component yield, whole genome bisulfite sequencing was employed to profile the global DNA methylation patterns of mammary gland tissues from 17 Canada Holstein cows with various milk fat and protein contents. A total of 706, 2420 and 1645 differentially methylated CpG sites (DMCs) were found between high vs. low milk fat (HMF vs. LMF), high vs. low milk protein (HMP vs. LMP), and high vs. low milk fat and protein (HMFP vs. LMFP) groups, respectively (q value < 0.1). Twenty-seven, 56 and 67 genes harboring DMCs in gene regions (denoted DMC genes) were identified for HMF vs. LMF, HMP vs. LMP and HMFP vs. LMFP, respectively. DMC genes from HMP vs. LMP and HMFP vs. LMFP comparisons were significantly overrepresented in GO terms related to aerobic electron transport chain and/or mitochondrial ATP (adenosine triphosphate) synthesis coupled electron transport. A total of 83 (HMF vs. LMF), 708 (HMP vs. LMP) and 408 (HMFP vs. LMFP) DMCs were co-located with 87, 147 and 158 quantitative trait loci (QTL) for milk component and yield traits, respectively. In conclusion, the identified methylation changes are potentially involved in the regulation of milk fat and protein yields, as well as the variation in reported co-located QTLs.

Combined approaches identify known and novel genes associated with sheep litter size and non-seasonal breeding

Improvement of ewe reproduction is considerable by appropriately increasing litter size and sustaining non-seasonal breeding. However, their genetic makeups have not been entirely elucidated. Genome-wide analyses of 821 individuals were performed by combining three genomic approaches (genome-wide association study, XP-nSL, and runs of homozygosity). Consequently, 35 candidate genes including three domestication genes (TSHR, GTF2A1, and KITLG) were identified. Other than the FecB mutation at BMPR1B, we described a significant association of a missense mutation rs406686139 at seasonal lambing-associated TSHR gene with litter size. Some promising novel genes may be relevant for sheep reproduction by multitude biological processes, such as FETUB functioning in fertilization, HNRNPA1 in oogenesis, DCUN1D1 in spermatogenesis, and HRG in fertility outcome. The present study suggests that improvement of ewe reproduction is attributed to selective breeding, and casts light on the genetic basis and improvement of sheep reproduction.

Global analysis of DNA methylation in hepatocellular carcinoma via a whole-genome bisulfite sequencing approach

Alterations in DNA methylation patterns are considered early events in hepatocellular carcinoma (HCC). However, their mechanism and significance remain to be elucidated. We studied the genome-wide DNA methylation landscape of HCC by applying whole-genome bisulfite sequencing (WGBS) techonlogy. Overall, HCC exhibits a genome-wide hypomethylation pattern. After further annotation, we obtained 590 differentially hypermethylated genes (hyper-DMGs) and 977 differentially hypomethylated genes (hypo-DMGs) from three groups. Hyper-DMGs were mainly involved in ascorbate and alternate metabolism pathways, while hypo-DMGs were mainly involved in focal adhesion. By integrating the DMGs with HCC-related differentially expressed genes (DEGs) and DMGs from the TCGA database, we constructed prognostic model based on thirteen aberrantly methylated DEGs, and verified our prognostic model in GSE14520 dataset. This study compares the patterns of global epigenomic DNA methylation during the development of HCC, focusing on the role of DNA methylation in the early occurrence and development of HCC, providing a direction for future research on its epigenetic mechanism.

Metabolism and Sex Differentiation in Animals from a Starvation Perspective

Animals determine their sex genetically (GSD: genetic sex determination) and/or environmentally (ESD: environmental sex determination). Medaka (Oryzias latipes) employ a XX/XY GSD system, however, they display female-to-male sex reversal in response to various environmental changes such as temperature, hypoxia, and green light. Interestingly, we found that 5 days of starvation during sex differentiation caused female-to-male sex reversal. In this situation, the metabolism of pantothenate and fatty acid synthesis plays an important role in sex reversal. Metabolism is associated with other biological factors such as germ cells, HPG axis, lipids, and epigenetics, and supplys substances and acts as signal transducers. In this review, we discuss the importance of metabolism during sex differentiation and how metabolism contributes to sex differentiation.

The ovarian follicle of ruminants: the path from conceptus to adult

This review resulted from an international workshop and presents a consensus view of critical advances over the past decade in our understanding of follicle function in ruminants. The major concepts covered include: (1) the value of major genes; (2) the dynamics of fetal ovarian development and its sensitivity to nutritional and environmental influences; (3) the concept of an ovarian follicle reserve, aligned with the rise of anti-Müllerian hormone as a controller of ovarian processes; (4) renewed recognition of the diverse and important roles of theca cells; (5) the importance of follicular fluid as a microenvironment that determines oocyte quality; (6) the 'adipokinome' as a key concept linking metabolic inputs with follicle development; and (7) the contribution of follicle development to the success of conception. These concepts are important because, in sheep and cattle, ovulation rate is tightly regulated and, as the primary determinant of litter size, it is a major component of reproductive efficiency and therefore productivity. Nowadays, reproductive efficiency is also a target for improving the 'methane efficiency' of livestock enterprises, increasing the need to understand the processes of ovarian development and folliculogenesis, while avoiding detrimental trade-offs as greater performance is sought.

Genome-Wide Analyses Reveal Genetic Convergence of Prolificacy between Goats and Sheep

The litter size of domestic goats and sheep is an economically important trait that shows variation within breeds. Strenuous efforts have been made to understand the genetic mechanisms underlying prolificacy in goats and sheep. However, there has been a paucity of research on the genetic convergence of prolificacy between goats and sheep, which likely arose because of similar natural and artificial selection forces. Here, we performed comparative genomic and transcriptomic analyses to identify the genetic convergence of prolificacy between goats and sheep. By combining genomic and transcriptomic data for the first time, we identified this genetic convergence in (1) positively selected genes (CHST11 and SDCCAG8), (2) differentially expressed genes (SERPINA14, RSAD2, and PPIG at follicular phase, and IGF1, GPRIN3, LIPG, SLC7A11, and CHST15 at luteal phase), and (3) biological pathways (genomic level: osteoclast differentiation, ErbB signaling pathway, and relaxin signaling pathway; transcriptomic level: the regulation of viral genome replication at follicular phase, and protein kinase B signaling and antigen processing and presentation at luteal phase). These results indicated the potential physiological convergence and enhanced our understanding of the overlapping genetic makeup underlying litter size in goats and sheep.

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.

Multiparametric Study of Antioxidant Effect on Ram Sperm Cryopreservation-From Field Trials to Research Bench

Simple Summary

The use of antioxidant compounds could be a successful tool to improve sperm cryopreservation protocols in ovine species. These molecules have been widely employed in different mammalian species with this purpose. It is important to consider the existence of a species-specific antioxidant effect discarding the extrapolations from other animal species. To corroborate the real effectiveness of these compounds is important to combine two approaches: in vitro sperm quality analyses and in vivo field trials based on fertility. In the first scenario, a multiparametric analyses and novel tests based on spermatozoa redox balance, as the main target of antioxidants, could improve the accuracy on antioxidant effectiveness on sperm quality. Moreover, an extensive field insemination study provides the definitive tool to select the best antioxidant treatment. All these aspects have been applied and extensively discussed throughout this manuscript. Novel approaches have been incorporated, such as RedoxSYS, to provide more accuracy in the integrative studies of Redox status in spermatozoa. The effectiveness of an antioxidant treatment, as trolox in our study, should be demonstrated in an integrative way, from in vivo (fertility trials) to in vitro analyses (sperm quality assays), especially when the final aim is to reach AI implementation.

Abstract

The optimization of sperm cryopreservation protocols in ram is a feasible tool to reinforce artificial insemination technologies considering the desirable application of sperm by vaginal/cervical or transcervical deposition. Cryopreservation provokes different types of damage on spermatozoa and many of these detrimental effects are triggered by redox deregulation. For this reason, the antioxidant supplementation in sperm cryopreservation protocols to decrease reactive oxygen species (ROS) levels and to equilibrate redox status has been widely employed in different species. Despite this, more fertility trials are necessary to provide the definitive tool to ensure the antioxidant effectiveness on sperm quality. For this reason, in this work, we performed a multiparametric analysis of some previously tested antioxidants (crocin, GSH and Trolox) on ram sperm cryopreservation from field trials to sperm quality analyses focused on new strategies to measure redox balance. Attending to fertility trial, Trolox supplementation registered an improvement concerning to fertility (when we considered high fertility males) and multiple lambing frequency and other complementary and descriptive data related to lambing performance such as prolificacy and fecundity. This positive effect was more evident in multiple lambing frequency when we considered low fertility males than in global male analysis. In vitro analyses of sperm quality confirmed in vivo trials registering a positive effect on sperm viability and redox balance. In this study, we provided the definitive evidence that the role of trolox on redox balance maintenance has a direct effect on fertility parameters, such as prolificacy. The effectiveness of antioxidant treatments was tested, for the first time in ovine species, using an integrative and multiparametric approach combining in vivo and in vitro analyses and novel approaches, such as RedoxSYS. These types of strategies should be applied to improve sperm conservation methods and optimize AI technologies upgrading the correlation between in vitro and in vivo analyses.

Analysis of Genome DNA Methylation at Inherited Coat Color Dilutions of Rex Rabbits

Background: The dilution of color in rabbits is associated with many different genetic mechanisms that form different color groups. A number of previous studies have revealed potential regulatory mechanisms by which epigenetics regulate pigmentation. However, the genome-wide DNA methylation involved in animal coat color dilution remains unknown. Results: We compared genome-wide DNA methylation profiles in Rex rabbit hair follicles in a Chinchilla group (Ch) and a diluted Chinchilla group (DCh) through whole-genome bisulfite sequencing (WGBS). Approximately 3.5% of the cytosine sites were methylated in both groups, of which the CG methylation type was in greatest abundance. In total, we identified 126,405 differentially methylated regions (DMRs) between the two groups, corresponding to 11,459 DMR-associated genes (DMGs). Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed that these DMGs were principally involved in developmental pigmentation and Wnt signaling pathways. In addition, two DMRs were randomly selected to verify that the WGBS data were reliable using bisulfite sequencing PCR, and seven DMGs were analyzed to establish the relationship between the level of DNA methylation and mRNA expression using qRT-PCR. Due to the limitation of small sample size, replication of the results with a larger sample size would be important in future studies. Conclusion: These findings provide evidence that there is an association between inherited color dilution and DNA methylation alterations in hair follicles, greatly contributing to our understanding of the epigenetic regulation of rabbit pigmentation.

Genetic Determinism Exists for the Global DNA Methylation Rate in Sheep

Recent studies showed that epigenetic marks, including DNA methylation, influence production and adaptive traits in plants and animals. So far, most studies dealing with genetics and epigenetics considered DNA methylation sites independently. However, the genetic basis of the global DNA methylation rate (GDMR) remains unknown. The main objective of the present study was to investigate genetic determinism of GDMR in sheep. The experiment was conducted on 1,047 Romane sheep allocated into 10 half-sib families. After weaning, all the lambs were phenotyped for global GDMR in blood as well as for production and adaptive traits. GDMR was measured by LUminometric Methylation Analysis (LUMA) using a pyrosequencing approach. Association analyses were conducted on some of the lambs (n = 775) genotyped by using the Illumina OvineSNP50 BeadChip. Blood GDMR varied among the animals (average 70.7 ± 6.0%). Female lambs had significantly higher GDMR than male lambs. Inter-individual variability of blood GDMR had an additive genetic component and heritability was moderate (h² = 0.20 ± 0.05). No significant genetic correlation was found between GDMR and growth or carcass traits, birthcoat, or social behaviors. Association analyses revealed 28 QTLs associated with blood GDMR. Seven genomic regions on chromosomes 1, 5, 11, 17, 24, and 26 were of most interest due to either high significant associations with GDMR or to the relevance of genes located close to the QTLs. QTL effects were moderate. Genomic regions associated with GDMR harbored several genes not yet described as being involved in DNA methylation, but some are already known to play an active role in gene expression. In addition, some candidate genes, CHD1, NCO3A, KDM8, KAT7, and KAT6A have previously been described to be involved in epigenetic modifications. In conclusion, the results of the present study indicate that blood GDMR in domestic sheep is under polygenic influence and provide new insights into DNA methylation genetic determinism.

Comparative Transcriptome Profiling of Ovary Tissue between Black Muscovy Duck and White Muscovy Duck with High- and Low-Egg Production

The egg-laying rate is an important indicator for evaluating fertility of poultry. In order to better understand the laying mechanism of Muscovy ducks, gene expression profiles and pathways of ovarian tissues in high- and low-laying black (BH and BL) and white Muscovy ducks (WH and WL) during the peak production period were performed by using RNA-seq. The total number of reads produced for each ovarian sample ranged from 44,344,070 to 47,963,328. A total of 113, 619 and 87 differentially expressed genes (DEGs) were identified in BH-vs-WH, BL-vs-BH and BL-vs-WL, respectively. Among them, 54, 356 and 49 genes were up regulated and 59, 263 and 38 genes were down regulated. In addition, there were only 10 up-regulated genes in WL-vs-WH. In the comparison of DEGs in black and white Muscovy ducks, two co-expressed DEG genes were detected between BH-vs-WH and BL-vs-WL and seven DEGs were co-expressed between BL-vs-BH and WL-vs-WH. The RNA-Seq data were confirmed to be reliable by qPCR. Numerous DEGs known to be involved in ovarian development were identified, including TGFβ2, NGFR, CEBPD, CPEB2, POSTN, SMOC1, FGF18, EFNA5 and SDC4. Gene Ontology (GO) annotations indicated that DEGs related to ovarian development were mainly enriched in biological processes of "circadian sleep/wake cycle process," "negative regulation of transforming growth factor-β secretion," "positive regulation of calcium ion transport" in BH-vs-WH and "cell surface receptor signaling pathway," "Notch signaling pathway" and "calcium ion transport" in BL-vs-BH. Besides, "steroid biosynthetic process," "granulosa cell development" and "egg coat formation" were mainly enriched in BL-vs-WL and "reproduction," "MAPK cascade" and "mitotic cell cycle" were mainly enriched in WL-vs-WH. KEGG pathway analysis showed that the PI3K-Akt signaling pathway and ovarian steroidogenesis were the most enriched in Muscovy duck ovary transcriptome data. This work highlights potential genes and pathways that may affect ovarian development in Muscovy duck.

Luzhong mutton sheep: inbreeding and selection signatures

Intense artificial selection has been imposed to Luzhong mutton sheep population in the past years. Improvements on growth and reproductive performance are two breeding goals in the present herd. Although some progresses were phenotypically observed possibly due to inbreeding induced by strong selection in terms of these traits, the genomic evaluation was poorly understood. Therefore, a high-density SNP array was used to characterize the pattern of runs of homozygosity (ROH), estimate inbreeding and inbreeding depressions on early growth performance and litter size based upon ROH, and scan positive selection signatures of recent population. Consequently, a low inbreeding level was observed which had negative effects on litter size, but not on early growth performance. And 160 genes were under selection, of which some were reported to be linked to several traits of sheep including body weight, litter size, carcass and meat quality, milk yield and composition, fiber quality and health, and the top genes were associated with growth (growth hormone [GH]- growth hormone receptor [GHR]- Insulin-like growth factor 1 [IGF1] axis) and litter size (bone morphogenic proteins [BMPs]-associated). The effectiveness of previous breeding measures was highlighted, but purging selection was proposed to alleviate the inbreeding depression on litter size, providing some genomic insights to breeding management of Luzhong mutton sheep.

Multi-omics analysis reveals that natural hibernation is crucial for oocyte maturation in the female Chinese alligator

BACKGROUND

Hibernation in an appropriate environment not only is important for the survival of hibernators in winter, but also is crucial for breeding in the following season for many hibernating species. However, the genetic and epigenetic mechanism underlying this process remain unclear. In the current study, we performed an integrative multi-omics analysis of gonads collected from Chinese alligators that overwintered in wild cave and artificial warmroom to explore transcriptomic and epigenomic alternations in these organs.

RESULTS

The data revealed that in the breeding season, female alligators were more strongly affected in terms of gene expression than males by non-hibernation because of overwintering in a warm room, especially for genes related to oocyte maturation, and this effect commenced in winter with the downregulation of STAR, which is the rate limiting factor of steroid biosynthesis. Further, miRNAs were found to play essential roles in this negative effect of overwintering in the warm room on hibernation. The upregulated miRNAs likely were responsible for the suppression of oocyte maturation in the breeding season. Finally, DNA methylome changes, especially hypomethylation, were found to play an important role in the alterations in ovarian function-related gene expression induced by non-hibernation.

CONCLUSIONS

Our study revealed the crucial role of hibernation quality for oocyte maturation in the Chinese alligator and the underlying genetic and epigenetic mechanisms, and highlights the importance of habitat, and especially, the overwintering site, in the conservation of not only the Chinese alligator, but also other endangered hibernators.

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

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

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

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

BMP7 is a candidate gene for reproductive traits in Yorkshire sows

Bone morphogenetic protein 7 (BMP7) is of the BMP subfamily, and has effects on female fertility by regulating steroidogenesis, granulosa cell states, and follicular development. In the present study, there was assessment of the combined genotypes formed by the three variants within the 3'-UTR of BMP7 gene as associations with sow reproductive functions. The 3'-UTR of the BMP7 gene of pigs was identified using the 3' RACE assay, and its full-length sequence was found to be 1538 bp in length. Multiple RNA regulatory elements were detected in this region, luciferase activity assays were performed and results indicated miR-22-3p affects BMP7 by directly binding to the miRNA response element in the 3'-UTR (c.2358-2382). In addition, two novel complete linkage variants, c.2256 G > C and a 7-bp indel (c.2259-2265), were identified within the 3'-UTR of the BMP7 gene of pigs. Importantly, combined genotypes with these two novel variants and c.1569A > G, a variant previously identified in the BMP7 3'-UTR of pigs, were associated with sow reproductive traits, including the total number of piglets born, number of dead piglets at birth, and litter weight in the Yorkshire pig population studies. Results from the present study confirm that BMP7 is a candidate gene for the reproductive traits in Yorkshire sows.

Comparative analysis of the transcriptome, methylome, and metabolome during pollen abortion of a seedless citrus mutant

Pollen abortion could be mainly attributed to abnormal meiosis in the mutant. Multiomics analysis uncovered significant epigenetic variations between the mutant and its wild type during the pollen abortion process. Male sterility caused by aborted pollen can result in seedless fruit. A seedless Ponkan mandarin mutant (bud sport) was used to compare the transcriptome, methylome, and metabolome with its progenitor to understand the mechanism of citrus pollen abortion. Cytological observations showed that the anther of the mutant could form microspore mother cells, although the microspores failed to develop fertile pollen at the anther dehiscence stage. Based on pollen phenotypic analysis, pollen abortion could be mainly attributed to abnormal meiosis in the mutant. A transcriptome analysis uncovered the molecular mechanisms underlying pollen abortion between the mutant and its wild type. A total of 5421 differentially expressed genes were identified, and some of these genes were involved in the meiosis, hormone biosynthesis and signaling, carbohydrate, and flavonoid pathways. A total of 50,845 differentially methylated regions corresponding to 15,426 differentially methylated genes in the genic region were found between the mutant and its wild type by the methylome analysis. The expression level of these genes was negatively correlated with their methylation level, especially in the promoter regions. In addition, 197 differential metabolites were identified between the mutant and its wild type based on the metabolome analysis. The transcription and metabolome analysis further indicated that the expression of genes in the flavonoid, carbohydrate, and hormone metabolic pathways was significantly modulated in the pollen of the mutant. These results indicated that demethylation may alleviate the silencing of carbohydrate genes in the mutant, resulting in excessive starch and sugar hydrolysis and thereby causing pollen abortion in the mutant.

Analysis of DNA methylation profiles during sheep skeletal muscle development using whole-genome bisulfite sequencing

BACKGROUND

DNA methylation is an epigenetic regulatory form that plays an important role in regulating the gene expression and the tissues development.. However, DNA methylation regulators involved in sheep muscle development remain unclear. To explore the functional importance of genome-scale DNA methylation during sheep muscle growth, this study systematically investigated the genome-wide DNA methylation profiles at key stages of Hu sheep developmental (fetus and adult) using deep whole-genome bisulfite sequencing (WGBS).

RESULTS

Our study found that the expression levels of DNA methyltransferase (DNMT)-related genes were lower in fetal muscle than in the muscle of adults. The methylation levels in the CG context were higher than those in the CHG and CHH contexts, and methylation levels were highest in introns, followed by exons and downstream regions. Subsequently, we identified 48,491, 17, and 135 differentially methylated regions (DMRs) in the CG, CHG, and CHH sequence contexts and 11,522 differentially methylated genes (DMGs). The results of bisulfite sequencing PCR (BSP) correlated well with the WGBS-Seq data. Moreover, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional annotation analysis revealed that some DMGs were involved in regulating skeletal muscle development and fatty acid metabolism. By combining the WGBS-Seq and previous RNA-Seq data, a total of 159 overlap genes were obtained between differentially expressed genes (DEGs) and DMGs (FPKM > 10 and fold change > 4). Finally, we found that 9 DMGs were likely to be involved in muscle growth and metabolism of Hu sheep.

CONCLUSIONS

We systemically studied the global DNA methylation patterns of fetal and adult muscle development in Hu sheep, which provided new insights into a better understanding of the epigenetic regulation of sheep muscle development.

Whole-Genome Methylation Analysis Reveals Epigenetic Variation in Cloned and Donor Pigs

Somatic cloning has had a significant impact on the life sciences and is important in a variety of processes, including medical research and animal production. However, the application of somatic cloning has been limited due to its low success rate. Therefore, potential epigenetic variations between cloned and donor animals are still unclear. DNA methylation, one of the factors which is responsible for phenotypic differences in animals, is a commonly researched topic in epigenetic studies of mammals. To investigate the epigenetic variations between cloned and donor animals, we selected blood and ear fibroblasts of a donor pig and a cloned pig to perform whole-genome bisulfite sequencing (WGBS). A total of 215 and 707 differential methylation genes (DMGs) were identified in blood and ear fibroblasts, respectively. Functional annotation revealed that DMGs are enriched in many pathways, including T/B or natural killer (NK) cell differentiation, oocyte maturation, embryonic development, and reproductive hormone secretion. Moreover, 22 DMGs in the blood and 75 in the ear were associated with immune responses (e.g., CD244, CDK6, CD5, CD2, CD83, and CDC7). We also found that 18 DMGs in blood and 53 in ear fibroblasts were involved in reproduction. Understanding the expression patterns of DMGs, especially in relation to immune responses and reproduction, will reveal insights that will aid the advancement of future somatic cloning techniques in swine.

Gene Expression Profiling in Ovaries and Association Analyses Reveal HEP21 as a Candidate Gene for Sexual Maturity in Chickens

Simple Summary

Chicken meat and egg productions are essential for human beings. Sexual maturity is important for both egg production and meat flavor. It is necessary to elucidate the genetic mechanism of chicken sexual maturity. In current study, we used digital gene expression (DGE) RNA-sequencing analysis to investigate differential expression of genes in pre-pubertal and post-pubertal ovaries in two different sub-breeds of chicken with different onsets of sexual maturity. After the analysis of RNA-sequencing data, numerous differentially expressed genes were found in both comparisons (32 day old, early-sexual-maturity pre-laying hens (P-F-O1) vs. 103 day old early-sexual-maturity laying hens (P-F-O2), and 32 day old late-sexual-maturity pre-laying hens (L-F-O1) vs. 153 day old late-sexual-maturity pre-laying hens (L-F-O2)). With the bioinformatic analysis, hen egg protein 21 kDa (HEP21) was chosen as the candidate gene to conduct following experiment. The variations in HEP21 were screened and association analyses between rs315156783 and reproductive traits were investigated in fifth-generation Ningdu Yellow chickens from a closely bred population. These results demonstrated that HEP21 is a candidate gene for sexual maturity and ovary development in chickens. However, the underlying mechanism of how HEP21 regulates chicken sexual maturity needs further focused studies.

Abstract

The age of onset of sexual maturity is an important reproductive trait in chickens. In this study, we explored candidate genes associated with sexual maturity and ovary development in chickens. We performed DGE RNA-sequencing analyses of ovaries of pre-laying (P-F-O1, L-F-O1) and laying (P-F-O2, L-F-O2) hens of two sub-breeds of Ningdu Yellow chicken. A total of 3197 genes were identified in the two comparisons, and 966 and 1860 genes were detected exclusively in comparisons of P-F-O1 vs. P-F-O2 and L-F-O1 vs. L-F-O2, respectively. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed that genes involved in transmembrane signaling receptor activity, cell adhesion, developmental processes, the neuroactive ligand–receptor interaction pathway, and the calcium signaling pathway were enriched in both comparisons. Genes on these pathways, including growth hormone (GH), integrin subunit beta 3 (ITGB3), thyroid stimulating hormone subunit beta (TSHB), prolactin (PRL), and transforming growth factor beta 3 (TGFB3), play indispensable roles in sexual maturity. As a gene unique to poultry, hen egg protein 21 kDa (HEP21) was chosen as the candidate gene. Differential expression and association analyses were performed. RNA-seq data and qPCR showed that HEP21 was significantly differentially expressed in pre-pubertal and pubertal ovaries. A total of 23 variations were detected in HEP21. Association analyses of single nucleotide polymorphisms (SNPs) in HEP21 and reproductive traits showed that rs315156783 was significantly related to comb height at 84 and 91 days. These results indicate that HEP21 is a candidate gene for sexual maturity in chickens. Our results contribute to a more comprehensive understanding of sexual maturity and reproduction in chickens.

Variants in BMP7 and BMP15 3'-UTRs Associated with Reproductive Traits in a Large White Pig Population

Bone morphogenetic protein 7 (BMP7) and BMP15, which encode members of the BMP family, have been identified by whole-genome resequencing as breeding-related genes that overlap with a known quantitative trait locus for reproductive traits. In this study, we investigated the effects of variants at the BMP7 and BMP15 gene loci on sow reproductive traits. We isolated 669 and 1213 bp sequences of the 3’-untranslated region (3’-UTR) of the porcine BMP7 and BMP15 genes, respectively, and detected several RNA regulatory elements, such as miRNA response elements and AU-rich elements. Pooled DNA sequencing identified two novel point mutations (viz., BMP7 c.1569A>G and BMP15 c.2366G>A) in the 3’-UTR. Association analysis showed that the c.1569A>G polymorphism was associated with the litter weight trait in a Large White pig population. Furthermore, analysis of the combined genetic effects revealed that AA/GA and AG/GG were the favorable combined genotypes for the total number of piglets born (TNB) and the total number of piglets born alive (NBA), whereas. Together, our findings confirm that BMP7 and BMP15 are candidate genes for porcine reproductive performance.

Single-base-resolution methylome of giant panda's brain, liver and pancreatic tissue

The giant panda (Ailuropoda melanoleuca) is one of the most endangered mammals, and its conservation has significant ecosystem and cultural service value. Cytosine DNA methylation (5mC) is a stable epigenetic modification to the genome and has multiple functions such as gene regulation. However, DNA methylome of giant panda and its function have not been reported as of yet. Bisulfite sequencing was performed on a 4-day-old male giant panda's brain, liver and pancreatic tissues. We found that the whole genome methylation level was about 0.05% based on reads normalization and mitochondrial DNA was not methylated. Three tissues showed similar methylation tendency in the protein-coding genes of their genomes, but the brain genome had a higher count of methylated genes. We obtained 467 and 1,013 different methylation regions (DMR) genes in brain vs. pancreas and liver, while only 260 DMR genes were obtained in liver vs pancreas. Some lncRNA were also DMR genes, indicating that methylation may affect biological processes by regulating other epigenetic factors. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analysis indicated that low methylated promoter, high methylated promoter and DMR genes were enriched at some important and tissue-specific items and pathways, like neurogenesis, metabolism and immunity. DNA methylation may drive or maintain tissue specificity and organic functions and it could be a crucial regulating factor for the development of newborn cubs. Our study offers the first insight into giant panda's DNA methylome, laying a foundation for further exploration of the giant panda's epigenetics.

Developmental Programming: Contribution of Epigenetic Enzymes to Antral Follicular Defects in the Sheep Model of PCOS

Prenatal testosterone (T)-treated sheep, similar to women with polycystic ovary syndrome (PCOS), manifest oligo-/anovulation, hyperandrogenism, and polyfollicular ovary. The polyfollicular ovarian morphology, a result of persistence of antral follicles, arises, in part, by transcriptional changes in key mediators of follicular development that, in turn, are driven by epigenetic mechanisms. We hypothesized that prenatal T excess induces, in a cell-specific manner, transcriptional changes in key mediators of follicular development associated with relevant changes in epigenetic machinery. Expression levels of key mediators of follicular development, DNA methyltransferases (DNMTs), and histone de-/methylases and de-/acetylases were determined in laser-capture microdissection-isolated antral follicular granulosa and theca and ovarian stromal cells from 21 months of age control and prenatal T-treated sheep (100 mg IM twice weekly from gestational day 30 to 90; term: 147 days). Changes in histone methylation were determined by immunofluorescence. Prenatal T treatment induced the following: (i) cell-specific changes in gene expression of key mediators of follicular development and steroidogenesis; (ii) granulosa, theca, and stromal cell-specific changes in DNMTs and histone de-/methylases and deacetylases, and (iii) increases in histone 3 trimethylation at lysine 9 in granulosa and histone 3 dimethylation at lysine 4 in theca cells. The pattern of histone methylation was consistent with the expression profile of histone de-/methylases in the respective cells. These findings suggest that changes in expression of key genes involved in the development of the polyfollicular phenotype in prenatal T-treated sheep are mediated, at least in part, by cell-specific changes in epigenetic-modifying enzymes.

Smad4 Feedback Enhances BMPR1B Transcription in Ovine Granulosa Cells

BMPR1B is a type 1B receptor of the canonical bone morphogenetic protein (BMP)/Sma- and mad-related protein (Smad) signaling pathway and is well known as the first major gene associated with sheep prolificacy. However, little is known about the transcriptional regulation of the ovine BMPR1B gene. In this study, we identified the ovine BMPR1B gene promoter and demonstrated that its transcription was regulated by Smad4. In sheep ovarian follicles, three transcriptional variants of BMPR1B gene with distinct transcription start sites were identified using 5′ RACE assay while variants II and III were more strongly expressed. Luciferase assay showed that the region −405 to −200 nt is the PII promoter region of variant II. Interestingly, two putative Smad4-binding elements (SBEs) were detected in this region. Luciferase and ChIP assay revealed that Smad4 enhances PII promoter activity of the ovine BMPR1B gene by directly interacting with SBE1 motif. Furthermore, in the ovine granulosa cells, Smad4 regulated BMPRIB expression, and BMPRIB-mediated granulosa cells apoptosis. Overall, our findings not only characterized the 5’ regulatory region of the ovine BMPR1B gene, but also uncovered a feedback regulatory mechanism of the canonical BMP/Smad signaling pathway and provided an insight into the transcriptional regulation of BMPR1B gene and sheep prolificacy.

Effect of CREB1 promoter non-CpG island methylation on its differential expression profile on sheep ovaries associated with prolificacy

This study aimed to investigate the effect of different methylated regions of cyclic-AMP response element binding protein 1 (CREB1) by comparing the high prolificacy (HP) group and low prolificacy (LP) group, which was detected in our previous study. The expression level of CREB1 mRNA in the ovaries of the HP group was higher than in the LP group (P <  0.05). The differential methylated region (DMR) had 4 methylated CG dinucleotides(CGs): -1546, -1544, -1494 and -1464. The DNA methylation levels of -1546 CGs and -1464 CGs were significantly higher in the HP group than in the LP group (P <  0.05). The activity from -1296 to +26 (without DMR) was significantly higher than the activity from -1598 to +26 (with DMR) (P <  0.05). The result of 5-aza-2'-deoxycytidine treatment indicated that the inhibition DNA methylation of DMR reduced the transcription of CREB1. The bioinformatics predictive analysis were found that the -1546 CG site was located in the CCAAT/enhancer-binding protein alpha (CEBPA) binding site and the -1464 CG site was located in the Sp1 binding site. Finally, this study revealed the relationship between the methylation of non-CpG sites of the promoter and transcription of CREB1. This study will provide a theoretical basis of the Hu sheep ovaries associated with DNA methylation.

Single nucleotide polymorphisms in BMP2 and BMP7 and the association with litter size in Small Tail Han sheep

Although it has been investigated for many years, the physiological processes regulating prolificacy in sheep remains unclear because of regulation by many genes. To better understand the effects of three single nucleotide polymorphisms (SNPs) comprising g.48462350C>T in BMP2, g.58171856C>G and g.58171886A>C in BMP7, a population genetic analysis was conducted using data obtained from genotyping in 768 sheep from six breeds (three polytocous and three monotocous). The results indicate that all the sheep breeds were considered to conform to the Hardy-Weinberg equilibrium (P > 0.05). The associations of these three SNPs with litter size in 384 Small Tail Han sheep were analyzed, therefore, and found to be correlated with fecundity as assessed by mean litter size (P < 0.05). Bioinformatic analysis indicated there was a transmembrane domain change that occurred after a mutation in BMP2 at g.48462350C>T, and changes involving transcription factors such as USF1, USF2 and INMS1 in the BMP7 promoter region might be involved in greater sheep prolificacy.

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

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

Non-CpG Methylation Revised

Textbook and scientific papers addressing DNA methylation usually still cite “DNA methylation occurs at CpG cytosines”. Methylation at cytosines outside the CpG nucleotide, the so-called “non-CpG methylation”, is usually considered a minor and not biologically relevant process. However, the technical improvements and additional studies in epigenetics have demonstrated that non-CpG methylation is present with frequency higher than previously thought and retains biological activity, potentially relevant to the understanding and the treatment of human diseases.

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.