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.

Identification of Selection Signatures and Candidate Genes Related to Environmental Adaptation and Economic Traits in Tibetan Pigs.

Tibetan pigs are indigenous to the Qinghai-Tibet Plateau and have been the subject of extensive genomic research primarily focused on their adaptation to high altitudes. However, genetic modifications associated with their response to low-altitude acclimation have not been thoroughly explored. To investigate the genetic basis underlying the low-altitude acclimation of Tibetan pigs, we generated and analyzed genotyping data of Tibetan pigs that inhabit high-altitude regions (average altitude 4000 m) and Tibetan pigs that have inhabited nearby low-altitude regions (average altitude 500 m) for approximately 20 generations. We found that the highland and lowland Tibetan pigs have distinguishable genotype and phenotype variations. We identified 46 and 126 potentially selected SNPs associated with 29 and 56 candidate genes in highland and lowland Tibetan pigs, respectively. Candidate genes in the highland Tibetan pigs were involved in immune response (NFYC and STAT1) and radiation (NABP1), whereas candidate genes in the lowland Tibetan pigs were related to reproduction (ESR2, DMRTA1, and ZNF366), growth and development (NTRK3, FGF18, and MAP1B), and blood pressure regulation (CARTPT). These findings will help to understand the mechanisms of environmental adaptation in Tibetan pigs and offer valuable information into the genetic improvement of Tibetan pigs pertaining to low-altitude acclimation and economic traits.

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.

A cell transcriptomic profile provides insights into adipocytes of porcine mammary gland across development.

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.

Expression Improvement of Recombinant Plasmids of the Interleukin-7 Gene in Chitosan-Derived Nanoparticles and Their Elevation of Mice Immunity.

To investigate a safe and effective approach for enhancing the in vivo expression of recombinant genes and improving the systemic immunity of animals against infectious diseases, we employed the interleukin-7 (IL-7) gene from Tibetan pigs to construct a recombinant eukaryotic plasmid (VRTPIL-7). We first examined VRTPIL-7's bioactivity on porcine lymphocytes in vitro and then encapsulated it with polyethylenimine (PEI), chitosan copolymer (CS), PEG-modified galactosylated chitosan (CS-PEG-GAL) and methoxy poly (ethylene glycol) (PEG) and PEI-modified CS (CS-PEG-PEI) nanoparticles using the ionotropic gelation technique. Next, we intramuscularly or intraperitoneally injected mice with various nanoparticles containing VRTPIL-7 to evaluate their immunoregulatory effects in vivo. We observed a significant increase in neutralizing antibodies and specific IgG levels in response to the rabies vaccine in the treated mice compared to the controls. Treated mice also exhibited increased leukocytes, CD8+ and CD4+ T lymphocytes, and elevated mRNA levels of toll-like receptors (TLR1/4/6/9), IL-1, IL-2, IL-4, IL-6, IL-7, IL-23, and transforming growth factor-beta (TGF-ß). Notably, the recombinant IL-7 gene encapsulated in CS-PEG-PEI induced the highest levels of immunoglobulins, CD4+ and CD8+ T cells, TLRs, and cytokines in the mice's blood, suggesting that chitosan-PEG-PEI may be a promising carrier for in vivo IL-7 gene expression and enhanced innate and adaptive immunity for the prevention of animal diseases.

Screening of candidate genes related to differences in growth and development between Chinese indigenous and Western pig breeds.

Neijiang (NJ) and Yacha (YC) are two indigenous pig breeds in the Sichuan basin of China, displaying higher resistance to diseases, lower lean ratio, and slower growth rate than the commercial Western pig breed Yorkshire (YS). The molecular mechanisms underlying the differences in growth and development between these pig breeds are still unknown. In the present study, five pigs from NJ, YC, and YS breeds were subjected to the whole genome resequencing, and then the differential single-nucleotide polymorphisms (SNPs) were screened using a 10-kb window sliding in 1-kb step using the Fst method. Finally, 48,924, 48,543, and 46,228 nonsynonymous single-nucleotide polymorphism loci (nsSNPs) were identified between NJ and YS, NJ and YC, and YC and YS, which highly or moderately affected 2,490, 800, and 444 genes, respectively. Moreover, three nsSNPs were detected in the genes of acetyl-CoA acetyltransferase 1 (ACAT1) insulin-like growth factor 2 receptor (IGF2R), insulin-like growth factor 2 and mRNA-binding protein 3 (IGF2BP3), which potentially affected the transformation of acetyl-CoA to acetoacetyl-CoA and the normal functions of the insulin signaling pathways. Moreover, serous determinations revealed significantly lower acetyl-CoA content in YC than in YS, supporting that ACAT1 might be a reason explaining the differences in growth and development between YC and YS breeds. Contents of phosphatidylcholine (PC) and phosphatidic acid (PA) significantly differed between the pig breeds, suggesting that glycerophospholipid metabolism might be another reason for the differences between Chinese and Western pig breeds. Overall, these results might contribute basic information to understand the genetic differences determining the phenotypical traits in pigs.

