Genome-Wide Identification of Histone Modifications Involved in Placental Development in Pigs

Analysis of histone H3K4me3 modifications in bovine placenta derived from different calf-production methods

In ruminants, the overgrowth of offspring produced by in vitro fertilization (IVF) is a common problem. Abnormal epigenetic modifications caused by environmental factors during the early embryonic period are suspected as an aetiology of overgrowth. In this study, we investigated the genome-wide histone H3K4me3 profiles of bovine placentae that play a pivotal role in foetal development and compared their characteristics between artificial insemination (AI)- and IVF-derived samples. Cotyledons were harvested from the placentae obtained at parturition of 5 AI- and 13 IVF-derived calves, and chromatin immunoprecipitation sequencing was performed for H3K4me3. We confirmed no significant maternal tissue contamination in the samples we used. The revealed H3K4me3 profiles reflected the general characteristics of the H3K4me3 modification, which is abundantly distributed in the promoter region of active genes. By extracting common modifications from multiple samples, the genes involved in placenta-specific biological processes could be enriched. Comparison with the H3K4me3 modifications of blastocyst samples was also effective for enriching the placenta-specific features. Principal component analysis suggested the presence of differential H3K4me3 modifications in AI- and IVF-derived samples. The genes contributing to the difference were related to the developmental biological processes. Imprinted genes such as BEGAIN, ZNF215 and DLX5 were among the extracted genes. Principal component and discriminant analyses using only male samples categorized the samples into three groups based on foetal weight and calf-production methods. To our knowledge, this is the first study to profile the genome-wide histone modifications of bovine foetal placentae and reveal their differential characteristics between different calf-production methods.

Dynamic chromatin architecture of the porcine adipose tissues with weight gain and loss

Using an adult female miniature pig model with diet-induced weight gain/weight loss, we investigated the regulatory mechanisms of three-dimensional (3D) genome architecture in adipose tissues (ATs) associated with obesity. We generated 249 high-resolution in situ Hi-C chromatin contact maps of subcutaneous AT and three visceral ATs, analyzing transcriptomic and chromatin architectural changes under different nutritional treatments. We find that chromatin architecture remodeling underpins transcriptomic divergence in ATs, potentially linked to metabolic risks in obesity development. Analysis of chromatin architecture among subcutaneous ATs of different mammals suggests the presence of transcriptional regulatory divergence that could explain phenotypic, physiological, and functional differences in ATs. Regulatory element conservation analysis in pigs and humans reveals similarities in the regulatory circuitry of genes responsible for the obesity phenotype and identified non-conserved elements in species-specific gene sets that underpin AT specialization. This work provides a data-rich tool for discovering obesity-related regulatory elements in humans and pigs.

Spatio-temporal transcriptome dynamics coordinate rapid transition of core crop functions in 'lactating' pigeon

Pigeons (Columba livia) are among a select few avian species that have developed a specialized reproductive mode wherein the parents produce a 'milk' in their crop to feed newborn squabs. Nonetheless, the transcriptomic dynamics and role in the rapid transition of core crop functions during 'lactation' remain largely unexplored. Here, we generated a de novo pigeon genome assembly to construct a high resolution spatio-temporal transcriptomic landscape of the crop epithelium across the entire breeding stage. This multi-omics analysis identified a set of 'lactation'-related genes involved in lipid and protein metabolism, which contribute to the rapid functional transitions in the crop. Analysis of in situ high-throughput chromatin conformation capture (Hi-C) sequencing revealed extensive reorganization of promoter-enhancer interactions linked to the dynamic expression of these 'lactation'-related genes between stages. Moreover, their expression is spatially localized in specific epithelial layers, and can be correlated with phenotypic changes in the crop. These results illustrate the preferential de novo synthesis of 'milk' lipids and proteins in the crop, and provides candidate enhancer loci for further investigation of the regulatory elements controlling pigeon 'lactation'.

IGF2 May Enhance Placental Fatty Acid Metabolism by Regulating Expression of Fatty Acid Carriers in the Growth of Fetus and Placenta during Late Pregnancy in Pigs

