Phoenixin as a New Target in the Development of Strategies for Endometriosis Diagnosis and Treatment

Small integral membrane protein 20/phoenixin (SMIM20/PNX) and its receptor GPR173 (G Protein-Coupled Receptor 173) play a role in the regulation of the hypothalamic–pituitary–gonadal axis (HPG). The aim of the study was to determine PNX, FSH, LH, and 17β-estradiol association in women with endometriosis, and the expression of SMIM20/PNX signaling via GPR173. Serum PNX, FSH, LH, and 17β-estradiol concentrations were measured by enzyme and electrochemiluminescence immunoassay. SMIM20/PNX and GPR173 expression in the eutopic and ectopic endometrium was assessed by qPCR and immunohistochemistry. Reduced PNX level, increased LH/FSH ratio and elevated 17β-estradiol concentration were found in patients with endometriosis. No differences in SMIM20 expression were observed between the studied endometria. GPR173 expression was lower in ectopic than in eutopic endometria. SMIM20 expression was mainly restricted to stroma. GPR173 was detected in some eutopic and ectopic stromal cells and in eutopic glandular epithelial cells. Discriminant analysis indicates the diagnostic relevance of PNX and LH/FSH ratio in patients with endometriosis. In women with endometriosis, reduced PNX levels and GPR173 expression may be responsible for HPG axis dysregulation. These new insights may contribute to a better understanding of the pathophysiology of endometriosis and provide the basis for a new strategy for diagnosis and treatment of endometriosis.


Introduction
Endometriosis is one of the most common benign gynecological disorders. It is defined as the presence of stroma and glandular cells of the endometrium outside the uterine cavity [1]. It affects approximately 6-10% of women of reproductive age, increasing to 20-30% in women diagnosed with infertility [2]. Despite numerous studies, the etiology of endometriosis is still unknown. The most widely accepted theory of its formation is Sampson's transplantation, and implantation theory, published a century ago. According to this theory, endometriosis results from retrograde menstrual reflux of endometrial fragments into the peritoneal cavity [3,4]. Noteworthy, most women have retrograde menstruation, and only a few develop endometriosis. Therefore, imbalance of the immune and endocrine systems, and the influence of environmental factors that ensure the implantation, growth, and persistence of the endometrium in ectopic sites are currently indicated [1][2][3][4][5].
The unclear etiology of endometriosis is the reason that a causal endometriosis treatment is essentially non-existent. Currently used pharmacotherapy masks the symptoms of the disease, often causing the occurrence of numerous side effects. This is sufficient to justify the search for new therapies, especially those that act causally.
In women, it has been shown that PNX can regulate the sexual cycle and ovarian follicle maturation and affect HPG axis by increasing hypothalamic (kisspeptin, GnRH) and pituitary (luteinizing hormone LH, follicle-stimulating hormone FSH) hormone production [7,9]. Decreased PNX-14 levels resulted in a delayed estrus phase in rats and a decreased GnRH receptor in the pituitary. Therefore, expression of PNX-14 is necessary for normal estrus cyclicity [7]. A positive correlation was observed between PNX-14 expression and LH, FSH, and testosterone and a negative correlation with serum estradiol and insulin levels in women with polycystic ovary syndrome (PCOS) [10]. In a rat model of PCOS, there was an increase in SMIM20 expression and PNX-14 synthesis in serum [11]. Moreover, GnRH analogs are often used to treat endometriosis and have been found to increase SMIM20 expression and decrease GPR173 receptor expression in the rat hypothalamus, pituitary, and ovaries [12]. It is worth noting that both SMIM20 and GPR173 are expressed in the human endometrium [13,14]. Currently, nothing is known about the role of PNX-14 in endometriosis.
The present study aimed to investigate whether patients with ovarian-localized endometriosis relative to patients without endometriosis have altered serum PNX and FSH, LH and 17β-estradiol levels. Expression of the SMIM20 precursor protein and the GPR173 receptor was also evaluated. We assume that the action of PNX outside the HPG axis is GnRH-dependent. In view of the reduced expression of the GnRH receptor, and GnRH itself, in the eutopic endometrium of women with endometriosis, altered expression of SMIM20 and the GPR173 receptor in the eutopic and ectopic endometrium of these women is expected.

