Comparative Analysis of the Ovary Transcriptome among Wanyue Black and Yorkshire Gilts Using RNA-Seq.

Pubertal genetic variations between the indigenous Chinese Wanyue Black pig breed and the imported Yorkshire breed significantly impact their reproductive capacity. In order to identify the differentially expressed genes, gene networks, and metabolic pathways in ovary transcriptome of gilts, the serum hormone levels were analyzed by ELISA, and RNA-seq was performed to analyze ovarian genes. Our results reveal higher estradiol (E2) levels in Wanyue black gilts compared to Yorkshire gilts, while Yorkshire gilts exhibit elevated progesterone (P4) and GnRH levels. We identified a total of 154 differentially expressed genes (DEGs), with 87 up-regulated and 67 down-regulated genes in the Wanyue black gilts ovaries compared to the Yorkshire gilts. GO enrichment analysis unveiled the participation of DEGs in processes such as "Reproduction", "Reproductive system development", and "Ovarian follicle development". Moreover, KEGG enrichment analysis revealed the involvement of DEGs in multiple signaling pathways associated with hormone biosynthesis and puberty, encompassing "Steroid hormone biosynthesis", "Estrogen signaling pathway", and "Prolactin signaling pathway". The subsequent bioinformatics analysis identified nine functional genes that potentially contribute to the disparity in ovaries between Wanyue black gilts and Yorkshire gilts. This study offers significant insights into the endocrine and genetic aspects of pubertal development in gilts.

Fecal Microbial Structure and Metabolic Profile in Post-Weaning Diarrheic Piglets.

(1) Background: Piglet diarrhea is one of the most serious diseases in pigs and has brought great economic losses to the pig industry. Alteration of the gut microbiota is an important factor in the etiology of piglet diarrhea. Therefore, this study aimed to analyze the differences in the gut microbial structures and fecal metabolic profile between post-weaning diarrhea and healthy Chinese Wannan Black pigs. (2) Methods: An integrated approach of 16S rRNA gene sequencing combined with LC/MS-based metabolomics was employed in this study. (3) Results: We found an increase in the relative abundance of the bacterial genus Campylobacter and a decrease in phylum Bacteroidetes and the species Streptococcus gallolyticus subsp. macedonicus. (S. macedonicus) in piglet diarrhea. Meanwhile, obvious changes in the fecal metabolic profile of diarrheic piglets were also detected, particularly higher levels of polyamines (spermine and spermidine). Moreover, there were substantial associations between the disturbed gut microbiota and the altered fecal metabolites, especially a strong positive relationship between spermidine and Campylobacter. (4) Conclusions: These observations may provide novel insights into potential etiologies related to post-weaning diarrhea and further enhance our understanding of the role of gut microbiota in host homeostasis and in modulating gut microbial structure.

Identification of microRNAs implicated in modulating resveratrol-induced apoptosis in porcine granulosa cells.

MicroRNAs (miRNAs) are small, noncoding RNAs that play a crucial role in the complex and dynamic network that regulates the apoptosis of porcine ovarian granulosa cells (POGCs). Resveratrol (RSV) is a nonflavonoid polyphenol compound that is involved in follicular development and ovulation. In previous study, we established a model of RSV treatment of POGCs, confirming the regulatory effect of RSV in POGCs. To investigate the miRNA-level effects of RSV on POGCs to reveal differentially expressed miRNAs, a control group (n = 3, 0 µM RSV group), a low RSV group (n = 3, 50 µM RSV group), and a high RSV group (n = 3, 100 µM RSV group) were created for small RNA-seq. In total, 113 differentially expressed miRNAs (DE-miRNAs) were identified, and a RT-qPCR analysis showed a correlation with the sequencing data. Functional annotation analysis revealed that DE-miRNAs in the LOW vs. CON group may be involved in cell development, proliferation, and apoptosis. In the HIGH vs. CON group, RSV functions were associated with metabolic processes and responses to stimuli, while the pathways were related to PI3K24, Akt, Wnt, and apoptosis. In addition, we constructed miRNA-mRNA networks related to Apoptosis and Metabolism. Then, ssc-miR-34a and ssc-miR-143-5p were selected as key miRNAs. In conclusion, this study provided an improved understanding of effects of RSV on POGCs apoptosis through the miRNA modulations. The results suggest that RSV may promote POGCs apoptosis by stimulating the miRNA expressions and provided a better understanding of the role of miRNAs combined with RSV in ovarian granulosa cell development in pigs.

