Research

The assembly and activation of the PANoptosome promote porcine granulosa cell programmed cell death during follicular atresia

1 ,1 ,1 ,1 ,1 ,1 ,1 ,2 ,3

Abstract

Background

Follicular atresia significantly impairs female fertility and hastens reproductive senescence. Apoptosis of granulosa cells is the primary cause of follicular atresia. Pyroptosis and necroptosis, as additional forms of programmed cell death, have been reported in mammalian cells. However, the understanding of pyroptosis and necroptosis pathways in granulosa cells during follicular atresia remains unclear. This study explored the effects of programmed cell death in granulosa cells on follicular atresia and the underlying mechanisms.

Results

The results revealed that granulosa cells undergo programmed cell death including apoptosis, pyroptosis, and necroptosis during follicular atresia. For the first time, we identified the formation of a PANoptosome complex in porcine granulosa cells. This complex was initially identified as being composed of ZBP1, RIPK3, and RIPK1, and is recruited through the RHIM domain. Additionally, we demonstrated that caspase-6 is activated and cleaved, interacting with RIPK3 as a component of the PANoptosome. Heat stress may exacerbate the activation of the PANoptosome, leading to programmed cell death in granulosa cells.

Conclusions

Our data identified the formation of a PANoptosome complex that promoted programmed cell death in granulosa cells during the process of follicular atresia. These findings provide new insights into the molecular mechanisms underlying follicular atresia.

Introduction

Folliculogenesis, an ongoing and intricately regulated process throughout an individual’s reproductive lifespan, is marked by significant redundancy, where the number of activated follicles is considerably higher than that of the follicles that ultimately ovulate. Notably, over 99% of these activated follicles naturally undergo atresia [1, 2], which is a complex natural process that acts as the intrinsic regulator of the follicle pool, thus influencing follicular development and encompassing diverse mechanisms of regulated cell death [3].

Changes in the follicular environment are critical in influencing the process of follicular atresia. Variations in hormonal levels, such as reduced follicle stimulating hormone (FSH) and a decrease in luteinizing hormone (LH) pulse amplitude, can hinder the development and survival of follicles [4,5,6]. Elevated oxidative stress within the follicle often leads to apoptosis of granulosa cells (GCs) and oocyte [7, 8]. Furthermore, disruptions in nutrient availability and blood supply can weaken the follicle’s microenvironment, diminishing its capacity to sustain vital cellular processes [9,10,11]. These environmental changes collectively accelerate the atresia process, leading to the degeneration of non-dominant follicles and a consequent decline in ovarian reserve. Among these mechanisms, apoptosis of GCs is recognized as the principal initiator of follicular atresia [12,13,14,15,16]. However, recent studies indicate that other forms of GC death, such as necroptosis and pyroptosis, could significantly influence follicular function during follicular atresia [17,18,19,20]. These three cell death forms are collectively called programmed cell death (PCD).

PCD in mature organisms is crucial for controlling infection through the elimination intracellular pathogen replication sites and preventing cancer by removing emerging neoplastic cells [21, 22]. Unlike caspase-3-driven apoptosis, necroptosis and pyroptosis are orchestrated by mixed lineage kinase domain-like protein (MLKL) and gasdermin proteins, such as gasdermin D (GSDMD), that form pores in cellular membranes upon activation [23]. Specifically, RIPK3 phosphorylation of MLKL induces significant structural changes and therefore promotes its migration to the plasma membrane and consequent membrane permeabilization [24]. In addition, caspase 1 cleaves GSDMD into N-terminal and C-terminal fragments, with the N-terminal segments assembling into transmembrane pores on the inner cell membrane, disrupting the cytosol-extracellular space boundary [25, 26]. Recent research has characterized the complex interplay and regulatory mechanisms among apoptosis, necroptosis, and pyroptosis.

Recent research advances have introduced the concept of the PANoptosome, a complex that amalgamates critical elements from various cell death pathways [22, 27,28,29,30]. The PANoptosome—which was initially identified in response to influenza A virus (IAV) infection—comprises key components such as RIPK1, RIPK3, CASP6, and ZBP1 [22, 31, 32]. Notably, PANoptosome facilitates the simultaneous activation of its constituent pathways and promotes pro-inflammatory cell death. However, the specific role of PANoptosome mediated cell death in GCs and its contribution to follicular atresia remains unclear. In light of this, here we provide evidence that the assembly and activation of the PANoptosome are crucial for the PCD of GCs, contributing to follicular atresia.

Materials and methods

Classification of healthy, slightly atretic, and atretic follicles

Ovaries were harvested from 5-month-old, commercially raised Large White gilts at a local abattoir and transported to the laboratory within 2 h in vacuum flasks containing sterile physiological saline at 30–35 °C. Upon arrival, the ovaries were rinsed twice with sterile physiological saline at 37 °C, enriched with 100 IU/L penicillin and 50 mg/L streptomycin.

Gilts alternate between ovaries characterized by small follicles (< 3 mm), known as honeycomb type (HT), and those with medium to large follicles (3–7 mm), referred to as grape type (GT), which are observable on the ovarian surface. Therefore, we focused on GT ovaries during the follicular phase (d 17–21) of the estrous cycle for our study, as follicle selection is believed to originate from the proliferative pool of 1–6 mm follicles present between d 14 and 16 of the estrous cycle. Additionally, during this phase, the atresia rate of small follicles increases, making it easier to assess follicular atresia through morphological observation. The follicles were categorized into healthy (H), slightly atretic (SA), and atretic (A) based on established morphological criteria: Healthy follicles exhibited a vascularized theca interna and clear, amber follicular fluid without debris. In contrast, follicles not meeting these criteria were deemed atretic, with slightly atretic and atretic follicles displaying a gray theca interna and flocculent follicular fluid to varying extents [33, 34].

Cell culture and experimental design

GCs, harvested from the ovarian follicles of gilts with diameters of 3–5 mm, were pooled and washed twice with phosphate-buffered saline (PBS) (Gibco; 1 mmol/L KH2PO4, 155 mmol/L NaCl, 3 mmol/L Na2HPO4·7H2O) containing penicillin–streptomycin. The washed GCs were then evenly distributed cultured in 35 mm dishes using Dulbecco’s Modified Eagle Medium (DMEM) supplemented with 10% fetal bovine serum (FBS) and incubated for 24 h at 37 °C in a 5% CO2 atmosphere.

Pigs are polyovulatory animals with a high number of ovarian follicles, which facilitates the assessment of follicular health and atresia. Consequently, we selected the porcine ovary as the model system for our study (Fig. 1).

Western blotting

Samples were incubated in lysis buffer on ice for 10 min. Next, the lysates were transferred into 1.5-mL centrifuge tubes and centrifuged at 12,000 × g for 10 min at 4 °C. The protein concentration was determined using a BCA assay kit (TransGen Biotech, Beijing, China). The samples were then prepared for SDS-PAGE by mixing with protein loading buffer (TransGen Biotech, Beijing, China) and boiling at 100 °C for 10 min. Following this, the 30 µg protein samples were subjected to gel electrophoresis at 150 V using a 10% Bis-Tris gel (Epizyme, Shanghai, China) for 1 h. Following electrophoresis, the proteins were transferred onto a nitrocellulose membrane (BioTrace™ NT; Pall Corp., FL, USA) using a wet transfer system and incubated with primary antibodies at 4 °C overnight. The primary antibodies (Table S1) used were diluted 1:1,000 with primary antibody diluent (Huaxing Bio, Beijing, China). After three washes with TBST, the membranes were incubated for 1 h at 37 °C with HRP-bound anti-rabbit IgG and HRP-bound anti-mouse IgG (ZSGB-BIO, Beijing, China). Each antibody was diluted to 1:5,000 with a secondary antibody diluent (Huaxing Bio, Beijing, China).

