Building Haplotype-Resolved 3D Genome Maps of Chicken Skeletal Muscle.

Haplotype-resolved 3D chromatin architecture related to allelic differences in avian skeletal muscle development has not been addressed so far, although chicken husbandry for meat consumption has been prevalent feature of cultures on every continent for more than thousands of years. Here, high-resolution Hi-C diploid maps (1.2-kb maximum resolution) are generated for skeletal muscle tissues in chicken across three developmental stages (embryonic day 15 to day 30 post-hatching). The sequence features governing spatial arrangement of chromosomes and characterize homolog pairing in the nucleus, are identified. Multi-scale characterization of chromatin reorganization between stages from myogenesis in the fetus to myofiber hypertrophy after hatching show concordant changes in transcriptional regulation by relevant signaling pathways. Further interrogation of parent-of-origin-specific chromatin conformation supported that genomic imprinting is absent in birds. This study also reveals promoter-enhancer interaction (PEI) differences between broiler and layer haplotypes in skeletal muscle development-related genes are related to genetic variation between breeds, however, only a minority of breed-specific variations likely contribute to phenotypic divergence in skeletal muscle potentially via allelic PEI rewiring. Beyond defining the haplotype-specific 3D chromatin architecture in chicken, this study provides a rich resource for investigating allelic regulatory divergence among chicken breeds.

Haplotype-resolved 3D chromatin architecture of the hybrid pig.

In diploid mammals, allele-specific three-dimensional (3D) genome architecture may lead to imbalanced gene expression. Through ultradeep in situ Hi-C sequencing of three representative somatic tissues (liver, skeletal muscle, and brain) from hybrid pigs generated by reciprocal crosses of phenotypically and physiologically divergent Berkshire and Tibetan pigs, we uncover extensive chromatin reorganization between homologous chromosomes across multiple scales. Haplotype-based interrogation of multi-omic data revealed the tissue dependence of 3D chromatin conformation, suggesting that parent-of-origin-specific conformation may drive gene imprinting. We quantify the effects of genetic variations and histone modifications on allelic differences of long-range promoter-enhancer contacts, which likely contribute to the phenotypic differences between the parental pig breeds. We also observe the fine structure of somatically paired homologous chromosomes in the pig genome, which has a functional implication genome-wide. This work illustrates how allele-specific chromatin architecture facilitates concomitant shifts in allele-biased gene expression, as well as the possible consequential phenotypic changes in mammals.

The ESCRT protein CHMP5 promotes T cell leukemia by controlling BRD4-p300-dependent transcription.

Oncogene activity rewires cellular transcription, creating new transcription networks to which cancer cells become addicted, by mechanisms that are still poorly understood. Using human and mouse models of T cell acute lymphoblastic leukemia (T-ALL), we identify an essential nuclear role for CHMP5, a cytoplasmic endosomal sorting complex required for transport (ESCRT) protein, in establishing and maintaining the T-ALL transcriptional program. Nuclear CHMP5 promoted the T-ALL gene program by augmenting recruitment of the co-activator BRD4 by the histone acetyl transferase p300 selectively at enhancers and super-enhancers, an interaction that potentiated H3K27 acetylation at these regulatory enhancers. Consequently, loss of CHMP5 diminished BRD4 occupancy at enhancers and super-enhancers and impaired RNA polymerase II pause release, which resulted in downregulation of key T-ALL genes, notably MYC. Reinforcing its importance in T-ALL pathogenesis, CHMP5 deficiency mitigated chemoresistance in human T-ALL cells and abrogated T-ALL induction by oncogenic NOTCH1 in vivo. Thus, the ESCRT protein CHMP5 is an essential positive regulator of the transcriptional machinery promoting T-ALL disease.

Comparative evaluation of the fecal microbiota of adult hybrid pigs and Tibetan pigs, and dynamic changes in the fecal microbiota of hybrid pigs.