Metabolomic and Proteomic Profiling of Porcine Intestinal Epithelial Cells Infected with Porcine Epidemic Diarrhea Virus.

Porcine epidemic diarrhea virus (PEDV) infection results in severe epidemic diarrhea and the death of suckling pigs. Although new knowledge about the pathogenesis of PEDV has been improved, alterations in metabolic processes and the functional regulators involved in PEDV infection with host cells remain largely unknow. To identify cellular metabolites and proteins related to PEDV pathogenesis, we synergistically investigated the metabolome and proteome profiles of PEDV-infected porcine intestinal epithelial cells by liquid chromatography tandem mass spectrometry and isobaric tags for relative and absolute quantification techniques. We identified 522 differential metabolites in positive and negative ion modes and 295 differentially expressed proteins after PEDV infection. Pathways of cysteine and methionine metabolism, glycine, serine and threonine metabolism, and mineral absorption were significantly enriched by differential metabolites and differentially expressed proteins. The betaine-homocysteine S-methyltransferase (BHMT) was indicated as a potential regulator involved in these metabolic processes. We then knocked down the BHMT gene and observed that down-expression of BHMT obviously decreased copy numbers of PEDV and virus titers (p < 0.01). Our findings provide new insights into the metabolic and proteomic profiles in PEDV-infected host cells and contribute to our further understanding of PEDV pathogenesis.

2,5-Hexanedione Affects Ovarian Granulosa Cells in Swine by Regulating the CDKN1A Gene: A Transcriptome Analysis.

N-hexane, a common industrial organic solvent, causes multiple organ damage owing to its metabolite, 2,5-hexanedione (2,5-HD). To identify and evaluate the effects of 2,5-HD on sows' reproductive performance, we used porcine ovarian granulosa cells (pGCs) as a vehicle and carried out cell morphology and transcriptome analyses. 2,5-HD has the potential to inhibit the proliferation of pGCs and induce morphological changes and apoptosis depending on the dose. RNA-seq analyses identified 4817 differentially expressed genes (DEGs), with 2394 down-regulated and 2423 up-regulated following 2,5-HD exposure treatment. The DEG, cyclin-dependent kinase inhibitor 1A (CDKN1A), according to the Kyoto Encyclopedia of Genes and Genomes enrichment analysis, was significantly enriched in the p53 signaling pathway. Thus, we evaluated its function in pGC apoptosis in vitro. Then, we knocked down the CDKN1A gene in the pGCs to identify its effects on pGCs. Its knockdown decreased pGC apoptosis, with significantly fewer cells in the G1 phase (p < 0.05) and very significantly more cells in the S phase (p < 0.01). Herein, we revealed novel candidate genes that influence pGCs apoptosis and cell cycle and provided new insights into the role of CDKN1A in pGCs during apoptosis and cell cycle arrest.

Copy Number Variation Analysis Revealed the Evolutionary Difference between Chinese Indigenous Pigs and Asian Wild Boars.

Copy number variation (CNV) has been widely used to study the evolution of different species. We first discovered different CNVs in 24 Anqingliubai pigs and 6 Asian wild boars using next-generation sequencing at the whole-genome level with 10× depth to understand the relationship between genetic evolution and production traits in wild boars and domestic pigs. A total of 97,489 CNVs were identified and divided into 10,429 copy number variation regions (CNVRs), occupying 32.06% of the porcine genome. Chromosome 1 had the most CNVRs, and chromosome 18 had the least. Ninety-six CNVRs were selected using VST 1% based on the signatures of all CNVRs, and sixty-five genes were identified in the selected regions. These genes were strongly correlated with traits distinguishing groups by enrichment in Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways, such as growth (CD36), reproduction (CIT, RLN), detoxification (CYP3A29), and fatty acid metabolism (ELOVL6). The QTL overlapping regions were associated with meat traits, growth, and immunity, which was consistent with CNV analysis. Our findings increase the understanding of evolved genome structural variations between wild boars and domestic pigs, and provide new molecular biomarkers to guide breeding and the efficient use of available genetic resources.

Detection of Selection Signatures in Anqing Six-End-White Pigs Based on Resequencing Data.