LncPLAAT3-AS Regulates PLAAT3-Mediated Adipocyte Differentiation and Lipogenesis in Pigs through miR-503-5p.

Animal fat deposition has a significant impact on meat flavor and texture. However, the molecular mechanisms of fat deposition are not well understood. LncPLAAT3-AS is a naturally occurring transcript that is abundant in porcine adipose tissue. Here, we focus on the regulatory role of lncPLAAT3-AS in promoting preadipocyte proliferation and adipocyte differentiation. By overexpressing or repressing lncPLAAT3 expression, we found that lncPLAAT3-AS promoted the transcription of its host gene PLAAT3, a regulator of adipocyte differentiation. In addition, we predicted the region of lncPLAAT3-AS that binds to miR-503-5p and showed by dual luciferase assay that lncPLAAT3-AS acts as a sponge to absorb miR-503-5p. Interestingly, miR-503-5p also targets and represses PLAAT3 expression and helps regulate porcine preadipocyte proliferation and differentiation. Taken together, these results show that lncPLAAT3-AS upregulates PLAAT3 expression by absorbing miR-503-5p, suggesting a potential regulatory mechanism based on competing endogenous RNAs. Finally, we explored lncPLAAT3-AS and PLAAT3 expression in adipose tissue and found that both molecules were expressed at significantly higher levels in fatty pig breeds compared to lean pig breeds. In summary, we identified the mechanism by which lncPLAAT3-AS regulates porcine preadipocyte proliferation and differentiation, contributing to our understanding of the molecular mechanisms of lipid deposition in pigs.

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.

Genome-wide association studies identified loci associated with both feed conversion ratio and residual feed intake in Yorkshire pigs.

Feed occupies a significant proportion in the production cost of pigs, and the feed efficiency (FE) in pigs is of utmost economic importance. Hence, the objective of this study is to identify single-nucleotide polymorphisms (SNPs) and candidate genes associated with FE-related traits, including feed conversion ratio (FCR) and residual feed intake (RFI). A genome-wide association study was conducted for FCR and RFI in 169 Yorkshire pigs using whole-genome sequencing data. A total of 23 and 33 suggestive significant SNPs (P < 1 × 10-6) were detected for FCR and RFI, respectively. However, none of SNPs achieved the genome-wide significance threshold (P < 5 × 10-8). Importantly, three common SNPs (SSC7:7987268, SSC13:42350250, and SSC13:42551718) were associated with both FCR and RFI. Additionally, the NEDD9 gene related to FCR and RFI traits was overlapped. This study detected novel SNPs on SSC7 and SSC13 common for FCR and RFI. These results provide new insights into the genetic mechanisms and candidate genes of FE-related traits in pigs.

Temporal microRNA expression profile of pig peripheral blood during postnatal development.

Gene expression profiles of blood can reflect the physiopathologic status of the immune system. The dynamic microRNA (miRNA) expression profiles of peripheral blood from pigs at different developmental stages, and how differential expression of miRNAs might relate to immune system development, are unknown. In this study, peripheral blood samples taken at five developmental stages were used to construct 15 miRNA libraries (three biological replicates/stage): 0 days (newborn), 30 days (weaning), 60 days (weaned), and 180 and 360 days (puberty). We identified 295 known mature miRNAs. Hierarchical clustering of the miRNA expression profile showed significant differences between individuals at the neonatal and postnatal stages. Functional enrichment analysis revealed that miRNAs differentially expressed between pairwise comparisons of the developmental stages were over-represented in immune-related pathways such as toll-like receptor signaling. The time-course of expression of the over-representated miRNAs exhibited a pattern of steady decline over time, for both the complete miRNA compendium and immune-related miRNAs. We identified six marker miRNAs that were highly negatively correlated with chronologic age and enriched for genes involved in immune-related pathways. This study of a peripheral blood miRNA transcriptome offers insight into immune system development in swine and provides a resource for pig genome annotation.

Allele-specific regulatory effects on the pig transcriptome.

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.

Differential DNA methylation analysis reveals key genes in Chinese Qingyu and Landrace pigs.

