Transcriptomic analysis of expression of genes regulating cell cycle progression in porcine ovarian granulosa cells during short-term in vitro primary culture

Preliminary construction of non-coding RNAs and ceRNA regulatory networks mediated by exosomes in porcine follicular fluid

BACKGROUND

Follicles are fundamental units of the ovary, regulated intricately during development. Exosomes and ovarian granulosa cells (OGCs) play pivotal roles in follicular development, yet the regulatory mechanisms governing exosomes remain elusive.

RESULTS

High-throughput sequencing was employed to evaluate the complete transcript expression profiles of six samples (three porcine ovarian granulosa cells-exosome co-culture samples (GCE) and three porcine ovarian granulosa cells (POGCs) samples). Differential expression analysis revealed 924 lncRNAs, 35 circRNAs, 49 miRNAs, and 9823 mRNAs in the GCE group. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses indicated enrichment of differentially expressed transcripts in pathways related to cell proliferation and apoptosis. Furthermore, a ceRNA regulatory network comprising 43 lncRNAs, 6 circRNAs, 11 miRNAs, and 126 mRNAs was constructed based on intergene co-expression correlations. Seven miRNAs associated with cell proliferation and apoptosis regulation were identified within this network, encompassing 92 subnet pairs as candidate genes for further exploration of exosome regulatory mechanisms. Additionally, preliminary verification at the cellular level demonstrated that exosomal miR-200b enhances the viability of POGCs.

CONCLUSIONS

Transcriptome analysis unveiled a pivotal candidate ceRNA network potentially implicated in exosome-mediated regulation of granulosa cell proliferation and apoptosis, thereby influencing porcine follicular development. These findings offer insights into the molecular mechanisms of follicular fluid exosome regulation, encompassing both coding and non-coding RNA perspectives.

Genome-Wide Association Analysis and Genetic Parameters for Egg Production Traits in Peking Ducks

Egg production traits are crucial in the poultry industry, including age at first egg (AFE), egg number (EN) at different stages, and laying rate (LR). Ducks exhibit higher egg production capacity than other poultry species, but the genetic mechanisms are still poorly understood. In this study, we collected egg-laying data of 618 Peking ducks from 22 to 66 weeks of age and genotyped them by whole-genome resequencing. Genetic parameters were calculated based on SNPs, and a genome-wide association study (GWAS) was performed for these traits. The SNP-based heritability of egg production traits ranged from 0.09 to 0.54. The GWAS identified nine significant SNP loci associated with AFE and egg number from 22 to 66 weeks. These loci showed that the corresponding alleles were positively correlated with a decrease in the traits. Moreover, three potential candidate genes (ENSAPLG00020011445, ENSAPLG00020012564, TMEM260) were identified. Functional enrichment analyses suggest that specific immune responses may have a critical impact on egg production capacity by influencing ovarian function and oocyte maturation processes. In conclusion, this study deepens the understanding of egg-laying genetics in Peking duck and provides a sound theoretical basis for future genetic improvement and genomic selection strategies in poultry.

Heterozygous gain of function variant in GUCY1A2 may cause autonomous ovarian hyperfunction

PURPOSE

The purpose of this study was to characterize the phenotype associated with a de novo gain-of-function variant in the GUCY1A2 gene.

METHODS

An individual carrying the de novo heterozygous variant c.1458G>T p.(E486D) in GUCY1A2 was identified by exome sequencing. The effect of the corresponding enzyme variant α2E486D/β1 was evaluated using concentration-response measurements with wild-type enzyme and the variant in cytosolic fractions of HEK293 cells, UV-vis absorbance spectra of the corresponding purified enzymes, and examination of overexpressed fluorescent protein-tagged constructs by confocal laser scanning microscopy.

RESULTS

The patient presented with precocious peripheral puberty resembling the autonomous ovarian puberty seen in McCune-Albright syndrome. Additionally, the patient displayed severe intellectual disability. In vitro activity assays revealed an increased nitric oxide affinity for the mutant enzyme. The response to carbon monoxide was unchanged, while thermostability was decreased compared to wild type. Heme content, susceptibility to oxidation, and subcellular localization upon overexpression were unchanged.