As a distinguished Chinese indigenous pig breed that exhibits disease resistance and high meat quality, the Anqing six-end-white (AQ) pig represents a valuable germplasm resource for improving the quality of the pig breeding industry. In this study, 24 AQ pigs that were distantly blood-related and 6 Asian Wild Boar (AWB) were selected for 10× deep-genome resequencing. The signatures of the selection were analyzed to explore the genetic basis of their germplasm characteristics and to identify excellent germplasm-related functional genes based on NGS data. A total of 49,289,052 SNPs and 6,186,123 indels were detected across the genome in 30 pigs. Most of the genetic variations were synonym mutations and existed in the intergenic region. We identified 275 selected regions (top 1%) harboring 85 genes by applying a crossover approach based on genetic differentiation (FST) and polymorphism levels (π ratio). Some genes were found to be positively selected in AQ pigs' breeding. The SMPD4 and DDX18 genes were involved in the immune response to pseudorabies virus (PRV) and porcine reproductive and respiratory syndrome virus (PRRSV). The BCL6 and P2RX6 genes were involved in biological regulation of immune T cells and phagocytes. The SLC7A4 and SPACA4 genes were related to reproductive performance. The MSTN and HIF1A genes were related to fat deposition and muscle development. Moreover, 138 overlapping regions were detected in selected regions and ROH islands of AQ pigs. Additionally, we found that the QTLs with the most overlapping regions were related to back fat thickness, meat color, pH value, fatty acid content, immune cells, parasitic immunity, and bacterial immunity. Based on functional enrichment analysis and QTLs mapping, we conducted further research on the molecular genetic basis of germplasm traits (disease resistance and excellent meat quality). These results are a reliable resource for conserving germplasm resources and exploiting molecular markers of AQ pigs.

Identification of Candidate Genes and Regulatory Competitive Endogenous RNA (ceRNA) Networks Underlying Intramuscular Fat Content in Yorkshire Pigs with Extreme Fat Deposition Phenotypes.

Intramuscular fat (IMF) content is vital for pork quality, serving an important role in economic performance in pig industry. Non-coding RNAs, with mRNAs, are involved in IMF deposition; however, their functions and regulatory mechanisms in porcine IMF remain elusive. This study assessed the whole transcriptome expression profiles of the Longissimus dorsi muscle of pigs with high (H) and low (L) IMF content to identify genes implicated in porcine IMF adipogenesis and their regulatory functions. Hundreds of differentially expressed RNAs were found to be involved in fatty acid metabolic processes, lipid metabolism, and fat cell differentiation. Furthermore, combing co-differential expression analyses, we constructed competing endogenous RNAs (ceRNA) regulatory networks, showing crosstalk among 30 lncRNAs and 61 mRNAs through 20 miRNAs, five circRNAs and 11 mRNAs through four miRNAs, and potential IMF deposition-related ceRNA subnetworks. Functional lncRNAs and circRNAs (such as MSTRG.12440.1, ENSSSCT00000066779, novel_circ_011355, novel_circ_011355) were found to act as ceRNAs of important lipid metabolism-related mRNAs (LEP, IP6K1, FFAR4, CEBPA, etc.) by sponging functional miRNAs (such as ssc-miR-196a, ssc-miR-200b, ssc-miR10391, miR486-y). These findings provide potential regulators and molecular regulatory networks that can be utilized for research on IMF traits in pigs, which would aid in marker-assisted selection to improve pork quality.

Genome-wide scan for runs of homozygosity in Asian wild boars and Anqing six-end-white pigs.

Inbreeding and loss of genetic diversity are serious problems in local Chinese pig breeds. Runs of homozygosity (ROHs) are contiguous lengths of homozygous genotypes that can provide information on inbreeding levels, mating schemes, selection pressure, and genetic events. The Anqing six-end-white (AQ) pig, an important autochthonous pig breed bred in the Anhui province of China, plays a key role in the local pig industry. In recent decades, AQ pigs have become a closed breed with a small population size in conservation farms, raising the issue of inbreeding decline. In this study, we used 10× resequencing to detect the whole genome of 24 AQ pigs and six Asian wild boars (AWBs). We used the Plink software to analyze the homozygosity of the genome and the distribution of ROHs in the genome based on the resequencing data. We obtained 935.04 Gb of raw data from the resequencing results and more than 45 822 239 SNPs. Additionally, we identified 28 702 ROHs. Compared with AWB, AQ pigs had higher ROH numbers, longer ROH rates, and higher FROH values, which revealed that artificial selection influenced ROH distribution and genome inbreeding level. A total of 307 and 205 ROH islands were identified in the AQ pigs and AWBs respectively. The genes of ROH islands in AQ pigs were mainly enriched in immune biological processes. Our findings provide a useful reference for developing breeding programs to maintain genetic diversity and germplasm resources in AQ pigs.