Fatty acids (FAs) are essential substances for the growth and development of the fetus and placenta. The growing fetus and placenta must obtain adequate FAs received from the maternal circulation and facilitated by various placental FA carriers, including FA transport proteins (FATPs), FA translocase (FAT/CD36), and cytoplasmic FA binding proteins (FABPs). Placental nutrition transport was regulated by imprinted genes H19 and insulin-like growth factor 2 (IGF2). Nevertheless, the relationship between the expression patterns of H19/IGF2 and placental fatty acid metabolism throughout pig pregnancy remains poorly studied and unclear. We investigated the placental fatty acid profile, expression patterns of FA carriers, and H19/IGF2 in the placentae on Days 40 (D40), 65 (D65), and 95 (D95) of pregnancy. The results showed that the width of the placental folds and the number of trophoblast cells of D65 placentae were significantly increased than those of D40 placentae. Several important long-chain FAs (LCFAs), including oleic acid, linoleic acid, arachidonatic acid, eicosapentaenoic acid, and docosatetraenoic acid, in the pig placenta showed dramatically increased levels throughout pregnancy. The pig placenta possessed higher expression levels of CD36, FATP4, and FABP5 compared with other FA carriers, and their expression levels had significantly upregulated 2.8-, 5.6-, and 12.0-fold from D40 to D95, respectively. The transcription level of IGF2 was dramatically upregulated and there were corresponding lower DNA methylation levels in the IGF2 DMR2 in D95 placentae relative to D65 placentae. Moreover, in vitro experimentation revealed that the overexpression of IGF2 resulted in a significant increase in fatty acid uptake and expression levels of CD36, FATP4, and FABP5 in PTr2 cells. In conclusion, our results indicate that CD36, FATP4, and FABP5 may be important regulators that enhance the transport of LCFAs in the pig placenta and that IGF2 may be involved in FA metabolism by affecting the FA carriers expression to support the growth of the fetus and placenta during late pregnancy in pigs.

Detailed molecular and epigenetic characterization of the pig IPEC-J2 and chicken SL-29 cell lines

The pig IPEC-J2 and chicken SL-29 cell lines are of interest because of their untransformed nature and wide use in functional studies. Molecular characterization of these cell lines is important to gain insight into possible molecular aberrations. The aim of this paper is to provide a molecular and epigenetic characterization of the IPEC-J2 and SL-29 cell lines, a cell-line reference for the FAANG community, and future biomedical research. Whole genome sequencing, gene expression, DNA methylation, chromatin accessibility, and ChIP-seq of four histone marks (H3K4me1, H3K4me3, H3K27ac, H3K27me3) and an insulator (CTCF) are used to achieve these aims. Heteroploidy (aneuploidy) of various chromosomes was observed from whole genome sequencing analysis in both cell lines. Furthermore, higher gene expression for genes located on chromosomes with aneuploidy in comparison to diploid chromosomes was observed. Regulatory complexity of gene expression, DNA methylation, and chromatin accessibility was investigated through an integrative approach.

Integrated Analysis of DNA Methylation and Gene Expression in Porcine Placental Development

Proper placental development is crucial for the conceptus to grow and survive, because the placenta is responsible for transporting nutrients and oxygen from the pregnant female to the developing fetus. However, the processes of placental morphogenesis and fold formation remain to be fully elucidated. In this study, we used whole-genome bisulfite sequencing and RNA sequencing to produce a global map of DNA methylation and gene expression changes in placentas from Tibetan pig fetuses 21, 28, and 35 days post-coitus. Substantial changes in morphology and histological structures at the uterine-placental interface were revealed via hematoxylin-eosin staining. Transcriptome analysis identified 3959 differentially expressed genes (DEGs) and revealed the key transcriptional properties in three stages. The DNA methylation level in the gene promoter was negatively correlated with gene expression. We identified a set of differentially methylated regions associated with placental developmental genes and transcription factors. The decrease in DNA methylation level in the promoter was associated with the transcriptional activation of 699 DEGs that were functionally enriched in cell adhesion and migration, extracellular matrix remodeling, and angiogenesis. Our analysis provides a valuable resource for understanding the mechanisms of DNA methylation in placental development. The methylation status of different genomic regions plays a key role in establishing transcriptional patterns from placental morphogenesis to fold formation.

Integrating Transcriptomic and ChIP-Seq Reveals Important Regulatory Regions Modulating Gene Expression in Myometrium during Implantation in Pigs