The Chinese Qingyu pig is a typical domestic fatty pig breed and an invaluable indigenous genetic resource in China. Compared with the Landrace pig, the Qingyu pig has unique meat characteristics, including muscle development, intramuscular fat, and other meat quality traits. At present, few studies have explored epigenetic differences due to DNA methylation between the Qingyu pig and the Landrace pig. In this study, 30 Qingyu pigs and 31 Landrace pigs were subjected to reduced representation bisulfite sequencing (RRBS). Genome-wide differential DNA methylation analysis was conducted. Six genomic regions, including regions on Sus scrofa chromosome (SSC) 1: 266.09-274.23 Mb, SSC5: 0.88-10.68 Mb, SSC8: 41.23-48.51 Mb, SSC12: 45.43-54.38 Mb, SSC13: 202.15-207.95 Mb, and SSC14: 126.43-139.85 Mb, were regarded as key regions that may be associated with phenotypic differences between the Qingyu pig and the Landrace pig. Furthermore, according to further analysis, five differentially methylated genes (ADCY1, FUBP3, GRIN2B, KIT, and PIK3R6) were identified as key candidate genes that might be associated with meat characteristics. Our findings provide new insights into the differences in DNA methylation between the Qingyu pig and the Landrace pig. These results enrich the epigenetic research of the Chinese Qingyu pig.

Genome-wide DNA methylation analysis in Chinese Chenghua and Yorkshire pigs.

BACKGROUND: The Chinese Chenghua pig (CHP) is a typical Chinese domestic fatty pig breed with superior meat quality characteristics, while the Yorkshire pig (YP) has the characteristics of fast growth and a high rate of lean meat. Long term natural selection and artificial selection resulted in great phenotypic differences between the two breeds, including growth, development, production performance, meat quality, and coat color. However, genome-wide DNA methylation differences between CHP and YP remain unclear. RESULTS: DNA methylation data were generated for muscle tissues of CHP and YP using reduced representation bisulfite sequencing (RRBS). In this study, a total of 2,416,211 CpG sites were identified. Besides, the genome-wide DNA methylation analysis revealed 722 differentially methylated regions (DMRs) and 466 differentially methylated genes (DMGs) in pairwise CHP vs. YP comparison. Six key genomic regions (Sus scrofa chromosome (SSC)1:253.47-274.23 Mb, SSC6:148.71-169.49 Mb, SSC7:0.25-9.86 Mb, SSC12:43.06-61.49 Mb, SSC14:126.43-140.95 Mb, and SSC18:49.17-54.54 Mb) containing multiple DMRs were identified, and differences of methylation patterns in these regions may be related to phenotypic differences between CHP and YP. Based on the functional analysis of DMGs, 8 DMGs (ADCY1, AGBL4, EXOC2, FUBP3, PAPPA2, PIK3R1, MGMT and MYH8) were considered as important candidate genes associated with muscle development and meat quality traits in pigs. CONCLUSIONS: This study explored the difference in meat quality between CHP and YP from the epigenetic point of view, which has important reference significance for the local pork industry and pork food processing.

A pig BodyMap transcriptome reveals diverse tissue physiologies and evolutionary dynamics of transcription.

A comprehensive transcriptomic survey of pigs can provide a mechanistic understanding of tissue specialization processes underlying economically valuable traits and accelerate their use as a biomedical model. Here we characterize four transcript types (lncRNAs, TUCPs, miRNAs, and circRNAs) and protein-coding genes in 31 adult pig tissues and two cell lines. We uncover the transcriptomic variability among 47 skeletal muscles, and six adipose depots linked to their different origins, metabolism, cell composition, physical activity, and mitochondrial pathways. We perform comparative analysis of the transcriptomes of seven tissues from pigs and nine other vertebrates to reveal that evolutionary divergence in transcription potentially contributes to lineage-specific biology. Long-range promoter-enhancer interaction analysis in subcutaneous adipose tissues across species suggests evolutionarily stable transcription patterns likely attributable to redundant enhancers buffering gene expression patterns against perturbations, thereby conferring robustness during speciation. This study can facilitate adoption of the pig as a biomedical model for human biology and disease and uncovers the molecular bases of valuable traits.

Temporal expression profiling of long noncoding RNA and mRNA in the peripheral blood during porcine development.

