Identification and Conservation Analysis of Cis-Regulatory Elements in Pig Liver

Multi-omics analysis reveals critical cis-regulatory roles of transposable elements in livestock genomes

Transposable elements (TEs) are important sources of genetic and regulatory variation, yet their functional roles in domesticated animals remain insufficiently explored. To address this gap, we comprehensively annotated TE types, ages, and distributions in the genomes of pig (Sus scrofa), cattle (Bos taurus), and chicken (Gallus gallus). Our analysis revealed species-specific patterns in TE abundance, amplification, and activity in modern genomes. By integrating transcriptomic and epigenomic data, we explored the impact of specific TE types on cis-regulatory elements (CREs) and constructed a TE expression atlas across five tissues in all three species. Our findings underscored the critical roles of tissue-specific TE expression and chromatin accessibility in regulating tissue-specific biological processes. Most notably, we developed a computational framework to uncover TE-mediated gene regulatory networks (TE-GRNs). Our findings provide valuable insights into the regulatory functions of TEs in livestock and offer a robust approach for studying TE-GRNs in diverse biological contexts.

Differential Regulation of Male-Hormones-Related Enhancers Revealed by Chromatin Accessibility and Transcriptional Profiles in Pig Liver

Surgical castration can effectively avoid boar taint and improve pork quality by removing the synthesis of androstenone in the testis, thereby reducing its deposition in adipose tissue. The expression of genes involved in testis-derived hormone metabolism was altered following surgical castration, but the upstream regulatory factors and underlying mechanism remain unclear. In this study, we systematically profiled chromatin accessibility and transcriptional dynamics in liver tissue of castrated and intact full-sibling Yorkshire pigs. First, we identified 897 differentially expressed genes and 6864 differential accessible regions (DARs) using RNA- and ATAC-seq. By integrating the RNA- and ATAC-seq results, 227 genes were identified, and a significant positive correlation was revealed between differential gene expression and the ATAC-seq signal. We constructed a transcription factor regulatory network after motif analysis of DARs and identified a candidate transcription factor (TF) SP1 that targeted the HSD3B1 gene, which was responsible for the metabolism of androstenone. Subsequently, we annotated DARs by incorporating H3K27ac ChIP-seq data, marking 2234 typical enhancers and 245 super enhancers involved in the regulation of all testis-derived hormones. Among these, four typical enhancers associated with HSD3B1 were identified. Furthermore, an in-depth investigation was conducted on the androstenone-related enhancers, and an androstenone-related mutation was identified in a newfound candidatetypical enhancer (andEN) with dual-luciferase assays. These findings provide further insights into how enhancers function as links between phenotypic and non-coding area variations. The discovery of upstream TF and enhancers of HSD3B1 contributes to understanding the regulatory networks of androstenone metabolism and provides an important foundation for improving pork quality.

The landscape of super-enhancer regulates remote target gene transcription through loop domains in adipose tissue of pig

BACKGROUND

A super-enhancer (SE) is a huge cluster of multiple enhancers that control the key genes for cell identity and function. The rise of advanced chromatin immunoprecipitation sequencing (ChIP-seq) technology such as Cleavage Under Targets and Tagmentation (CUT&Tag) allows more SEs to be discovered. However, SE studies in Luchuan and Duroc pigs are very rare in animal husbandry.

RESULTS

We used the CUT&Tag technique to identify 145 and 378 SEs from the adipose tissues of Luchuan and Duroc pigs, respectively. There were significant differences in the peak coverage ratio of SE peaks in the gene promoter region between the two breeds. Not only that, peak signals at the start and end point of the SE peak profile showed obvious spikes. The proximal target genes of SE were highly expressed compared with the background genes and the typical enhancer target genes. Subsequently, in conjoint analysis with high-throughput chromosome conformation capture sequencing (Hi-C seq) data, we predicted the remote regulatory genes of SE and found that their expression level was related to the distance of SE extended to the loop's anchor, but not the length of loops. According to our prediction model, SEs can maintain promoter accessibility of partial remote target genes through loop domains. Finally, a batch of SEs closely related to fat metabolism traits were obtained by performing a coalition analysis of quantitative trait loci and SE data.

CONCLUSIONS