Biological Materials
The study used archival material collected from 2008-2012: frozen sera, isolated RNA, and formalin-fixed paraffin-embedded endometrium. The age of patients included in the study ranged from 28 to 42 years and all had normal, regular menstrual cycles. The control group consisted of patients undergoing laparoscopy for non-endometrial pathologies (leiomyomata, cervical intraepithelial neoplasia, and benign ovarian cyst). The studied group consisted of women diagnosed with ovarian endometriosis (III/IV degree). Patients with endometriosis had not been treated pharmacologically for at least 6 months. Exclusion criteria for the studied group were gynecological malignant neoplasms and hormonal therapy. All women were operated on in the first phase of the menstrual cycle, between days 7th and 12th.
These preliminary classifications were then verified by histological analysis. Women with histological evidence of necrosis and active inflammation were excluded from the study.
Frozen sera were obtained from peripheral blood before surgery from the antecubital vein of 23 healthy volunteers and 32 patients diagnosed with endometriosis.
Archival RNA was isolated from eutopic and ectopic endometria obtained from 108 patients: 46 controls and 62 with ovarian endometriosis.
All samples were assayed without any dilution. The final fluorescence intensity values were measured using a microplate reader (Synergy 2 Multi-Mode Microplate Reader, BioTek, Winooski, VT, USA) with the excitation and emission wavelengths of 325 and 420 nm, respectively. The RNA template, primers, and water mixture were incubated for 10 minutes at 65 • C and subsequently chilled on ice. After adding the remaining reaction compounds, the reaction was processed as described previously [16] and obtained cDNA was used immediately or stored at −20 • C.

Real-Time PCR
mRNA expression of SMIM20 and GPR173 were assessed by quantitative PCR in a LightCycler 2.0 carousel glass capillary-based system (Roche, Manheim, Germany). Eva Green was used as a detection dye. The qPCR reaction mixture of 10 µL contained: 1X HotFire Pol Eva Green qPCR Mix Plus, 5 pm/µL of each genes' specific forward and reverse primers (Table 1), and 2 µL cDNA. The following thermal profile was applied: pre-incubation step (12 minutes, 95 • C), followed by 45 cycles of amplification with the end-point acquisition (15 s at 95 • C, 30 s at 60 • C and 20 s at 72 • C for SMIM20, B2M, GAPDH, HPRT1 or 15 s at 95 • C, 20 s at 55 • C and 6 s at 72 • C for GPR173). Amplification step was followed by a melting curve analysis at 95 • C for 0 sec and 65 • C for 15 s and 97 • C at the step acquisition mode (continuous fluorescence measurement; ramping rate, 0.1). Three reference genes: beta-2-microglobulin (B2M), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and hypoxanthine phosphoribosyltransferase 1 (HPRT1) were analyzed to choose the most stable. As a negative control, the reverse transcription reaction of mixed mRNAs without reverse transcriptase was used. qPCR reactions were made in duplicates. Standard curves of each gene were performed to calculate the efficiency of the PCR reaction using serial dilutions of cDNA.
mRNA expression levels were normalized to the most stable reference gene and calculated using the relative quantification method. The concentration ratio (Cr) was used in the further analyses.
2.5. Immunohistochemistry 2.5.1. SMIM20 Immunostaining Slides were deparaffinized and rehydrated as mentioned above. After antigen retrieval and blocking in 2.5% goat serum, sections were incubated overnight at 4 • C with 1:500 polyclonal anti-SMIM20 antibody (ThermoFisher Scientific, Carlsbad, CA, USA). Next, slides were washed in TBS-T and stained with a secondary anti-rabbit DyLight 594-conjugated antibody. DAPI was used to detect nuclei. Imaging was performed using a Zeiss LSM 780 confocal microscopy system (Carl Zeiss Meditec AG, Jena, Germany). In all immunohistochemical negative control reactions, the primary antibody incubation step was omitted.
In the IHC reactions, first, slides were incubated in TBS-T buffer with 2.5% horse serum at room temperature for one hour to block the non-specific binding of the antibody. Next, sections were incubated in a humid chamber overnight at 4 • C with rabbit polyclonal anti-GPR173 antibodies (1:500; ThermoFisher Scientific, Waltham, MA, USA). Slides were then washed twice in TBS-T buffer (5 minutes) and incubated in darkness at room temperature for 1 h with a secondary horse anti-rabbit DyLight 488-conjugated antibody (Vector Laboratories, Inc., Burlingame, CA, USA).
After rinsing the unbounded antibodies three times in TBS-T buffer for 5 min and blocking in TBS-T buffer supplemented with 2.5% goat serum (Vector Laboratories, Inc., Burlingame, CA, USA), subsequently, second immunohistochemistry staining was per-formed. The slides were incubated overnight at 4 • C with polyclonal anti-PNX-14 antibody (1:500; Phoenix Pharmaceuticals, Inc., Burlingame, CA, USA). After washing in TBS-T buffer (2 × 5 minutes), slides were incubated with a secondary goat anti-rabbit DyLight 594-conjugated antibody (in the dark, room temperature, 1 h Vector Laboratories, Inc., Burlingame, CA, USA). Next, slides were washed 3X in TBS-T buffer and stained with 1 µg/mL DAPI at room temperature for 5 minutes (ThermoFisher Scientific, Carlsbad, CA, USA) to visualize the nuclei. Imaging was performed using a Zeiss LSM 780 confocal microscopy system (Carl Zeiss Meditec AG, Jena, Germany).

