Genome-wide DNA methylation analysis of the porcine hypothalamus-pituitary-ovary axis

In search of epigenetic hallmarks of different tissues: an integrative omics study of horse liver, lung, and heart

DNA methylation and microRNA (miRNA) expression are epigenetic mechanisms essential for regulating tissue-specific gene expression and metabolic processes. However, high-resolution transcriptome, methylome, or miRNAome data is only available for a few model organisms and selected tissues. Up to date, only a few studies have reported on gene expression, DNA methylation, or miRNA expression in adult equine tissues at the genome-wide level. In the present study, we used RNA-Seq, miRNA-seq, and reduced representation bisulfite sequencing (RRBS) data from the heart, lung, and liver tissues of healthy cold-blooded horses to identify differentially expressed genes (DEGs), differentially expressed miRNA (DE miRNA) and differentially methylated sites (DMSs) between three types of horse tissues. Additionally, based on integrative omics analysis, we described the observed interactions of epigenetic mechanisms with tissue-specific gene expression alterations. The obtained data allowed identification from 4067 to 6143 DMSs, 9733 to 11,263 mRNAs, and 155 to 185 microRNAs, differentially expressed between various tissues. We pointed out specific genes whose expression level displayed a negative correlation with the level of CpG methylation and miRNA expression and revealed biological processes that they enrich. Furthermore, we confirmed and validated the accuracy of the Next-Generation Sequencing (NGS) results with bisulfite sequencing PCR (BSP) and quantitative PCR (qPCR). This comprehensive analysis forms a strong foundation for exploring the epigenetic mechanisms involved in tissue differentiation, especially the growth and development of the equine heart, lungs, and liver.

Epigenetic disruptions in the offspring hypothalamus in response to maternal infection

DNA methylation is an epigenetic modification that can alter gene expression, and the incidence can vary across developmental stages, inflammatory conditions, and sexes. The effects of viral maternal viral infection and sex on the DNA methylation patterns were studied in the hypothalamus of a pig model of immune activation during development. DNA methylation at single-base resolution in regions of high CpG density was measured on 24 individual hypothalamus samples using reduced representation bisulfite sequencing. Differential over- and under-methylated sites were identified and annotated to proximal genes and corresponding biological processes. A total of 120 sites were differentially methylated (FDR-adjusted p-value < 0.05) between maternal infection or sex groups. Among the 66 sites differentially methylated between groups exposed to inflammatory signals and control, most sites were over-methylated in the challenged group and included sites in the promoter regions of genes SIRT3 and NRBP1. Among the 54 differentially methylated sites between females and males, most sites were over-methylated in females and included sites in the promoter region of genes TNC and EIF4G1. The analysis of the genes proximal to the differentially methylated sites suggested that biological processes potentially impacted include immune response, neuron migration and ensheathment, peptide signaling, adaptive thermogenesis, and tissue development. These results suggest that translational studies should consider that the prolonged effect of maternal infection during gestation may be enacted through epigenetic regulatory mechanisms that may differ between sexes.

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.

How Epigenetics Can Enhance Pig Welfare?

Epigenetics works as an interface between the individual and its environment to provide phenotypic plasticity to increase individual adaptation capabilities. Recently, a wide variety of epi-genetic findings have indicated evidence for its application in the development of putative epi-biomarkers of stress in farm animals. The purpose of this study was to evaluate previously reported stress epi-biomarkers in swine and encourage researchers to investigate potential paths for the development of a robust molecular tool for animal welfare certification. In this literature review, we report on the scientific concerns in the swine production chain, the management carried out on the farms, and the potential implications of these practices for the animals' welfare and their epigenome. To assess reported epi-biomarkers, we identified, from previous studies, potentially stress-related genes surrounding epi-biomarkers. With those genes, we carried out a functional enrichment analysis of differentially methylated regions (DMRs) of the DNA of swine subjected to different stress-related conditions (e.g., heat stress, intrauterine insult, and sanitary challenges). We identified potential epi-biomarkers for target analysis, which could be added to the current guidelines and certification schemes to guarantee and certify animal welfare on farms. We believe that this technology may have the power to increase consumers' trust in animal welfare.

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 using next-generation sequencing to reveal candidate genes responsible for boar taint in pigs

