Pig-eRNAdb: a comprehensive enhancer and eRNA dataset of pigs.

Enhancers and the enhancer RNAs (eRNAs) have been strongly implicated in regulations of transcriptions. Based the multi-omics data (ATAC-seq, ChIP-seq and RNA-seq) from public databases, Pig-eRNAdb is a dataset that comprehensively integrates enhancers and eRNAs for pigs using the machine learning strategy, which incorporates 82,399 enhancers and 37,803 eRNAs from 607 samples across 15 tissues of pigs. This user-friendly dataset covers a comprehensive depth of enhancers and eRNAs annotation for pigs. The coordinates of enhancers and the expression patterns of eRNAs are downloadable. Besides, thousands of regulators on eRNAs, the target genes of eRNAs, the tissue-specific eRNAs, and the housekeeping eRNAs are also accessible as well as the sequence similarity of eRNAs with humans. Moreover, the tissue-specific eRNA-trait associations encompass 652 traits are also provided. It will crucially facilitate investigations on enhancers and eRNAs with Pig-eRNAdb as a reference dataset in pigs.

A compendium of genetic regulatory effects across pig tissues.

The Farm Animal Genotype-Tissue Expression (FarmGTEx) project has been established to develop a public resource of genetic regulatory variants in livestock, which is essential for linking genetic polymorphisms to variation in phenotypes, helping fundamental biological discovery and exploitation in animal breeding and human biomedicine. Here we show results from the pilot phase of PigGTEx by processing 5,457 RNA-sequencing and 1,602 whole-genome sequencing samples passing quality control from pigs. We build a pig genotype imputation panel and associate millions of genetic variants with five types of transcriptomic phenotypes in 34 tissues. We evaluate tissue specificity of regulatory effects and elucidate molecular mechanisms of their action using multi-omics data. Leveraging this resource, we decipher regulatory mechanisms underlying 207 pig complex phenotypes and demonstrate the similarity of pigs to humans in gene expression and the genetic regulation behind complex phenotypes, supporting the importance of pigs as a human biomedical model.

PigBiobank: a valuable resource for understanding genetic and biological mechanisms of diverse complex traits in pigs.

