ABSTRACT
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.
Subject(s)
Chromatin , Chromosomes , Animals , Swine/genetics , Chromatin/genetics , Haplotypes , Chromosomes/genetics , Genome , Mammals/geneticsABSTRACT
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'.
Subject(s)
Columbidae , Transcriptome , Animals , Female , Transcriptome/genetics , Columbidae/genetics , Columbidae/metabolism , Gene Expression Profiling , Milk , LactationABSTRACT
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.
Subject(s)
Adipose Tissue , Chromatin , Genome , Animals , Humans , Chromatin/genetics , Chromatin Assembly and Disassembly , Genomics , Homeostasis , Mammals , Adipose Tissue/metabolismABSTRACT
Lacking PTRF (polymerase I and transcript release factor), an essential caveolae component, causes a secondary deficiency of caveolins resulting in muscular dystrophy. The transcriptome responses of different types of muscle fibers and mononuclear cells in skeletal muscle to muscular dystrophy caused by Ptrf deletion have not been explored. Here, we created muscular dystrophy mice by Ptrf knockout and applied single-nucleus RNA sequencing (snRNA-seq) to unveil the transcriptional changes of the skeletal muscle at single-nucleus resolution. 11 613 muscle nuclei (WT, 5838; Ptrf KO, 5775) were classified into 12 clusters corresponding to 11 nuclear types. Trajectory analysis revealed the potential transition between type IIb_1 and IIb_2 myonuclei upon muscular dystrophy. Functional enrichment analysis indicated that apoptotic signaling and enzyme-linked receptor protein signaling pathway were significantly enriched in type IIb_1 and IIb_2 myonuclei of Ptrf KO, respectively. The muscle structure development and the PI3K-AKT signaling pathway were significantly enriched in type IIa and IIx myonuclei of Ptrf KO. Meanwhile, metabolic pathway analysis showed a decrease in overall metabolic pathway activity of myonuclei subtypes upon muscular dystrophy, with the most decrease in type IIb_1 myonuclei. Gene regulatory network analysis found that the activity of Mef2c, Mef2d, Myf5, and Pax3 regulons was enhanced in type II myonuclei of Ptrf KO, especially in type IIb_2 myonuclei. In addition, we investigated the transcriptome changes in adipocytes and found that muscular dystrophy enhanced the lipid metabolic capacity of adipocytes. Our findings provide a valuable resource for exploring the molecular mechanism of muscular dystrophy due to Ptrf deficiency.
Subject(s)
Muscular Dystrophies , Transcriptome , Mice , Animals , Phosphatidylinositol 3-Kinases/metabolism , Muscular Dystrophies/genetics , Muscle, Skeletal/metabolism , Muscle Fibers, Skeletal/metabolismABSTRACT
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.
Subject(s)
Chromatin , Heterochromatin , Animals , Chickens/genetics , DNA Transposable Elements , Euchromatin/genetics , Heterochromatin/genetics , Mammals/genetics , Vertebrates/geneticsABSTRACT
BACKGROUND: The iconic giant panda (Ailuropoda melanoleuca), as both a flagship and umbrella species endemic to China, is a world famous symbol for wildlife conservation. The giant panda has several specific biological traits and holds a relatively small place in evolution. A high-quality genome of the giant panda is key to understanding the biology of this vulnerable species. FINDINGS: We generated a 2.48-Gb chromosome-level genome (GPv1) of the giant panda named "Jing Jing" with a contig N50 of 28.56 Mb and scaffold N50 of 134.17 Mb, respectively. The total length of chromosomes (n = 21) was 2.39-Gb, accounting for 96.4% of the whole genome. Compared with the previously published four genomes of the giant panda, our genome is characterized by the highest completeness and the correct sequence orientation. A gap-free and 850 kb length of immunoglobulin heavy-chain gene cluster was manually annotated in close proximity to the telomere of chromosome 14. Additionally, we developed an algorithm to predict the centromere position of each chromosome. We also constructed a complete chromatin structure for "Jing Jing", which includes inter-chromosome interaction pattern, A/B compartment, topologically associated domain (TAD), TAD-clique and promoter-enhancer interaction (PEI). CONCLUSIONS: We presented an improved chromosome-level genome and complete chromatin structure for the giant panda. This is a valuable resource for the future genetic and genomic studies on giant panda.
Subject(s)
Ursidae , Animals , Ursidae/genetics , Genomics , China , Chromosomes/genetics , ChromatinABSTRACT
The study of chromatin accessibility across tissues and developmental stages is essential for elucidating the transcriptional regulation of various phenotypes and biological processes. However, the chromatin accessibility profiles of multiple tissues in newborn pigs and across porcine liver development remain poorly investigated. Here, we used ATAC-seq and rRNA-depleted RNA-seq to profile open chromatin maps and transcriptional features of heart, kidney, liver, lung, skeletal muscle, and spleen in newborn pigs and porcine liver tissue in the suckling and adult stages, respectively. Specifically, by analyzing a union set of protein-coding genes (PCGs) and two types of transcripts (lncRNAs and TUCPs), we obtained a comprehensive annotation of consensus ATAC-seq peaks for each tissue and developmental stage. As expected, the PCGs with tissue-specific accessible promoters had active transcription and were relevant to tissue-specific functions. In addition, other non-coding tissue-specific peaks were involved in both physical activity and the morphogenesis of neonatal tissues. We also characterized stage-specific peaks and observed a close association between dynamic chromatin accessibility and hepatic function transition during liver postnatal development. Overall, this study expands our current understanding of epigenetic regulation in mammalian tissues and organ development, which can benefit both economic trait improvement and improve the biomedical usage of pigs.
Subject(s)
Chromatin , Epigenesis, Genetic , Animals , Swine/genetics , Chromatin/genetics , Gene Expression Regulation , Chromatin Immunoprecipitation Sequencing , Promoter Regions, Genetic , Mammals/geneticsABSTRACT
Long non-coding RNAs (lncRNAs) have been extensively studied, and their crucial roles in adipogenesis, lipid metabolism, and gene expression have been revealed. However, the exact regulatory or other mechanisms by which lncRNAs influence the functioning of mesenteric adipose tissue (MAT) remain largely unknown. In this paper, we report the identification of a new lncRNA, named G8110, from the MAT of Bama pigs. The coordinated expression levels of lncRNA G8110 and NFE2L1 were significantly decreased in the MAT of obese Bama pigs compared with those in the MAT of lean pigs. Using a bone mesenchymal stem cell adipogenic differentiation model, we found that lncRNA G8110 played a role in adipocyte differentiation by positively regulating NFE2L1. We also found that lncRNA G8110 inhibited the formation of intracellular lipid synthesis, promoted lipid metabolism, and inhibited the expression of inflammatory cytokines. Our findings regarding lipid synthesis may further promote the role of lncRNAs in driving adipose tissue remodeling and maintaining metabolic health.
Subject(s)
Adipogenesis , RNA, Long Noncoding , Animals , Swine , Adipogenesis/genetics , RNA, Long Noncoding/genetics , RNA, Long Noncoding/metabolism , Adipocytes/metabolism , Cell Differentiation/genetics , LipidsABSTRACT
BACKGROUND: Goose, descendants of migratory ancestors, have undergone extensive selective breeding, resulting in their remarkable ability to accumulate fat in the liver and exhibit a high tolerance for significant energy intake. As a result, goose offers an excellent model for studying obesity, metabolic disorders, and liver diseases in mammals. Although the impact of the three-dimensional arrangement of chromatin within the cell nucleus on gene expression and transcriptional regulation is widely acknowledged, the precise functions of chromatin architecture reorganization during fat deposition in goose liver tissues still need to be fully comprehended. RESULTS: In this study, geese exhibited more pronounced changes in the liver index and triglyceride (TG) content following the consumption of the high-fat diet (HFD) than mice without significant signs of inflammation. Additionally, we performed comprehensive analyses on 10 goose liver tissues (5 HFD, 5 normal), including generating high-resolution maps of chromatin architecture, conducting whole-genome gene expression profiling, and identifying H3K27ac peaks in the livers of geese and mice subjected to the HFD. Our results unveiled a multiscale restructuring of chromatin architecture, encompassing Compartment A/B, topologically associated domains, and interactions between promoters and enhancers. The dynamism of the three-dimensional genome architecture, prompted by the HFD, assumed a pivotal role in the transcriptional regulation of crucial genes. Furthermore, we identified genes that regulate chromatin conformation changes, contributing to the metabolic adaptation process of lipid deposition and hepatic fat changes in geese in response to excessive energy intake. Moreover, we conducted a cross-species analysis comparing geese and mice exposed to the HFD, revealing unique characteristics specific to the goose liver compared to a mouse. These chromatin conformation changes help elucidate the observed characteristics of fat deposition and hepatic fat regulation in geese under conditions of excessive energy intake. CONCLUSIONS: We examined the dynamic modifications in three-dimensional chromatin architecture and gene expression induced by an HFD in goose liver tissues. We conducted a cross-species analysis comparing that of mice. Our results contribute significant insights into the chromatin architecture of goose liver tissues, offering a novel perspective for investigating mammal liver diseases.
ABSTRACT
Cebpa is a master transcription factor gene for adipogenesis. However, the mechanisms of enhancer-promoter chromatin interactions controlling Cebpa transcriptional regulation during adipogenic differentiation remain largely unknown. To reveal how the three-dimensional structure of Cebpa changes during adipogenesis, we generated high-resolution chromatin interactions of Cebpa in 3T3-L1 preadipocytes and 3T3-L1 adipocytes using circularized chromosome conformation capture sequencing (4C-seq). We revealed dramatic changes in chromatin interactions and chromatin status at interaction sites during adipogenic differentiation. Based on this, we identified five active enhancers of Cebpa in 3T3-L1 adipocytes through epigenomic data and luciferase reporter assays. Next, epigenetic repression of Cebpa-L1-AD-En2 or -En3 by the dCas9-KRAB system significantly down-regulated Cebpa expression and inhibited adipocyte differentiation. Furthermore, experimental depletion of cohesin decreased the interaction intensity between Cebpa-L1-AD-En2 and the Cebpa promoter and down-regulated Cebpa expression, indicating that long-range chromatin loop formation was mediated by cohesin. Two transcription factors, RXRA and PPARG, synergistically regulate the activity of Cebpa-L1-AD-En2. To test whether Cebpa-L1-AD-En2 plays a role in adipose tissue development, we injected dCas9-KRAB-En2 lentivirus into the inguinal white adipose tissue (iWAT) of mice to suppress the activity of Cebpa-L1-AD-En2. Repression of Cebpa-L1-AD-En2 significantly decreased Cebpa expression and adipocyte size, altered iWAT transcriptome, and affected iWAT development. We identified functional enhancers regulating Cebpa expression and clarified the crucial roles of Cebpa-L1-AD-En2 and Cebpa promoter interaction in adipocyte differentiation and adipose tissue development.
Subject(s)
Adipogenesis , Chromatin , Animals , Mice , Adipocytes , Adipogenesis/genetics , Adipose Tissue , Cell DifferentiationABSTRACT
Adipose tissue plays critical roles in an individual's aging process. In this research, we use single-nucleus RNA sequencing to create highly detailed transcriptional maps of subcutaneous adipose tissue and visceral adipose tissue in young and aged mice. We comprehensively identify the various cell types within the white adipose tissue of mice, our study has elucidated seven distinct cell types within this tissue. Further analyses focus on adipocytes, fibro-adipogenic progenitors, and immune cells, revealing age-related declines in the synthetic metabolic activity of adipocytes, diminished immune regulation, and reduced maturation or proliferation of fibroblasts in undifferentiated adipocytes. We confirm the presence of distinct subpopulations of adipocytes, highlighting decreases in adipogenesis subgroups due to aging. Additionally, we uncover a reduction in immune cell subpopulations, driven by age-associated immune system dysregulation. Furthermore, pseudo-time analyses indicate that Adipocyte1 represents the 'nascent' phase of adipocyte development, while Adipocyte2 represents the 'mature' phase. We use cell-cell interaction to explore the age-dependent complexities of the interactions between FAPs and adipocytes, and observed increased expression of the inflammation-related Retn-Tlr4 interaction in older mice, while the anti-inflammatory Angpt1-Tek interaction was only detected in young mice. These transcriptional profiles serve as a valuable resource for understanding the functional genomics underlying metabolic disorders associated with aging in human adipose tissue.
Subject(s)
Adipocytes , Aging , Gene Expression Profiling , Animals , Aging/genetics , Mice , Adipocytes/metabolism , Transcriptome , Adipogenesis/genetics , Adipose Tissue/metabolism , Intra-Abdominal Fat/metabolism , Male , Mice, Inbred C57BL , Adipose Tissue, White/metabolism , Single-Cell AnalysisABSTRACT
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.
Subject(s)
Chickens , Haplotypes , Muscle, Skeletal , Animals , Chickens/genetics , Haplotypes/genetics , Muscle, Skeletal/metabolism , Chromosome Mapping/methods , Genome/geneticsABSTRACT
The storage of lipids as energy in adipose tissue (AT) has been conserved over the course of evolution. However, substantial differences in ATs physiological activities were reported among species. Hence, establishing the mechanisms shaping evolutionarily divergence in ATs transcriptomes could provide a deeper understanding of AT regulation and its roles in obesity-related diseases. While previous studies performed anatomical, physiological and morphological comparisons between ATs across different species, little is currently understood at the molecular phenotypic levels. Here, we characterized transcriptional and lipidomic profiles of available subcutaneous and visceral ATs samples across 15 vertebrate species, spanning more than 300 million years of evolution, including placental mammals, birds and reptiles. We provide detailed descriptions of the datasets produced in this study and report gene expression and lipid profiles across samples. We demonstrate these data are robust and reveal the AT transcriptome and lipidome vary greater among species than within the same species. These datasets may serve as a resource for future studies on the functional differences among ATs in vertebrate species.
Subject(s)
Lipidomics , Transcriptome , Vertebrates , Animals , Gene Expression Profiling , Mammals , Vertebrates/geneticsABSTRACT
Sympathetic innervation of white adipose tissue (WAT) plays a key role in the regulation of lipid metabolism. Sympathetic activation promotes release of norepinephrine (NE), which binds to adrenergic receptors on adipocytes, promoting adipocyte lipolysis and enhanced oxidative metabolism. However, the mechanism by which sympathetic nerves regulate lipid metabolism in pig adipose tissue remains unclear. We used NE to simulate the process of sympathetic driving in pig adipocytes. RNA sequencing (RNA-seq) was used to determine the gene expression profile of pig adipocytes responding to NE stimulation. Our data suggests that the lipolytic signaling pathway is activated in pig adipocytes upon acute stimulation of NE, resulting in enhanced lipid metabolism and lipolysis, consistent with the phenomena found in humans and mice. Specifically, differentially expressed protein coding genes (PCGs) (SIRT4, SLC27A1) are mainly associated with functions that inhibit fatty acid oxidation and promote lipid synthesis. Similarly, we investigated the changes in regulatory transcripts such as long non-coding RNAs (lncRNAs) and transcripts of uncertain coding potential (TUCP) in response to NE and found that differentially expressed lncRNAs (lncG47338, lncG30660, lncG29516, lncG3790) and TUCP (TUCP_G38001) were co-expressed with target genes related to the promotion of fatty acid ß-oxidation, lipolysis and oxidative metabolism, thus acting as regulators. These results indicate a broad suite of gene expression alterations in response to NE stimulation and promote the understanding of the molecular mechanisms by which NE regulates lipid metabolism in pigs.
Subject(s)
Lipolysis , RNA, Long Noncoding , Mice , Humans , Animals , Swine , Lipolysis/genetics , Norepinephrine/pharmacology , Norepinephrine/metabolism , RNA, Long Noncoding/metabolism , Adipocytes/metabolism , Gene Expression Profiling , Fatty Acids/metabolismABSTRACT
BACKGROUND: Studying the composition and developmental mechanisms in mammary gland is crucial for healthy growth of newborns. The mammary gland is inherently heterogeneous, and its physiological function dependents on the gene expression of multiple cell types. Most studies focused on epithelial cells, disregarding the role of neighboring adipocytes. RESULTS: Here, we constructed the largest transcriptomic dataset of porcine mammary gland cells thus far. The dataset captured 126,829 high-quality nuclei from physiological mammary glands across five developmental stages (d 90 of gestation, G90; d 0 after lactation, L0; d 20 after lactation, L20; 2 d post natural involution, PI2; 7 d post natural involution, PI7). Seven cell types were identified, including epithelial cells, adipocytes, endothelial cells, fibroblasts cells, immune cells, myoepithelial cells and precursor cells. Our data indicate that mammary glands at different developmental stages have distinct phenotypic and transcriptional signatures. During late gestation (G90), the differentiation and proliferation of adipocytes were inhibited. Meanwhile, partly epithelial cells were completely differentiated. Pseudo-time analysis showed that epithelial cells undergo three stages to achieve lactation, including cellular differentiation, hormone sensing, and metabolic activation. During lactation (L0 and L20), adipocytes area accounts for less than 0.5% of mammary glands. To maintain their own survival, the adipocyte exhibited a poorly differentiated state and a proliferative capacity. Epithelial cells initiate lactation upon hormonal stimulation. After fulfilling lactation mission, their undergo physiological death under high intensity lactation. Interestingly, the physiological dead cells seem to be actively cleared by immune cells via CCL21-ACKR4 pathway. This biological process may be an important mechanism for maintaining homeostasis of the mammary gland. During natural involution (PI2 and PI7), epithelial cell populations dedifferentiate into mesenchymal stem cells to maintain the lactation potential of mammary glands for the next lactation cycle. CONCLUSION: The molecular mechanisms of dedifferentiation, proliferation and redifferentiation of adipocytes and epithelial cells were revealed from late pregnancy to natural involution. This cell transcriptomic profile constitutes an essential reference for future studies in the development and remodeling of the mammary gland at different stages.
ABSTRACT
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.
Subject(s)
Chromatin , Weight Gain , Adult , Humans , Female , Swine , Animals , Obesity , Adipose Tissue , Chromatin Assembly and Disassembly , Weight Loss , MammalsABSTRACT
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.
ABSTRACT
BACKGROUND: Allele-specific expression (ASE) refers to the preferential expression of one allele over the other and contributes to adaptive phenotypic plasticity. Here, we used a reciprocal cross-model between phenotypically divergent European Berkshire and Asian Tibetan pigs to characterize 2 ASE classes: imprinting (i.e., the unequal expression between parental alleles) and sequence dependent (i.e., unequal expression between breed-specific alleles). We examined 3 transcript types, including protein-coding genes (PCGs), long noncoding RNAs, and transcripts of unknown coding potential, across 7 representative somatic tissues from hybrid pigs generated by reciprocal crosses. RESULTS: We identified a total of 92 putative imprinted transcripts, 69 (75.00%) of which are described here for the first time. By combining the transcriptome from purebred Berkshire and Tibetan pigs, we found â¼6.59% of PCGs are differentially expressed between breeds that are regulated by trans-elements (e.g., transcriptional factors), while only â¼1.35% are attributable to cis (e.g., promoters). The higher prevalence of trans-PCGs indicates the dominated effects of trans-regulation in driving expression differences and shaping adaptive phenotypic plasticity between breeds, which were supported by functional enrichment analysis. We also found strong evidence that expression changes mediated by cis-effects were associated with accumulated variants in promoters. CONCLUSIONS: Our study provides a comprehensive map of expression regulation that constitutes a valuable resource for the agricultural improvement of pig breeds.
Subject(s)
RNA, Long Noncoding , Transcriptome , Animals , Swine/genetics , Alleles , RNA, Long Noncoding/geneticsABSTRACT
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.
Subject(s)
Avian Proteins/genetics , Chickens/genetics , Chromatin/ultrastructure , Granulosa Cells/metabolism , Oogenesis/genetics , Transcriptome , Animals , Avian Proteins/classification , Avian Proteins/metabolism , Chickens/growth & development , Chickens/metabolism , Chromatin/chemistry , Enhancer Elements, Genetic , Female , Gene Expression Profiling , Gene Expression Regulation, Developmental , Gene Ontology , Granulosa Cells/cytology , Molecular Sequence Annotation , Oocytes/cytology , Oocytes/metabolism , Promoter Regions, GeneticABSTRACT
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.