The myometrium is the outer layer of the uterus. Its contraction and steroidogenic activities are required for embryo implantation. However, the molecular mechanisms underlying its functions remain unknown in pigs. The myometrium includes the inner circular muscle (CM) and the outer longitudinal muscle (LM) layers. In this study, we collected the CM and LM samples from the mesometrial side (named M) of the uterus on days 12 (pre-implantation stage) and 15 (implantation stage) of pregnancy and day 15 of the estrous cycle. The transcriptomic results revealed distinct differences between the uterine CM and LM layers in early pregnancy: the genes expressed in the LM layer were mainly related to contraction pathways, whereas the transcriptional signatures in the CM layer on day 15 of pregnancy were primarily involved in the immune response processes. Subsequent comparisons in the CM layer between pregnant and cyclic gilts show that the transcriptional signatures of the CM layer are implantation-dependent. Next, we investigated the genome-wide profiling of histone H3 lysine 27 acetylation (H3K27ac) and histone H3 lysine 4 trimethylation (H3K4me3) in pig uterine CM and LM layers. The genomic regions that had transcriptional activity and were associated with the expression of genes in the two layers were characterized. Taken together, the regulatory regions identified in the study may contribute to modulating the gene expression in pig uterine CM and LM layers during implantation.

Mapping and analysis of a spatiotemporal H3K27ac and gene expression spectrum in pigs