To fully unlock the potential of pigs as both agricultural species for animal-based protein food and biomedical models for human biology and disease, a comprehensive understanding of molecular and cellular mechanisms underlying various complex phenotypes in pigs and how the findings can be translated to other species, especially humans, are urgently needed. Here, within the Farm animal Genotype-Tissue Expression (FarmGTEx) project, we build the PigBiobank (http://pigbiobank.farmgtex.org) to systematically investigate the relationships among genomic variants, regulatory elements, genes, molecular networks, tissues and complex traits in pigs. This first version of the PigBiobank curates 71 885 pigs with both genotypes and phenotypes from over 100 pig breeds worldwide, covering 264 distinct complex traits. The PigBiobank has the following functions: (i) imputed sequence-based genotype-phenotype associations via a standardized and uniform pipeline, (ii) molecular and cellular mechanisms underlying trait-associations via integrating multi-omics data, (iii) cross-species gene mapping of complex traits via transcriptome-wide association studies, and (iv) high-quality results display and visualization. The PigBiobank will be updated timely with the development of the FarmGTEx-PigGTEx project, serving as an open-access and easy-to-use resource for genetically and biologically dissecting complex traits in pigs and translating the findings to other species.

HDAC1-Mediated lncRNA Stimulatory Factor of Follicular Development to Inhibit the Apoptosis of Granulosa Cells and Regulate Sexual Maturity through miR-202-3p-COX1 Axis.

Abnormal sexual maturity exhibits significant detrimental effects on adult health outcomes, and previous studies have indicated that targeting histone acetylation might serve as a potential therapeutic approach to regulate sexual maturity. However, the mechanisms that account for it remain to be further elucidated. Using the mouse model, we showed that Trichostatin A (TSA), a histone deacetylase (HDAC) inhibitor, downregulated the protein level of Hdac1 in ovaries to promote the apoptosis of granulosa cells (GCs), and thus arrested follicular development and delayed sexual maturity. Using porcine GCs as a cell model, a novel sexual maturity-associated lncRNA, which was named as the stimulatory factor of follicular development (SFFD), transcribed from mitochondrion and mediated by HDAC1, was identified using RNA sequencing. Mechanistically, HDAC1 knockdown significantly reduced the H3K27ac level at the -953/-661 region of SFFD to epigenetically inhibit its transcription. SFFD knockdown released miR-202-3p to reduce the expression of cyclooxygenase 1 (COX1), an essential rate-limited enzyme involved in prostaglandin synthesis. This reduction inhibited the proliferation and secretion of 17ß-estradiol (E2) while promoting the apoptosis of GCs. Consequently, follicular development was arrested and sexual maturity was delayed. Taken together, HDAC1 knockdown-mediated SFFD downregulation promoted the apoptosis of GCs through the miR-202-3p-COX1 axis and lead to delayed sexual maturity. Our findings reveal a novel regulatory network modulated by HDAC1, and HDAC1-mediated SFFD may be a promising new therapeutic target to treat delayed sexual maturity.

Comprehensive Analysis of CircRNA Expression Profiles in Multiple Tissues of Pigs.

Circular RNAs (circRNAs) are a class of non-coding RNAs with diverse functions, and previous studies have reported that circRNAs are involved in the growth and development of pigs. However, studies about porcine circRNAs over the past few years have focused on a limited number of tissues. Based on 215 publicly available RNA sequencing (RNA-seq) samples, we conducted a comprehensive analysis of circRNAs in nine pig tissues, namely, the gallbladder, heart, liver, longissimus dorsi, lung, ovary, pituitary, skeletal muscle, and spleen. Here, we identified a total of 82,528 circRNAs and discovered 3818 novel circRNAs that were not reported in the CircAtlas database. Moreover, we obtained 492 housekeeping circRNAs and 3489 tissue-specific circRNAs. The housekeeping circRNAs were enriched in signaling pathways regulating basic biological tissue activities, such as chromatin remodeling, nuclear-transcribed mRNA catabolic process, and protein methylation. The tissue-specific circRNAs were enriched in signaling pathways related to tissue-specific functions, such as muscle system process in skeletal muscle, cilium organization in pituitary, and cortical cytoskeleton in ovary. Through weighted gene co-expression network analysis, we identified 14 modules comprising 1377 hub circRNAs. Additionally, we explored circRNA-miRNA-mRNA networks to elucidate the interaction relationships between tissue-specific circRNAs and tissue-specific genes. Furthermore, our conservation analysis revealed that 19.29% of circRNAs in pigs shared homologous positions with their counterparts in humans. In summary, this extensive profiling of housekeeping, tissue-specific, and co-expressed circRNAs provides valuable insights into understanding the molecular mechanisms of pig transcriptional expression, ultimately deepening our understanding of genetic and biological processes.

Single-cell multi-omics profiling reveals key regulatory mechanisms that poise germinal vesicle oocytes for maturation in pigs.

The molecular mechanisms controlling the transition from meiotic arrest to meiotic resumption in mammalian oocytes have not been fully elucidated. Single-cell omics technology provides a new opportunity to decipher the early molecular events of oocyte growth in mammals. Here we focused on analyzing oocytes that were collected from antral follicles in different diameters of porcine pubertal ovaries, and used single-cell M&T-seq technology to analyze the nuclear DNA methylome and cytoplasmic transcriptome in parallel for 62 oocytes. 10× Genomics single-cell transcriptomic analyses were also performed to explore the bi-directional cell-cell communications within antral follicles. A new pipeline, methyConcerto, was developed to specifically and comprehensively characterize the methylation profile and allele-specific methylation events for a single-cell methylome. We characterized the gene expressions and DNA methylations of individual oocyte in porcine antral follicle, and both active and inactive gene's bodies displayed high methylation levels, thereby enabled defining two distinct types of oocytes. Although the methylation levels of Type II were higher than that of Type I, Type II contained nearly two times more of cytoplasmic transcripts than Type I. Moreover, the imprinting methylation patterns of Type II were more dramatically divergent than Type I, and the gene expressions and DNA methylations of Type II were more similar with that of MII oocytes. The crosstalk between granulosa cells and Type II oocytes was active, and these observations revealed that Type II was more poised for maturation. We further confirmed Insulin Receptor Substrate-1 in insulin signaling pathway is a key regulator on maturation by in vitro maturation experiments. Our study provides new insights into the regulatory mechanisms between meiotic arrest and meiotic resumption in mammalian oocytes. We also provide a new analytical package for future single-cell methylomics study.

Local GHR roles in regulation of mitochondrial function through mitochondrial biogenesis during myoblast differentiation.

BACKGROUND: Myoblast differentiation requires metabolic reprogramming driven by increased mitochondrial biogenesis and oxidative phosphorylation. The canonical GH-GHR-IGFs axis in liver exhibits a great complexity in response to somatic growth. However, the underlying mechanism of whether local GHR acts as a control valve to regulate mitochondrial function through mitochondrial biogenesis during myoblast differentiation remains unknown. METHODS: We manipulated the GHR expression in chicken primary myoblast to investigate its roles in mitochondrial biogenesis and function during myoblast differentiation. RESULTS: We reported that GHR is induced during myoblast differentiation. Local GHR promoted mitochondrial biogenesis during myoblast differentiation, as determined by the fluorescence intensity of Mito-Tracker Green staining and MitoTimer reporter system, the expression of mitochondrial biogenesis markers (PGC1α, NRF1, TFAM) and mtDNA encoded gene (ND1, CYTB, COX1, ATP6), as well as mtDNA content. Consistently, local GHR enhanced mitochondrial function during myoblast differentiation, as determined by the oxygen consumption rate, mitochondrial membrane potential, ATP level and ROS production. We next revealed that the regulation of mitochondrial biogenesis and function by GHR depends on IGF1. In terms of the underlying mechanism, we demonstrated that IGF1 regulates mitochondrial biogenesis via PI3K/AKT/CREB pathway. Additionally, GHR knockdown repressed myoblast differentiation. CONCLUSIONS: In conclusion, our data corroborate that local GHR acts as a control valve to enhance mitochondrial function by promoting mitochondrial biogenesis via IGF1-PI3K/AKT/CREB pathway during myoblast differentiation. Video Abstract.

Advances in Oocyte Maturation In Vivo and In Vitro in Mammals.

The quality and maturation of an oocyte not only play decisive roles in fertilization and embryo success, but also have long-term impacts on the later growth and development of the fetus. Female fertility declines with age, reflecting a decline in oocyte quantity. However, the meiosis of oocytes involves a complex and orderly regulatory process whose mechanisms have not yet been fully elucidated. This review therefore mainly focuses on the regulation mechanism of oocyte maturation, including folliculogenesis, oogenesis, and the interactions between granulosa cells and oocytes, plus in vitro technology and nuclear/cytoplasm maturation in oocytes. Additionally, we have reviewed advances made in the single-cell mRNA sequencing technology related to oocyte maturation in order to improve our understanding of the mechanism of oocyte maturation and to provide a theoretical basis for subsequent research into oocyte maturation.

Study on detection of CNVs using human whole genome bisulfite sequencing data.

It has been reported that the aberrant DNA methylation may result in copy number variations (CNVs), and the CNVs may alter the levels of DNA methylation. Whole genome bisulfite sequencing (WGBS) is able to generate the sequencing data of DNAs, and shows the potential ability to detect CNVs. However, the evaluations and performances on the detections of CNVs using WGBS data is still unclear. In this study, five software with different strategies for CNV detections, e.g., BreakDancer, cn.mops, CNVnator, DELLY and Pindel, were selected to explore and benchmark the performances of CNV detections with WGBS data. Based on the real (2.62 billion reads) and simulated (12.35 billion reads) WGBS data of humans, we calculated the number, precision, recall, relative ability, memory usage, and running time of CNV detections by 150 times, and tried to figure out the optimal strategy for CNV detections with WGBS data. Based on the real WGBS data, Pindel detected the most deletions and duplications, CNVnator detected the deletions with the highest precision, cn.mops detected the duplications with the highest precision, Pindel detected the deletions with the highest recall, and cn.mops detected the duplications with the highest recall. Based on the simulated WGBS data, BreakDancer detected the most deletions, and cn.mops detected the most duplications. The CNVnator showed the highest precision and recall for both deletions and duplications. In real and simulated WGBS data, the ability of CNVnator to detect CNVs was likely to overtake that in the whole genome sequencing data. Additionally, DELLY and BreakDancer displayed the lowest peak of memory usage and the lest CPU runtime, while CNVnator expressed the highest peak of memory usage and the most CPU runtime. Taken together, CNVnator and cn.mops showed the excellent performances of CNV detections with WGBS data. These results suggested that it was feasible to detect CNVs using WGBS data, and provided the useful information to further investigate both CNVs and DNA methylation using WGBS data alone.

DNMT1-mediated lncRNA IFFD controls the follicular development via targeting GLI1 by sponging miR-370.

DNA methylation and long noncoding RNAs (lncRNAs) exhibit an indispensable role in follicular development. However, the specific mechanisms regarding lncRNAs mediated by DNA methylation in follicular development remain unclearly. In this study, we found that inhibiting the expression of DNMT1 promoted granulosa cells (GCs) apoptosis to inhibit follicular development. A novel follicular development-associated lncRNA named inhibitory factor of follicular development (IFFD) was mediated by DNMT1 and showed to arrest follicular development by inhibiting GCs proliferation and estrogen (E2) secretion but promoting GCs apoptosis. Mechanistically, the deactivated Cas9-TET1 demonstrated that the hypomethylation in -1261/-1254 region of IFFD promoted the transcription of IFFD by recruiting SP1. IFFD induced the expression of GLI family zinc finger 1 through competitive binding miR-370, thereby up-regulating the expression of CASP3 to promote GCs apoptosis, as well as downregulating the expressions of PCNA and CYP19A1 to inhibit GCs proliferation and E2 secretion. Collectively, DNMT1-mediated IFFD might be a novel target for the regulation of follicular development.

Evaluation of six machine learning classification algorithms in pig breed identification using SNPs array data.

Breed identification utilizing multiple information sources and methods is widely applicated in the field of animal genetics and breeding. Simultaneously, with the development of artificial intelligence, the integration of high-throughput genomic data and machine learning techniques is increasingly used for breed identification. In this context, we used 654 individuals from 15 pig breeds, evaluating the performance of machine learning and stacking ensemble learning classifiers, as well as the function of feature selection and anomaly detection in different scenarios. Our results showed that, when using a training set of 16 individuals per breed and 32 features (SNPs), the accuracy of breed identification with feature selection (eXtreme Gradient Boosting, XGBoost) could exceed 95.00% (nine breeds), and was improved by 7.04% over the results with random selection. For stacking ensemble learning, feature selection methods (including random selection method) were used before different base learners. When these base learners' training set had 16 individuals per breed and 32 features, the accuracy of stacking ensemble learning improved by 9.24% over the best base learner (nine breeds), but did not significantly increase the advantage over the models with XGBoost feature selection. When using a training set of 16 individuals and 512 features per breed, breed identification with anomaly detection (local outlier factor, LOF) and random selection could achieve an accuracy of 89.06% (15 breeds). These results show that machine learning could be an effective tool for breed identification and this study will also provide useful information for the application of machine learning in animal genetics and breeding.

Expression Profile of Housekeeping Genes and Tissue-Specific Genes in Multiple Tissues of Pigs.

Pigs have become an ideal model system for human disease research and development and an important farm animal that provides a valuable source of nutrition. To profile the all-sided gene expression and their biological functions across multiple tissues, we conducted a comprehensive analysis of gene expression on a large scale around the side of housekeeping genes (HKGs), tissue specific genes (TSGs), and the co-expressed genes in 14 various tissues. In this study, we identified 2351 HKGs and 3018 TSGs across tissues, among which 4 HKGs (COX1, UBB, OAZ1/NPFF) exhibited low variation and high expression levels, and 31 particular TSGs (e.g., PDC, FKBP6, STAT2, and COL1A1) were exclusively expressed in several tissues, including endocrine brain, ovaries, livers, backfat, jejunum, kidneys, lungs, and longissimus dorsi muscles. We also obtained 17 modules with 230 hub genes (HUBGs) by weighted gene co-expression network analysis. On the other hand, HKGs functions were enriched in the signaling pathways of the ribosome, spliceosome, thermogenesis, oxidative phosphorylation, and nucleocytoplasmic transport, which have been highly suggested to involve in the basic biological tissue activities. While TSGs were highly enriched in the signaling pathways that were involved in specific physiological processes, such as the ovarian steroidogenesis pathway in ovaries and the renin-angiotensin system pathway in kidneys. Collectively, these stable, specifical, and co-expressed genes provided useful information for the investigation of the molecular mechanism for an understanding of the genetic and biological processes of complex traits in pigs.

Benchmarking DNA methylation analysis of 14 alignment algorithms for whole genome bisulfite sequencing in mammals.

Whole genome bisulfite sequencing (WGBS) is an essential technique for methylome studies. Although a series of tools have been developed to overcome the mapping challenges caused by bisulfite treatment, the latest available tools have not been evaluated on the performance of reads mapping as well as on biological insights in multiple mammals. Herein, based on the real and simulated WGBS data of 14.77 billion reads, we undertook 936 mappings to benchmark and evaluate 14 wildly utilized alignment algorithms from reads mapping to biological interpretation in humans, cattle and pigs: Bwa-meth, BSBolt, BSMAP, Walt, Abismal, Batmeth2, Hisat_3n, Hisat_3n_repeat, Bismark-bwt2-e2e, Bismark-his2, BSSeeker2-bwt, BSSeeker2-soap2, BSSeeker2-bwt2-e2e and BSSeeker2-bwt2-local. Specifically, Bwa-meth, BSBolt, BSMAP, Bismark-bwt2-e2e and Walt exhibited higher uniquely mapped reads, mapped precision, recall and F1 score than other nine alignment algorithms, and the influences of distinct alignment algorithms on the methylomes varied considerably at the numbers and methylation levels of CpG sites, the calling of differentially methylated CpGs (DMCs) and regions (DMRs). Moreover, we reported that BSMAP showed the highest accuracy at the detection of CpG coordinates and methylation levels, the calling of DMCs, DMRs, DMR-related genes and signaling pathways. These results suggested that careful selection of algorithms to profile the genome-wide DNA methylation is required, and our works provided investigators with useful information on the choice of alignment algorithms to effectively improve the DNA methylation detection accuracy in mammals.

Ovary-derived circular RNAs profile analysis during the onset of puberty in gilts.

BACKGROUND: In mammals, the ovary is the essential system of female reproduction for the onset of puberty, and the abnormal puberty has negative outcomes on health. CircRNA is a non-coding RNA produced by non-canonical alternative splicing (AS). Several studies have reported that circRNA is involved in the gene regulation and plays an important role in some human diseases. However, the contribution of circRNA has received little known within the onset of puberty in ovary. RESULTS: Here, the profiles of ovarian circRNAs across pre-, in- and post-pubertal stages were established by RNA-sEq. In total, 972 circRNAs were identified, including 631 stage-specific circRNAs and 8 tissue-specific circRNAs. The biological functions of parental genes of circRNAs were enriched in steroid biosynthesis, autophagy-animal, MAPK signaling pathway, progesterone-mediated oocyte maturation and ras signaling pathway. Moreover, 5 circRNAs derived from 4 puberty-related genes (ESR1, JAK2, NF1 and ARNT) were found in this study. The A3SS events were the most alternative splicing, but IR events were likely to be arose in post-pubertal ovaries. Besides, the circRNA-miRNA-gene networks were explored for 10 differentially expressed circRNAs. Furthermore, the head-to-tail exon as well as the expressions of 10 circRNAs were validated by the divergent RT-qPCR and sanger sequencing. CONCLUSIONS: In summary, the profiles of ovarian circRNAs were provided during pubertal transition in gilts, and these results provided useful information for the investigation on the onset of puberty at the ovarian-circRNAs-level in mammals.

DNA methylation mediated RSPO2 to promote follicular development in mammals.

In female mammals, the proliferation, apoptosis, and estradiol-17ß (E2) secretion of granulosa cells (GCs) have come to decide the fate of follicles. DNA methylation and RSPO2 gene of Wnt signaling pathway have been reported to involve in the survival of GCs and follicular development. However, the molecular mechanisms for how DNA methylation regulates the expression of RSPO2 and participates in the follicular development are not clear. In this study, we found that the mRNA and protein levels of RSPO2 significantly increased during follicular development, but the DNA methylation level of RSPO2 promoter decreased gradually. Inhibition of DNA methylation or DNMT1 knockdown could decrease the methylation level of CpG island (CGI) in RSPO2 promoter and upregulate the expression level of RSPO2 in porcine GCs. The hypomethylation of -758/-749 and -563/-553 regions in RSPO2 promoter facilitated the occupancy of transcription factor E2F1 and promoted the transcriptional activity of RSPO2. Moreover, RSPO2 promoted the proliferation of GCs with increasing the expression level of PCNA, CDK1, and CCND1 and promoted the E2 secretion of GCs with increasing the expression level of CYP19A1 and HSD17B1 and inhibited the apoptosis of GCs with decreasing the expression level of Caspase3, cleaved Caspase3, cleaved Caspase8, cleaved Caspase9, cleaved PARP, and BAX. In addition, RSPO2 knockdown promoted the apoptosis of GCs, blocked the development of follicles, and delayed the onset of puberty with decreasing the expression level of Wnt signaling pathway-related genes (LGR4 and CTNNB1) in vivo. Taken together, the hypomethylation of -758/-749 and -563/-553 regions in RSPO2 promoter facilitated the occupancy of E2F1 and enhanced the transcription of RSPO2, which further promoted the proliferation and E2 secretion of GCs, inhibited the apoptosis of GCs, and ultimately ameliorated the development of follicles through Wnt signaling pathway. This study will provide useful information for further exploration on DNA-methylation-mediated RSPO2 pathway during follicular development.

Alternative Splicing Dynamics of the Hypothalamus-Pituitary-Ovary Axis During Pubertal Transition in Gilts.

The timing of puberty in mammals marks the point at which reproduction becomes possible. Abnormalities in the timing of puberty may exert a series of negative effects on subsequent health outcomes. Alternative splicing (AS) has not only emerged as a significant factor in the transcription of genes but it is also reported to play a role in the timing of puberty. However, to date, the changes and dynamics of AS during the onset of puberty is extremely seldom explored. In the present study, we used gilts as a research model to investigated the dynamics of AS and differentially expressed AS (DEAS) events within the hypothalamus-pituitary-ovary (HPO) axis across pre-, in-, and post-puberty. We detected 3,390, 6,098, and 9,085 DEAS events in the hypothalamus, pituitary, and ovary when compared across pre-, in-, and post-pubertal stages, respectively. Within the entire HPO axis, we also identified 22,889, 22,857, and 21,055 DEAS events in the pre-, in-, and post-pubertal stages, respectively. Further analysis revealed that the differentially spliced genes (DSGs) associated with staged DEAS events were likely to be enriched in the oxytocin signaling pathway, thyroid hormone signaling pathway, GnRH signaling pathway, and oocyte meiosis signaling pathway. The DSGs associated with DEAS events across the entire HPO axis were enriched in endocytosis signaling pathway, the MAPK signaling pathway, and the Rap1 signaling pathway. Moreover. the ASs of TAC1, TACR3, CYP19A1, ESR1, ESRRA, and FSHR were likely to regulate the functions of the certain HPO tissues during the onset of puberty. Collectively, the AS dynamics and DEAS events were comprehensively profiled in hypothalamus, pituitary, and ovary across the pre-, in-, and post-pubertal stages in pigs. These findings may enhance our knowledge of how puberty is regulated by AS and shed new light on the molecular mechanisms underlying the timing of puberty in mammals.

Identification of Circular RNAs in Hypothalamus of Gilts during the Onset of Puberty.

The disorders of puberty have shown negative outcomes on health of mammals, and the hypothalamus is thought to be the main regulator of puberty by releasing GnRH. Many studies show that the circular RNAs (circRNAs) might be implicated in the timing of puberty in mammals. However, the circRNAs in the hypothalamus of gilts have not been explored. To profile the changes and biological functions of circRNAs in the hypothalamus during the onset of puberty, RNA-seq was utilized to establish pre-, in-, and post-pubertal hypothalamic circRNAs profiles. In this study, the functions of hypothalamic circRNAs were enriched in the signaling pathway of neurotrophin, progesterone-mediated oocyte maturation, oocyte meiosis, insulin, ErbB, and mTOR, which have been highly suggested to be involved in the timing of puberty. Furthermore, 53 circRNAs were identified to be putative hypothalamus-specific expressed circRNAs, and some of them were exclusively expressed in the one of three pubertal stages. Moreover, 22 differentially expressed circRNAs were identified and chosen to construct the circRNA-miRNA-gene network. Moreover, 10 circRNAs were found to be driven by six puberty-related genes (ESR1, NF1, APP, ENPP2, ARNT, and DICER1). Subsequently, the expression changes of several circRNAs were confirmed by RT-qPCR. Collectively, the preliminary results of hypothalamic circRNAs provided useful information for the investigation of the molecular mechanism for the timing of puberty in gilts.

STAT4 targets KISS1 to promote the apoptosis of ovarian granulosa cells.

BACKGROUND: In mammals, it is known that the estradiol-17ß (E2) is mainly synthetized in ovarian granulosa cells (GCs), and the excessive apoptosis of GCs induces the follicular atresia. Many studies have implicated the essential role of KISS1, with the pro-synthetic effect of E2 and the anti-apoptotic effect on GCs, in the mammalian folliculogenesis, and several STAT4 potential binding sites were previously predicted on the promoter of KISS1 in pigs. However, the biological effects of STAT4 on GCs and the molecular regulation between STAT4 and KISS1 remained largely unknown. METHODS: Using the porcine GCs as the cellular model, the overexpression plasmid, small interfering RNA, 5'-deletion and luciferase assay were applied to investigate the molecular mechanisms for STAT4 regulating the expression of KISS1. RESULTS: In this study, the STAT4 negatively regulated the mRNA and protein levels of KISS1 in porcine GCs, and the mRNA level of STAT4 was observed to significantly decrease from immature to mature follicles, which was inversed with that of KISS1. The relative luciferase activity of KISS1 promoter was significantly increased with deletion of the fourth potential binding site (- 305/- 295), and ChIP further confirmed that the STAT4 bound at - 305/- 295 region of KISS1. Besides, the STAT4 significantly regulated the mRNA levels of PDK1, FOXO3 and TSC2 of PI3K signaling pathway to promote the cell apoptosis and the percentage of cells at G0/G1 phase of cell cycle in GCs. Alternatively, the STAT4 significantly decreased the mRNA levels of CYP17, 3B-HSD, 17B-33 HSD, ESR1, and ESR2, as well as the concentration of E2 in GCs. Furthermore, interfering with the expression of STAT4 was observed to significantly stimulate the pro-synthetic effect of E2 and anti-apoptotic effect of KISS1 in GCs. CONCLUSIONS: Collectively, the STAT4 might directly target at - 305/- 295 region of KISS1 to negatively regulate the transcription of KISS1, promote the cell apoptosis via PI3K signaling pathway, suppress the synthesis of E2 through the estrogen signaling pathway in porcine GCs. These proposed works could provide useful insight in further investigations on the molecular functionalities of STAT4 and KISS1 in the folliculogenesis of mammals.

Pituitary-Derived Circular RNAs Expression and Regulatory Network Prediction During the Onset of Puberty in Landrace × Yorkshire Crossbred Pigs.

Being the center of the hypothalamus-pituitary-ovary (HPO) axis, the pituitary plays a key role in the onset of puberty. Recent studies show that circular RNAs (circRNAs) can perform as miRNA sponges to regulate development in animals. However, the function of pituitary-derived circRNAs in first estrus remains unclear in pigs. In this study, we performed a genome-wide identification and characterization of circRNAs using pituitaries from Landrace × Yorkshire crossbred pigs at three stages: pre-, in-, and post-puberty, to describe such pituitary-derived circRNAs in pigs. A total of 5148 circRNAs were found in the gilts' pituitaries, averaging 18 682 bp in genomic distance, which consisted of approximately 91% exonic, 6% intergenic, and 3% intronic circRNAs. Furthermore, 158 novel circRNAs were identified for the first time and classified as putative pituitary-specific circRNAs. Their expression levels during the onset of puberty, significantly exceeded those of the other circRNAs, and the parental genes of these putative pituitary-specific circRNAs were enriched in "ssc04917: prolactin signaling pathway," "ssc04080: neuroactive ligand-receptor interaction," and "ssc04728: dopaminergic synapse" pathways, all of which were consistent with pituitary functioning. Additionally, 17 differentially regulated circRNAs were found and investigated for their potential interaction with miRNAs, along with genes, by constructing a circRNA-targeted miRNA-gene network. Taken together, these results provide new insight into the circRNA-mediated timing of puberty in gilts at the pituitary level.

Dynamic DNA methylation of ovaries during pubertal transition in gilts.

BACKGROUND: In female mammals, the initiation of puberty, coupling with the dramatically morphological changes in ovaries, indicates the sexual and follicular maturation. Previous studies have suggested that the disrupted DNA methylation results in the delayed puberty. However, to date, the changes in ovarian methylomes during pubertal transition have not been investigated. In this study, using gilts as a pubertal model, the genome-wide DNA methylation were profiled to explore their dynamics during pubertal transition across Pre-, In- and Post-puberty. RESULTS: During pubertal transition, the follicles underwent maturation and luteinization, coupled with the significant changes in the mRNA expression of DNMT1 and DNMT3a. DNA methylation levels of In-puberty were higher than that of Pre- and Post-puberty at the locations of genes and CpG islands (CGIs). Analysis of the DNA methylation changes identified 12,313, 20,960 and 17,694 differentially methylated CpGs (DMCs) for the comparisons of Pre- vs. In-, In vs. Post-, and Pre- vs. Post-puberty, respectively. Moreover, the CGIs, upstream and exonic regions showed a significant underrepresentation of DMCs, but the CGI shores, CGI shelves, intronic, downstream and intergenic regions showed a significant overrepresentation of DMCs. Furthermore, biological functions of these methylation changes enriched in PI3K-Akt signaling pathway, GnRH signaling pathway, and Insulin secretion, and the mRNA expressions of several genes of these signaling pathway, including MMP2, ESR1, GSK3B, FGF21, IGF1R, and TAC3, were significantly changed across Pre-, In- and Post-puberty in ovaries. CONCLUSIONS: During pubertal transition in gilts, the DNA methylation changes of ovaries were likely to affect the transcription of genes related to PI3K-Akt signaling pathway, GnRH signaling pathway, and Insulin secretion. These observations can provide new insight into the epigenetic mechanism of follicular and sexual maturation during pubertal transition in mammals.