Integrated analysis of lncRNA and mRNA for the apoptosis of porcine ovarian granulosa cells after polyphenol resveratrol treatment.

Resveratrol (RES) is a non-flavonoid polyphenol compound that can be involved in follicular development and ovulation. However, the mechanism by which resveratrol regulates the apoptosis of porcine ovarian granulosa cells (POGCs) through long non-coding RNA (lncRNA) is poorly understood. We generated POGCs models of different doses of RES (0, 25, 50, 75, and 100 µM). It was observed that the cell viability was the highest in the 50 µM group, and the highest apoptosis rates were recorded in the 100 µM group. Therefore, a control group (n = 3, 0 µM RES group), a low RES group (n = 3, 50 µM RES group), and a high RES group (n = 3, 100 µM RES group) of POGCs were created for next RNA sequencing. Gene Ontology (GO) indicated that differentially expressed lncRNAs associated with apoptotic process were highly enriched. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of lncRNA target genes found that the Wnt signaling pathway and PI3K-Akt signaling pathway were both enriched. Furthermore, we constructed lncRNA-mRNA networks related to Metabolic and Cell Apoptosis, respectively. In the networks, five key-lncRNAs were screened, which may play a significant role in the process of POGCs metabolism and apoptosis. Furthermore, we focused on the function of a lnc-GAM (lncRNA associated with Granulosa cells Apoptosis and Metabolism) and verified that lnc-GAM could influence cell apoptosis in POGCs development by affecting the mRNA expression of apoptosis-related markers, and also affects the secretion of steroid hormones and related genes expression in POGCs cultured in vitro. Our study provides seminal data and important new insights into the regulation of reproductive mechanisms in porcine and other female mammals.

Genome-wide scan for runs of homozygosity identifies candidate genes in Wannan Black pigs.

OBJECTIVE: Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes that can reveal inbreeding levels, selection pressure, and mating schemes. In this study, ROHs were evaluated in Wannan Black pigs to assess the inbreeding levels and the genome regions with high ROH frequency. METHODS: In a previous study, we obtained 501.52 GB of raw data from resequencing (10×) of the genome and identified 21,316,754 single-nucleotide variants in 20 Wannan Black pig samples. We investigated the number, length, and frequency of ROH using resequencing data to characterize the homozygosity in Wannan Black pigs and identified genomic regions with high ROH frequencies. RESULTS: In this work, 1,813 ROHs (837 ROHs in 100 to 500 kb, 449 ROHs in 500 to 1,000 kb, 527 ROHs in >1,000 kb) were identified in all samples, and the average genomic inbreeding coefficient (FROH) in Wannan Black pigs was 0.5234. Sixty-one regions on chromosomes 2, 3, 7, 8, 13, 15, and 16 harbored ROH islands. In total, 105 genes were identified in 42 ROH islands, among which some genes were related to production traits. CONCLUSION: This is the first study to identify ROH across the genome of Wannan Black pigs, the Chinese native breed of the Anhui province. Overall, Wannan Black pigs have high levels of inbreeding due to the influence of ancient and recent inbreeding due to the genome. These findings are a reliable resource for future studies and contribute to save and use the germplasm resources of Wannan Black pigs.

MiR-21-5p actions at the Smad7 gene during pig ovarian granulosa cell apoptosis.

MicroRNAs (miRNAs) are endogenous non-coding RNAs in eukaryotic cells that modulate apoptosis of ovarian granulosa cells (GCs), which is an important cause of mammalian follicular atresia. In the present study, associations were evaluated between miR-21-5p and the extent of Smad7 protein production in regulation of ovarian granulosa cell (pGC) apoptosis. There was detection of miR-21-5p and Smad7 primarily in the cytoplasm and nucleus of pGCs, respectively. When there was an enhanced abundance of miR-21-5p and decreased abundance of Smad7 there were similar effects in pGCs, including inducing proliferation, inhibiting apoptosis, increasing the number of cells in S and G2/M phases, increasing serum estradiol, and decreasing serum progesterone concentrations. Furthermore, the Smad7 mRNA transcript was identified as a target for miR-21-5p actions, with enhanced abundances of miR-21-5p being associated with a lesser abundance of Smad7 mRNA transcript and protein in pGCs. Overall, results from the present study indicate that miR-21-5p has actions on the Smad7 mRNA transcript during the process of ovarian granulosa cell apoptosis in pigs.

