1. Introduction
Human reproduction is a complex sequential process that includes the long-lasting journey of motile sperm through diverse compartments of the female genital tract and the final fusion with mature oocytes in the ampullary part of the oviduct. In recent years, the physiology of heterotopic olfactory receptors has become increasingly important in human reproduction, in particular, related to sperm chemotaxis [
1]. Olfactory receptors (ORs) are G-protein-coupled receptors, capable of detecting volatile chemicals that act as odorants [
1]. The activation of ORs initiates an intracellular signal transduction pathway triggered by the stimulation of adenylate cyclase III, which leads to Ca
2+ influx and membrane depolarization thereafter [
2].
Firstly, detected in the olfactory epithelium, ORs are subsequently identified in several tissues that include, but are not limited to, the brain, bladder, prostate, and testis [
1,
3,
4,
5,
6,
7,
8,
9]. However, the physiological consequences of the activation of ORs outside of the nose is still uncertain due to the scarce knowledge of the odorant ligands involved.
As recent studies have shown, the testis possess the richest OR expression outside of the nose [
6,
9,
10]. Late spermatids and spermatozoa have been shown to have a selectively distinct anatomical compartment expression of ORs [
10,
11,
12] suggesting their possible involvement in processes, such as maturation, migration, and fertilization [
13].
Spermatozoa chemotaxis is mostly dependent on intracellular calcium release. ORs seem to be among the receptors involved in activating this signaling pathway [
10,
11,
12,
14,
15]. However, parallel Ca
2+-independent signaling has been described, as evidenced by a signal transduction cascade independently activated by adenylate cyclase’s second messenger cAMP [
12,
16]. Of note, the CatSper calcium channel in human spermatozoa seems to be stimulated by odorant ligands too [
16]. Neuhaus and colleagues reported ERK phosphorylation as a further molecular signature related to OR-mediated sperm activation. Indeed, they treated sperm cells with the odorant bourgeonal, and they showed an accumulation of pERK upon the odorant stimulus, concomitant with the activation and acrosomal relocalization of the cognate receptor OR17-4 [
11]. As a general sign of sperm-cell activation, it is tempting to speculate that this molecular signature underlying the OR activity is related to sperm chemotaxis.
In fact, growing evidence of the close relationship between ORs and spermatozoa chemotaxis is spreading in the current literature. OR17-4 (also known as OR1D2) is localized in the spermatozoon midpiece [
11]. In 2003, Spehr and colleagues discovered the aforementioned bourgeonal, a ligand for human OR17-4, capable of increasing Ca
2+ as well as inducing chemotaxis in human spermatozoa [
14]. Although bourgeonal is most likely not an endogenous compound in the human body, the study supported the idea that human ORs act as chemoreceptors for small molecules in sperm. Interestingly, a reduced olfactory perception of burgeonal has been correlated with idiopathic infertility [
17].
Other odorants, such as Myrac and PI-23472, ligands of OR7A5 and OR4D1, respectively, have been demonstrated to influence the activation of spermatozoa motility. These same receptors seem to be activated by other ligands found in vaginal secretions and follicular fluid, in particular, 5α-androst-16-en-3-one for OR4D1 and 4-hydroxy-2,5-dimethyl-3[2H]-furanone for OR7A5 [
18].
Among the evaluated odorants, great importance has recently been attributed to volatile short-chain fatty acids (SCFAs), the metabolic by-products of anaerobic bacterial fermentation. OR51E2 and its murine ortholog (Olfr78) have been reported to respond to lactate [
19] and SCFAs, such as propionate (hereafter PA) and acetate (hereafter AA), in the kidneys [
20]. However, OR51E2 is not exclusively limited to the kidneys. Indeed, Aisenberg and collaborators reported OR51E2 expression in airway smooth muscle cells, in which they found a cytoskeleton remodeling mediated by OR-SCFAs binding [
21]. Moreover, Flegel and colleagues showed that OR51E2 is expressed on the extracellular membrane of the flagella, midpiece, and, albeit less frequently, on the acrosomal cap of human spermatozoa [
10]. Intriguingly, Hartmann and colleagues detected SCFAs in human follicular fluid and vaginal secretion by gas chromatography–olfactometry. In particular, AA was observed in both biological fluids while PA was only in follicular fluid [
18].
As a further demonstration of PA presence in the female reproductive tract,
vaginimicrobium propionicum (VP), a Gram-stain positive anaerobic rod-shaped bacterium, was isolated using culturomics from vaginal discharge [
22]. VP belongs to the Propionibacteria family, which are producers of PA by definition. Despite the fact that these results provide evidence that SCFAs reside within gynecological fluids, the data in the literature is still missing about their specific role in human reproduction. Since OR51E2 has been observed on the surface of spermatozoa and it has been shown to bind PA in other cell types, we wondered if OR51E2/PA has a functional role in sperm activation. Additionally, OR51E2 could be involved in the interplay between microbiota and human reproduction as mediators of spermatozoa chemotaxis. The aim of the present study is to gain an insight into the role of OR51E2 as a player of the sperm chemotactic process through its binding to specific ligands in SCFAs.
3. Discussion
The observation that olfactory receptors are expressed in male and female reproductive systems could uncover how environmental components, represented by active ligands for the identified chemosensory receptors, affect successful reproduction.
The non-orthotopic expression of ORs, meant as the expression outside the olfactory epithelium, has been reported to date in many different organs and in the testis, too [
4]. Previous research from our group identified more than 20 ORs in sperm cells, playing roles in the guidance of spermatozoa towards the oocyte, at multiple levels [
9]. Furthermore, we showed the presence of different ORs in the ducts of the testis and epididymis [
9]. The presence of these receptors on mature cells could be related to sperm motility, while their presence in testicular tubules and epididymis might suggest specific roles in the regulation of sperm maturation and migration. Moreover, DNA sequencing of human testicular biopsies showed the expression of a human olfactory receptor, OR17-4, in both the testis and surface of mature spermatozoa. This OR was reported to influence the chemotactic movement of sperm cells in response to the odorant bourgeonal [
14]. A major limit in the proof of concept for the pro-chemotactic role of bourgeonal is that it has not been reported any endogenous corresponding molecule for this odorant [
24].
Indeed, although new information has emerged concerning the profile of distinct olfactory receptors and coupled molecular pathways in human reproductive systems, endogenous ligands of specific chemoreceptors have only been rudimentarily explored [
1]. In this study, we reported for the first time the relocalization of OR51E2 on a human sperm surface and the consequent molecular activation upon SCFAs treatment, in particular, PA. This in turn is related to sperm migration towards SCFA-enriched medium. Moreover, we evidenced the presence of PA and AA in cervical mucus and showed that an incubation of these samples with spermatozoa is able to trigger pERK accumulation, similar to PA and AA treatments. Consequently, due to the exclusive SCFA production by microbiota metabolism, the presence of SCFAs in the female genital tract could suggest a key role of microbiota in fertility and reproduction through, at least in part, the ORs. Similarly, the possible presence of SCFAs in the male reproductive system might also be involved in the regulation of sperm maturation and migration.
Previously, it was reported that ORs are specifically localized on distinct spermatozoon anatomical areas [
10,
11]. Flegel and colleagues showed that OR51E2 is expressed on the flagella, midpiece, and, occasionally, acrosomal cap [
10]. The authors concluded that different OR localizations can subtend different functions. We found that, upon PA stimulus, OR51E2 was able to relocalize from the flagella and midpiece to the acrosomal cap, just as it was previously reported for OR17-4 in response to the odorant bourgeonal [
11]. The OR51E2 acrosomal relocalization upon ligand stimulus suggests this receptor’s involvement in sperm chemotaxis towards the oocyte.
At a molecular level, it was reported that the process of chemotaxis requires the activation of several pathways. In particular, Teves and colleagues showed a temporal sequence of chemotactic signaling events in which the activation of adenylate cyclase was followed by a cascade of protein tyrosine phosphorylation and calcium influx [
24].
Within the protein-phosphorylation cascade, ERK phosphorylation (pERK) is a well-established signal for sperm activation [
25,
26]. We analyzed pERK levels in sperm cells upon PA stimulus and we found a dose and time-dependent accumulation of ERK protein phosphorylation, indicating that the OR51E2 acrosomal clusterization upon ligand binding could subtend a molecular cascade of sperm-cell activation.
We then assessed whether PA stimulus could lead to a Ca
2+ influx, which has been previously reported to be the key event in the activation of sperm chemotaxis [
24,
27]. In particular, its role as a downstream molecular effector of ORs–ligand binding should be noted [
14,
15]. Indeed, Spehr et al. demonstrated that the OR17-4 ligand, bourgeonal, can evoke calcium release and subsequently, sperm chemotaxis [
14]. We found that PA stimulus was able to induce a significant increase in intracellular Ca
2+ in comparison to the control, confirming the involvement of this central molecular signature in the SCFA-mediated activation of sperm cells.
Overall, these results support the hypothesis that PA, through OR51E2 acrosomal relocalization, is able to mediate molecular activation in sperm cells, involving ERK and calcium influx. We demonstrated that spermatozoa activation mediated by PA relies on OR51E2, because the adoption of OR51E2-specific antibody pre-treatment interferes with the previously described pathway. Both pERK accumulation and calcium influx promoted by PA stimulus were significantly reduced with α-OR51E2 antibody pretreatment, supporting the idea that the PA-mediated molecular activation of spermatozoa depends on the binding of PA to this olfactory receptor.
As stated above, the PA-mediated OR51E2 acrosomal clusterization on the outermost part of the sperm head suggests that this ligand–receptor binding could be involved in the sperm migratory response. To further address the involvement of PA-OR51E2 binding in sperm migration, we performed a modified direct-swim-up assay. Furthermore, we analyzed the differential spermatozoa migration towards a medium with or without PA or AA. In this model, higher quantities of sperm cells in both PA- and AA-enriched mediums were recovered in comparison to the control one. This result suggests a key role for SCFAs in sperm chemotaxis. The CASA analysis of kinematics parameters quantified the phenomenon. Indeed, we found higher VSL and directionality-related parameters, such as LIN and STR, upon SCFA addition. Moreover, we observed that VAP and VCL did not increase in the spermatozoa fraction migrated to PA- and AA-enriched mediums. VSL is the velocity (µm/s) recorded along the straight-line path of the spermatozoon, from the first to the last point of the recorded path itself; on the contrary, VCL is calculated as the velocity (µm/s) through the curvilinear path performed by the sperm cell. LIN is calculated as the percentage ratio VSL/VCL [
28]. Therefore, a higher percentage of straight-directed spermatozoa, instead of curvilinear ones, increased LIN. Indeed, a sperm cell characterized by a circling path would have low LIN due to a higher curvilinear path than a straight one [
29]. STR, which represents the straightness index, defined by the ratio VSL/VAP, increased in PA-treated spermatozoa as they were characterized by a straighter trajectory along the average recorded path. Therefore, the increase in LIN and STR, related to the increase in VSL, indicated that PA and AA not only promoted a general migration of sperm cells, but specifically a more linear orientation. These data fit with SCFA’s possible role as a chemoattractant through the activation of OR51E2. On the other hand, OR51E2 has been demonstrated to bind SCFA as a specific ligand able to promote downstream molecular activation in many different cellular models [
20,
30,
31]; intriguingly, SCFAs, such as PA and AA, are human endogenous compounds produced by gut microbiota metabolism, through dietary non-digestible carbohydrate fermentation [
32,
33]. Beyond gut microbiota, its presence has been reported in a vaginal swab from a woman without vaginosis of specific bacterial species producing propionic acid, the vaginomicrobium propionicum [
22]. The authors reported propionic and acetic acids as the major end products of bacterium metabolism. This novel microbial identification revealed the local presence of specific bacteria-producing SCFAs as part of the vaginal microbiota. Moreover, Hartmann and colleagues applied gas chromatography–olfactometry coupled to mass spectrometry analysis of both vaginal secretion and follicular fluid aiming to identify odorous compounds able to promote an in vitro activation of three different ORs (OR4D1, OR7A5, and OR1D2) [
18].
4. Materials and Methods
4.1. Sperm Samples: Ethical Approval
Sperm collection and analysis were approved by the local ethics committee of the Fondazione Policlinico Universitario A. Gemelli, IRCCS, Rome, Italy (protocol number ID 3943). Human sperm were freshly obtained by masturbation (3–5 days of sexual abstinence) from young, healthy donors who provided informed, signed consent. All samples fulfilled WHO’s 2021 semen-quality criteria for cell count (concentration > 15 million/mL), motility (motile cells > 40%), and viability (live cells > 58%).
4.2. Biological Preparations
The sperm were handled and prepared as previously described [
34,
35]. Briefly, liquefied semen was overlaid on a 50% Percoll (cell-culture tested, Sigma-Aldrich, MO, USA) density gradient and centrifuged at room temperature for 30 min at 400 g. The recovered sperm cells were then resuspended in 1 mL of phosphate-buffered saline (PBS; Sigma-Aldrich) and subjected to a residual leukocyte depletion using Dynabeads
TM CD45 (Invitrogen, Carlsbad, CA, USA). Then, the samples were washed twice and the pellet was collected, washed in Ringer’s solution (140 mM NaCl, 5 mM KCl, 2 mM CaCl
2, 2 mM MgCl
2, 10 mM Hepes, 10 mM glucose, pH 7.4), and again centrifuged for 15 min. Then, the motile spermatozoa pellet was resuspended in Ringer’s solution and used for further experiments.
4.3. Analysis of OR51E2 Expression and Localization
Immunofluorescence and confocal microscopy were used to detect the expression and localization of OR51E2 in human spermatozoa and HeLa cells. In HeLa cells (kindly provided by Dr. Fabiola Moretti), we overexpressed FlagRho-OR51E2 (kindly provided by Dr. Jennifer Pluznick) to verify the specificity of the α-OR51E2 antibody. Human sperm cells, purified as described above, were treated with PA 50 mM or Ringer’s solution (CTR) for 10 min at room temperature. Then, the sperm cells were washed and smeared on polylysine slides, and then allowed to air-dry. Then, HeLa and sperm cells were fixed in 3.7% formaldehyde for 15 min at room temperature (RT) and subsequently treated with 5% BSA to block nonspecific binding. Then, both cellular models were treated with polyclonal anti-OR51E2 (1:200) overnight (Aviva System Biology, San Diego, CA, USA) antibody, washed three times by PBS, and incubated for 45 min with goat Alexa Fluo-488 anti-rabbit IgG (A11034 Molecular Probes, Thermo Fisher, Waltham, MA, USA). For the HeLa cells, DNA was stained with Prolong Gold DAPI (P36935 Molecular Probes), while in sperm cells, acrosomal compartment was counterstained with lectin PNA (L32459 Molecular Probes) following the manufacturer’s instructions.
4.4. Confocal Microscopy
The cells were imaged with an inverted confocal microscope (Nikon A1-MP). Fluorescence images (excitation: 402 nm for the blue channel, 488 nm for green channel, and 561 nm for red channel) were collected in three separated channels (emission filter: 450/50 nm for the blue channel, 525/50 nm for the green channel, and 595/50 nm for the red channel) using a 60× immersion-oil objective with 1.4 Numerical Aperture (NA). Internal photon multiplier tubes collected 2048 × 2048-pixel images in 16 bit at a 0.063 ms dwell time.
4.5. Colocalization Analysis
To quantify the compartmentalization of OR51E2 within the acrosome of different spermatozoa, we applied a colocalization approach. Colocalization analysis was performed through the Colocalization Threshold plugin available via the open source software ImageJ (NIH). The analysis of fluorescence colocalization was graphically represented in scatterplots, where the intensity of the green channel was plotted versus the intensity of the red channel for each pixel. A proportional fluorescence intensity of the two probes resulted in the distribution of points along a straight line, with the slope reflecting the ratio of the fluorescence of the two probes. To quantify the intranuclear compartmentalization of OR51E2, we evaluated the Manders’ Colocalization Coefficients. In particular, for two probes, denoted as
B and
G, respectively, the coefficient
M1 provides the fraction of
B, i.e., the fraction of pixels expressing the OR51E2 protein, in compartments containing
G, which is PNA (a marker of the acrosome), according to the following expression:
The threshold value was automatically identified by applying the Costes method [
36] and the coefficient was only evaluated for pixels above threshold (
tM1), with values ranging from 1, when the protein was expressed in the whole acrosome, to 0, in case of no expression of the protein.
4.6. Western Blot
For the Western blot analysis, after PA treatment (5 and 50 mM for 5, 10, and 15 min) the sperm was lysed in RIPA buffer (50 mM Tris–Cl, pH 7.5, 150 mM NaCl, 1% Nonidet P-40, 0.5% Na deoxycholate, 0.1% SDS, 1 mM EDTA) and sonicated for 15 s at 20% intensity and 4 °C, and placed on ice for 30 min. The samples were then centrifuged for 15 min at 15,000× g and 4 °C, and supernatant containing proteins were then heat-denatured at 100 °C for 5 min in reducing Laemmli sample buffer 4× (BioRad, Philadelphia, PA, USA). Proteins were resolved by SDS-PAGE and then transferred onto PVDF membranes (Millipore, Burlington, MA, USA). All buffers contained a cocktail of protease inhibitors (Boehringer, Ingelheim am Rhein, Germany). Membranes were developed using enhanced chemiluminescence (ECL westar, Cynagen, Bologna, Italy). Bands were analyzed by a chemiluminescence imaging system, Alliance 2.7 (UVITEC, Cambridge, UK), and quantified by the software Alliance V_1607. The following primary antibodies were used: rabbit α-pERK (cell signaling 9101S) and rabbit α-ERK tot (cell signaling 9102S). In some experiments, to assess the specificity of ERK phosphorylation in response to PA binding to its specific receptor, the sperm cells were pre-treated with polyclonal antibody anti-OR51E2 (OAAF-05032 Aviva System Biology, San Diego, CA, USA) or with non-related antibody polyclonal anti-OR6B2 (NBP2-13711 Novus Biologicals, Milan, Italy) for 10 min at room temperature.
4.7. Measurement of Changes in [Ca2+]
The changes in the free intracellular Ca2+ concentration [Ca2+] in human sperm cells were measured in 96 multi-well plates in a fluorescence plate reader (Varioskan Lux, ThermoFisher, Walthman, MA, USA). Briefly, an aliquot of seminal fluids was incubated with the fluorescent Ca2+ indicator Fluo-4, AM (10 μM, Fluo-4 Direct Calcium Assay kit, ThermoFisher, Walthman, MA, USA) for 30 min at 37°C. Fluorescence was excited at 480 nm and emission was recorded at 520 nm with bottom optics. Fluorescence was recorded before and after, with or without the addition of compounds (PA 50 mM or bourgeonal 500 µM) and controls to duplicate wells. The changes in Fluo-4 fluorescence were shown as relative fluorescent signals (RFSs). The change in fluorescence signals was evaluated with respect to the basal fluorescence before the addition of compounds. In a set of experiments, to evaluate the specificity of the [Ca2+] increase upon PA binding to its specific receptor, the sperm cells were pre-treated with polyclonal antibody anti-OR51E2 or with non-related polyclonal antibody polyclonal anti-OR6B2 for 10 min at room temperature. The changes in Fluo-4 fluorescence were shown as relative fluorescent signals (RFSs). The change in fluorescence signals was evaluated with respect to the basal fluorescence before the addition of compounds.
4.8. Adenosine 3′5′-Cyclic Monophosphate (cAMP) Assay
The determination of cAMP levels was performed using the Parametertm cAMP assay (R&D Systems®, Inc., Minneapolis, MN, USA) according to the manufacturer’s instructions. Briefly, purified sperm cells were treated with 5 or 50 mM of PA for 15 min, and were then washed 3 times in cold PBS and resuspended in cell lysis buffer to a concentration of 1 × 107 cells/mL. Cells were subjected to three repeated cycles of freeze/thaw at −20 °C. Lysates were centrifuged at 600× g for 10 min at 4 °C to remove cellular debris. The supernatants were immediately assayed for cAMP determination. The absorbance was read at 450 nm. The intensity of the color was inversely proportional to the concentration of cAMP in the sample.
4.9. Modified Direct-Swim-Up Assay
After the fluidification of the seminal fluid, the entire volume was divided in fractions of 0.5 mL into round-bottom tubes, and 0.5 mL of culture medium was placed over the semen in each tube, adding or not PA 50 mM, AA 5 mM, or 0.6% BSA, as a known chemoattractant positive control. The tubes were put in the thermo block in a vertical position and incubated at 37 °C for 15 min. After the incubation, 0.2 mL of the upper part of each tube (distal fraction) was gently removed, aspirating the sperm from the upper meniscus downwards [
37]. Highly motile sperm cells recovered were counted using a Makler Counting Chamber.
4.10. Evaluation of Sperm-Motility Parameters by CASA
The computer-assisted sperm analysis (CASA) system (Microptic S.L., Barcelona, Spain) was used to analyze sperm motility. Motile spermatozoa obtained from modified swim-up, as described above, were evaluated for each treatment for the following parameters: average path velocity (VAP, µm/s), straight-line velocity (VSL, µm/s), curvilinear velocity (VCL, µm/s), straightness (STR, %), and linearity (LIN, %). Briefly, 3 μL of the sperm sample was pipetted and loaded into a pre-warmed (37 °C) standard-count four-chamber Leja slide (SC 20-01-04-B). To analyze sperm motility parameters, a minimum of 200 sperm at 5 different fields were examined in each group.
4.11. UPLC-MS/MS
Two short-chain fatty acids (SCFAs), acetic acid and propionic acid, were measured in cervical mucus samples by ultraperformance liquid-chromatography mass spectrometry (UPLC-MS/MS). After obtaining informed consent, cervical mucus samples were collected from healthy, fertile women attending routine gynecological visits. Cervical mucus was collected by a disposable swab at the ectocervix under vaginal speculum examination before any other operation. All participants voluntarily agreed to participate in this study approved by the local ethics committee of Fondazione Policlinico Universitario A. Gemelli, IRCCS, Rome, Italy—Protocol number ID 3943. These freshly collected cervical mucus samples were kept in sterile tubes, transported to the laboratory on ice, and stored at −80 °C until utilization. The UPLC-MS/MS system consisted of an UPLC and autosampler ExionLC AD system (ABSciex, Framingham, MA, USA) and a Qtrap 6500+ (ABSciex, Framingham, MA, USA) equipped with an electrospray ion source. Analyses were conducted in positive-ion mode. 4-Acetamido-7-mercaptobenzofurazan (AABD-SH) was purchased from Tokyo Chemical Industry (TCI, Nihonbashi-honcho, Tokyo, Japan). Short-chain fatty acids–Mixture 2 (acetic, propionic, isobutyric, n-butyric, 2-methylbutyric, isovaleric, and n-valeric acids) was purchased from Cayman Chemical (Cayman Chemical, East Ellesworth, MI, USA). Dichloromethane, triphenylhphosphine (TPP), 2,2′-Dipyridyl disulfide (DPDS), [
2H
4] acetic acid, and [
2H
3] propionic acid, were purchased from Sigma (Sigma-Aldrich, St. Louis, MO, USA). Water, acetonitrile, methanol, and formic acid (LC-MS grade) were purchased from Biosolve (Biosolve Chimie, Dieuze, France). The analytical procedure required a deproteinization step, performed by adding 300 µL of acetonitrile containing 15 µg/mL of [
2H
4] acetic acid and 15 µg/mL of [
2H
3] propionic acid to 100 µL of standard solution or 0.1 g of cervical mucus sample. After vigorous agitation, the sample was centrifuged at 14,000 rpm for 5 min at room temperature; then the supernatant was derivatized. For derivatization, 20 µL of AABD-SH (20 mM in dichloromethane), 20 µL of supernatant, 20 µL of TPP (20 mM in acetonitrile), and 20 µL of DPDS (20 mM in acetonitrile) were added to a glass tube and vigorously vortexed. Derivatization was performed for 15 min at room temperature. The reaction solution was dried under a nitrogen flow and reconstituted with 1 mL of methanol. Samples were loaded onto an ACQUITY PRM BEH C18 1.7 µ 2.1 × 100 mm (Waters Corporation, Milford, MA, USA). The chromatographic separation was performed with a gradient of mobile phases A (water containing 0.1% formic acid) and B (methanol containing 0.1% formic acid), with a flow rate of 0.500 mL/min. The gradient followed this pattern: 0–0.5 min 25% B, 0.5–12.5 min 50% B, 12.51–13.50 min 90% B, and 13.51–15.0 min 25% B. The oven temperature was set at 40 °C. The injection volume was 2 µL, and the total analysis time was 15 min. The optimized parameters for the ion source included: temperature at 450 °C, curtain gas at 20, nebulizing gas (GS1) at 60, drying gas (GS2) at 60, collision-activated-dissociation gas at medium, and ion-spray voltage at 5500 V. The multiple-reaction-monitoring (MRM) transitions for each analyte, their respective collision energy, and cone-voltage values are reported in
Table 1.
Data acquisition was performed using mass spectrometer software (Analyst Software 1.7.1, ABSciex, Framingham, MA, USA), while for quantitative analysis we used processing software (Sciex OS, ABSciex, Framingham, MA, USA).
4.12. Statistical Analysis
The results are expressed as values of mean ± standard deviation (SD). The statistical tests used were paired two-tailed Student’s or one-sample t-tests when the mean value of the control group was arbitrarily set to one. Significant difference was defined as p < 0.05. The software used was GraphPad Prism 7.04.