INTRODUCTION
Oocyte quality and maturation not only affect the embryo and fertilization success but also have long-term impacts on fetal growth and development. Notably, mature oocytes arrested at metaphase II (MII) can be successfully fertilized within a restricted time window after ovulation. When fertilization does not occur within an optimal time frame, MII oocytes experience a process of degradation known as “post-ovulatory aging”(POA) [
1]. The POA refers to the oocytes released from the ovary. Over time, the cytoplasm ages and eventually dies, which is inevitable[
1]. In mice, the quality of oocytes rapidly deteriorates six hours after ovulation if they are not fertilized, with them subsequently being lost after 12 hours [
2]. For human oocytes, the time window for fertilization is 4–12 hours after ovulation [
3]. However, women do not show natural visual signs of ovulation; fertilization can occur hours later and can involve aged oocytes and freshly ejaculated spermatozoa.
POA induces many abnormal effects on cell biology, including partial exocytosis of cortical granules [4, 5], hardening of the zona pellucida [
5], decline in maturation-promoting factor and MAPK levels [
6], abnormalities in the cytoskeleton, and condensation of the chromosome [
7]. It can also induce mitochondrial dysfunction, leading to apoptosis [8-10]; perturbation of Ca
2+ homeostasis; oxidative damage to lipids, proteins, and DNA components of the cell [
11]; and epigenetic changes [
12]. Furthermore, POA reduces the fertilization rate and embryo quality and increases the likelihood of abnormalities in the offspring. Therefore, the timely fertilization of oocytes is crucial to ensure embryonic development. Although oocytes fertilized via intracytoplasmic sperm injection (ICSI) can develop into viable embryos, the procedure takes 6–12 hours and, hence, induces POA [
13]. Aging-induced defects such as oxidative stress, mitochondrial dysfunction, and chromosomal abnormalities have been detected in oocytes during assisted reproductive technology (ART) procedures, which can result in poor oocyte quality, lower fertilization rates, embryo development aberrancy, and even unhealthy offsprings [14-16]. Therefore, developing strategies to improve in-vitro fertilization (IVF) outcomes remains imperative.
Research comparing the anti-aging ability of in vivo- and vitro-matured porcine oocytes has not yet been reported. Compared with other species, porcine oocytes are more similar to human oocytes in many aspects, including the closed oocyte volume (120–125 μm in diameter), average time for oocyte maturation (40–44 h for pig, 40 h for human) [
17], similar core transcriptional network required to maintain pluripotency [
18], and similar developmental stage of embryonic genome activation [17, 19]. Moreover, they both contain a large quantity of endogenous lipids [20, 21]. Therefore, porcine oocytes are a good model for human reproductive research and clinical-assisted reproductive technology applications. Therefore, we collected in vivo- and in vitro-matured and aged porcine oocytes for RNA sequencing to explore the effects of POA on porcine oocytes and the key pathways and genes that provide in vivo-matured oocytes with stronger anti-aging ability.
In this study, we used transcriptomic analysis and experimental validation to study the similarities and differences between in vivo- and vitro-matured porcine oocytes during POA. Addressing this question will help us understand the detailed mechanisms of POA and provide strategies for researching oocyte aging.
DISCUSSION
The exact mechanisms underlying ovarian aging are poorly understood. Oocyte senescence is thought to play a central role in ovarian aging. Previous studies have compared the differences in gene expression profiles in GV- or MII-stage oocytes from young and old individuals [24, 25], which were limited by the differences in in vivo and in vitro culture environments that are hard to ignore. Herein, we used RNA sequencing to determine the gene expression profiles of oocytes in the following four groups: fresh in vivo, fresh in vitro, aged in vivo, and aged in vitro. The similar and different effects of POA on in vivo- and in vitro-matured oocytes were explored.
Most DEGs appeared between fresh
in vivo and
in vitro oocytes, with a total of 5880 (
Figure 2B). Different culture environments had a profound influence on the quality and anti-aging abilities of porcine oocytes. To exclude the influence of
in vitro conditions, DEGs between F_vivo and A_vivo and DEGs between F_vitro and A_vitro were intersected to explore the co-effects of POA.
We performed enrichment analyses of common DEGs caused by POA to explore the effects of POA on both
in vivo and
in vitro-matured oocytes. The results revealed that many DEGs were related to important cellular function. These DEGs were enriched for intracellular, organelle, and membrane-bounded organelle. Among them, we are concerned that mitochondria-related terms are enriched (
Figure 3B).
Our research showed that the dysfunction of mitochondrial function, including electron transport chain, cell redox homeostasis and oxidative phosphorylation, contributed to the aging of oocytes. As oocyte maturation requires a large amount of ATP for continuous transcription and translation, the availability of the right number of functional mitochondria is crucial[
26]. It has been speculated that the higher ATP content in human and mouse MII oocytes was associated with better embryo potential for development and implantation[
27]. The normal function of mitochondria is the guarantee of oocyte quality and embryo developmental potentiality, while oocyte aging is related to mitochondrial dysfunction and disturbance of energy metabolism[
28]. The protein machinery to control the health of the mitochondria via unfolded protein response and mitophagy may be compromised in oocytes from aged females, which may result in inefficient clearance of dysfunctional mitochondria and reduced oocyte quality[
29]. This highlights the importance of mitochondrial function in aging.
Abnormal distribution of mitochondria as well as mitochondrial dysfunction, resulting in severely impaired germinal vesicle breakdown (GVBD) of mouse oocytes. This conclusion was also confirmed by our data, as KEGG enrichment analysis showed that the meiosis pathway was enriched, and the expression of related genes was reduced in aging oocytes[
30](
Figure 4A, 4C).
Another enriched pathway in both
in vivo and
in vitro-matured oocytes was protein export (
Figure 4A, 4B). This pathway contains DEGs named signal recognition particle (SRP), which was a ribonucleoprotein particle crucial for co-translational targeting of secretory and membrane proteins to the endoplasmic reticulum [
31]. The eukaryotic SPC is composed of five subunits including two isoforms of catalytic subunits SEC11A and SEC11C and three regulatory subunits including signal peptidase complex subunit 1 (SPCS1), SPCS, and SPCS3 [32, 33]. Among the SPC subunits, SEC11 and SPCS3 are essential for signal peptidase activity and cell survival[34, 35]. SPCS2 interacts with the b subunit of SEC61 translocon complex likely facilitating co-translational signal peptidase processing[
36] and, at high-temperature conditions, regulates the catalytic activity of the SPC and viability of yeast cells[
37]. A previous study has shown that nuclear protein export is a common hallmark of pathological and physiological aging in the Hutchinson Gilford syndrome cellular phenotype of normal fibroblasts [
38].
The co-effects of POA on in vivo and in vitro oocytes were observed through enrichment analysis, and then we performed WGCNA to explore the unique gene expression patterns of
in vivo-matured oocytes during POA. 22 co-expression network modules were determined, and successfully mined the specific module related to A_vivo. WGCNA divides modules into soft thresholds, which reflect the effectiveness of biological networks more effectively than hard thresholds[
39]. Functional enrichment analysis using GO/KEGG identified RNA binding and spliceosome as key pathways against POA in vivo (
Figure 7A, 7B).
RNA-binding proteins (RBPs) achieve their biological function essentially by post-transcriptional gene regulation[
40]. In eukaryotic cells, following RNA polymerase II-mediated transcript synthesis, RBPs dictate extensive pre-mRNA processing by interacting with the target RNA and partner proteins[
41]. Alternative splicing generates diverse transcripts by removing introns and splicing exons. They are responsible for adding modifications to the transcript that affect stability and translation efficiency, including 5’-end capping and 3’-end polyadenylations. RBPs play a crucial role in the transport of transcripts from the nucleus to the cytoplasm, where protein synthesis occurs. Additionally, RBPs are involved in modulating the localization and half-life of mRNA within the cell by either promoting or delaying transcript degradation. Hence, RBPs have a significant impact on almost every aspect of RNA biology[42, 43].
Decades of work on aging have shed light into the fundamental role played within this context by a class of proteins termed RBPs[44, 45]. Loss of intracellular RBP AU-rich-element factor-1 will alter post-transcriptional regulation of targets particularly relevant for protection of genomic integrity and gene regulation, thus concurring to responses related to oxidative stress and accelerated aging[
46]. RBP PUM1 overexpression protected MSCs against H
2O
2-induced cellular senescence by suppressing TLR4-mediated NF-κB activity[
47].
Spliceosome induced by pharmacological and genetic inhibition of spliceosome genes, have been reported to trigger cell senescence, suggesting a key role of spliceosome genes as a gatekeeper [
48]. Pre-mRNA splicing is fundamental for gene expression and regulation. Similar results were observed in transcriptome sequencing data from human mature defective oocytes, with spliceosomes being the most abundant pathway [
49].
The hub genes
DNAJC7,
DDX18, and
IK were screened out by constructing a gene regulatory relationship network. These hub genes are involved to a certain extent in maintaining genome stability, DNA damage repair and other life processes.
DNAJC7 participates in
p53/MDM2 negative feedback regulatory pathway, and thereby enhancing the stability and activity of tumor suppressor
p53 which promotes apoptosis to eliminate cells with seriously damaged DNA to maintain genomic integrity[
50].
DDX18 depletion leads to γH2AX accumulation and genome instability[
51].
IK is a splicing factor that promotes spliceosome activation and contributes to pre-mRNA splicing[
52]. And all the three genes had higher expression levels in A_vivo than A_vitro group.
These results demonstrate the impact of post-transcriptional regulation on cellular senescence. In vivo-matured oocytes have higher expression levels on these pathways than in vitro-matured oocytes. We predict that in vivo-matured oocytes can resist POA through post-transcriptional regulation of gene expression.
To summarize our findings, we have observed that POA has an impact on the quality of porcine oocytes. This is likely due to its effects on mitochondrial function and protein export. We have also noted that there are variations in the expression patterns of in vivo and in vitro oocytes during POA, particularly in pathways related to post-transcriptional regulation such as RNA-binding and spliceosome pathways (
Figure 9). Our data serves as a foundation for understanding the mechanisms that underlie POA, but more research is needed to fully explore the relationship between these factors.
Figure 1.
Experimental flowchart. In brief, GV-stage oocytes were collected from pig ovaries and cultured for 44 and 92 h to obtain two groups of oocytes, fresh in vitro and aged in vitro, respectively. Fresh in vivo and aged in vivo oocytes were obtained by in vitro culture of MII-stage oocytes collected from porcine fallopian tubes for 0 and 48 h, respectively.
Figure 1.
Experimental flowchart. In brief, GV-stage oocytes were collected from pig ovaries and cultured for 44 and 92 h to obtain two groups of oocytes, fresh in vitro and aged in vitro, respectively. Fresh in vivo and aged in vivo oocytes were obtained by in vitro culture of MII-stage oocytes collected from porcine fallopian tubes for 0 and 48 h, respectively.
Figure 2.
Global transcriptome analysis of fresh and aged in vivo and vitro-matured porcine oocytes. (A) Spearman’s correlation heatmap shows the reproducibility between repetitive samples and the difference between groups. The color gradient indicates the magnitude of the correlation coefficient. (B) Scatterplot of differentially expressed genes between each two groups. The criteria for differentially expressed genes were FDR<0.05, and fold change>2. (C) Heatmap of mRNA expression profiles in fresh and aging oocytes, showing changes in a subset of genes in response to POA. The color key (from blue to red) of the Z-score value indicates low to high expression levels. (D) Relative mRNA expression levels of BTF3, CEP95, CMC1, DBI, and NOUFB5 in four groups.
Figure 2.
Global transcriptome analysis of fresh and aged in vivo and vitro-matured porcine oocytes. (A) Spearman’s correlation heatmap shows the reproducibility between repetitive samples and the difference between groups. The color gradient indicates the magnitude of the correlation coefficient. (B) Scatterplot of differentially expressed genes between each two groups. The criteria for differentially expressed genes were FDR<0.05, and fold change>2. (C) Heatmap of mRNA expression profiles in fresh and aging oocytes, showing changes in a subset of genes in response to POA. The color key (from blue to red) of the Z-score value indicates low to high expression levels. (D) Relative mRNA expression levels of BTF3, CEP95, CMC1, DBI, and NOUFB5 in four groups.
Figure 3.
GO enrichment analysis of common DEGs caused by POA in both in vivo and in vitro matured oocytes (A) Differential gene Venn diagram. The intersection of DEGs between F_vivo and A_vivo and DEGs between F_vitro and A_vitro were defined as common DEGs. (B) GO enrichment analysis on common DEGs showing the top 15 enriched GO items. (C) The heat map for gene expression levels associated with electron transport chain in enriched GO terms. (D) The heat map for gene expression levels associated with Cell REDOX homeostasis in enriched GO terms.
Figure 3.
GO enrichment analysis of common DEGs caused by POA in both in vivo and in vitro matured oocytes (A) Differential gene Venn diagram. The intersection of DEGs between F_vivo and A_vivo and DEGs between F_vitro and A_vitro were defined as common DEGs. (B) GO enrichment analysis on common DEGs showing the top 15 enriched GO items. (C) The heat map for gene expression levels associated with electron transport chain in enriched GO terms. (D) The heat map for gene expression levels associated with Cell REDOX homeostasis in enriched GO terms.
Figure 4.
KEGG enrichment analysis of common DEGs caused by POA in both in vivo and vitro-matured oocytes (A) KEGG enrichment analysis on common DEGs showing the top 15 enriched KEGG pathways. (B) The heat map for gene expression levels associated with protein export in enriched KEGG pathway. (C) The heat map for gene expression levels associated with oocyte meiosis in enriched KEGG pathway. (D) The heat map for gene expression levels associated with oxidative phosphorylation in enriched KEGG pathway.
Figure 4.
KEGG enrichment analysis of common DEGs caused by POA in both in vivo and vitro-matured oocytes (A) KEGG enrichment analysis on common DEGs showing the top 15 enriched KEGG pathways. (B) The heat map for gene expression levels associated with protein export in enriched KEGG pathway. (C) The heat map for gene expression levels associated with oocyte meiosis in enriched KEGG pathway. (D) The heat map for gene expression levels associated with oxidative phosphorylation in enriched KEGG pathway.
Figure 5.
Construction of weighted gene co-expression network analysis (WGCNA). (A) Module hierarchical clustering tree. Gene modules are divided according to the clustering relationship between genes. Genes with similar expression patterns will be classified into the same module. The branches of the cluster tree are cut and distinguished to generate different modules. Each color represents a module, gray Indicates genes that cannot be assigned to any module. After the preliminary module division, the preliminary divided “Dynamic Tree Cut” is obtained. Modules with similar expression patterns are then merged based on the similarity of module feature values to obtain the final dividend “Merged dynamic”. (B) The bar plot shows the number of genes in every module.
Figure 5.
Construction of weighted gene co-expression network analysis (WGCNA). (A) Module hierarchical clustering tree. Gene modules are divided according to the clustering relationship between genes. Genes with similar expression patterns will be classified into the same module. The branches of the cluster tree are cut and distinguished to generate different modules. Each color represents a module, gray Indicates genes that cannot be assigned to any module. After the preliminary module division, the preliminary divided “Dynamic Tree Cut” is obtained. Modules with similar expression patterns are then merged based on the similarity of module feature values to obtain the final dividend “Merged dynamic”. (B) The bar plot shows the number of genes in every module.
Figure 6.
(A)Module–trait relationships. Each row presents a module eigengene, each column presents a trait. Each cell contains the corresponding correlation and p-value (“*” means p-value < 0.05, “**” means p-value < 0.01). The table is color-coded by correlation according to the color legend. (B) The heat map for gene expression levels of genes in the red module among the four groups. (C) The (module membership) MM- gene significance (GS) correlation of genes in the red module. Use the GS and MM values of genes to analyze the correlation between each trait and the module. Modules with high correlation play an important biological role in the trait. (D) Intramodular connectivity (K.in) and GS correlation analysis of genes in the red module. Analyze the association between modules, genes, and traits using the connectivity of genes within the module (K.in value) and the correlation value between genes and traits (GS).
Figure 6.
(A)Module–trait relationships. Each row presents a module eigengene, each column presents a trait. Each cell contains the corresponding correlation and p-value (“*” means p-value < 0.05, “**” means p-value < 0.01). The table is color-coded by correlation according to the color legend. (B) The heat map for gene expression levels of genes in the red module among the four groups. (C) The (module membership) MM- gene significance (GS) correlation of genes in the red module. Use the GS and MM values of genes to analyze the correlation between each trait and the module. Modules with high correlation play an important biological role in the trait. (D) Intramodular connectivity (K.in) and GS correlation analysis of genes in the red module. Analyze the association between modules, genes, and traits using the connectivity of genes within the module (K.in value) and the correlation value between genes and traits (GS).
Figure 7.
Enrichment analysis of genes in the red module. (A) GO enrichment analysis of genes in the red module showing the top 15 enriched GO terms (B) KEGG enrichment analysis of genes in the red module showing the top 15 enriched KEGG pathways.
Figure 7.
Enrichment analysis of genes in the red module. (A) GO enrichment analysis of genes in the red module showing the top 15 enriched GO terms (B) KEGG enrichment analysis of genes in the red module showing the top 15 enriched KEGG pathways.
Figure 8.
Hub genes in the red module from the gene regulatory relationship network. (A) Construction of gene regulatory network. Each node in the figure is a gene, and each line represents the regulatory relationship between the nodes. The darker and larger the node color is, the higher the abundance and the stronger the connectivity. The weight value defines the color and thickness of the line. The darker the color and the thicker the line, the stronger the regulatory relationship between genes. (B) The boxplot shows the expression of hub genes in the A_vivo and A_vitro groups.
Figure 8.
Hub genes in the red module from the gene regulatory relationship network. (A) Construction of gene regulatory network. Each node in the figure is a gene, and each line represents the regulatory relationship between the nodes. The darker and larger the node color is, the higher the abundance and the stronger the connectivity. The weight value defines the color and thickness of the line. The darker the color and the thicker the line, the stronger the regulatory relationship between genes. (B) The boxplot shows the expression of hub genes in the A_vivo and A_vitro groups.
Figure 9.
Schematic representation depicting the same and different effects of POA on in vivo and in vitro-matured porcine oocytes. POA adversely affects the quality of porcine oocytes, both in vivo and in vitro. This is most likely due to the impairment of mitochondrial function and protein export. RNA-binding and spliceosome pathways were the most differentially enriched pathways between oocytes that matured in vivo and those that matured in vitro.
Figure 9.
Schematic representation depicting the same and different effects of POA on in vivo and in vitro-matured porcine oocytes. POA adversely affects the quality of porcine oocytes, both in vivo and in vitro. This is most likely due to the impairment of mitochondrial function and protein export. RNA-binding and spliceosome pathways were the most differentially enriched pathways between oocytes that matured in vivo and those that matured in vitro.
Table 1.
Information of primers used for qPCR.
Table 1.
Information of primers used for qPCR.
Genes |
Primer sequences(5’-3’)F |
Primer sequences(5’-3’)R |
BTF3 |
GTGTGTGCGCCTTATCTCAG |
GTTTGGCGAGTTTCTCCTGG |
CEP95 |
AGAGGGCAGGAGAGAGGTTA |
ACATCCTCCTCTTCACAGCC |
CMC1 |
CGCAGAACAGCATCTCAGAC |
TCCAGAGTCCTTGCAGCATT |
DBI |
ACAGCCACTACAAACAAGCG |
ACGCTTTCATGGCATCTTCC |
NDUFB5 |
GCTTTGCCCTCAGTCAACAT |
CATGGCTACTATGGGCGAGA |
18s |
CGCGGTTCTATTTTGTTGGT |
AGTCGGCATCGTTTATGGTC |