Identification of Signatures of Selection by Whole-Genome Resequencing of a Chinese Native Pig.

Identification of genomic signatures of selection that help reveal genetic mechanisms underlying traits in domesticated pigs is of importance. Anqing six-end-white pig (ASP), a representative of the native breeds in China, has many distinguishing phenotypic characteristics. To identify the genomic signatures of selection of the ASP, whole-genome sequencing of 20 ASPs produced 469.01 Gb of sequence data and more than 26 million single-nucleotide polymorphisms. Combining these data with the available whole genomes of 13 Chinese wild boars, 157 selected regions harboring 48 protein-coding genes were identified by applying the polymorphism levels (θπ) and genetic differentiation (F ST ) based cross approaches. The genes found to be positively selected in ASP are involved in crucial biological processes such as coat color ( MC1R ), salivary secretion ( STATH ), reproduction ( SPIRE2 , OSBP2 , LIMK2 , FANCA , and CABS1 ), olfactory transduction ( OR5K4 ), and growth ( NPY1R , NPY5R , and SELENOM ). Our research increased the knowledge of ASP phenotype-related genes and help to improve our understanding of the underlying biological mechanisms and provide valuable genetic resources that enable effective use of pigs in agricultural production.

Transcriptomic comparison of liver tissue between Anqing six-end-white pigs and Yorkshire pigs based on RNA sequencing.

Chinese indigenous pig and Western commercial pig breeds show different patterns of lipid metabolism, fat deposition, and fatty acid composition; for these reasons, they have become vitally important models of energy metabolism and obesity in humans. To compare the mechanisms underlying lipid metabolism between Yorkshire pigs (lean type) and Anqing six-end-white pigs (obese type), the liver transcriptomes of six castrated boars with a body weight of approximately 100 kg (three Yorkshire and three Anqing) were analyzed by RNA-seq. The total number of reads produced for each liver sample ranged from 47.05 to 62.6 million. Among 362 differentially expressed genes, 142 were up-regulated and 220 were down-regulated in Anqing six-end-white pigs. Based on these data, 79 GO terms were significantly enriched. The top 10 (the 10 with lowest corrected P-value) significantly enriched GO terms were identified, including lipid metabolic process and carboxylic acid metabolic process. Pathway analysis revealed three significantly enriched KEGG pathways including PPAR signaling pathway, steroid hormone biosynthesis, and retinol metabolism. Based on protein-protein interaction networks, multiple genes responsible for lipid metabolism were identified, such as PCK1, PPARA, and CYP7A1, and these were considered promising candidate genes that could affect porcine liver lipid metabolism and fat deposition. Our results provide abundant transcriptomic information that will be useful for animal breeding and biomedical research.

Genomic analysis reveals selection signatures of the Wannan Black pig during domestication and breeding.

OBJECTIVE: The Wannan Black pig is a typical Chinese indigenous, disease-resistant pig breed with high fertility, and a crude-feed tolerance that has been bred by artificial selection in the south of Anhui province for a long time. However, genome variation, genetic relationships with other pig breeds, and domestication, remain poorly understood. Here, we focus on elucidating the genetic characteristics of the Wannan Black pig and identifying selection signatures during domestication and breeding. METHODS: We identified the whole-genome variation in the Wannan Black pig and performed population admixture analyses to determine genetic relationships with other domesticated pig breeds and wild boars. Then, we identified the selection signatures between the Wannan Black pig and Asian wild boars in 100-kb windows sliding in 10 kb steps by using two approaches: the fixation index (FST) and π ratios. RESULTS: Resequencing the Wannan Black pig genome yielded 501.52 G of raw data. After calling single-nucleotide variants (SNVs) and insertions/deletions (InDels), we identified 21,316,754 SNVs and 5,067,206 InDels (2,898,582 inserts and 2,168,624 deletions). Additionally, we found genes associated with growth, immunity, and digestive functions. CONCLUSION: Our findings help in explaining the unique genetic and phenotypic characteristics of Wannan Black pigs, which in turn can be informative for future breeding programs of Wannan Black pigs.