The breed of pig can affect the diversity and composition of fecal microbiota, but there is a lack of research on the fecal microbiota of hybrid pigs. In this study, feces samples from Chuanxiang black pigs (a hybrid of Tibetan and Duroc pigs) aged 3 days (n = 24), 70 days (n = 31), 10 months (n = 13) and 2 years (n = 30) and Tibetan pigs aged 10 months (n = 14) and 2 years (n = 15) were collected and sequenced by 16S rRNA gene sequencing technology. We also measured the weight of all the tested pigs and found that the 10-month-old and two-year-old Chuanxiang black pigs weighed about three times the weight of Tibetan pigs of the same age. After comparing the genus-level microbiota composition of Tibetan pigs and Chuanxiang black pigs at 10 months and two years of age, we found that Treponema and Streptococcus were the two most abundant bacteria in Chuanxiang black pigs, while Treponema and Chirstensenellaceae_R.7_group were the two most abundant bacteria in Tibetan pigs. Prediction of microbial community function in adult Chuanxiang black pigs and Tibetan pigs showed changes in nutrient absorption, disease resistance, and coarse feeding tolerance. In addition, we also studied the changes in fecal microbiota in Chuanxiang black pigs at 3 days, 70 days, 10 months, and 2 years of age. We found that the ecologically dominant bacteria in fecal microbiota of Chuanxiang black pigs changed across developmental stages. For example, the highest relative abundance of 70-day-old Chuanxiang black pigs at the genus level was Prevotella. We identified specific microbiota with high abundance at different ages for Chuanxiang black pigs, and revealed that the potential functions of these specific microbiota were related to the dominant phenotype such as fast growth rate and strong disease resistance. Our findings help to expand the understanding of the fecal microbiota of hybrid pigs and provide a reference for future breeding and management of hybrid pigs.

Revisiting Assessment of Computational Methods for Hi-C Data Analysis.

The performances of algorithms for Hi-C data preprocessing, the identification of topologically associating domains, and the detection of chromatin interactions and promoter-enhancer interactions have been mostly evaluated using semi-quantitative or synthetic data approaches, without utilizing the most recent methods, since 2017. In this study, we comprehensively evaluated 24 popular state-of-the-art methods for the complete end-to-end pipeline of Hi-C data analysis, using manually curated or experimentally validated benchmark datasets, including a CRISPR dataset for promoter-enhancer interaction validation. Our results indicate that, although no single method exhibited superior performance in all situations, HiC-Pro, DomainCaller, and Fit-Hi-C2 showed relatively balanced performances of most evaluation metrics for preprocessing, topologically associating domain identification, and chromatin interaction/promoter-enhancer interaction detection, respectively. The comprehensive comparison presented in this manuscript provides a reference for researchers to choose Hi-C analysis tools that best suit their needs.

Milletdb: a multi-omics database to accelerate the research of functional genomics and molecular breeding of millets.