Boar taint (BT) is an offensive flavor observed in non-castrated male pigs that reduces the carcass price. Surgical castration effectively avoids the taint but is associated with animal welfare concerns. The functional annotation of farm animal genomes for understanding the biology of complex traits can be used in the selection of breeding animals to achieve favorable phenotypic outcomes. The characterization of pig epigenomes/methylation changes between animals with high and low BT and genome-wide epigenetic markers that can predict BT are lacking. Reduced representation bisulfite sequencing of DNA methylation patterns based on next-generation sequencing is an efficient technology to identify candidate epigenetic biomarkers associated with BT. Three different BT levels were analyzed using reduced representation bisulfite sequencing data to calculate the methylation levels of cytosine and guanine dinucleotide (CpG) sites. The co-analysis of differentially methylated CpG sites identified by this study and differentially expressed genes identified by a previous study found 32 significant co-located genes. The joint analysis of GO terms and pathways revealed that methylation and gene expression of seven candidate genes were associated with BT; in particular, FASN plays a key role in fatty acid biosynthesis, and PEMT might be involved in estrogen regulation and the development of BT. This study is the first to report the genome-wide DNA methylation profiles of BT in pigs using next-generation sequencing and summarize candidate genes associated with epigenetic markers of BT, which could contribute to the understanding of the functional biology of BT traits and selective breeding of pigs against BT based on epigenetic biomarkers.

Dynamic DNA methylation of ovaries during pubertal transition in gilts

Background

In female mammals, the initiation of puberty, coupling with the dramatically morphological changes in ovaries, indicates the sexual and follicular maturation. Previous studies have suggested that the disrupted DNA methylation results in the delayed puberty. However, to date, the changes in ovarian methylomes during pubertal transition have not been investigated. In this study, using gilts as a pubertal model, the genome-wide DNA methylation were profiled to explore their dynamics during pubertal transition across Pre-, In- and Post-puberty.

Results

During pubertal transition, the follicles underwent maturation and luteinization, coupled with the significant changes in the mRNA expression of DNMT1 and DNMT3a. DNA methylation levels of In-puberty were higher than that of Pre- and Post-puberty at the locations of genes and CpG islands (CGIs). Analysis of the DNA methylation changes identified 12,313, 20,960 and 17,694 differentially methylated CpGs (DMCs) for the comparisons of Pre- vs. In-, In vs. Post-, and Pre- vs. Post-puberty, respectively. Moreover, the CGIs, upstream and exonic regions showed a significant underrepresentation of DMCs, but the CGI shores, CGI shelves, intronic, downstream and intergenic regions showed a significant overrepresentation of DMCs. Furthermore, biological functions of these methylation changes enriched in PI3K-Akt signaling pathway, GnRH signaling pathway, and Insulin secretion, and the mRNA expressions of several genes of these signaling pathway, including MMP2, ESR1, GSK3B, FGF21, IGF1R, and TAC3, were significantly changed across Pre-, In- and Post-puberty in ovaries.

Conclusions

During pubertal transition in gilts, the DNA methylation changes of ovaries were likely to affect the transcription of genes related to PI3K-Akt signaling pathway, GnRH signaling pathway, and Insulin secretion. These observations can provide new insight into the epigenetic mechanism of follicular and sexual maturation during pubertal transition in mammals.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5884-x) contains supplementary material, which is available to authorized users.

An Epigenome-Wide DNA Methylation Map of Testis in Pigs for Study of Complex Traits

