Pig epiblast-derived pluripotent stem cells are considered to have great potential and broad prospects for human therapeutic model development and livestock breeding. Despite ongoing attempts since the 1990s, no stably defined pig epiblast-derived stem cell line has been established. Here, guided by insights from a large-scale single-cell transcriptome analysis of pig embryos from embryonic day (E) 0 to E14, specifically, the tracing of pluripotency changes during epiblast development, we developed an in vitro culture medium for establishing and maintaining stable pluripotent stem cell lines from pig E10 pregastrulation epiblasts (pgEpiSCs). Enabled by chemical inhibition of WNT-related signaling in combination with growth factors in the FGF/ERK, JAK/STAT3, and Activin/Nodal pathways, pgEpiSCs maintain their pluripotency transcriptome features, similar to those of E10 epiblast cells, and normal karyotypes after more than 240 passages and have the potential to differentiate into three germ layers. Strikingly, ultradeep in situ Hi-C analysis revealed functional impacts of chromatin 3D-spatial associations on the transcriptional regulation of pluripotency marker genes in pgEpiSCs. In practice, we confirmed that pgEpiSCs readily tolerate at least three rounds of successive gene editing and generated cloned gene-edited live piglets. Our findings deliver on the long-anticipated promise of pig pluripotent stem cells and open new avenues for biological research, animal husbandry, and regenerative biomedicine.

To date, research on poultry egg production performance has only been conducted within inter or intra-breed groups, while those combining both inter- and intra-breed groups are lacking. Egg production performance is known to differ markedly between Sichuan white goose (Anser cygnoides) and Landes goose (Anser anser). In order to understand the mechanism of egg production performance in geese, we undertook this study. Here, 18 ovarian stromal samples from both Sichuan white goose and Landes goose at the age of 145 days (3 individuals before egg production initiation for each breed) and 730 days (3 high- and low egg production individuals during non-laying periods for each breed) were collected to reveal the genome-wide expression profiles of ovarian mRNAs and lncRNAs between these two geese breeds at different physiological stages. Briefly, 58, 347, 797, 777, and 881 differentially expressed genes (DEGs) and 56, 24, 154, 105, and 224 differentially expressed long non-coding RNAs (DElncRNAs) were found in LLD vs. HLD (low egg production Landes goose vs. high egg production Landes goose), LSC vs. HSC (low egg production Sichuan White goose vs. high egg production Sichuan white goose), YLD vs. YSC (young Landes goose vs. young Sichuan white goose), HLD vs. HSC (high egg production Landes goose vs. high egg production Sichuan white goose), and LLD vs. LSC (low egg production Landes goose vs. low egg production Sichuan white goose) groups, respectively. Functional enrichment analysis of these DEGs and DElncRNAs suggest that the “neuroactive ligand–receptor interaction pathway” is crucial for egg production, and particularly, members of the 5-hydroxytryptamine receptor (HTR) family affect egg production by regulating ovarian metabolic function. Furthermore, the big differences in the secondary structures among HTR1F and HTR1B, HTR2B, and HTR7 may lead to their different expression patterns in goose ovaries of both inter- and intra-breed groups. These results provide novel insights into the mechanisms regulating poultry egg production performance.

Myofibres (primary and secondary myofibre) are the basic structure of muscle and the determinant of muscle mass. To explore the skeletal muscle developmental processes from primary myofibres to secondary myofibres in pigs, we conducted an integrative three-dimensional structure of genome and transcriptomic characterization of longissimus dorsi muscle of pig from primary myofibre formation stage [embryonic Day 35 (E35)] to secondary myofibre formation stage (E80). In the hierarchical genomic structure, we found that 11.43% of genome switched compartment A/B status, 14.53% of topologically associating domains are changed intradomain interactions (D-scores) and 2,730 genes with differential promoter–enhancer interactions and (or) enhancer activity from E35 to E80. The alterations of genome architecture were found to correlate with expression of genes that play significant roles in neuromuscular junction, embryonic morphogenesis, skeletal muscle development or metabolism, typically, NEFL, MuSK, SLN, Mef2D and GCK. Significantly, Sox6 and MATN2 play important roles in the process of primary to secondary myofibres formation and increase the regulatory potential score and genes expression in it. In brief, we reveal the genomic reorganization from E35 to E80 and construct genome-wide high-resolution interaction maps that provide a resource for studying long-range control of gene expression from E35 to E80.