OBJECTIVE: We investigated the temporal expression profiles of long noncoding RNA (lncRNA) and mRNA in the peripheral blood of pigs during development and identified the lncRNAs that are related to the blood-based immune system. METHODS: Peripheral blood samples were obtained from the pigs at 0, 7, 28, and 180 days and 2 years of age. RNA sequencing was performed to survey the lncRNA and mRNA transcriptomes in the samples. Short time-series expression miner (STEM) was used to show temporal expression patterns in the mRNAs and lncRNAs. Gene ontology and Kyoto encyclopedia of genes and genomes analyses were performed to assess the genes' biological relevance. To predict the functions of the identified lncRNAs, we extracted mRNAs that were nearby loci and highly correlated with the lncRNAs. RESULTS: In total of 5,946 lncRNA and 12,354 mRNA transcripts were identified among the samples. STEM showed that most lncRNAs and mRNAs had similar temporal expression patterns during development, indicating the expressional correlation and functional relatedness between them. The five stages were divided into two classes: the suckling period and the late developmental stage. Most genes were expressed at low level during the suckling period, but at higher level during the late stages. Expression of several T-cell-related genes increased continuously during the suckling period, indicating that these genes are crucial for establishing the adaptive immune system in piglets at this stage. Notably, lncRNA TCONS-00086451may promote blood-based immune system development by upregulating nuclear factor of activated T-cells cytoplasmic 2 expression. CONCLUSION: This study provides a catalog of porcine peripheral blood-related lncRNAs and mRNAs and reveals the characteristics and temporal expression profiles of these lncRNAs and mRNAs during peripheral blood development from the newborn to adult stages in pigs.

Genomic analyses provide insights into breed-of-origin effects from purebreds on three-way crossbred pigs.

Crossbreeding is widely used aimed at improving crossbred performance for poultry and livestock. Alleles that are specific to different purebreds will yield a large number of heterozygous single-nucleotide polymorphisms (SNPs) in crossbred individuals, which are supposed to have the power to alter gene function or regulate gene expression. For pork production, a classic three-way crossbreeding system of Duroc × (Landrace × Yorkshire) (DLY) is generally used to produce terminal crossbred pigs with stable and prominent performance. Nonetheless, little is known about the breed-of-origin effects from purebreds on DLY pigs. In this study, we first estimated the distribution of heterozygous SNPs in three kinds of three-way crossbred pigs via whole genome sequencing data originated from three purebreds. The result suggested that DLY is a more effective strategy for three-way crossbreeding as it could yield more stably inherited heterozygous SNPs. We then sequenced a DLY pig family and identified 95, 79, 132 and 42 allele-specific expression (ASE) genes in adipose, heart, liver and skeletal muscle, respectively. Principal component analysis and unrestricted clustering analyses revealed the tissue-specific pattern of ASE genes, indicating the potential roles of ASE genes for development of DLY pigs. In summary, our findings provided a lot of candidate SNP markers and ASE genes for DLY three-way crossbreeding system, which may be valuable for pig breeding and production in the future.

Genome-Wide Chromatin Structure Changes During Adipogenesis and Myogenesis.

The recently developed high-throughput chromatin conformation capture (Hi-C) technology enables us to explore the spatial architecture of genomes, which is increasingly considered an important regulator of gene expression. To investigate the changes in three-dimensional (3D) chromatin structure and its mediated gene expression during adipogenesis and myogenesis, we comprehensively mapped 3D chromatin organization for four cell types (3T3-L1 pre-adipocytes, 3T3-L1-D adipocytes, C2C12 myoblasts, and C2C12-D myotubes). We demonstrate that the dynamic spatial genome architecture affected gene expression during cell differentiation. A considerable proportion (~22%) of the mouse genome underwent compartment A/B rearrangement during adipogenic and myogenic differentiation, and most (~80%) upregulated marker genes exhibited an active chromatin state with B to A switch or stable A compartment. More than half (65.4%-73.2%) of the topologically associating domains (TADs) are dynamic. The newly formed TAD and intensified local interactions in the Fabp gene cluster indicated more precise structural regulation of the expression of pro-differentiation genes during adipogenesis. About half (32.39%-59.04%) of the differential chromatin interactions (DCIs) during differentiation are promoter interactions, although these DCIs only account for a small proportion of genome-wide interactions (~9.67% in adipogenesis and ~4.24% in myogenesis). These differential promoter interactions were enriched with promoter-enhancer interactions (PEIs), which were mediated by typical adipogenic and myogenic transcription factors. Differential promoter interactions also included more differentially expressed genes than nonpromoter interactions. Our results provide a global view of dynamic chromatin interactions during adipogenesis and myogenesis and are a resource for studying long-range chromatin interactions mediating the expression of pro-differentiation genes.

Comparison of stomach microRNA transcriptomes of Tibetan and Yorkshire pigs by deep sequencing.