Millets are a class of nutrient-rich coarse cereals with high resistance to abiotic stress; thus, they guarantee food security for people living in areas with extreme climatic conditions and provide stress-related genetic resources for other crops. However, no platform is available to provide a comprehensive and systematic multi-omics analysis for millets, which seriously hinders the mining of stress-related genes and the molecular breeding of millets. Here, a free, web-accessible, user-friendly millets multi-omics database platform (Milletdb, http://milletdb.novogene.com) has been developed. The Milletdb contains six millets and their one related species genomes, graph-based pan-genomics of pearl millet, and stress-related multi-omics data, which enable Milletdb to be the most complete millets multi-omics database available. We stored GWAS (genome-wide association study) results of 20 yield-related trait data obtained under three environmental conditions [field (no stress), early drought and late drought] for 2 years in the database, allowing users to identify stress-related genes that support yield improvement. Milletdb can simplify the functional genomics analysis of millets by providing users with 20 different tools (e.g., 'Gene mapping', 'Co-expression', 'KEGG/GO Enrichment' analysis, etc.). On the Milletdb platform, a gene PMA1G03779.1 was identified through 'GWAS', which has the potential to modulate yield and respond to different environmental stresses. Using the tools provided by Milletdb, we found that the stress-related PLATZs TFs (transcription factors) family expands in 87.5% of millet accessions and contributes to vegetative growth and abiotic stress responses. Milletdb can effectively serve researchers in the mining of key genes, genome editing and molecular breeding of millets.

A mutator-derived prognostic eRNA signature provides insight into the pathogenesis of breast cancer.

Abundant evidence suggests that enhancer RNA (eRNA) is closely related to tumorigenesis, and the role of eRNA transcription in promoting genomic instability in cancers is gradually unveiled. However, research on the evaluation of the prognostic value and molecular mechanisms of genomic instability associated eRNAs in breast cancer is long overdue. Here, we integratively analyzed eRNA expression and somatic mutation profiles in breast cancer genome. We identified genomic instability associated eRNAs and developed a prognostic signature based on these eRNAs with the area under the curve (AUC) around 0.8 at 9-year survival. We further found the prognostic value of this signature is independent of common clinical factors and is better than TP53 status. Higher expression of genomic instability associated genes in the high-risk group was observed, suggesting that this eRNA signature may serve as an indicator of genomic instability in breast cancer. We found prognostic eRNA co-expressed genes are mainly enriched in Gene set 'Breast Cancer 8P12-P11 Amplicon', Gene set 'Metabolism of lipids' and GO process 'Ubiquitin protein ligase binding'. Furthermore, 11 eRNA-signature regulated genes are identified by assessing promoter-enhancer interaction. Among these genes, F11R, BHLHE40, and NECTIN4 are previously reported oncogenes and EGOT is a tumor suppressor gene, indicating the direct roles of eRNAs in tumorigenesis.

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'.

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.

Comparative three-dimensional genome architectures of adipose tissues provide insight into human-specific regulation of metabolic homeostasis.

Elucidating the regulatory mechanisms of human adipose tissues (ATs) evolution is essential for understanding human-specific metabolic regulation, but the functional importance and evolutionary dynamics of three-dimensional (3D) genome organizations of ATs are not well defined. Here, we compared the 3D genome architectures of anatomically distinct ATs from humans and six representative mammalian models. We recognized evolutionarily conserved and human-specific chromatin conformation in ATs at multiple scales, including compartmentalization, topologically associating domain (TAD), and promoter-enhancer interactions (PEI), which have not been described previously. We found PEI are much more evolutionarily dynamic with respect to compartmentalization and topologically associating domain. Compared to conserved PEIs, human-specific PEIs are enriched for human-specific sequence, and the binding motifs of their potential mediators (transcription factors) are less conserved. Our data also demonstrated that genes involved in the evolutionary dynamics of chromatin organization have weaker transcriptional conservation than those associated with conserved chromatin organization. Furthermore, the genes involved in energy metabolism and the maintenance of metabolic homeostasis are enriched in human-specific chromatin organization, while housekeeping genes, health-related genes, and genetic variations are enriched in evolutionarily conserved compared to human-specific chromatin organization. Finally, we showed extensively divergent human-specific 3D genome organizations among one subcutaneous and three visceral ATs. Together, these findings provide a global overview of 3D genome architecture dynamics between ATs from human and mammalian models and new insights into understanding the regulatory evolution of human ATs.

Analysis of mRNA and lncRNA Expression Profiles of Breast Muscle during Pigeon (Columbalivia) Development.

The breast muscle is essential for flight and determines the meat yield and quality of the meat type in pigeons. At present, studies about long non-coding RNA (lncRNA) expression profiles in skeletal muscles across the postnatal development of pigeons have not been reported. Here, we used transcriptome sequencing to examine the White-King pigeon breast muscle at four different ages (1 day, 14 days, 28 days, and 2 years old). We identified 12,918 mRNAs and 9158 lncRNAs (5492 known lncRNAs and 3666 novel lncRNAs) in the breast muscle, and 7352 mRNAs and 4494 lncRNAs were differentially expressed in the process of development. We found that highly expressed mRNAs were mainly related to cell-basic and muscle-specific functions. Differential expression and time-series analysis showed that differentially expressed genes were primarily associated with muscle development and functions, blood vessel development, cell cycle, and energy metabolism. To further predict the possible role of lncRNAs, we also conducted the WGCNA and trans/cis analyses. We found that differentially expressed lncRNAs such as lncRNA-LOC102093252, lncRNA-G12653, lncRNA-LOC110357465, lncRNA-G14790, and lncRNA-LOC110360188 might respectively target UBE2B, Pax7, AGTR2, HDAC1, Sox8 and participate in the development of the muscle. Our study provides a valuable resource for studying the lncRNAs and mRNAs of pigeon muscles and for improving the understanding of molecular mechanisms in muscle development.

The Analysis of Transcriptomes and Microorganisms Reveals Differences between the Intestinal Segments of Guinea Pigs.

The intestine is a tubular organ with multiple functions such as digestion absorption and immunity, but the functions of each intestinal segments are different. Intestinal regionalization is necessary for normal physiological function, but it also means the research results obtained at specific sites may not be applicable to other intestinal segments. In order to comprehensively describe the functional changes in the intestine, different intestinal segments and their contents (duodenum, jejunum, ileum, cecum, colon, and rectum) of guinea pigs were collected for RNA seq and 16S rRNA seq, respectively. The results showed differential genes of each intestinal segment mainly involve mucosa, digestion, absorption, and immunity. The gene sets related to fat, bill salts, vitamins, aggregates, amino acids, and water absorption were highly expressed in the small intestine, and the gene sets related to metal ions, nucleotides, and SCFAs were highly expressed in the large intestine. In terms of immunity, the CD8+ T, Th1, eosinophils, pDCs, and natural killer (NK) T cells in the small intestine showed higher scores than those in the large intestine, while the pattern-recognition receptor signaling pathway-related genes are highly expressed in the large intestine. In terms of microbial composition, Proteobacteria and Actinobacteria are abundant in the small intestine, while Firmicutes and Spirochaete are abundant in large intestine. The correlation analysis showed a high correlation between intestinal microorganisms and gene modules related to digestion and absorption. In addition, cross-species analysis showed the SCFA metabolism gene expression trends in human and rodent intestine were different. In conclusion, we analyzed the changes in substance transport, immune and microbial composition between different intestinal segments of guinea pigs, and explored the relationship between intestinal transcriptome and microorganisms, our research will provides a reference for subsequent intestinal-related research.

Transcriptome analysis revealed the roles of long non-coding RNA and mRNA in the bursa of Fabricius during pigeon (Columba livia) development.

The bursa of Fabricius (BF) is the critical humoral immune organ to birds, playing an essential role in B lymphocyte differentiation. However, unlike other poultries, surgical removal of pigeon BF did not limit humoral immune responsiveness. To investigate the expression profiles and the potential role of mRNA and long non-coding RNA (LncRNA) in squab BFs, transcriptome analysis was performed by RNA-Sequencing (RNA-Seq) over three developmental stages (1-day, 13 and 26 days old). We identified 13,072 mRNAs and 19,129 lncRNAs, of which 2,752 mRNAs and 1,515 lncRNAs were differential expressed (DE) in pigeon BFs over three developmental stages. Cluster analysis presented different expression patterns in DE mRNAs and lncRNAs. Functional enrichment analysis revealed that DE lncRNAs and mRNAs with distinct expression patterns might play crucial roles in the immune system process and tissue morphogenesis. In particular, some DE genes and lncRNAs with higher expression levels in 13D or 26D are related to lymphocyte activation and differentiation, adaptive immune response, positive regulation of immune response, leukocyte migration, etc. Protein-protein interaction (PPI) network and Molecular Complex Detection (MCODE) analysis sreened six significant modules containing 37 genes from immune-related DE gene cluster, which is closely linked in B cell activation, lymphocyte differentiation, B cell receptor signaling pathway, etc. Our study characterizes mRNA and lncRNA transcriptomic variability in pigeon BFs over different developmental stages and enhances understanding of the mechanisms underlying physiological functions of pigeon BF.

Exploring high-resolution chromatin interaction changes and functional enhancers of myogenic marker genes during myogenic differentiation.

Skeletal muscle differentiation (myogenesis) is a complex and highly coordinated biological process regulated by a series of myogenic marker genes. Chromatin interactions between gene's promoters and their enhancers have an important role in transcriptional control. However, the high-resolution chromatin interactions of myogenic genes and their functional enhancers during myogenesis remain largely unclear. Here, we used circularized chromosome conformation capture coupled with next generation sequencing (4C-seq) to investigate eight myogenic marker genes in C2C12 myoblasts (C2C12-MBs) and C2C12 myotubes (C2C12-MTs). We revealed dynamic chromatin interactions of these marker genes during differentiation and identified 163 and 314 significant interaction sites (SISs) in C2C12-MBs and C2C12-MTs, respectively. The interacting genes of SISs in C2C12-MTs were mainly involved in muscle development, and histone modifications of the SISs changed during differentiation. Through functional genomic screening, we also identified 25 and 41 putative active enhancers in C2C12-MBs and C2C12-MTs, respectively. Using luciferase reporter assays for putative enhancers of Myog and Myh3, we identified eight activating enhancers. Furthermore, dCas9-KRAB epigenome editing and RNA-Seq revealed a role for Myog enhancers in the regulation of Myog expression and myogenic differentiation in the native genomic context. Taken together, this study lays the groundwork for understanding 3D chromatin interaction changes of myogenic genes during myogenesis and provides insights that contribute to our understanding of the role of enhancers in regulating myogenesis.

Identification of enhancers responsible for the coordinated expression of myosin heavy chain isoforms in skeletal muscle.

BACKGROUND: Skeletal muscles consist of fibers of differing contractility and metabolic properties, which are primarily determined by the content of myosin heavy chain (MYH) isoforms (MYH7, MYH2, MYH1, and MYH4). The regulation of Myh genes transcription depends on three-dimensional chromatin conformation interaction, but the mechanistic details remain to be determined. RESULTS: In this study, we characterized the interaction profiles of Myh genes using 4C-seq (circular chromosome conformation capture coupled to high-throughput sequencing). The interaction profile of Myh genes changed between fast quadriceps and slow soleus muscles. Combining chromatin immunoprecipitation-sequencing (ChIP-seq) and transposase accessible chromatin with high-throughput sequencing (ATAC-seq), we found that a 38 kb intergenic region interacting simultaneously with fast Myh genes promoters controlled the coordinated expression of fast Myh genes. We also identified four active enhancers of Myh7, and revealed that binding of MYOG and MYOD increased the activity of Myh7 enhancers. CONCLUSIONS: This study provides new insight into the chromatin interactions that regulate Myh genes expression.

Chromatin accessibility landscape of stromal subpopulations reveals distinct metabolic and inflammatory features of porcine subcutaneous and visceral adipose tissue.

Background: Fat accumulation in visceral adipose tissue (VAT) confers increased risk for metabolic disorders of obesity, whereas accumulation of subcutaneous adipose tissue (SAT) is associated with lower risk and may be protective. Previous studies have shed light on the gene expression profile differences between SAT and VAT; however, the chromatin accessibility landscape differences and how the cis-regulatory elements govern gene expression changes between SAT and VAT are unknown. Methods: Pig were used to characterize the differences in chromatin accessibility between the two adipose depots-derived stromal vascular fractions (SVFs) using DNase-sequencing (DNase-seq). Using integrated data from DNase-seq, H3K27ac ChIP-sequencing (ChIP-seq), and RNA-sequencing (RNA-seq), we investigated how the regulatory locus complexity regulated gene expression changes between SAT and VAT and the possible impact that these changes may have on the different biological functions of these two adipose depots. Results: SVFs form SAT and VAT (S-SVF and V-SVF) have differential chromatin accessibility landscapes. The differential DNase I hypersensitive site (DHS)-associated genes, which indicate dynamic chromatin accessibility, were mainly involved in metabolic processes and inflammatory responses. Additionally, the Krüppel-like factor family of transcription factors were enriched in the differential DHSs. Furthermore, the chromatin accessibility data were highly associated with differential gene expression as indicated using H3K27ac ChIP-seq and RNA-seq data, supporting the validity of the differential gene expression determined using DNase-seq. Moreover, by combining epigenetic and transcriptomic data, we identified two candidate genes, NR1D1 and CRYM, could be crucial to regulate distinct metabolic and inflammatory characteristics between SAT and VAT. Together, these results uncovered differences in the transcription regulatory network and enriched the mechanistic understanding of the different biological functions between SAT and VAT.

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.

Comparative 3D genome architecture in vertebrates.

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.

Dynamic transcriptome and chromatin architecture in granulosa cells during chicken folliculogenesis.

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.

Generation and characterization of stable pig pregastrulation epiblast stem cell lines.

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.