Statistical Analyses
Statistical analyses were performed using Statistica ® Version 13.5.0 software for Windows (TIBCO Software Inc., Palo Alto, CA, USA). The results were compared in groups: controls vs. cases. All continuous variables were checked for outliers and were winsorized if any were present using the equation (mean ± 2* standard deviations) [18]. The Shapiro-Wilk test was used for the normality of continuous variable distribution assessment. The median and interquartile range were used to describe experimental results. The differences in expression levels and serum concentration between the controls and cases were evaluated using the Mann-Whitney U test, and the results were described as median [Q1 (lower quartile)-Q3 (upper quartile)]. Associations between hormone levels were correlated using a nonparametric Spearman's rank test. Moreover, using a system of linear discriminant equations, discriminant analysis of the serum hormone concentrations profile was performed. The patient was classified into the group for which the value of the equation was the highest. Data were considered statistically significant at p < 0.05. Reference gene expression stability was analyzed using the Microsoft Excel plugin NormFinder v0.953 (available at: https://moma.dk/normfinder-software, last accessed: 8 October 2021).

Analysis of PNX, FSH, LH, 17β-Estradiol Serum Levels
The statistical analysis of PNX, FSH, LH and 17β-estradiol concentrations were presented in Table 2   The final discriminant analysis model used for the tested serum hormone levels includes the PNX concentration and the LH/FSH ratio. Table 3 shows the Wilks' Lambda values, F test, and significance levels. Based on the presented results, the proposed model including PNX level and LH to FSH ratio distinguishes sera of women with endometriosis from healthy women at a confidence level of up to 99.9%. The classification accuracy is 96% in the case of women with endometriosis and 67% in the case of healthy patients.