Epigenetic changes are important for understanding complex trait variation and inheritance in pigs that are also a valuable biomedical model for human health research. Testis is the main organ for reproduction and boar taint in pigs; however, there have been no studies to-date on adult pig testis epigenome. The main objective of this study was to establish a genome-wide DNA methylation map of pig testis that would help identify candidate epigenetic biomarkers and methylated genes for complex traits such as male reproduction, fertility or boar taint. Reduced Representation Bisulfite Sequencing (RRBS) was used to study methylation levels of cytosine in nine pig testis samples. The results showed that genome-wide methylation status of nine samples overlapped greatly and their variation among pigs were low. The methylation levels of promoter, exon, intron, cytosine and guanine dinucleotide (CpG) islands and CpG island shores regions were 0.15, 0.47, 0.55, 0.39, and 0.53, respectively. Cytosines binding to CpG islands showed different methylation levels between exon and intron regions. All methylation levels of CpG islands were lower than CpG island shores in different genic features. The distribution of 12,738 differentially methylated cytosines (DMCs) within CpG islands, CpG island shores and other regions was 36.86, 21.65, and 41.49%, respectively, and was 0.33, 1.71, 5.95, and 92.01% in promoter, exon, intron and intergenic regions, respectively. Methylation levels of DMCs in promoter, exon and intron regions were significantly different between CpG islands and CpG island shores (P < 0.05). A total of 898 genes with 2089 DMCs were enriched in 112 Gene Ontology (GO) terms. Fifteen methylated genes from our study were associated with fertility or boar taint traits. Our analysis revealed the methylation patterns in different genic features and CpG island regions of testis in pigs, and summarized several candidate genes associated with DMCs and the involved GO terms. These findings are helpful to understand the relationship between DNA methylation and genic CpG islands, to provide candidate epigenetic regions or biomarkers for pig production and welfare and for translational epigenomic studies that use pigs as an animal model for human research.

Genome-Wide DNA Methylation Analysis of Hypothalamus During the Onset of Puberty in Gilts

Although selection of the early age at puberty in gilts will make for a favorable effect on the reproductivity of sow, a large proportion of phenotypic variation in age at puberty of gilts cannot be explained by genetics. Previous studies have implicated hypothalamic DNA methylation in the onset of puberty in mammals. However, the underlying molecular mechanism regarding the regulation of the onset of puberty has remained largely unexplored in gilts. Herein, the genome-scale DNA methylation of hypothalamus was acquired, using the reduced representation bisulfite sequencing, to compare and describe the changes of DNA methylation across Pre-, In- and Post-pubertal gilts. In this study, the average methylation levels of CpGs and CpHs (where H = C, T, or A) in CpG islands- and gene-related regions were gradually decreased in hypothalamic methylomes during the pubertal transition. Comparisons of Pre- vs. In-, In- vs. Post-, and Pre- vs. Post-pubertal stage revealed that there were 85726, 92914, and 100421 differentially methylated CpGs and 5940, 14804, and 16893 differentially methylated CpHs (where H = C, T, or A) in the hypothalamic methylomes. The methylation changes of CpHs were more dynamic than that of CpGs, and methylation changes of CpGs and CpHs were likely to be, respectively, involved in the developmental processes of reproduction and the molecular processes of cellular communications in the hypothalamus. Moreover, methylation changes of CpHs were observed to overrepresent in the quantitative trait loci of age at puberty, and the biological function of these CpH methylation changes was enriched in the pancreas development in gilts. Furthermore, the mRNA levels of several differentially CpG or CpH methylated genes related to the transcription of RNA II polymerase, GnRH signaling pathway, Estrogen signaling pathway, PI3K-AKt signaling pathway, and Insulin signaling pathway, including MAX, MMP2, FGF11, IGF1R, FGF21, and GSK3B, were significantly changed across these pubertal stages in the hypothalamus. These results will help our understanding of how DNA methylation contributes to phenotypic variation of age at puberty.

Genome-wide DNA methylation analysis of pituitaries during the initiation of puberty in gilts