CONCLUSION

Our data define a syndromic autonomous ovarian puberty likely due to the activating allele p.(E486D) in GUCY1A2 leading to an increase in cGMP. The overlap with the ovarian symptoms of McCune-Albright syndrome suggests an impact of this cGMP increase on the cAMP pathway in the ovary. Additional cases will be needed to ensure a causal link.

Cumulus Cell and Oocyte Gene Expression in Prepubertal Gilts and Sows Identifies Cumulus Cells as a Prime Informative Parameter of Oocyte Quality

Cumulus cells (CCs) are pivotal during oocyte development. This study aimed to identify novel marker genes for porcine oocyte quality by examining the expression of selected genes in CCs and oocytes, employing the model of oocytes from prepubertal animals being of reduced quality compared to those from adult animals. Total RNA was extracted either directly after follicle aspiration or after in vitro maturation, followed by RT-qPCR. Immature gilt CCs accumulated BBOX1 transcripts, involved in L-carnitine biosynthesis, to a 14.8-fold higher level (p < 0.05) relative to sows, while for CPT2, participating in fatty acid oxidation, the level was 0.48 (p < 0.05). While showing no differences between gilt and sow CCs after maturation, CPT2 and BBOX1 levels in oocytes were higher in gilts at both time points. The apparent delayed lipid metabolism and reduced accumulation of ALDOA and G6PD transcripts in gilt CCs after maturation, implying downregulation of glycolysis and the pentose phosphate pathway, suggest gilt cumulus-oocyte complexes have inadequate ATP stores and oxidative stress balance compared to sows at the end of maturation. Reduced expression of BBOX1 and higher expression of CPT2 in CCs before maturation and higher expression of G6PD and ALDOA after maturation are new potential markers of oocyte quality.

Mutation of the ETS1 3'UTR interacts with miR-216a-3p to regulate granulosa cell apoptosis in sheep

ETS1, an important member of the ETS transcription factor family, is involved in a variety of physiological processes in living organisms, such as cell development, differentiation, proliferation and apoptosis, and is thought to be associated with embryonic development and reproduction. However, the polymorphism of ETS1 has been rarely studied, and its potential impact on the formation of reproductive traits in sheep remains unclear. Here, we first analyzed polymorphisms of ETS1 in a population of 382 small-tailed Han sheep with a lambing number record using the Kompetitive Allele Specific PCR (KASP) technique. The results showed the presence of a SNP locus rs161611767 (T > C) in the 3'UTR of ETS1. The association analysis showed the lambing number of first, second and third parity in the individuals with the CC genotype (2.51 ± 0.108, 2.51 ± 0.179, 1.27 ± 0.196) was higher than that of individuals with the TT genotype (1.79 ± 0.086, 1.56 ± 0.102, 0.56 ± 0.100) (P < 0.05). Then, molecular biotechnologies were used to investigate the effects of the EST1 rs161611767 mutant locus on host gene expression in sheep and the underlying mechanism of its effect on sheep reproduction. The RT‒qPCR results showed that the expression of ETS1 was higher in individuals with the CC genotype than in those with the TT genotype (P < 0.05). The dual luciferase reporter assay showed that the luciferase activity of ETS1 in sheep with the TT genotype was decreased compared to CC genotype (P < 0.05), confirming the existence of EST1 rs161611767 in the 3'UTR as a functional SNP. Given that the 3'UTR is an important regulatory region of gene transcription and translation, we performed bioinformatics prediction and confirmed that the SNP rs161611767 of ETS1 was a direct functional target of miR-216a-3p using dual luciferase activity assay, and the binding capacity of allele T was stronger than that of allele C. Subsequently, the cell transfection results showed that miR-216a-3p suppressed the endogenous expression of ETS1 in sheep primary granulosa cells (GCs). Finally, CCK-8, EdU, WB detection of marker proteins and flow cytometry were used to detect the effects of miR-216a-3p on GCs viability and proliferation/apoptosis, respectively. The results showed that miR-216a-3p inhibited the proliferation of GCs while promoting apoptosis of GCs. In conclusion, these results demonstrate that the SNP rs161611767 of ETS1 is associated with lambing number in small-tailed Han sheep, and miR-216a-3p can act as a regulatory element binding to the T mutation in rs161611767 to regulate ETS1 expression and affect GCs development, which may indirectly affect the number of lambs in sheep. These studies provide evidence for the involvement of ETS1 polymorphisms in sheep reproduction and are expected to provide new insights to elucidate the molecular genetic mechanisms of lambing traits in sheep.