MiRNAs regulate the expression of target genes in diverse cellular processes and hence play important roles in different physiological processes, yet little is known about the stomach microRNAome (miRNAome) of the Tibetan pig. The objective of this experiment was to investigate differentially expressed stomach miRNAs participating in digestion. Firstly, we isolated total RNA by Trizol reagent from three Tibetan and three Yorkshire purebred pigs stomach samples at 90-day-old. Secondly, a comprehensive analysis of Tibetan and Yorkshire pig stomach miRNAomes was performed by small RNA sequencing in the Illumina HiSeq 2000 system. Finally, SYBR Green Real-time RT-PCR was performed to validate the differentially expressed miRNAs. We identified 318 unique miRNAs, 260 were co-expressed in both libraries, 17 and 31 miRNAs were specifically expressed in Tibetan and Yorkshire pigs respectively. Fifty six differentially expressed miRNAs were identified by the identifying differentially expressed genes 6 (IDEG6). Kyoto encyclopedia of genes and genomes analysis revealed that some of the differentially expressed miRNAs were associated with protein and fat digestion. Two differentially expressed miRNAs (miR-214-3p and ssc-un39) participating in the digestion of lipid were identified. Additionally, qRT-PCR results suggested that a higher expression of miR-214-3p in the Tibetan pig stomach could lead to relatively lower expression of calcium-dependent phospholipase A2, which is an enzyme important for the digestion of glycerol phospholipid. This study has delineated the different stomach miRNAs expression patterns of Tibetan and Yorkshire pigs, which would help explain the regulatory mechanisms of miRNAs in digestion of Tibetan pigs, and contribute to utilize a the unique digestion merits of Tibetan pig in future porcine hybridization breeding.

Testosterone-Dependent miR-26a-5p and let-7g-5p Act as Signaling Mediators to Regulate Sperm Apoptosis via Targeting PTEN and PMAIP1.

Recent evidence suggests that testosterone deficiency can dramatically decrease the quality of sperm. MicroRNAs (miRNAs) are conserved mediators of post-transcriptional gene regulation in eukaryotes. However, the systemic regulation and function of miRNAs in sperm quality decline induced by testosterone deficiency has not been investigated. Here, we found that the sperm apoptosis was significantly enhanced and the sperm motility was dramatically decreased in hemicastrated pigs. We then used small RNA sequencing to detect miRNA profiles of sperm from pigs with prepubertal hemicastration (HC) and compared them with control libraries. We identified 16 differentially expressed (DE) miRNAs between the sperm of prepubertal HC and control (CT) pigs. Functional enrichment analysis indicated that the target genes of these DE miRNAs were mainly enriched in apoptosis-related pathways including the p53, mitogen-activated protein kinase (MAPK), and mammalian target of rapamycin (mTOR) pathways. Furthermore, gain- and loss-of-function analyses demonstrated potential anti-apoptotic effects of the DE miRNAs miR-26a-5p and let-7g-5p on sperm cells. The luciferase reporter assay confirmed that PTEN and PMAIP1 are targets of miR-26a-5p and let-7g-5p, respectively. Spearman’s correlation analysis revealed significantly positive correlations between the sperm and its corresponding seminal plasma exosomes regarding the miRNA expression levels. In conclusion, testosterone deficiency-induced changes in the miRNA components of seminal plasma exosomes secreted by the genital tract may partially elucidate sperm miRNAome alterations, which are further responsible for the decline of sperm motility.

Genome-wide profiling of gene expression and DNA methylation provides insight into low-altitude acclimation in Tibetan pigs.

Efforts have been made to characterize the high-altitude adaption in Tibetan pigs and identified vast of genes or genomic regions undergone natural selection. Nonetheless, information concerning gene expression and DNA methylation changes response to low-altitude acclimation in Tibetan pigs is long overdue. To explore the exceptional mechanisms of gene expression and DNA methylation that are induced by low altitude environments in Tibetan pigs, we performed a comparative transcriptomic and DNA methylation analysis of skeletal muscle between indigenous Tibetan pigs that reside in high altitude regions (~4000m) and their counterparts that migrated to the geographically neighboring low-altitude regions (~500m) for nearly ten generations. Many genes that related to hypoxia response (EGLN3 and FLT1) and energy metabolism (TFB2M) were differentially expressed, but without significant DNA methylation changes. We also found genes embedded in differentially methylated regions were mainly involved in 'Starch and sucrose metabolism', 'glucuronosyltransferase activity' processes. Specifically, our results showed increased SIN3A mRNA expression, with hypomethylation status of its promoter, in longissimus dorsi muscle of low-altitude Tibetan pig. Another gene, CACNG6, showed decreasing expression level with an elevated methylation in its intron 3. These results indicated DNA-methylation-mediated expression alterations in low-altitude acclimation. We envision that this study will serve as a valuable resource for mammal acclimation research and agricultural food industry.