The limited knowledge of genomic noncoding and regulatory regions has restricted our ability to decipher the genetic mechanisms underlying complex traits in pigs. In this study, we characterized the spatiotemporal landscape of putative enhancers and promoters and their target genes by combining H3K27ac-targeted ChIP-Seq and RNA-Seq in fetal (prenatal days 74-75) and adult (postnatal days 132-150) tissues (brain, liver, heart, muscle and small intestine) sampled from Asian aboriginal Bama Xiang and European highly selected Large White pigs of both sexes. We identified 101,290 H3K27ac peaks, marking 18,521 promoters and 82,769 enhancers, including peaks that were active across all tissues and developmental stages (which could indicate safe harbor locus for exogenous gene insertion) and tissue- and developmental stage-specific peaks (which regulate gene pathways matching tissue- and developmental stage-specific physiological functions). We found that H3K27ac and DNA methylation in the promoter region of the XIST gene may be involved in X chromosome inactivation and demonstrated the utility of the present resource for revealing the regulatory patterns of known causal genes and prioritizing candidate causal variants for complex traits in pigs. In addition, we identified an average of 1,124 super-enhancers per sample and found that they were more likely to show tissue-specific activity than ordinary peaks. We have developed a web browser to improve the accessibility of the results ( http://segtp.jxau.edu.cn/pencode/?genome=susScr11 ).

Dynamic 3D genome reorganization during development and metabolic stress of the porcine liver

Liver development is a complex process that is regulated by a series of signaling pathways. Three-dimensional (3D) chromatin architecture plays an important role in transcriptional regulation; nonetheless, its dynamics and role in the rapid transition of core liver functions during development and obesity-induced metabolic stress remain largely unexplored. To investigate the dynamic chromatin architecture during liver development and under metabolic stress, we generated high-resolution maps of chromatin architecture for porcine livers across six major developmental stages (from embryonic day 38 to the adult stage) and under a high-fat diet-induced obesity. The characteristically loose chromatin architecture supports a highly plastic genome organization during early liver development, which fundamentally contributes to the rapid functional transitions in the liver after birth. We reveal the multi-scale reorganization of chromatin architecture and its influence on transcriptional regulation of critical signaling processes during liver development, and show its close association with transition in hepatic functions (i.e., from hematopoiesis in the fetus to metabolism and immunity after birth). The limited changes in chromatin structure help explain the observed metabolic adaptation to excessive energy intake in pigs. These results provide a global overview of chromatin architecture dynamics associated with the transition of physiological liver functions between prenatal development and postnatal maturation, and a foundational resource that allows for future in-depth functional characterization.

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

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

Unraveling three-dimensional chromatin structural dynamics during spermatogonial differentiation

Spermatogonial stem cells (SSCs) are able to undergo both self-renewal and differentiation. Unlike self-renewal, which replenishes the SSC and progenitor pool, differentiation is an irreversible process committing cells to meiosis. Although the preparations for meiotic events in differentiating spermatogonia (Di-SG) are likely to be accompanied by alterations in chromatin structure, the three-dimensional chromatin architectural differences between SSCs and Di-SG, and the higher-order chromatin dynamics during spermatogonial differentiation, have not been systematically investigated. Here, we performed in situ high-throughput chromosome conformation capture, RNA-seq, and chromatin immunoprecipitation-sequencing analyses on porcine undifferentiated spermatogonia (which consist of SSCs and progenitors) and Di-SG. We identified that Di-SG exhibited less compact chromatin structural organization, weakened compartmentalization, and diminished topologically associating domains in comparison with undifferentiated spermatogonia, suggesting that diminished higher-order chromatin architecture in meiotic cells, as shown by recent reports, might be preprogrammed in Di-SG. Our data also revealed that A/B compartments, representing open or closed chromatin regions respectively, and topologically associating domains were related to dynamic gene expression during spermatogonial differentiation. Furthermore, we unraveled the contribution of promoter-enhancer interactions to premeiotic transcriptional regulation, which has not been accomplished in previous studies due to limited cell input and resolution. Together, our study uncovered the three-dimensional chromatin structure of SSCs/progenitors and Di-SG, as well as the interplay between higher-order chromatin architecture and dynamic gene expression during spermatogonial differentiation. These findings provide novel insights into the mechanisms for SSC self-renewal and differentiation and have implications for diagnosis and treatment of male sub-/infertility.

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.

Genes Associated With Chromatin Modification Within the Swine Placenta Are Differentially Expressed Due to Factors Associated With Season

Seasonal reproductive inefficiency is still observed in modern swine facilities. We previously reported global placental methylation activity was reduced from summer breedings and tended to be less from semen collected during cooler periods. The objective of the current study was to evaluate chromatin modification marks within swine placenta in relationship to breeding season, semen collection season, and semen storage. White composite gilts were artificially inseminated in August or January using single-sire semen that was collected during warm or cool periods and stored as either cryopreserved or cooled-extended. Gilts were harvested 45 days post-breeding, and placental samples from the smallest, average, and largest fetus in each litter were collected and stored at −80°C until RNA extraction. An RT² Profiler assay featuring 84 known chromatin modification enzyme targets was performed using placental RNA pooled by litter. Real-time quantitative polymerase chain reaction results were analyzed using the MIXED procedure, and P-values were Hochberg corrected using the MULTTEST procedure in SAS. The complete model included the fixed effects of breeding season (winter or summer), semen collection season (cool or warm), semen storage (cooled-extended or cryopreserved), interactions; boar as repeated effect; and plate as random effect. If interactions were not significant, only the main effects were tested. The genes, ATF2, AURKA, and KDM5B, were different (P < 0.05) by interaction of breeding season, semen collection season, and semen storage. In general, the greatest (P < 0.05) expression was in placentas derived from summer breedings. Expression of AURKA was also influenced by semen collection and storage. Expression of placental KDM5B from winter breedings was also greater (P < 0.05) from semen collected during cool periods. Placental expressions of ASH2L, DNMT3B, ESCO1, HDAC2, ING3, KDM6B, MYSM1, and SMYD3 were greater (P < 0.05) from summer breedings. Increased expressions of known chromatin modification genes, from placentas derived from summer breedings, are likely responsible for differences in gene transcription between summer- or winter-derived placentas.

Transcriptomic and ChIP-seq Integrative Analysis Reveals Important Roles of Epigenetically Regulated lncRNAs in Placental Development in Meishan Pigs

The development of the placental fold, which increases the maternal–fetal interacting surface area, is of primary importance for the growth of the fetus throughout the whole pregnancy. However, the mechanisms involved remain to be fully elucidated. Increasing evidence has revealed that long non-coding RNAs (lncRNAs) are a new class of RNAs with regulatory functions and could be epigenetically regulated by histone modifications. In this study, 141 lncRNAs (including 73 up-regulated and 68 down-regulated lncRNAs) were identified to be differentially expressed in the placentas of pigs during the establishment and expanding stages of placental fold development. The differentially expressed lncRNAs and genes (DElncRNA-DEgene) co-expression network analysis revealed that these differentially expressed lncRNAs (DElncRNAs) were mainly enriched in pathways of cell adhesion, cytoskeleton organization, epithelial cell differentiation and angiogenesis, indicating that the DElncRNAs are related to the major events that occur during placental fold development. In addition, we integrated the RNA-seq (RNA sequencing) data with the ChIP-seq (chromatin immunoprecipitation sequencing) data of H3K4me3/H3K27ac produced from the placental samples of pigs from the two stages (gestational days 50 and 95). The analysis revealed that the changes in H3K4me3 and/or H3K27ac levels were significantly associated with the changes in the expression levels of 37 DElncRNAs. Furthermore, several H3K4me3/H3K27ac-lncRNAs were characterized to be significantly correlated with genes functionally related to placental development. Thus, this study provides new insights into understanding the mechanisms for the placental development of pigs.