Background The three-dimensional (3D) architecture of the genome has a highly ordered and hierarchical nature, which influences the regulation of essential nuclear processes at the basis of gene expression, such as gene transcription. While the hierarchical organization of heterochromatin and euchromatin can underlie differences in gene expression that determine evolutionary differences among species, the way 3D genome architecture is affected by evolutionary forces within major lineages remains unclear. Here, we report a comprehensive comparison of 3D genomes, using high resolution Hi-C data in fibroblast cells of fish, chickens, and 10 mammalian species. Results This analysis shows a correlation between genome size and chromosome length that affects chromosome territory (CT) organization in the upper hierarchy of genome architecture, whereas lower hierarchical features, including local transcriptional availability of DNA, are selected through the evolution of vertebrates. Furthermore, conservation of topologically associating domains (TADs) appears strongly associated with the modularity of expression profiles across species. Additionally, LINE and SINE transposable elements likely contribute to heterochromatin and euchromatin organization, respectively, during the evolution of genome architecture. Conclusions Our analysis uncovers organizational features that appear to determine the conservation and transcriptional regulation of functional genes across species. These findings can guide ongoing investigations of genome evolution by extending our understanding of the mechanisms shaping genome architecture.

Folliculogenesis is a complex biological process involving a central oocyte and its surrounding somatic cells. Three-dimensional chromatin architecture is an important transcription regulator; however, little is known about its dynamics and role in transcriptional regulation of granulosa cells during chicken folliculogenesis. We investigate the transcriptomic dynamics of chicken granulosa cells over ten follicular stages and assess the chromatin architecture dynamics and how it influences gene expression in granulosa cells at three key stages: the prehierarchical small white follicles, the first largest preovulatory follicles, and the postovulatory follicles. Our results demonstrate the consistency between the global reprogramming of chromatin architecture and the transcriptomic divergence during folliculogenesis, providing ample evidence for compartmentalization rearrangement, variable organization of topologically associating domains, and rewiring of the long-range interaction between promoter and enhancers. These results provide key insights into avian reproductive biology and provide a foundational dataset for the future in-depth functional characterization of granulosa cells.

Cryopreservation induces sperm cryoinjuries, including physiological and functional changes. However, the molecular mechanisms of sperm cryoinjury and cryoresistance are still unknown. Cryoresistance or the freeze tolerance of sperm varies across species, and boar sperm is more susceptible to cold stress. Contrary to boar sperm, giant panda sperm appears to be strongly freeze-tolerant and is capable of surviving repeated cycles of freeze-thawing. In this study, differentially expressed (DE) PIWI-interacting RNAs (piRNAs) of fresh and frozen-thawed sperm with different freeze tolerance capacity from giant panda and boar were evaluated. The results showed that 1,160 (22 downregulated and 1,138 upregulated) and 384 (110 upregulated and 274 downregulated) DE piRNAs were identified in giant panda and boar sperm, respectively. Gene ontology (GO) enrichment analysis revealed that the target DE messenger RNAs (mRNAs) of DE piRNAs were mainly enriched in biological regulation, cellular, and metabolic processes in giant panda and boar sperm. Moreover, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that the target DE mRNAs of DE piRNAs were only distributed in DNA replication and the cyclic adenosine monophosphate (cAMP) signaling pathway in giant panda, but the cAMP, cyclic guanosine monophosphate (cGMP), and mitogen-activated protein kinase (MAPK) signaling pathways in boar sperm were considered as part of the olfactory transduction pathway. In conclusion, we speculated that the difference in the piRNA profiles and the DE piRNAs involved in the cAMP signaling pathway in boar and giant panda may have contributed to the different freeze tolerance capacities between giant panda and boar sperm, which helps to elucidate the molecular mechanism behind sperm cryoinjury and cryoresistance.

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.

Background A growing body of evidence has revealed that the mammalian genome is organized into hierarchical layers that are closely correlated with and may even be causally linked with variations in gene expression. Recent studies have characterized chromatin organization in various porcine tissues and cell types and compared them among species and during the early development of pigs. However, how chromatin organization differs among pig breeds is poorly understood. Results In this study, we investigated the 3D genome organization and performed transcriptome characterization of two adipose depots (upper layer of backfat [ULB] and greater omentum [GOM]) in wild boars and Bama pigs; the latter is a typical indigenous pig in China. We found that over 95% of the A/B compartments and topologically associating domains (TADs) are stable between wild boars and Bama pigs. In contrast, more than 70% of promoter-enhancer interactions (PEIs) are dynamic and widespread, involving over a thousand genes. Alterations in chromatin structure are associated with changes in the expression of genes that are involved in widespread biological functions such as basic cellular functions, endocrine function, energy metabolism and the immune response. Approximately 95% and 97% of the genes associated with reorganized A/B compartments and PEIs in the two pig breeds differed between GOM and ULB, respectively. Conclusions We reported 3D genome organization in adipose depots from different pig breeds. In a comparison of Bama pigs and wild boar, large-scale compartments and TADs were mostly conserved, while fine-scale PEIs were extensively reorganized. The chromatin architecture in these two pig breeds was reorganized in an adipose depot-specific manner. These results contribute to determining the regulatory mechanism of phenotypic differences between Bama pigs and wild boar.