Relationship among porcine lncRNA TCONS_00010987, miR-323, and leptin receptor based on dual luciferase reporter gene assays and expression patterns.

OBJECTIVE: Considering the physiological and clinical importance of leptin receptor (LEPR) in regulating obesity and the fact that porcine LEPR expression is not known to be controlled by lncRNAs and miRNAs, we aim to characterize this gene as a potential target of SSC-miR-323 and the lncRNA TCONS_00010987. METHODS: Bioinformatics analyses revealed that lncRNA TCONS_00010987 and LEPR have SSC-miR-323-binding sites and that LEPR might be a target of lncRNA TCONS_00010987 based on cis prediction. Wild-type and mutant TCONS_00010987-target sequence fragments and wild-type and mutant LEPR 3'-UTR fragments were generated and cloned into pmiR-RB-REPORTTM-Control vectors to construct respective recombinant plasmids. HEK293T cells were co-transfected with the SSC-miR-323 mimics or a negative control with constructs harboring the corresponding binding sites and relative luciferase activities were determined. Tissue expression patterns of lncRNA TCONS_00010987, SSC-miR-323, and LEPR in Anqing six-end-white (AQ, the obese breed) and Large White (LW, the lean breed) pigs were detected by real-time quantitative polymerase chain reaction; backfat expression of LEPR protein was detected by western blotting. RESULTS: Target gene fragments were successfully cloned, and the four recombinant vectors were constructed. Compared to the negative control, SSC-miR-323 mimics significantly inhibited luciferase activity from the wild-type TCONS_00010987-target sequence and wild-type LEPR-3'-UTR (p<0.01 for both) but not from the mutant TCONS_00010987-target sequence and mutant LEPR-3'-UTR (p>0.05 for both). Backfat expression levels of TCONS_ 00010987 and LEPR in AQ pigs were significantly higher than those in LW pigs (p<0.01), whereas levels of SSC-miR-323 in AQ pigs were significantly lower than those in LW pigs (p<0.05). LEPR protein levels in the backfat tissues of AQ pigs were markedly higher than those in LW pigs (p<0.01). CONCLUSION: LEPR is a potential target of SSC-miR-323, and TCONS_00010987 might act as a sponge for SSC-miR-323 to regulate LEPR expression.

MiR-26a promotes apoptosis of porcine granulosa cells by targeting the 3ß-hydroxysteroid-Δ24-reductase gene.

OBJECTIVE: Apoptosis of ovarian granulosa cells (GCs) affects mammalian follicular development and fecundity. This study aimed to explore the regulatory relationship between microRNA-26a (miR-26a) and the 3ß-hydroxysteroid-Δ24-reductase gene (DHCR24) gene in porcine follicular granular cells (pGCs), and to provide empirical data for the development of methods to improve the reproductive capacity of pigs. METHODS: The pGCs were transfected with miR-26a mimic, miR-26a inhibitor and DHCR24-siRNA in vitro. The cell apoptosis rate of pGCs was detected by the flow cytometry. The secretion levels of estradiol (E2) and progesterone (P) in pGCs were detected by enzymelinked immunosorbent assay. Double luciferase validation system was used to detect the binding sites between miR-26a and DHCR24 3'-UTR region. Qualitative real-time polymerase chain reaction and Western blotting were used to verify the DHCR24 mRNA and protein expression in pGCs, respectively, after transfecting with miR-26a mimic and miR-26a inhibitor. RESULTS: Results showed that enhancement of miR-26a promoted apoptosis, and inhibited E2 and P secretion in pGCs. Meanwhile, inhibition of DHCR24 also upregulated the Caspase-3 expression, reduced the BCL-2 expression, promoted pGCs apoptosis, and inhibited E2 and P secretion in pGCs. There were the binding sites of miR-26a located within DHCR24 3'-UTR. Up-regulation of miR-26a inhibited DHCR24 mRNA and protein expression in pGCs. CONCLUSION: This study demonstrates that miR-26a can promote cell apoptosis and inhibit E2 and P secretion by inhibiting the expression of DHCR24 in pGCs.

Transcriptional analysis of deoxynivalenol-induced apoptosis of sow ovarian granulosa cell.