Expression of genes regulating cell division in porcine follicular granulosa cells

BACKGROUND

Cell cycle regulation influences the proliferation of granulosa cells and affects many processes related to ovarian folliclular growth and ovulation. Abnormal regulation of the cell cycle can lead to many diseases within the ovary. The aim of this study was to describe the expression profile of genes within granulosa cells, which are related to the formation of the cytoskeleton, organization of cell organelles inside the cell, and regulation of cell division. Established in vitro primary cultures from porcine ovarian follicle granulosa cells were maintained for 48, 96, 144 h and evaluated via microarray expression analysis.

RESULTS

Analyzed genes were assigned to 12 gene ontology groups "actin cytoskeleton organization", "actin filament organization", "actin filament-based process", "cell-matrix adhesion", "cell-substrate adhesion", "chromosome segregation", "chromosome separation", "cytoskeleton organization", "DNA integrity checkpoint", "DNA replication initiation", "organelle fision", "organelle organization". Among the genes with significantly changed expression, those whose role in processes within the ovary are selected for consideration. Genes with increased expression include (ITGA11, CNN1, CCl2, TPM2, ACTN1, VCAM-1, COL3A1, GSN, FRMD6, PLK2). Genes with reduced expression inlcude (KIF14, TACC3, ESPL1, CDC45, TTK, CDC20, CDK1, FBXO5, NEK2-NIMA, CCNE2). For the results obtained by microarray expressions, quantitative validation by RT-qPCR was performed.

CONCLUSIONS

The results indicated expression profile of genes, which can be considered as new molecular markers of cellular processes involved in signaling, cell structure organization. The expression profile of selected genes brings new insight into regulation of physiological processes in porcine follicular granulosa cells during primary in vitro culture.

Exploration of the Antioxidant Effect of Spermidine on the Ovary and Screening and Identification of Differentially Expressed Proteins

Spermidine is a naturally occurring polyamine compound that has many biological functions, such as inducing autophagy and anti-inflammatory and anti-aging effects. Spermidine can affect follicular development and thus protect ovarian function. In this study, ICR mice were fed exogenous spermidine drinking water for three months to explore the regulation of ovarian function by spermidine. The results showed that the number of atretic follicles in the ovaries of spermidine-treated mice was significantly lower than that in the control group. Antioxidant enzyme activities (SOD, CAT, T-AOC) significantly increased, and MDA levels significantly decreased. The expression of autophagy protein (Beclin 1 and microtubule-associated protein 1 light chain 3 LC3 II/I) significantly increased, and the expression of the polyubiquitin-binding protein p62/SQSTM 1 significantly decreased. Moreover, we found 424 differentially expressed proteins (DEPs) were upregulated, and 257 were downregulated using proteomic sequencing. Gene Ontology and KEGG analyses showed that these DEPs were mainly involved in lipid metabolism, oxidative metabolism and hormone production pathways. In conclusion, spermidine protects ovarian function by reducing the number of atresia follicles and regulating the level of autophagy protein, antioxidant enzyme activity, and polyamine metabolism in mice.

Bisphenol A and bisphenol S both disrupt ovine granulosa cell steroidogenesis but through different molecular pathways

BACKGROUND