Immunofluorescence

Porcine follicles at various stages measuring 3–5 mm in diameter and at various stages, were preserved in 4% phosphate-buffered formaldehyde before being embedded in paraffin. The follicle sections with a thickness of 5 μm were chosen at random and subjected to staining. Immunofluorescence was performed following previously established protocols [35]. For analysis via cellular immunofluorescence, 293T cells were cultured on glass coverslips and fixed with 4% paraformaldehyde for 1 h, then washed with PBS. The cells were permeabilized using 0.5% Triton X-100 for 10 min at 4 °C and blocked with 1% BSA at room temperature for 1 h. Subsequently, the cells were stained with the appropriate primary antibody at 4 °C overnight. Following 3 washes with PBS, the cells were incubated with a secondary antibody diluted in PBS for 1 h at 37 °C, and the nuclei were counterstained with Hoechst 33342 for 5 min. Finally, immunofluorescence samples were analyzed using a confocal microscope (A1HD25, Nikon, Japan), with the images being duly recorded.

ROS measurement

ROS (reactive oxygen species) levels were measured using the Reactive Oxygen Species Assay Kit (S0033S, Beyotime, Shanghai, China) following the manufacturer’s protocol. This fluorescent method relies on the oxidation of dichlorodihydrofluorescein (DCFH) to dichlorofluorescein (DCF), which emits green fluorescence when excited at 488 nm. The cells were visualized using a confocal microscope (A1HD25, Nikon, Japan).

Plasmid construction and transfection

The coding sequences of the ZBP1, RIPK3, and RIPK1 genes, as well as those of their truncation mutants, were amplified from GC cDNA and cloned into the pcDNA3.1 expression vector (Invitrogen). Accordingly, site-directed mutagenesis was performed as previously described [33]. GFP-Casp6 variants S76A, S76D, S254A, S254D, S76AS254A, and S76DS254D were generated using the Fast Mutagenesis System (FM111; TransGen Biotech). In porcine caspase6, Ser76 and Ser254 were mutated to alanine (caspase6-S76A, caspase6-S254A) and aspartic acid (caspase6-S76D, caspase6-S254D), respectively. 293T cells for transfection were cultured in 60 mm dishes using Dulbecco’s DMEM supplemented with 10% FBS and incubated at 37 °C with 5% CO2 until reaching the desired confluency for transfection. 293T cells were transfected with these plasmids using Lipo8000™ (Beyotime, Shanghai, China) in Opti-MEM.

Molecular docking analysis