This work enabled us to obtain hundreds of SEs from Luchuan and Duroc pigs. Our model provides a new method for predicting the SE remote target genes based on loop domains, and to further explore the potential role of super-enhancer in the regulation of fat metabolism.

Experimental Validation and Prediction of Super-Enhancers: Advances and Challenges

Super-enhancers (SEs) are cis-regulatory elements of the human genome that have been widely discussed since the discovery and origin of the term. Super-enhancers have been shown to be strongly associated with the expression of genes crucial for cell differentiation, cell stability maintenance, and tumorigenesis. Our goal was to systematize research studies dedicated to the investigation of structure and functions of super-enhancers as well as to define further perspectives of the field in various applications, such as drug development and clinical use. We overviewed the fundamental studies which provided experimental data on various pathologies and their associations with particular super-enhancers. The analysis of mainstream approaches for SE search and prediction allowed us to accumulate existing data and propose directions for further algorithmic improvements of SEs' reliability levels and efficiency. Thus, here we provide the description of the most robust algorithms such as ROSE, imPROSE, and DEEPSEN and suggest their further use for various research and development tasks. The most promising research direction, which is based on topic and number of published studies, are cancer-associated super-enhancers and prospective SE-targeted therapy strategies, most of which are discussed in this review.

Chromatin Accessibility and Transcriptional Landscape during Inhibition of Salmonella enterica by Lactobacillus reuteri in IPEC-J2 Cells

Lactobacillus reuteri is a probiotic with bacteriostatic effects, which can effectively inhibit the activity of pathogens. However, the molecular mechanism underlying the inhibition of pathogens by L. reuteri in intestinal cells remains unclear. Using the porcine intestinal cell line IPEC-J2 as a model, we combined RNA-seq and ATAC-seq methods to delineate the porcine genome-wide changes in biological processes and chromatin accessibility in IPEC-J2 cells stimulated by Salmonella enterica BNCC186354, as well as L. reuteri ATCC 53608. Overall, we found that many porcine transcripts were altered after S. enterica BNCC186354 treatment, while L. reuteri ATCC 53608 treatment partially restored this alteration, such as salmonella infection and PI3K/AKT and MAPK pathways. Combined analysis of these two datasets revealed that 26 genes with similar trends overlapped between gene expression and chromatin accessibility. In addition, we identified potential host functional transcription factors (TFs), such as GATA1, TAL1, TBP, RUNX1, Gmeb1, Gfi1b, RARA, and RXRG, in IPEC-J2 cells that might play a critical role and are targeted by L. reuteri ATCC 53608. Moreover, we verified that PI3K/AKT, MAPK, and apoptosis pathways are potentially regulated by S. enterica BNCC186354 but restored by L. reuteri ATCC 53608. The PI3K/AKT pathway was activated by L. reuteri ATCC 53608, thereby potentially inhibiting S. enterica BNCC186354 infection. In conclusion, our data provide new insights into the expression pattern of functional genes and the epigenetic alterations in IPEC-J2 cells underlying the bacteriostatic action of L. reuteri ATCC 53608.

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.

A Comparative Analysis of Super-Enhancers and Broad H3K4me3 Domains in Pig, Human, and Mouse Tissues

Super-enhancers (SEs) and broad H3K4me3 domains (BDs) are crucial regulators in the control of tissue identity in human and mouse. However, their features in pig remain largely unknown. In this study, by integrative computational analyses of epigenomic and transcriptomic data, we have characterized SEs and BDs in six pig tissues and analyzed their conservation in comparison with human and mouse tissues. Similar to human and mouse, pig SEs and BDs display higher tissue specificity than their typical counterparts. Genes proximal to SEs and BDs are associated with tissue identity in most tissues. About 55-182 SEs (5-17% in total) and 99-309 BDs (8-16% in total) across pig tissues are considered as functionally conserved elements because they have orthologous SEs and BDs in human and mouse. However, these elements do not necessarily exhibit sequence conservation. The functionally conserved SEs are correlated to tissue identity in majority of pig tissues, while those conserved BDs are linked to tissue identity in a few tissues. Our study provides resources for future gene regulatory studies in pig. It highlights that SEs are more effective in defining tissue identity than BDs, which is contrasting to a previous study. It also provides novel insights on understanding the sequence features of functionally conserved elements.