Real-Time PCR: SMIM20 and GPR173 Expression
Three different reference genes were analyzed: B2M, GAPDH, and HPRT1. Accordingly, using NormFinder software, the most stable was B2M apart from tissue origin, and thus this gene was used for normalization and concentration ratio calculations. GPR173 and SMIM20 mRNA expression was confirmed in 104 and 108 analyzed samples, respectively. (p = 0.0215). There was no difference observed between the groups in the case of the SMIM20 normalized expression (p = 0.7702). The Mann-Whitney U test results are presented in Table 4. Analyzed gene expression Cr values, normalized to B2M, are presented in Figure 1.   Figure 2 shows SMIM20 (red color) immunolocalization in eutopic (a-d) and ectopic (e-h) endometrium. SMIM20 is a precursor protein of PNX-14, and its expression was restricted mainly to the stromal cells (Figure 2).  Figure 2 shows SMIM20 (red color) immunolocalization in eutopic (a-d) and ectopic (e-h) endometrium. SMIM20 is a precursor protein of PNX-14, and its expression was restricted mainly to the stromal cells ( Figure 2).  Figure 2 shows SMIM20 (red color) immunolocalization in eutopic (a-d) and ectopic (e-h) endometrium. SMIM20 is a precursor protein of PNX-14, and its expression was restricted mainly to the stromal cells ( Figure 2).  in some eutopic and ectopic stromal cells and eutopic glandular epithelial cells. Interestingly, some stromal cells showed a positive signal of PNX-14 and GPR173, suggesting a possible autocrine effect of PNX-14.

Discussion
An investigation of PNX in the proper functioning of the HPG axis was initiated by Yosten et al. in 2013 [7]. Researchers hypothesized that PNX might affect the glandular portion of the pituitary gland. Further studies have shown administration of PNX-14 or