To explore intermolecular relationships, we performed rigid protein–protein docking analysis using GRAMM-X (http://gramm.compbio.ku.edu/). The structural domains of RIPK3, RIPK1, ZBP1, and CASP6 were obtained from the AlphaFold Protein Structure Database (https://alphafold.ebi.ac.uk/). We utilized PyMOL (Version 2.4) and PDBePISA (https://www.ebi.ac.uk/pdbe/pisa/) to study protein–protein interactions and conduct visual examinations.

Co-immunoprecipitation assay

293T cells were lysed in immunoprecipitation buffer (Beyotime, Shanghai, China). The lysates were incubated with 2 μg of specific primary antibodies and gently rotated at 4 °C overnight and then combined with Protein A/G magnetic beads and rotated for 2 h. After 3 washes with IP buffer, the immunoprecipitated samples were eluted and boiled for 10 min in 1% (w/v) SDS sample buffer. Finally, the immunoprecipitated samples were analyzed using Western blot.

Gel filtration chromatography

GC proteins from porcine atretic follicles were extracted at 4 °C using lysis buffer (20 mmol/L Tris, 150 mmol/L NaCl, 1% Triton X-100), with protease inhibitors added at a ratio of 1/100 of the final volume. Following two rounds of centrifugation at 13,000 × g at 4 °C, protein concentration was measured using a BCA protein assay kit. Gel filtration was carried out at 4 °C on a Superose 6 Increase 10/300 GL chromatography column (Cytiva, USA), equilibrated with PBS buffer (10 mmol/L phosphate buffer, 140 mmol/L NaCl, pH 7.4) at a flow rate of 0.5 mL/min. Fractions of 0.5 mL were collected, and 30 μL from each fraction was analyzed by Western blotting.

Heat stress exposure

GCs were cultured in an incubator (Thermo Fisher Scientific) at 43 °C with humidified 5% CO2 for 4 h, followed by incubation at 37 °C for the specified durations. HSP70 expression increases in response to stress, serving as a reliable biomarker for heat stress [36, 37].

Data availability

The proteome data for porcine GCs from healthy, slightly atretic, and atretic follicles were obtained from DOI: https://doi.org/10.3389/fcell.2020.624985.

Statistical analysis

All statistical analyses were performed using SPSS 16.0 software (IBM, USA). Data are presented as mean ± SD. Differences between groups were assessed using ANOVA, followed by an LSD post hoc test. P < 0.05 was considered statistically significant.

Fig. 1
figure 1

The workflow schematic of this study. Firstly, the levels of programmed cell death marker proteins in the granulosa cells from healthy (H), slightly atretic (SA), and atretic (A) follicles were compared. Next, the interactions of PANoptosome molecules and their interaction domains will be directly studied through vector construction and co-immunoprecipitation techniques. Finally, a heat stress model will be used to study the process of programmed cell death in granulosa cells

Results

PCD is activated in GCs during porcine follicular atresia

To explore the involvement of PCD in porcine GCs during follicular atresia, we analyzed differential gene expression using an existing proteome dataset from GCs comprising healthy, slightly atresic, and atretic follicles. Figure 2A illustrates a heatmap that highlights the specific gene regulation within the PCD pathway, showing changes in protein levels of key genes in GCs from atretic follicles. Further investigation through Western blot analysis identified the presence of PCD markers in GCs from healthy, slightly atretic, and atretic follicle (Fig. 2B). Compared to healthy cells, those from atretic follicles demonstrated notable increases in the expression of necroptosis markers (RIPK1, RIPK3, and MLKL) and enhanced phosphorylation (P < 0.01, Fig. 2C). Additionally, the expression levels of pyroptosis markers (NLRP3, CASP1, and GSDMD) were increased, with CASP1 and GSDMD showing signs of cleavage in the atretic categories (P < 0.01, Fig. 2C). Immunofluorescence imaging (Fig. 2D) revealed the cytoplasmic localization and expression patterns of GSDMD and MLKL in GCs, where the intensity of the staining was correlated with the progression of follicular atresia. These findings indicate that GCs may undergo both pyroptosis and necroptosis during follicular atresia.

Fig. 2
figure 2

PCD is activated in GCs during porcine follicular atresia. A Heatmap of differentially regulated transcripts enriched in PCD-related genes between healthy and atretic follicular GCs. B Immunoblots of GC lysates from healthy (H), slight atretic (S), and atretic (A) in porcine follicles. C Quantitative analysis of protein levels in B. D Immunofluorescence staining of GSDMD and MLKL in H, S and A porcine follicles. The scale bar represents 100 μm. Nuclei were stained with DAPI (blue). The results are presented as mean ± SD. *P < 0.05; **P < 0.01; n.s., not significant

PANoptosome molecules directly interact

We conducted gel filtration experiments to confirm that the core PANoptosome proteins, RIPK3, RIPK1, and ZBP1, form a complex (Fig. 3A). To evaluate the interaction between the components of the cell death complex directly, we constructed plasmids encoding porcine RIPK3, RIPK1, and ZBP1 for transfection into 293T cells, as the endogenous proteins in GCs cannot be pulled down with currently available antibodies. Subsequent experiments demonstrated reciprocal co-immunoprecipitation between RIPK3 and both RIPK1 and ZBP1, as illustrated in Fig. 3B and C. Additionally, RIPK1 and ZBP1 exhibited mutual co-immunoprecipitation (Fig. 3D). Immunofluorescence analysis revealed co-localization of RIPK3, RIPK1, and ZBP1 (Fig. 3E–G). Additionally, bioinformatic analysis demonstrated the formation of hydrogen bonds between specific amino acid residues of these proteins, suggesting a stable docking model (Fig. 3H–I). Collectively, the results support the hypothesis that these porcine proteins interact directly to form a PANoptosome complex.

Fig. 3
figure 3

PANoptosome molecules directly interact. A Gel filtration analysis of the PANoptosome complex. The Western blots of column fractions were performed using the indicated antibodies. B HA-RIPK1 and Myc-RIPK3 were co-transfected into 293T cells for 48 h for immunoprecipitation. C Flag-ZBP1 and Myc-RIPK3 were co-transfected into 293T cells for 48 h for immunoprecipitation. D Flag-ZBP1 and HA-RIPK1 were co-transfected into 293T cells for 48 h for immunoprecipitation. E Co-transfection of HA-RIPK1 and Myc-RIPK3 in 293T cells, showing DNA stained with DAPI. Different tag proteins were stained with mouse- and rabbit-derived primary and secondary antibodies. Scale bar: 5 = μm. F Co-transfection of Myc-RIPK3 and Flag-ZBP1 in 293T cells, showing DNA stained with DAPI. Different tag proteins were stained with mouse- and rabbit-derived primary and secondary antibodies. Scale bar: 5 = μm. G Co-transfection of HA-RIPK1 and Flag-ZBP1 in 293T cells, showing DNA stained with DAPI. Different tag proteins were stained with mouse- and rabbit-derived primary and secondary antibodies. Scale bar: 5 = μm. H Interfacing residues between porcine RIPK3 and RIPK1 proteins (RIPK3, red; RIPK1, yellow; hydrogen bond interaction indicated by dotted lines). I Interfacing residues between porcine RIPK3 and ZBP1 proteins (RIPK3, red; ZBP1, purple; hydrogen bond interaction are indicated by dotted lines). J Interfacing residues between porcine RIPK1 and ZBP1 proteins (RIPK1, yellow; ZBP1, purple; hydrogen bond interactions are indicated by dotted lines)

RIPK3, RIPK1, and ZBP1 interact via the RHIM domain

We examined the interaction domains of RIPK3, RIPK1, and ZBP1. A comparative analysis of RIPK3, RIPK1, and ZBP1 amino acid sequences across porcine, mouse, and human species revealed the presence of a (V/I)-Q-(V/I/L/C)-G motif—characteristic of the RIP homotypic interaction motif (RHIM)—at the C-termini of all three proteins (Fig. 4A). To further detemine whether the interactions between these proteins depended on the RHIM domain, truncated mutants of RIPK3, RIPK1, and ZBP1 were generated in 293T cells to map the essential interaction domains (Fig. 4B). Both the amino-terminal and carboxy-terminal domains of RIPK1 interacted with the corresponding domains of RIPK3 (Fig. 4C). The carboxyl termini of these two proteins exhibited strong binding (Fig. 4C, indicated by the red arrows). Interestingly, deficiency in RHIM prevented ZBP1–RIPK3 binding as well as RIPK1 binding (Fig. 4D), implying the crucial role of the RHIM domain in mediating protein interactions.

Fig. 4
figure 4

RIPK3, RIPK1, and ZBP1 interact via the RHIM domain. A Amino acid sequences of RHIM-containing regions of human, pig, and mouse proteins RIPK1, RIPK3, and ZBP1. B Schematic representation of domains in porcine CASP6, RIPK1, and RIPK3 in full-length (F), N-terminal (N), and C-terminal (C) constructs. C Immunoprecipitates and total lysates were obtained from 293T cells following co-transfection with the specified variants of RIPK3 and RIPK1 for a duration of 48 h. D Immunoprecipitates and total lysates were obtained from 293T cells following co-transfection with the specified variants of RIPK3 and ZBP1 for a duration of 48 h. E Immunoprecipitates and total lysates were obtained from 293T cells following co-transfection with the specified variants of RIPK1 and ZBP1 for a duration of 48 h. F Diagram of the interaction mode of RIPK3, RIPK1, and ZBP1

Caspase-6 is a component of the PANoptosome in GCs

To explore the mechanistic role of caspase-6 in GCs during the process of follicular atresia, we initially assessed the variations in caspase-6 expression levels among GCs from follicles in different states. Caspase-6 undergoes cleavage in the GCs of atretic follicles, accompanied by a decrease in phosphorylation trends at its Ser76 and Ser254 sites (P < 0.01, Fig. 5A). Both Ser254 and Ser76, along with their respective adjacent sequences in caspase-6, demonstrated conservation across diverse species (Fig. 5B). Moreover, wild-type caspase-6, the non-phospho-mimetic Ser76, S254A mutant, and the phospho-mimetic S76D, S254D mutants were overexpressed in 293T cells. Following this, the cells underwent treatment with low doses of TNFα and CHX to induce the cleavage of caspase-6. The S254A, S76A, and S254AS76A mutants exhibited increased sensitivity to cleavage, whereas the S254D, S76D, and S76DS254D mutants displayed complete resistance (Fig. 5D).

Fig. 5
figure 5

Caspase-6 as a component of the PANoptosome in GCs. A Immunoblot of GC lysates from healthy (H), slight atretic (S), and artretic (A) follicles. B A quantitative analysis of the protein levels shown in A. C Alignment of the amino acid sequences surrounding Ser257 and Ser76 of caspase-6 in humans, mice, and pigs. D 293T cells overexpressing GFP-caspase-6 wild type, S76A, S76D, S254A, S254D, S76AS254A, or S76DS254D mutants were treated with 10 μg/ml CHX and 25 ng/ml TNFα for 2 h. The immunoblot analysis of cell lysates is shown. (Single-letter abbreviations for amino acid residues are as follows: A for Ala, D for Asp, and S for Ser. In the case of the mutants, specific amino acid substitutions occurred at defined positions. For examples, S76A signifies the replacement of serine at position 76 with alanine.) E Surface diagram of the docking model and the interfacing residues between CASP6 and RIPK3 proteins (CASP6, green; RIPK3, red; hydrogen bond interactions are indicated by dotted lines). F Co-transfection of GFP-CASP6 and Myc-RIPK3 in 293T cells is shown. DNA was stained with DAPI. Scale bar, 5 μm. G Co-immunoprecipitation assay of the interactions between RIPK3 and caspase-6 or its mutants at site 254. H Co-immunoprecipitation assay to assess the interactions between RIPK3 and caspase-6 or its mutants at site 76. I Caspase-6 as a schematic diagram of the PANoptosome structure. The results are presented as the means ± SD. *P < 0.05; **P < 0.01; n.s., not significant

To investigate whether caspase-6 contributes to the formation of the PANoptosome complex, we performed a protein–protein docking analysis between caspase-6 and RIPK3. The analysis revealed the formation of hydrogen bonds at specific amino acid residue sites, which indicated the establishment of a robust and stable protein docking model between caspase-6 and RIPK3 (Fig. 5E). Subsequently, we co-overexpressed caspase-6 along with RIPK1, ZBP1, and RIPK3 in 293T cells, enabling us to evaluate the interaction between caspase-6 and the components of the complex. Additionally, immunoprecipitation assays demonstrated the co-immunoprecipitation of caspase-6 and its mutant with RIPK3. Notably, both the full-length and cleaved forms of caspase-6 could interact with RIPK3 (Fig. 5G and H). Immunofluorescence co-localization signals further supported the interaction between caspase-6 and RIPK3 (Fig. 5F). However, we did not observe the binding of caspase-6 to either ZBP1 or RIPK1 (Fig. S1), suggesting a potential recruitment mechanism of caspase-6 to the complex via RIPK3 (Fig. 5I and S2).

Heat stress triggers GC pyroptosis and necroptosis

To investigate the effect of heat stress on GC PCD, cultured GCs were exposed to 43 °C for 4 h. This exposure significantly increased the phosphorylation levels of RIPK3, RIPK1, and MLKL proteins, as well as the cleavage of GSDMD, within 12 h post-heat stress. (P < 0.01, Fig. 6A, C–K). Furthermore, these cells demonstrated an elevated accumulation of intracellular reactive oxygen species (ROS) compared to the controls following heat stress (Fig. 6B).

Fig. 6
figure 6

Heat stress triggers granulosa cell pyroptosis and necroptosis. A Western blot analysis of the quantity of indicated proteins in GCs at the indicated time points after heat stress (43 °C for 4 h). B Fluorescent images showing ROS staining in control and heat stress-exposed GCs. Scale bar: 10 μm. C–K Quantitative analysis of protein levels in A. The results are presented as mean ± SD. *P < 0.05; **P < 0.01; n.s., not significant

Discussion

Follicular development and ovulation are key determinants of the reproductive potential of sows. However, the onset of follicular atresia and degeneration, triggered by the death of GCs, is a limiting factor for sow fertility. These intricate processes of follicular atresia and GC death are tightly regulated by a diverse array of factors. In the present study, we discovered that GCs undergo not only undergo apoptosis but also pyroptosis and necroptosis during the process of follicular atresia. Furthermore, we identified the assembly and activation of a PANoptosome during follicular atresia that mediates the occurrence of these three types of cell death (Fig. 7). Our results provide new insights into the mechanisms of follicular atresia and GC PCD.

Fig. 7
figure 7

The mechanism of PANoptosome formation and regulation of granulosa cell programmed cell death during follicular atresia

Apoptosis is the primary mechanism underlying GC death processes during follicle atresia [3]. Beyond apoptosis, necroptosis has been observed in GCs within the ovary as a response to the accumulation of reactive oxygen species (ROS). In vitro studies have shown that subjecting human GCs to serum starvation leads to the generation of ROS, triggering both necroptosis and apoptosis [38]. Notably, mediators of necroptosis, such as phosphorylated MLKL, RIPK1, and RIPK3 proteins were identifiable within the GCs of both pre-antral and antral follicles in macaque ovaries, as well as in the corpora lutea of human ovaries [39, 40]. Furthermore, recent research has demonstrated the upregulation of gene expression of RIPK1 and RIPK3 in atretic bovine follicles, in contrast to healthy follicles [18]. Notably, oxidative stress can prime the NLRP3 inflammasome, leading to caspase-1 activation and pyroptosis [41, 42]. Meanwhile, a recent study involving dairy cows revealed that culturing GCs with non-esterified fatty acids induced oxidative stress, pyroptosis, and inflammation [19]. Consistent with these results, we observed an upregulation of genes related to pyroptosis and necroptosis, including ZBP1, RIPK3, RIPK1, and MLKL, as well as their phosphorylated forms. Additionally, cleavage of GSDMD and caspase-1 was occurred, indicating the occurrence of pyroptosis and necroptosis in GCs during follicle atresia. This suggests a potential contribution of GC PCD to follicular atresia.

Interestingly, it has been found that the above proteins can act as integral components of larger cell death-inducing complexes known as PANoptosomes, which integrate components from multiple cell death pathways and drive a form of lytic, innate immune inflammatory cell death termed PANoptosis [21, 43, 44]. Previous studies have suggested that a singular cell death complex may regulate the activation of these pathways by orchestrating interactions among their key components [30, 45]. Accordingly, to verify PANoptosome formation in GCs during follicular atresia, we generated porcine RIPK3, RIPK1, and ZBP1 plasmids and conducted immunoprecipitation experiments. Our findings demonstrated interactions among RIPK1, RIPK3, and ZBP1, suggesting the assemble of a PANoptosome within the GCs. This aligns with previous observations demonstrating that macrophages assemble and activate the PANoptosome to induce PANoptosis in response to invasion by pathogenic microorganisms [30]. Nevertheless, the composition of PANoptosomes has varied among studies, and the phenotypes appear to depend on the stimulus provided [45].

Homotypic interactions are essential in various signaling complexes. While the compositions of individual PANoptosome may differ, each contains components with key domains such as RHIMs that are crucial for assembly [29]. The RHIM domains in these proteins include a core tetrad characterized by the consensus sequence (V/I)-Q-(V/I/L/C)-G that is essential for intermolecular protein interactions and signal transduction [22, 43, 46]. Specifically, RIPK3-dependent signal transduction relies heavily on the precise assembly of RHIM complexes. Thus, any substitution mutants in the RHIM domain residues can significantly disrupt downstream signaling events [47]. For example, altering a single valine residue within the RIPK3 core tetrad can prevent necroptosis in vitro [48]. In mouse bone marrow-derived macrophages, a lack of RHIM impeded the binding between ZBP1 and RIPK1 [49]. Consistent with these findings, our study demonstrated that ZBP1 engages RIPK3 and RIPK1 through the receptor interacting protein RHIM, and in the absence of the RHIM structure of ZBP1 cannot interact with RIPK3 or RIPK1, suggesting that RHIM is essential for the assembly of a PANoptosome during porcine follicle atresia.

Caspase-6 has long been regarded as a mysterious protein. In colorectal adenocarcinoma cells, phosphorylation of the Ser254 of caspase-6 by AMPK related protein kinase 5 has been identified [50]. This phosphorylation event has been shown to suppress both the cleavage and activation of caspase-6 [51]. Our previous proteomic phosphorylation studies have shown that the Ser76 site is significantly phosphorylated in GCs of healthy follicles, with dephosphorylation occurring during the atresia process. To validate this phosphorylation event, we generated a polyclonal antibody against the phospho-S76 peptide. Consistent with previous studies, we observed that the phosphorylation of caspase-6 at the Ser76 and Ser254 site in GCs gradually decreases, thereby leading to its cleavage during follicular atresia (Fig. 4A). Ser76 serves as another active site of caspase-6. Dephosphorylation at this site facilitates the cleavage and activation of caspase-6, ultimately leading to apoptosis [33]. Despite being classified as an apoptotic executor caspase and being implicated in various diseases such as nonalcoholic steatohepatitis, the full extent of caspase’s functional significance and its involvement in other biological processes remain unclear [52]. Recent studies have suggested that caspase-6 plays a key role in innate immunity and the activation of the ZBP1-NLRP3 inflammasome [22, 53, 54]. Notably, observations suggest a scaffolding function for caspase-6, as evidenced by the enhanced interaction between RIPK3 and ZBP1 in the presence of catalytically inactive or uncleavable caspase-6 mutants, particularly during infections with pathogens such as IAV [22]. In the present study, we observed that caspase-6 was a component of the PANoptosome. Interestingly, both the full-length and cleaved forms of caspase-6 interacted with RIPK3 but showed no interaction with RIPK1 or ZBP1. This observation suggests that the function of caspase-6 depends on the RIPK3 pathway.

The activation of the PANoptosome integrates multiple programmed cell death pathways, providing GCs with a flexible mechanism to adapt to stress, efficiently clear damaged cells, and protect the follicular microenvironment. Pyroptosis and necroptosis promote inflammatory responses and coordinate cell death with immune responses, making the PANoptosome more effective at eliminating GCs during follicular atresia compared to a single cell death pathway, thereby ensuring proper ovarian function. Different from apoptotic bodies formed by apoptosis, the most prominent feature of pyrodeath and necroptosis is the formation of pores in the cell membrane [26, 55]. This process results in the release of pro-inflammatory cytokines, particularly IL-1β and IL-18. The release of these cytokines, along with other intracellular contents like DAMPs, into the extracellular space triggers a strong inflammatory response. This inflammation recruits immune cells to the site of infection or damage, amplifies the immune response, and promotes the clearance of pathogens or damaged cells. It also provides evidence for the increase in macrophages in follicles during follicular atresia [56,57,58]. Therefore, further experiments are needed to prove the role of pyroptosis and necroptosis in follicular atresia.

The activation mechanisms of the PANoptosome in follicular atresia of GCs are currently unclear. In previous studies, the activation and assembly of the PANoptosome depended on various pathogen-induced innate immune responses [21, 22, 27]. Recent research has demonstrated that heat stress can induce PANoptosis [59]. Furthermore, under heat stress, the generated ROS can also induce necroptosis and pyroptosis [59,60,61,62,63]. A recent study has shown that heat stress can trigger necroptosis, apoptosis, and pyroptosis via the ZBP1-RIPK3 signaling pathway. In addition, the genetic deletion or targeted inhibition of ZBP1, RIPK3, or both caspase-8 and MLKL can prevent these pathways of cell death, thereby improving cell survival [59]. GCs exposed to heat stress in vitro demonstrated impaired physiological functions, including increasing intracellular accumulation of ROS-inducing apoptosis, and reduced synthesis of estradiol and progesterone [36, 37]. The above results indicate that cells may undergo apoptosis, necroptosis, and pyroptosis under heat stress and oxidative stress. Our previous studies have shown that heat stress induces apoptosis in porcine GCs [13]. In the present study, we found that heat stress could activate the PANoptosome, promoting necroptosis and pyroptosis in porcine GCs. Additionally, intracellular ROS levels also increased, which suggested a mechanistic connection between ROS production and PANoptosis. These results provide new insights into how heat waves during the summer cause follicular atresia, leading to reduced reproductive capacity in livestock [36, 37, 64].

As a critical determinant of follicular fate, FSH promotes the proliferation and differentiation of GCs. However, whether it also inhibits the activation of the PANoptosome in GCs as a protective mechanism for follicles requires further investigation. The ovarian follicular vascular system plays a vital role in supplying the essential nutrients, oxygen, and hormones necessary for follicle development [65]. Angiogenesis was active and increasing in selected follicles, while those undergoing atresia exhibited minimal blood vessel development [10, 66]. In atretic follicles, the reduction in blood vessel density within the basal membrane creates a hypoxic environment in GCs, leading to apoptosis [67]. Recent studies have shown that hypoxic conditions can simultaneously induce pyroptosis and necroptosis in mammalian cells [68,69,70]. The hypoxic environment within the follicle is a potential detrimental factor that triggers the PANoptosome in GCs, leading to PANoptosis. Furthermore, ovarian aging is often accompanied by an increase in follicular atresia and a decrease in follicle count [71]. In addition to apoptosis, pyroptosis and necrosis also play a role in the process of ovarian aging [72, 73]. Recent studies have shown that nicotinamide mononucleotide and spermidine can improve the quality of mouse oocytes, thereby delaying aging [74, 75]. However, these studies did not investigate the mechanisms underlying their effects on GCs. It is possible that they inhibit PANoptosis in GCs, thereby enhancing the communication between GCs and oocytes. This will be a key focus of future research.

Conclusions

Our study revealed that GCs undergo PANoptosis during follicular atresia. By demonstrating the interaction among RIPK1, RIPK3, ZBP1, and CASP6, we provide evidence for the existence of a cell death complex, the PANoptosome, which orchestrates various pathways in GCs. This PANoptosome likely plays a crucial role in regulating follicular atresia (Fig. 7). As such, understanding the molecular events underlying follicular atresia is crucial for elucidating the physiological processes involved in ovarian follicle development and dynamics.

Data Availability

The proteome data of porcine GCs from healthy, slightly atresia and atretic follicles was downloaded from doi: https://doi.org/10.3389/fcell.2020.624985 .

Abbreviations

  • A:: Atretic
  • FSH:: Follicle stimulating hormone
  • GC:: Granulosa cell
  • GSDMD:: Gasdermin D
  • H:: Healthy
  • IAV:: Influenza A virus
  • MLKL:: Mixed lineage kinase domain-like protein
  • PBS:: Phosphate buffered saline
  • PCD:: Programmed cell death
  • RHIM:: RIP homotypic interaction motif
  • ROS:: Reactive oxygen species
  • SA:: Slightly atretic

References

  1. 1.Pelosi E, Forabosco A, Schlessinger D. Genetics of the ovarian reserve. Front Genet. 2015;6.(2015)308. https://doi. org/10.3389/fgene.: 308.
  2. 2.Wallace WH, Kelsey TW. Human ovarian reserve from conception to the menopause. PLoS One. 2010;5(1).(2010)pone.0008772.
  3. 3.Regan SLP, Knight PG, Yovich JL, Leung Y, Arfuso F, Dharmarajan A. Granulosa cell apoptosis in the ovarian follicle—a changing view. Front Endocrinol (Lausanne). 2018;9.(2018)61. https://doi. org/10.3389/fendo.: 61.
  4. 4.Devorshak-Harvey E, Peluso JJ, Bona-Gallo A, Gallo RV. Effect of alterations in pulsatile luteinizing hormone release on ovarian follicular atresia and steroid secretion on diestrus 1 in the rat estrous cycle. Biol Reprod. 1985;33(1).(1985)1.103.: 103.
  5. 5.Wu G, Li C, Tao J, Liu Z, Li X, Zang Z, et al. FSH mediates estradiol synthesis in hypoxic granulosa cells by activating glycolytic metabolism through the HIF-1α-AMPK-GLUT1 signaling pathway. J Biol Chem. 2022;298(5).(2022)101830. https://doi. org/10. 1016/j.jbc.: 101830.
  6. 6.Chu YL, Xu YR, Yang WX, Sun Y. The role of FSH and TGF-β superfamily in follicle atresia. Aging (Albany NY). 2018;10(3).(2018)18632/aging.101391.: 305.
  7. 7.Liu S, Jia Y, Meng S, Luo Y, Yang Q, Pan Z. Mechanisms of and potential medications for oxidative stress in ovarian granulosa cells.(2023)A review.Int J Mol Sci.: 9205.
  8. 8.Shen M, Jiang Y, Guan Z, Cao Y, Sun SC, Liu H. FSH protects mouse granulosa cells from oxidative damage by repressing mitophagy. Sci Rep. 2016;6.(2016)org/10.1038/srep38090.: 38090.
  9. 9.Li C, Liu Z, Wu G, Zang Z, Zhang JQ, Li X, et al. FOXO1 mediates hypoxia-induced G0/G1 arrest in ovarian somatic granulosa cells by activating the TP53INP1-p53-CDKN1A pathway. Development. 2021;148(14).(2021)1242/dev.199453.
  10. 10.Li SY, Bhandary B, Gu X, DeFalco T. Perivascular cells support folliculogenesis in the developing ovary. Proc Natl Acad Sci U S A. 2022;119(41).(2022)1073/pnas.2213026119.
  11. 11.Fraser HM. Regulation of the ovarian follicular vasculature. Reprod Biol Endocrinol. 2006;4.(2006)org/10.1186/1477-7827-4-18.: 18.
  12. 12.Li M. The role of P53 up-regulated modulator of apoptosis (PUMA) in ovarian development, cardiovascular and neurodegenerative diseases. Apoptosis. 2021;26(5–6).(2021)org/10.1007/s10495-021-01667-z.: 235.
  13. 13.Yang F, Chen Y, Liu Q, Dai S, Zeng S. Dynamics and regulations of BimEL Ser65 and Thr112 phosphorylation in porcine granulosa cells during follicular atresia. Cells. 2020;9(2).(2020)org/10.3390/cells9020402.: 402.
  14. 14.Zhou X, He Y, Li N, Bai G, Pan X, Zhang Z, et al. DNA methylation mediated RSPO2 to promote follicular development in mammals. Cell Death Dis. 2021;12(7).(2021)org/10.1038/s41419-021-03941-z.: 653.
  15. 15.Manabe N, Goto Y, Matsuda-Minehata F, Inoue N, Maeda A, Sakamaki K, et al. Regulation mechanism of selective atresia in porcine follicles.(2004)regulation of granulosa cell apoptosis during atresia.J Reprod Dev.: 493.
  16. 16.Li Y, Zhang Z, Wang S, Du X, Li Q. miR-423 sponged by lncRNA NORHA inhibits granulosa cell apoptosis. J Anim Sci Biotechnol. 2023;14.(2023)org/10.1186/s40104-023-00960-y.: 154.
  17. 17.Sapuleni J, Szymanska M, Meidan R. Diverse actions of sirtuin-1 on ovulatory genes and cell death pathways in human granulosa cells. Reprod Biol Endocrinol. 2022;20.(2022)org/10.1186/s12958-022-00970-x.: 104.
  18. 18.McEvoy MJ, Sinderewicz E, Creedon L, McAfee M, Jonczyk AW, Piotrowska-Tomala KK, et al. Death processes in Bovine theca and granulosa cells modelled and analysed using a systems biology approach. Int J Mol Sci. 2021;22(9).(2021)org/10.3390/ijms22094888.: 4888.
  19. 19.Wang Y, Li C, Ali I, Li L, Wang G. N-acetylcysteine modulates non-esterified fatty acid-induced pyroptosis and inflammation in granulosa cells. Mol Immunol. 2020;127.(2020)157–63. https://doi. org/10. 1016/j.molimm.: 157.
  20. 20.Hou J, Lei Z, Cui L, Hou Y, Yang L, An R, et al. Polystyrene microplastics lead to pyroptosis and apoptosis of ovarian granulosa cells via NLRP3/Caspase-1 signaling pathway in rats. Ecotoxicol Environ Saf. 2021;212.(2021)112012. https://doi. org/10. 1016/j.ecoenv.: 112012.
  21. 21.Lee S, Karki R, Wang Y, Nguyen LN, Kalathur RC, Kanneganti TD. AIM2 forms a complex with pyrin and ZBP1 to drive PANoptosis and host defence. Nature. 2021;597(7876).(2021)org/10.1038/s41586-021-03875-8.: 415.
  22. 22.Zheng M, Karki R, Vogel P, Kanneganti TD. Caspase-6 Is a key regulator of innate immunity, inflammasome activation, and host defense. Cell. 2020;181(3).(2020)674–87. https://doi. org/10. 1016/j.cell.: 674.
  23. 23.Frank D, Vince JE. Pyroptosis versus necroptosis.(2019)similarities, differences, and crosstalk.Cell Death Differ.: 99.
  24. 24.Khan N, Lawlor KE, Murphy JM, Vince JE. More to life than death.(2014)molecular determinants of necroptotic and non-necroptotic RIP3 kinase signaling.Curr Opin Immunol.: 76.
  25. 25.Kayagaki N, Stowe IB, Lee BL, O’Rourke K, Anderson K, Warming S, et al. Caspase-11 cleaves gasdermin D for non-canonical inflammasome signalling. Nature. 2015;526(7575).(2015)org/10.1038/nature15541.: 666.
  26. 26.Shi J, Zhao Y, Wang K, Shi X, Wang Y, Huang H, et al. Cleavage of GSDMD by inflammatory caspases determines pyroptotic cell death. Nature. 2015;526(7575).(2015)org/10.1038/nature15514.: 660.
  27. 27.Malireddi RKS, Gurung P, Mavuluri J, Dasari TK, Klco JM, Chi H, et al. TAK1 restricts spontaneous NLRP3 activation and cell death to control myeloid proliferation. J Exp Med. 2018;215(4).(2018)1084/jem.20171922.: 1023.
  28. 28.Malireddi RKS, Gurung P, Kesavardhana S, Samir P, Burton A, Mummareddy H, et al. Innate immune priming in the absence of TAK1 drives RIPK1 kinase activity-independent pyroptosis, apoptosis, necroptosis, and inflammatory disease. J Exp Med. 2020;217(3).(2020)1084/jem.20191644.
  29. 29.Samir P, Malireddi RKS, Kanneganti TD. The PANoptosome.(2020)a deadly protein complex driving pyroptosis, apoptosis, and necroptosis (PANoptosis).Front Cell Infect Microbiol.: 238.
  30. 30.Christgen S, Zheng M, Kesavardhana S, Karki R, Malireddi RKS, Banoth B, et al. Identification of the PANoptosome.(2020)a molecular platform triggering pyroptosis, apoptosis, and necroptosis (PANoptosis).Front Cell Infect Microbiol.: 237.
  31. 31.Zhang T, Yin C, Boyd DF, Quarato G, Ingram JP, Shubina M, et al. Influenza virus Z-RNAs induce ZBP1-mediated necroptosis. Cell. 2020;180(6).(2020)1115–29. https://doi. org/10. 1016/j.cell.: 1115.
  32. 32.Kuriakose T, Man SM, Malireddi RK, Karki R, Kesavardhana S, Place DE, et al. ZBP1/DAI is an innate sensor of influenza virus triggering the NLRP3 inflammasome and programmed cell death pathways. Sci Immunol. 2016;1(2).(2016)1126/sciimmunol.aag2045.
  33. 33.Yang F, Liu Q, Chen Y, Ye H, Wang H, Zeng S. Integrative proteomic and phosphoproteomic analyses of granulosa cells during follicular atresia in porcine. Front Cell Dev Biol. 2020;8.(2020)624985. https://doi. org/10.3389/fcell.: 624985.
  34. 34.Wang XL, Wu Y, Tan LB, Tian Z, Liu JH, Zhu DS, et al. Follicle-stimulating hormone regulates pro-apoptotic protein Bcl-2-interacting mediator of cell death-extra long (BimEL)-induced porcine granulosa cell apoptosis. J Biol Chem. 2012;287(13).(2012)M111.293274.: 10166.
  35. 35.Wu Y, Wang Y, Liu Q, Zhu LJ, Gao H, Cui M, et al. Conserved microRNA mediates heating tolerance in germ cells versus surrounding somatic cells. RNA Biol. 2019;16(10).(2019)1494–503. https://doi. org/10.1080/15476286.: 1494.
  36. 36.Khan A, Dou J, Wang Y, Jiang X, Khan MZ, Luo H, et al. Evaluation of heat stress effects on cellular and transcriptional adaptation of bovine granulosa cells. J Anim Sci Biotechnol. 2020;11.(2020)org/10.1186/s40104-019-0408-8.: 25.
  37. 37.Li J, Gao H, Tian Z, Wu Y, Wang Y, Fang Y, et al. Effects of chronic heat stress on granulosa cell apoptosis and follicular atresia in mouse ovary. J Anim Sci Biotechnol. 2016;7.(2016)org/10.1186/s40104-016-0116-6.: 57.
  38. 38.Tsui KH, Wang PH, Lin LT, Li CJ. DHEA protects mitochondria against dual modes of apoptosis and necroptosis in human granulosa HO23 cells. Reproduction. 2017;154(2).(2017)org/10.1530/rep-17-0016.: 101.
  39. 39.Blohberger J, Kunz L, Einwang D, Berg U, Berg D, Ojeda SR, et al. Readthrough acetylcholinesterase (AChE-R) and regulated necrosis.(2015)pharmacological targets for the regulation of ovarian functions? Cell Death Dis.
  40. 40.Du Y, Bagnjuk K, Lawson MS, Xu J, Mayerhofer A. Acetylcholine and necroptosis are players in follicular development in primates. Sci Rep. 2018;8.(2018)org/10.1038/s41598-018-24661-z.: 6166.
  41. 41.Bauernfeind F, Bartok E, Rieger A, Franchi L, Núñez G, Hornung V. Cutting edge.(2011)reactive oxygen species inhibitors block priming, but not activation, of the NLRP3 inflammasome.J Immunol.: 613.
  42. 42.Won JH, Park S, Hong S, Son S, Yu JW. Rotenone-induced impairment of mitochondrial electron transport chain confers a selective priming signal for NLRP3 inflammasome activation. J Biol Chem. 2015;290(45).(2015)M115.667063.: 27425.
  43. 43.Sun X, Yin J, Starovasnik MA, Fairbrother WJ, Dixit VM. Identification of a novel homotypic interaction motif required for the phosphorylation of receptor-interacting protein (RIP) by RIP3. J Biol Chem. 2002;277(11).(2002)1074/jbc.M109488200.: 9505.
  44. 44.Malireddi RKS, Kesavardhana S, Karki R, Kancharana B, Burton AR, Kanneganti TD. RIPK1 distinctly regulates Yersinia-induced inflammatory cell death. PANoptosis Immunohorizons. 2020;4(12).(2020)4049/immunohorizons.2000097.: 789.
  45. 45.Wang L, Zhu Y, Zhang L, Guo L, Wang X, Pan Z, et al. Mechanisms of PANoptosis and relevant small-molecule compounds for fighting diseases. Cell Death Dis. 2023;14(12).(2023)org/10.1038/s41419-023-06370-2.: 851.
  46. 46.Rebsamen M, Heinz LX, Meylan E, Michallet MC, Schroder K, Hofmann K, et al. DAI/ZBP1 recruits RIP1 and RIP3 through RIP homotypic interaction motifs to activate NF-kappaB. EMBO Rep. 2009;10(8).(2009)916–22. https://doi. org/10.1038/embor.: 916.
  47. 47.Hu H, Wu X, Wu G, Nan N, Zhang J, Zhu X, et al. RIP3-mediated necroptosis is regulated by inter-filament assembly of RIP homotypic interaction motif. Cell Death Differ. 2021;28(1).(2021)org/10.1038/s41418-020-0598-9.: 251.
  48. 48.Zhang H, Wu X, Li X, Li M, Li F, Wang L, et al. Crucial roles of the RIP homotypic interaction motifs of RIPK3 in RIPK1-dependent cell death and lymphoproliferative disease. Cell Rep. 2020;31(7).(2020)107650. https://doi. org/10. 1016/j.celrep.: 107650.
  49. 49.Muendlein HI, Connolly WM, Magri Z, Smirnova I, Ilyukha V, Gautam A, et al. ZBP1 promotes LPS-induced cell death and IL-1β release via RHIM-mediated interactions with RIPK1. Nat Commun. 2021;12.(2021)org/10.1038/s41467-020-20357-z.: 86.
  50. 50.Suzuki A, Kusakai G, Kishimoto A, Shimojo Y, Miyamoto S, Ogura T, et al. Regulation of caspase-6 and FLIP by the AMPK family member ARK5. Oncogene. 2004;23(42).(2004)onc.1207963.: 7067.
  51. 51.Velázquez-Delgado EM, Hardy JA. Phosphorylation regulates assembly of the caspase-6 substrate-binding groove. Structure. 2012;20(4).(2012)742–51. https://doi. org/10. 1016/j.str.: 742.
  52. 52.Zhao P, Sun X, Chaggan C, Liao Z, In Wong K, He F, et al. An AMPK-caspase-6 axis controls liver damage in nonalcoholic steatohepatitis. Science. 2020;367(6478).(2020)1126/science.aay0542.: 652.
  53. 53.Zheng M, Karki R, Kancharana B, Berns H, Pruett-Miller SM, Kanneganti TD. Caspase-6 promotes activation of the caspase-11-NLRP3 inflammasome during gram-negative bacterial infections. J Biol Chem. 2021;297(6).(2021)101379. https://doi. org/10. 1016/j.jbc.: 101379.
  54. 54.Lin Y, Sheng M, Qin H, Zhang P, Wang C, Fu W, et al. Caspase 6 promotes innate immune activation by functional crosstalk between RIPK1-IκBα axis in liver inflammation. Cell Commun Signal. 2023;21(1).(2023)org/10.1186/s12964-023-01287-x.: 282.
  55. 55.Samson AL, Zhang Y, Geoghegan ND, Gavin XJ, Davies KA, Mlodzianoski MJ, et al. MLKL trafficking and accumulation at the plasma membrane control the kinetics and threshold for necroptosis. Nat Commun. 2020;11.(2020)org/10.1038/s41467-020-16887-1.: 3151.
  56. 56.Gaytán F, Morales C, Bellido C, Aguilar E, Sánchez-Criado JE. Ovarian follicle macrophages.(1998)is follicular atresia in the immature rat a macrophage-mediated event? Biol Reprod.: 52.
  57. 57.Miyake Y, Sakurai M, Tanaka S, Tunjung WA, Yokoo M, Matsumoto H, et al. Expression of hyaluronan synthase 1 and distribution of hyaluronan during follicular atresia in pig ovaries. Biol Reprod. 2009;80(2).(2009)108.067694.: 249.
  58. 58.Miyake Y, Matsumoto H, Yokoo M, Miyazawa K, Kimura N, Tunjung WA, et al. Expression and glycosylation with polylactosamine of CD44 antigen on macrophages during follicular atresia in pig ovaries. Biol Reprod. 2006;74(3).(2006)105.045641.: 501.
  59. 59.Yuan F, Cai J, Wu J, Tang Y, Zhao K, Liang F, et al. Z-DNA binding protein 1 promotes heatstroke-induced cell death. Science. 2022;376(6593).(2022)1126/science.abg5251.: 609.
  60. 60.Geng Y, Ma Q, Liu YN, Peng N, Yuan FF, Li XG, et al. Heatstroke induces liver injury via IL-1β and HMGB1-induced pyroptosis. J Hepatol. 2015;63(3).(2015)622–33. https://doi. org/10. 1016/j.jhep.: 622.
  61. 61.Zhou J, Qin X, Li L, Tian D, Zou Z, Gu Z, et al. Heat stress-induced intestinal epithelial cells necroptosis via TLR3-TRIF-RIP3 pathway was dependent on p53. Int Immunopharmacol. 2023;122.(2023)110574. https://doi. org/10. 1016/j.intimp.: 110574.
  62. 62.Huang W, Xie W, Gong J, Wang W, Cai S, Huang Q, et al. Heat stress induces RIP1/RIP3-dependent necroptosis through the MAPK, NF-κB, and c-Jun signaling pathways in pulmonary vascular endothelial cells. Biochem Biophys Res Commun. 2020;528(1).(2020)206–12. https://doi. org/10. 1016/j.bbrc.: 206.
  63. 63.Sundaram B, Pandian N, Mall R, Wang Y, Sarkar R, Kim HJ, et al. NLRP12-PANoptosome activates PANoptosis and pathology in response to heme and PAMPs. Cell. 2023;186(13).(2023)2783–801. e20. https://doi. org/10. 1016/j.cell.: 2783.
  64. 64.Cheng CY, Tu WL, Chen CJ, Chan HL, Chen CF, Chen HH, et al. Functional genomics study of acute heat stress response in the small yellow follicles of layer-type chickens. Sci Rep. 2018;8.(2018)org/10.1038/s41598-017-18335-5.: 1320.
  65. 65.Jiang JY, Shimizu T, Sasada H, Tsang BK, Sato E. Increased ovarian follicular angiogenesis and dynamic changes of follicular vascular plexuses induced by equine chorionic gonadotropin in the gilt. Cell Tissue Res. 2004;316(3).(2004)org/10.1007/s00441-004-0887-5.: 349.
  66. 66.Brown HM, Russell DL. Blood and lymphatic vasculature in the ovary.(2014)development, function and disease.Hum Reprod Update.: 29.
  67. 67.Liu Z, Feng C, Li C, He T, Wu G, Fu C, et al. IGF-I protects porcine granulosa cells from hypoxia-induced apoptosis by promoting homologous recombination repair through the PI3K/AKT/E2F8/RAD51 pathway. FASEB J. 2024;38(1).(2024)1096/fj.202301464R.
  68. 68.Xin X, Yang K, Liu H, Li Y. Hypobaric hypoxia triggers pyroptosis in the retina via NLRP3 inflammasome activation. Apoptosis. 2022;27(3–4).(2022)org/10.1007/s10495-022-01710-7.: 222.
  69. 69.Cheng SB, Nakashima A, Huber WJ, Davis S, Banerjee S, Huang Z, et al. Pyroptosis is a critical inflammatory pathway in the placenta from early onset preeclampsia and in human trophoblasts exposed to hypoxia and endoplasmic reticulum stressors. Cell Death Dis. 2019;10(12).(2019)org/10.1038/s41419-019-2162-4.: 927.
  70. 70.Huang Y, Feng Y, Cui L, Yang L, Zhang Q, Zhang J, et al. Autophagy-related LC3 accumulation interacted directly with LIR containing RIPK1 and RIPK3, stimulating necroptosis in hypoxic cardiomyocytes. Front Cell Dev Biol. 2021;9.(2021)679637. https://doi. org/10.3389/fcell.: 679637.
  71. 71.Wang S, Zheng Y, Li J, Yu Y, Zhang W, Song M, et al. Single-cell transcriptomic atlas of primate ovarian aging. Cell. 2020;180(3).(2020)585–600. e19. https://doi. org/10. 1016/j.cell.: 585.
  72. 72.Navarro-Pando JM, Alcocer-Gómez E, Castejón-Vega B, Navarro-Villarán E, Condés-Hervás M, Mundi-Roldan M, et al. Inhibition of the NLRP3 inflammasome prevents ovarian aging. Sci Adv. 2021;7(1).(2021)1126/sciadv.abc7409.
  73. 73.Li D, Chen J, Guo J, Li L, Cai G, Chen S, et al. A phosphorylation of RIPK3 kinase initiates an intracellular apoptotic pathway that promotes prostaglandin(2α)-induced corpus luteum regression. Elife. 2021;10.(2021)7554/eLife.67409.
  74. 74.Miao Y, Cui Z, Gao Q, Rui R, Xiong B. Nicotinamide mononucleotide supplementation reverses the declining quality of maternally aged oocytes. Cell Rep. 2020;32(5).(2020)107987. https://doi. org/10. 1016/j.celrep.: 107987.
  75. 75.Zhang Y, Bai J, Cui Z, Li Y, Gao Q, Miao Y, et al. Polyamine metabolite spermidine rejuvenates oocyte quality by enhancing mitophagy during female reproductive aging. Nat Aging. 2023;3(11).(2023)org/10.1038/s43587-023-00498-8.: 1372.

Acknowledgements

We thank Professor Junyi Jiang and Xueqing Liu from the National Protein Center for providing technical support with gel filtration.

Funding

This work was supported by the National Key Research and Development Program of China (2022YFD1300405), the National Natural Science Foundation of China (32202671), and the Open project of State Key Laboratory of Animal Biotech Breeding (Grant No. 2024SKLAB6-5).

Ethics Declaration

Ethics approval and consent to participate

Animal management in this study was approved by Ethical Committee of China Agricultural University with the permit number AW21404202-1-1.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Rights and Permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/ . The Creative Commons Public Domain Dedication waiver ( http://creativecommons.org/publicdomain/zero/1.0/ ) applies to the data made available in this article, unless otherwise stated in a credit line to the data. Reprints and permissions