Ovarian granulosa cells (GC) are essential for the development and maturation of a proper oocyte. GC are sensitive to endocrine disruptors, including bisphenol A (BPA) and its analogue bisphenol S (BPS), plasticisers present in everyday consumer products. BPA exhibits greater binding affinity for the membrane oestrogen receptor (GPER) than for the nuclear oestrogen receptors (ERα and ERβ). Here, we analysed the effects of BPA and BPS on the steroidogenesis of ovine GC in vitro, as well as their early mechanisms of action, the ovine being a relevant model to study human reproductive impairment. Disruption of GC steroidogenesis might alter oocyte quality and consequently fertility rate. In addition, we compared the effects of a specific GPER agonist (G-1) and antagonist (G-15) to those of BPA and BPS. Ewe GC were cultured with BPA or BPS (10 or 50 µM) or G-1 (1 µM) and/or G-15 (10 µM) for 48 h to study steroidogenesis.

RESULTS

Both BPA and BPS (10 µM) altered the secretion of progesterone, however, only BPS (10 µM) affected oestradiol secretion. RNA-seq was performed on GC after 1 h of culture with BPA or BPS (50 µM) or G-1 (10 µM), followed by real-time PCR analyses of differentially expressed genes after 12, 24 and 48 h of culture. The absence of induced GPER target genes showed that BPA and BPS did not activate GPER in GC after 1 h of treatment. These molecules exhibited mainly independent early mechanisms of action. Gene ontology analysis showed that after 1 h of treatment, BPA mainly disrupted the expression of the genes involved in metabolism and transcription, while BPS had a smaller effect and impaired cellular communications. BPA had a transient effect on the expression of CHAC1 (NOTCH signalling and oxidative balance), JUN (linked to MAPK pathway), NR4A1 (oestradiol secretion inhibition), ARRDC4 (endocytose of GPCR) and KLF10 (cell growth, differentiation and apoptosis), while expression changes were maintained over time for the genes LSMEM1 (linked to MAPK pathway), TXNIP (oxidative stress) and LIF (cell cycle regulation) after 12 and 48 h, respectively.

CONCLUSION

In conclusion, although they exhibited similar effects, BPA and BPS impaired different molecular pathways in GC in vitro. New investigations will be necessary to follow the temporal changes of these genes over time, as well as the biological processes involved.

New Gene Markers Expressed in Porcine Oviductal Epithelial Cells Cultured Primary In Vitro Are Involved in Ontological Groups Representing Physiological Processes of Porcine Oocytes

Changes that occur within oviducts after fertilization are dependent on post-ovulation events, including oocyte-oviduct interactions. Although general processes are well-defined, the molecular basis are poorly understood. Recently, new marker genes involved in ‘cell development’, ‘cell growth’, ‘cell differentiation’ and ‘cell maturation’ processes have been identified in porcine oocytes. The aim of the study was to assess the expression profile of genes in primary in vitro cultured oviductal epithelial cells (OECs), clustered in Gene Ontology groups which enveloped markers also identified in porcine oocytes. OECs (from 45 gilts) were surgically removed and cultured in vitro for ≤ 30 days, and then subjected to molecular analyses. The transcriptomic and proteomic profiles of cells cultured during 7, 15 and 30 days were investigated. Additionally, morphological/histochemical analyzes were performed. The results of genes expression profiles were validated after using RT-qPCR. The results showed a significant upregulation of UNC45B, NOX4, VLDLR, ITGB3, FMOD, SGCE, COL1A2, LOX, LIPG, THY1 and downregulation of SERPINB2, CD274, TXNIP, CELA1, DDX60, CRABP2, SLC5A1, IDO1, ANPEP, FST. Detailed knowledge of the molecular pathways occurring in the OECs and the gametes that contact them may contribute both to developments of basic science of physiology, and new possibilities in advanced biotechnology of assisted reproduction.

Effect of selected bisphenol derivatives on nuclear receptor expression in ovarian cell line COV434