It has been widely recognized that the early or delayed puberty appears to display harmful effects on adult health outcomes. During the timing of puberty, pituitaries responds to the hypothalamus and then introduce the following response of ovaries in hypothalamic-pituitary-gonadal axis. DNA methylation has been recently suggested to regulate the onset of puberty in female mammals. However, to date, the changes of DNA methylation in pituitaries have not been investigated during pubertal transition. In this study, using gilts as the pubertal model, the genome-scale DNA methylation of pituitaries was profiled and compared across Pre-, In- and Post-puberty by using the reduced representation bisulfite sequencing. We found that average methylation levels of each genomic feature in Post- were lower than Pre- and In-pubertal stage in CpG context, but they were higher in In- than that in Pre- and Post-pubertal stage in CpH (where H = A, T, or C) context. The methylation patterns of CpHs were more dynamic than that of CpGs at the location of high CpG content, low CpG content promoter genes, and differently genomic CGIs. Furthermore, the differently genomic CGIs were likely to show in a similar manner in CpG context but display in a stage-specific manner in the CpH context across the Pre-, In- and Post-pubertal stage. Among these pubertal stages, 5 kb upstream regions of the transcription start sites were protected from both CpG and CpH methylation changes. 12.65% of detected CpGs were identified as the differentially methylated CpGs, regarding 4301 genes which were involved in the fundamental functions of pituitaries. 0.35% of detected CpHs were identified as differentially methylated CpHs, regarding 3691 genes which were involved in the biological functions of releasing gonadotropin hormones. These observations and analyses would provide valuable insights into epigenetic mechanism of the initiation of puberty in pituitary level.

Methylation status of KLF4 and HS3ST2 genes as predictors of endometrial cancer and hyperplastic endometrial lesions

Endometrial carcinoma is one of the most common tumours in developed countries. In addition to the active role of genetic factors, epigenetic changes also have an important effect. The present study analysed the methylation status of kruppel like factor 4 (KLF4) and heparan sulfate‑glucosamine 3‑sulfotransferase 2 (HS3ST2) genes in three endometrial tissue types for carcinoma prediction. The sample comprised 91 women with histologically‑confirmed endometrial carcinoma (64.16±9.64 years old), 36 women with hyperplasia (53.39±9.64 years old) and 45 with no signs or symptoms of malignancy (48.53±11.11 years old). The CpG dinucleotide methylation levels were examined by quantitative pyrosequencing, and the discrimination accuracy of the model was calculated using the Random Forest classification algorithm of the area under the ROC curve (AUC). The mean values of KLF4 and HS3ST2 methylation indices were 23.83±11.39 and 8.52±2.57 in the control samples; 30.40±8.52 and 33.76±20.66 in hyperplasia and 34.72±10.79 and 34.49±18.39 in the cancerous tissues. Multinomial logistic regression indicated that the HS3ST2 CpG1 methylation status is a predictor of hyperplasia (P<0.05) and that the KLF4 CpG2 dinucleotide can predict carcinoma formation (P<0.001). The AUC value of 0.95 indicates high discrimination accuracy of the CpG nucleotides methylation status model between the controls and the two other diagnoses. The results of the present study establish the likelihood that aberrations in KLF4 and HS3ST2 gene methylation levels are important in the development of endometrial hyperplasia and carcinoma, with hyperplasia an intermediate step between healthy and tumour tissues.

Reduced representation bisulfite sequencing (RRBS) of dairy goat mammary glands reveals DNA methylation profiles of integrated genome-wide and critical milk-related genes

DNA methylation (DNAm), a major element of epigenetics, plays critical roles in individual development. Reduced representation bisulfite sequencing (RRBS) is an effective and economical method for analyzing the DNA methylation of a single base. The aims of this study were to determine the DNAm profiles of the methylation contexts (CGs and non-CGs) of lactation and dry periods of goat mammary glands using the RRBS, and to identify potential milk-related genes. The proportion of CG was the highest among all the sequence contexts. The highest CG levels (72.44% to 75.24%) occurred in the 3′ UTR region, followed by the gene body region (61.14% to 65.45%). The non-CG levels were low compared to the CG levels. Bioinformatic analysis demonstrated that the CGs were mainly enriched at high methylation levels (>90%), while non-CGs were enriched at low methylation levels. Methylation levels of 95 and 54 genes in the lactation period were up- or downregulated, respectively, relative to the dry period, such as PPARα, RXRα and NPY genes. The bisulfite sequencing PCR results showed that the methylation level of goat PPARα gene during the lactation period was significant lower than in the dry period, while the methylation level of the RXRα gene was lower in the dry period than in the lactation period. Meanwhile, the methylation levels of human PPARα and NPY genes were significantly higher in MCF-7 than in MCF-10A cells. These findings provide essential information for DNA methylation profiles of goat mammary gland and detect some potential milk-related genes in dairy goats.