Litter size is one of the most important economic traits in pig production. Recent studies identified that deoxynivalenol (DON), a widespread toxin in fodder, was associated with animal prolificacy. However, the underlying mechanisms have not yet been completely elucidated. Here, we used porcine ovary granulosa cells (pGCs) as a vector to establish DON concentration-time models and performed cell morphology and transcriptome analysis to identify and analyse the effects of DON on reproductive performance in swine. The results showed that DON can induce morphological changes and apoptosis of pGCs, while inhibiting cell proliferation. Moreover, these effects of DON on pGCs were dose-dependent. After treatment of pGCs with different concentrations of DON, the percentage of cells in S phase and G2/M phase increased. RNA-seq analyses revealed 5,937 differentially expressed genes, of which 1995 were down-regulated and 3,942 were up-regulated after DON treatment. KEGG enrichment analysis indicated important metabolic pathways such as IL-17 signalling pathway, eukaryotic ribosome synthesis pathway, RNA transport pathway and RNA degradation. Based on our results, we speculate that the effects of DON are related to the DNA damage process. Our study provides novel insights and a foundation to further understand the effect of DON on swine prolificacy.

Combined microRNAome and transcriptome analysis of follicular phase and luteal phase in porcine ovaries.

The aim of this study was to explore the expression difference of miRNAs and mRNAs between the follicular phase (FP) and luteal phase (LP) in porcine ovaries and provide a theoretical basis for the research on mammalian reproductive regulation. RNA-Seq and miRNA-Seq were used to identify differentially expressed genes (DEGs) and miRNAs (DEMs) between the FP and LP in ovaries of six sows (3-year-old Yorkshire pigs with similar weights and same parities). Bioinformatic analysis was used to screen potential genes and miRNAs related to porcine ovarian function. Real-time qualitative PCR was used to validate the sequencing results. RNA-Seq results showed that 3,078 genes were up-regulated, and 1,444 genes were down-regulated in the LP compared with the FP, and DEGs were significantly enriched in 242 Gene Ontology (GO) terms and 33 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. miRNA-Seq identified 112 DEMs, of which 25 were up-regulated and 87 were down-regulated in the LP compared with the FP. We obtained 186 intersection genes (IGs) between the 4,522 DEGs and 2,444 target genes predicted from the 112 DEMs. After constructing a miRNA-gene-pathway network, we identified key miRNAs and genes including miR-17-3p, miR-214, miR-221-5p, miR-125b, FGF1, YWHAG, YWHAZ, FDFT1 and DHCR24, which are enriched in Hippo and PI3K-Akt signalling pathways, and various metabolic pathways. These results indicate that these key genes and miRNAs may play important roles in the developmental transition from FP to LP in porcine ovaries and represent candidate targets for further study.

MiR-222 inhibits apoptosis in porcine follicular granulosa cells by targeting the THBS1 gene.

Apoptosis of granulosa cells affects follicular atresia and reproduction and is regulated by miRNAs and the expression of certain genes. For the present study, we investigated the regulatory relationship between microRNA-222 (miR-222) and THBS1 in porcine follicular granulosa cells (pGCs) and its effects on apoptosis to provide empirical data for developing methods to improve pig fecundity. Results revealed that miR-222 promotes the proliferation of pGCs. MiRNA mimics and luciferase reporter assays revealed that miR-222 functions as an anti-apoptotic factor in pGCs. MiR-222 mimics in pGCs result in the upregulation of the anti-apoptotic BCL-2 gene, down-regulation of the proapoptotic caspase-3 gene, and inhibition of apoptosis. MiR-222 inhibitors reduced BCL-2 and had no significant effect on caspase-3. MiR-222 mimics promoted estrogen levels. Inhibition of THBS1 inhibited pGC apoptosis. Transfection of THBS1-siRNA reduced the proapoptotic BAX gene. MiR-222 can directly target the 3'-untranslated region of the THBS1 gene. MiR-222 mimics suppressed THBS1 mRNA and proteins, but these were upregulated by the miR-222 inhibitor. Transfection of THBS1-siRNA resulted in the inhibition of the miR-222 inhibitor, which suggests that miR-222 inhibits pGC apoptosis by targeting THBS1. These findings suggest that miR-222 and THBS1 play important roles in follicular atresia, ovarian development, and female reproduction.