Objectives. Bisphenol A (BPA), as an indispensable plastic additive, has also been proven as an endocrine disruptor associated with adverse health effects including impaired ovarian function and cancer. Due to the restrictions of its usage, several analogs have been employed to replace BPA. Although many studies revealed a harmfulness in the biological effects of BPA analogs, their specific targets remain largely unknown. Nuclear receptors (NRs) may be one of the most important targets of bisphenols. Therefore, in this study, our attention was directed to explore the effect of BPA and its analogs, AF and S, on the mRNA expression of selected NRs involved in the steroidogenic and carcinogenic pathways in the human granulosa cell line COV434. The NRs investigated included: thyroid hormone receptor α (THRA), peroxisome proliferator activating receptor β/δ (PPARD), retinoid X receptor α (RXRA), chicken ovalbumin upstream promoter-transcription factor II (COUPTFII), nuclear receptor-related protein 1 (NURR1), and liver receptor homolog-1 (LRH1).Methods. COV434 cells were treated with the bisphenols at the concentrations of 10^-9 M, 10^-7 M, and 10^-5 M, and after 24 and 48 h, cell viability was monitored by the MTS assay and gene expressions were analyzed using RT-qPCR.Results. Bisphenol treatment did not alter the COV434 cell viability. After 24 h, the expression of neither of the NRs was changed. Likewise, after 48 h, the expression of the selected genes was not altered. However, both BPAF and BPS increased, at the highest concentration (10^-5 M) used, the mRNA levels of both PPARD and NURR1 NRs after 48 h of the treatment. In the BPA-treated groups, no significant upregulation was observed.Conclusions. In the present study, the effect of bisphenols on COUP-TFII, Nurr1, and LRH-1 NRs was investigated for the first time. Although generally we did not observe that BPs provoked any alterations in the expression of the selected NRs in COV434 cells, at specific concentrations and time points they might alter mRNA expression of certain NRs (NURR1, PPARD).

Identification of Crucial lncRNAs, miRNAs, mRNAs, and Potential Therapeutic Compounds for Polycystic Ovary Syndrome by Bioinformatics Analysis

BACKGROUND

This study was aimed at mining crucial long noncoding RNAs (lncRNAs), microRNAs (miRNAs), and messenger RNAs (mRNAs) for the development of polycystic ovary syndrome (PCOS) based on the coexpression and the competitive endogenous RNA (ceRNA) theories and investigating the underlying therapeutic drugs that may function by reversing the expression of lncRNAs, miRNAs, and mRNAs.

METHODS

RNA (GSE106724, GSE114419, GSE137684, and GSE138518) or miRNA (GSE84376 and GSE138572) expression profile datasets of PCOS patients were downloaded from the Gene Expression Omnibus database. The weighted gene coexpression network analysis (WGCNA) using four RNA datasets was conducted to construct the lncRNA-mRNA coexpression networks, while the common differentially expressed miRNAs in two miRNA datasets and module RNAs were used to establish the ceRNA network. A protein-protein interaction (PPI) network was created to explore the potential interactions between genes. Gene Ontology and KEGG pathway enrichment analyses were performed to explore the functions of genes in networks. Connectivity Map (CMap) and Comparative Toxicogenomics Database (CTD) analyses were performed to identify potential therapeutic agents for PCOS.

RESULTS

Three modules (black, magenta, and yellow) were identified to be PCOS-related after WGCNA analysis, in which KLF3-AS1-PLCG2, MAPKAPK5-AS1-MAP3K14, and WWC2-AS2-TXNIP were important coexpression relationship pairs. WWC2-AS2-hsa-miR-382-PLCG2 was a crucial ceRNA loop in the ceRNA network. The PPI network showed that MAP3K14 and TXNIP could interact with hub genes PLK1 (degree = 21) and TLR1 (degree = 18), respectively. These genes were enriched into mitosis (PLK1), immune response (PLCG2 and TLR1), and cell cycle (TXNIP and PLK1) biological processes. Ten small molecule drugs (especially quercetin) were considered to be therapeutical for PCOS.

CONCLUSION

Our study may provide a novel insight into the mechanisms and therapy for PCOS.