Discussion
An investigation of PNX in the proper functioning of the HPG axis was initiated by Yosten et al. in 2013 [7]. Researchers hypothesized that PNX might affect the glandular portion of the pituitary gland. Further studies have shown administration of PNX-14 or PNX-20 increases the release of GnRH due to a regulatory effect on the GnRH receptor (GnRH-R). Furthermore, intrathecal administration of phoenixin significantly reduced pain perception in mice injected intraperitoneally with acetic acid. As a corollary, phoenixin may preferentially suppress visceral pain [19].
Reduction of PNX expression in the hypothalamus using small interfering RNA in vivo has been shown to inhibit the onset of estrus in female rats. By reducing the expression of the pituitary GnRH receptor, this suggests a role of PNX in normal ovarian cyclicity [12]. Since endometriosis is an estrogen-dependent disease, GnRH analogs are used in its treatment and work by inhibiting the pituitary-ovarian axis. Moreover, it has been shown that the direct action of gonadotropin-releasing hormone 1 (GnRH I) on endometriotic cells inhibits proliferation and induction of apoptosis [20]. These observations prompted us to investigate the expression of SMIM20, the precursor protein of PNX, and GPR173, the putative PNX receptor, in eutopic and ectopic endometrium. Additionally, we also assessed the level of PNX, FSH, LH and 17β-estradiol in the serum of women diagnosed with ovarian endometriosis and compare it to the level of normal controls.
The receptor for PNX is considered to be GPR173 based on a ligand binding study by Stein et al. 2016 [21]. Interestingly, PNX regulates Oct-1 and C/EBPβ transcription factor activity, essential for regulating GnRH expression. Subsequent research into the role of GPR173 in the hypothalamus has shown that PNX is essential for a healthy estrus cycle. Injection of siRNA targeting GPR173 has been shown to double estrus cycle length in female rats and eliminate the PNX-induced increase in plasma LH concentration [8]. In summary, it is assumed that PNX potentially has two methods of action on the HPG axis. First is a hypothalamic mode of action as PNX may directly trigger the release of GnRH from hypothalamic neurons. Alternatively, the pituitary may directly activate pituitary gonadotrophs via the pituitary portal vessels. In our studies, we found no expression of SMIM20 differences between eutopic and ectopic endometria and a decrease in the expression of GPR173 in the ectopic endometrium. Surprisingly, the PNX serum level was significantly lower in women with endometriosis compared to non-endometriosis cases. Based on these results, we speculate that reduced serum PNX levels in women with endometriosis may be responsible for the intensified pain perception and hormonal imbalance maintained by the HPG axis. A correlation of phoenixin levels in sera of endometriosis patients with the level of pain would have been performed to confirm this hypothesis. Sadly, this was not possible in the present study since the study was performed on archival material.
Due to the limited numbers of studies concerning the central and peripheral effects of PNX on diminished fertility, these results are difficult to interpret. Due to decreased PNX concentration in the serum, we expect a disturbance of the central control of the sexual cycle. This likely leads to more frequent anovulatory cycles, forming luteinized unruptured follicle functional cysts. This disorder creates improper progesterone levels in phase II of the cycle therefore promoting the development of endometriosis. At the same time, it does not harm the development of endometrial implants residing in the peritoneum or the ovary. In these foci, there is a self-propelling mechanism of estrogen production through the activity of COX-2 in the foci of the ectopic endometrium. The negative influence on the cycle diminishes the protective effect of progesterone, and at the same time, the ectopic focus of the endometrium retains its autonomy. Furthermore, it is likely that in women with stage III or IV endometriosis, the control axis of GnRH release under PNX is simply disabled. This led us to perform an evaluation of the association of the hormonal profile of FSH, LH and 17β-estradiol with PNX in the blood serum of the studied patients. The analysis revealed an increased LH to FSH ratio and 17β-estradiol levels in the serum of women with endometriosis. Discriminant analysis showed that the LH/FSH ratio and the level of PNX can be used as an algorithm for the non-invasive procedure for detection of ectopic endometrial foci. As mentioned previously, PNX is necessary to maintain cyclicity of both ovaries and thus endometrial changes. The absence of PNX effects results in the impaired release of GnRH. Presumably, in the absence of GnRH in cells forming the ectopic focus, a mechanism initiating intracellular signaling events, such as modulation of transcription factors regulating SMIM20 and GPR173 gene expression, is triggered. Another hypothesis that may explain the decreased PNX expression in serum of women with endometriosis is that GPR173 is not the unique receptor for PNX binding. PNX may be a ligand for an unidentified membrane receptor, not excluding GnRH-R itself. This assumption is the most appropriate and consistent with the studies performed in rats. The authors explain the down-regulation of GPR173 and the increased level of SMIM20 by the existence of molecular interactions between GnRH receptors and PNX signaling in the HPG axis of female rats during the reproductive cycle [12].
Based on immunohistochemical studies, no difference in staining intensity was found between the eutopic and ectopic endometrium. Positive staining for GPR173 was found in eutopic endometrial glands and several stromal cells. In ectopic endometrium, the localization of the examined receptor was restricted to only some of the stromal cells.
Remarkably, some fibroblasts in the studied endometria showed a positive signal not only from GPR173 but also from PNX, suggesting the possibility of an autocrine mechanism of PNX action. On the other hand, in the case of SMIM20, expression was mainly confined to stromal cells. This is the first publication presenting data on tissue-specific localization and expression of SMIM20, the precursor protein of PNX-14, and its receptor, GPR173, in the eutopic and ectopic endometrium. The specificity of transcript localization requires further research. Moreover, decreased serum PNX-14 concentration in patients with endometriosis suggests the role of PNX-14 in disease pathogenesis as well as in enhancing pelvic pain associated with cyclic changes within the ectopic endometrium.
These new insights may provide not only a better understanding of endometriosis pathophysiology but also lay the potential groundwork for the diagnosis and treatment of endometriosis by targeting SMIM20/PNX. It is worth mentioning that GnRH analogs are often used in the treatment of endometriosis, inducing a menopause-like condition, thus creating side effects that discourage patients from therapy continuation. Perhaps agonists and/or antagonists of PNX or the GPR173 receptor could be used in the future to alter GnRH neuronal function to completely cure endometriosis.

Informed Consent Statement:
Written informed consent has been obtained from all subjects involved in the study.

Data Availability Statement:
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.