Research | Open | Published:
Spatio-temporal expression profile of matrix metalloproteinase (Mmp) modulators Reck and Sparc during the rat ovarian dynamics
Reproductive Biology and Endocrinologyvolume 16, Article number: 116 (2018)
Matrix metalloproteinases (Mmps) and their tissue inhibitors (Timps) are widely recognized as crucial factors for extracellular matrix remodeling in the ovary and are involved in follicular growth, ovulation, luteinization, and luteolysis during the estrous cycle. Recently, several genes have been associated to the modulation of Mmps activity, including Basigin (Bsg), which induces the expression of Mmps in rat ovaries; Sparc, a TGF-β modulator that is related to increased expression of Mmps in cancer; and Reck, which is associated with Mmps inhibition. However, the expression pattern of Mmp modulators in ovary dynamics is still largely uncharacterized.
To characterize the expression pattern of Mmps network members in ovary dynamics, we analyzed the spatio-temporal expression pattern of Reck and Sparc, as well as of Mmp2, Mmp9 and Mmp14 proteins, by immunohistochemistry (IHC), in pre-pubertal rat ovaries obtained from an artificial cycle induced by eCG/hCG, in the different phases of the hormone-induced estrous cycle. We also determined the gene expression profiles of Mmps (2, 9, 13 14), Timps (1, 2, 3), Sparc, Bsg, and Reck to complement this panel.
IHC analysis revealed that Mmp protein expression peaks at the early stages of folliculogenesis and ovulation, decreases during ovulation-luteogenesis transition and luteogenesis, increasing again during corpus luteum maintenance and luteolysis. The protein expression patterns of these metalloproteinases and Sparc were inverse relative to the pattern displayed by Reck. We observed that the gene expression peaks of Mmps inhibitors Reck and Timp2 were closely paraleled by Mmp2 and Mmp9 suppression. The opposite was also true: increased Mmp2 and Mmp9 expression was concomitant to reduced Reck and Timp2 levels.
Therefore, our results generate a spatio-temporal expression profile panel of Mmps and their regulators, suggesting that Reck and Sparc seem to play a role during ovarian dynamics: Reck as a possible inhibitor and Sparc as an inducer of Mmps.
The ovarian dynamics comprehended by different phases of the estrous cycle involves extensive tissue remodeling , and major extracellular matrix (ECM) reorganization is a key part of this process. This tissue remodeling is dependent upon cyclic hormonal variation, proteins and peptide growth factors and is a pre-requisite for expansion of the growing ovulatory follicle, breakdown of follicular walls during ovulation, luteinization of the postovulatory follicles, and regression of the corpus luteum (CL) [2, 3]. Matrix metalloproteinases (Mmps), a large family of zinc-dependent proteases, and their also numerous inhibitors (Timps) play a major role in ovarian ECM reorganization [2,3,4,5]. Together, Mmps and Timps form a network that regulates tissue remodelling and is modulated not only by gonadotropins and steroid hormones, but, also, by the expression of other factors, such as TNF-α, IGF-1 [5,6,7,8,9], and Bsg, which is an inducer of Mmps during the luteinization process . A complete map of the components of this network, as well as a precise characterization of their expression pattern, may allow us to better understand this tight regulation involved in ovarian tissue remodelling that occur during the estrous cycle.
It is already known from the literature that Mmp2, Mmp9, Mmp14 (Mt1-Mmp), and Mmp19 expression and activity in granulosa and/or theca cells are extensively linked to the disruption of ovarian ECM in rats, which is rich in collagen, laminin, and fibronectin, leading to follicular growth and release of the oocyte during the ovulatory process . MMP1, MMP2 e MMP9 have been detected in cultures of luteinized granulosa cells [12,13,14].
The expression and activity of Mmps have been known to be upregulated by gonadotropins during these phases [15,16,17,18,19]. However, hCG also was shown to reduce the expression of MMP2 and MMP9 [12, 13]. Following ovulation, the ruptured follicle is transformed into a CL by extensive cellular reorganization, migration, and neovascularization [20,21,22], which is associated with notable gelatinolytic activity, whereas the regression and absorption of CL is marked by an increase in Mmp13 levels and activity [23, 24].
Tissue inhibitors of metalloproteinases (Timps) have also been extensively studied in ovaries during the estrous cycle. Timp1 is upregulated by gonadotropins during folliculogenesis and ovulation in rats [25,26,27], providing proteolytic homeostasis to the ECM, through inhibition of metalloproteinases in the corpus luteum, while Timp3 levels decrease slightly .
Timp1 mRNA levels increase during luteal formation and regression [24, 29], while Timp3 levels are elevated during luteal maintenance . TIMP-2 has a role in the local regulation of MMP1 and MMP2 in the corpus luteum . Moreover, in addition to these inhibitors, there are also proteins, which have already been described as being responsible for inducing the expression of Mmps and ensuring system homeostasis in the ovary. In this context, it is worth mentioning Basigin (Bsg), a protein associated with Mmps induction and also with tumor progression and endometriosis, being expressed in granulosa and theca cells of pre-ovulatory follicles, and also in CL [10, 31]. However, the expression pattern of Mmp modulators in ovary dynamics is still largely uncharacterized, and there is evidence that members of the network are still to be identified. In particular, two proteins implied in tumor aggressiveness and invasive potential, through modulation of ECM integrity, namely, Sparc and Reck, may have an important role in ovarian tissue remodeling, since they are expressed in ovarian cells and are involved in the regulation of the Mmp system in other tissues [32,33,34,35,36].
Sparc (secreted protein acid and rich in cysteine), also known as Osteonectin, a 43 K protein, and basement-membrane protein 40 (BM-40), is a calcium-binding matrix cellular protein that plays a role in matrix mineralization , modulates TGF-β , and is expressed in the internal theca and follicular basal lamina of ovine follicles and theca-derived small luteal cells . Even though SPARC is highly expressed in several tumors, it may be able to inhibit tumorigenesis or tumor progression in human ovarian cancer  and its expression is reduced in ovarian cancer cell lines . Additionally, Sparc is implicated in indirect Mmp modulation and turnover of many physiological processes, but it has not yet been described in whole ovary dynamics. Despite the fact that Bagavandoss et al. described Sparc in the ovary of the late stages of estrous cycle, its expression in the initial phases (folliculogenesis, ovulation and ovulation-luteinization transition) of this process is still unknown . On the other hand, RECK (REversion-inducing-Cystein-rich protein with Kazal motifs) is a key regulator of extracellular matrix integrity and angiogenesis that is linked to the cell membrane by a GPI-anchor  and inhibits the activity of MMP2, MMP9, and MMP14 at different steps of their activation cascades in humans [36, 39]. Reck gene is expressed in several mouse tissues during development, including the uterus and ovaries, but at lower levels in most tissues, when compared to those of Timps 1, 2, 3, and 4 . Reck expression is inhibited by estrogen within the mouse uterus, but this inhibition is partially blocked by progesterone . To date, no information is available regarding the expression profile and function of Reck in rat ovaries development and maturation.
To analyse the ovarian extracellular matrix remodeling, we analyzed an artificial estral cycle, with a combination of equine chorionic gonadotropin (eCG) and human CG (hCG) in pre-pubertal rats , that have not yet started their estrous cycle, to induce multiple follicular development and ovulation. The rats only start their estrus cycles after the vaginal orifice opens, which tends to occur between postnatal 32 and 36 days . The literature describes that the administration of 10UI eCG stimulates the folliculogenesis and after 48 h the injection of 10UI of hCG can induce ovulation [42, 43]. In sequence, it is observed formation of corpus luteum, and 24 and 48 h post-hCG injection, the luteal formation in the ovaries. After 4 and 8 days, a peak of progesterone production is maximal resulting in the luteal maintenance, and corpus luteum regression can be analyzed at 14 days.
Therefore, understanding Reck and Sparc expression during the estrous cycle is essential to understand which role, if any, they play in the regulation of ovarian tissue remodeling, which allows all of the complex physiological processes that take place in the ovary. In this context, our hypothesis is that the Mmps regulators Reck and Sparc are modulated in ovarian dynamics possibly controlling the expression of Mmps, Reck being a possible inhibitor and Sparc a possible inducer of Mmps in this process.
To address these questions, we characterized the spatio-temporal expression pattern of Reck and Sparc proteins, as well as of Mmp2, Mmp9 and Mmp14 by immunohistochemistry in pre-pubertal rat ovaries during different phases of the hCG/eCG induced estrous cycle. To complement these data, we used qRT-PCR to investigate Mmps 2, 9, 13, 14, and 19, Timps 1, 2, and 3, Reck, Sparc, and Bsg mRNA levels and generate a complete expression profile panel of these targets during induced folliculogenesis, ovulation, luteogenesis, and luteolysis in prepubertal rat ovaries.
All animal experimentation was approved by the Ethics Committee for Animal Use of the Chemistry Institute (CEUA) on 08/28/2010, University of São Paulo, São Paulo, Brazil, in accordance with the National Council for the Control of Animal Experimentation (CONCEA).
To create an artificial estral cycle, we carried out the protocol to induce multiple follicular development and ovulation described by Jo and Curry and Ying and Meyer [42, 43]. First, 21 day-old pre-pubertal Sprague-Dawley female rats were induced by subcutaneous injection of 10 IU of equine chorionic gonadotropin (eCG) (Folligon, MSD Saúde Animal, São Paulo, SP, Brazil), resulting in animals in folliculogenesis phase. In sequence, ovulation and CL formation were induced by subcutaneous injection of 10 IU of human chorionic gonadotropin (hCG) (Choragon, Ferring GmbH, Kiel, Alemanha) 48 h after eCG injection. To analyze the relative expression pattern of the genes of interest throughout the hormonally induced estrous cycle, 10 experimental groups were established and different hormonal treatments were adopted to recreate the various stages of the estrous cycle. Treatment groups varied depending on the hormones applied (eCG only or both eCG and hCG) and induction and euthanasia periods of time. Animals were euthanized at 0, 24, or 48 h after eCG administration or 12, 24, or 48 h or 4, 8, or 14 days after hCG administration (Table 1). Two independent experiments were carried out, using three animals per experimental group. Animals injected with phosphate-buffered saline (PBS [137 mM NaCl, 2.7 mM KCl, 8 mM anhydrous Na2HPO4, and 1.4 mM KH2PO4; pH = 7.2]) were used as controls.
Both ovaries were collected from each animal. One ovary was stored in RNAholder solution (BioAgency Laboratories, São Paulo, SP, Brazil) at − 80 °C for total RNA extraction and the contralateral ovary was fixed in 4% paraformaldehyde solution (Sigma-Aldrich, St. Louis, MO, USA), stored at 4 °C, and subsequently included in paraffin for immunohistochemical analysis.
The paraffin-embedded tissues were sectioned into 4 μm-thick sections using a microtome (Leica Biosystems, Wetzlar, HE, Germany) and placed onto polylysine-treated slides. The sections were then dewaxed in xylene (Synth, São Paulo, SP, Brazil), incubated with decreasing concentrations of ethanol (Synth), and rehydrated in PBS-T (0.1% Triton X-100 in PBS [123 mM NaCl, 9.7 mM Na2HPO4, and 51.1 mM KH2PO4; pH = 7.4]). Next, tissue antigens were retrieved using phosphate-citrate buffer (pH = 6.0; Sigma-Aldrich) at 98 °C for 20 min, slowly cooled down to room temperature, and washed in PBS-T. To block endogenous peroxidase, the tissues were incubated with 3% hydrogen peroxide at room temperature for 45 min, washed in PBS-T, and incubated for 20 min in 3% nonfat milk solution in PBS for blocking of non-specific sites. The tissues were individually incubated overnight with rabbit polyclonal antibodies against Mmp2 (1:1000 dilution, product no. ab79781; Abcam, Inc., Cambridge, MA, USA), Mmp9 (1:100 dilution, product no. ab38898; Abcam, Inc.), and Sparc (1:100 dilution, product no. #5420; Cell Signaling Technology, Danvers, MA, USA) or rabbit monoclonal antibodies anti-Mmp14 (1:100 dilution, product no. ab51074; Abcam, Inc.) and anti-Reck (1:200 dilution, product no. #3433; Cell Signaling Technology) at 4 °C. Information on the antibodies and dilutions used in the IHC experiments are summarized in Additional file 1: Table S1. Antibody dilutions were made in PBS with 0.1% bovine serum albumin (BSA). The signal was amplified using the labeled streptavidin-biotin (LSAB) method with the Universal Dako LSAB + peroxidase kit (Dako North America, Inc., Via Real Carpinteria, CA, USA). Tissues were incubated for 30 min with the biotinylated link, washed twice with PBS-T, and incubated with streptavidin conjugated to horseradish peroxidase. For the IHC-negative control, the primary antibody was omitted. Tissues were revealed using the DAB kit (Dako North America) and counterstained with Harris hematoxylin (Easypath, São Paulo, SP, Brazil). Lastly, the sections were dehydrated, xylene embedded, and assembled with VectaMount permanent mounting medium (Vector Laboratories, Burlingame, CA, USA). Immunohistochemically stained sections were scanned using the Pannoramic P250 Flash II Slide Scanner (3DHISTECH, Budapest, Hungary), and the images were selected using the CaseViewer software (3DHISTECH, version 2.0). Representative areas of the tissue were selected for each condition to compose Fig. 1 and IHC scores were generated after analysis by three independent researchers (Table 2). IHC kits were provided by Dako North America.
RNA isolation, cDNA synthesis, and reverse transcriptase PCR
Total RNA was isolated using the TRIzol (Invitrogen, Carlsbad, CA, USA) protocol. RNA quality was evaluated by the spectrophotometric absorbance ratios at 280/260 nm and 230/260 nm (Nanodrop, Thermo Scientific, Waltham, MA, USA). Equal amounts of RNA treated with DNaseI (Thermo Scientific) from each of the three ovaries were mixed into a single pool for reverse transcription. Complementary DNA was synthetized from 2μg of total RNA pool using ImProm-II Reverse Transcriptase (Promega, Madison, WI, USA) according to the manufacturer’s instructions.
The relative mRNA expressions of the target genes were assessed using the Maxima SYBR Green/ROX qPCR Master Mix (Thermo Scientific) in a GeneAmp 7300 Sequence Detection System (Applied Biosystems, Foster City, CA, USA), according to the manufacturer’s instructions. The primers used to amplify the Mmps (2, 9, 13, 14, and 19), Timps (1, 2, and 3), Sparc, Bsg, and Reck genes are listed in Additional file 1: Table S2. A dissociation cycle was performed after each run to check for non-specific amplification. Relative mRNA expression levels were estimated using the method described in , which generates a normalization value in geNorm Software using β-Actin and Hprt as housekeeping genes.
Relative mRNA expression levels across treatment groups were normalized relative to the average expression of control animals treated with saline solution and compared using One-way analysis of variance (ANOVA), followed by Dunnett’s post hoc comparison test. Spearman correlations were calculated for each gene pair (Additional file 1: Table S3). All analyses were carried out using the SPSS Software (IBM SPSS Statistics; version 20). A p value of p < 0.05 was considered significant (* p < 0.05, ** p < 0.01, and *** p < 0.001).
Expression profile of Mmps during ovary ECM remodeling across hormone-induced estrous cycle
We studied the spatio-temporal patterns of expression of Mmps 2, 9, and 14 proteins within the ovarian tissue using immunohistochemistry. Figure 1a shows representative images for each condition, with the target staining intensities and respective localizations being summarized in Table 2. All Mmps analyzed were mainly restricted to the theca and stromal cells, except during the late stages of the estrous cycle (G7 and G8), when expression was also detected in the corpus luteum. Oocytes and granulosa cells were immune-reactive to gelatinases only. Ovaries showed a global increase in Mmp protein expression during the early stages of folliculogenesis (G2), followed by another expression peak during ovulation (G4). Nevertheless, the expression levels of gelatinases were lower during the ovulation-luteogenesis transition and CL formation and increased during CL maintenance, more specifically, at luteolysis.
Figure 2 shows the relative expression pattern of Mmps, at the mRNA level, by qRT-PCR. Mmp2 mRNA levels were significantly lower in ovaries from the control group, when compared to the early folliculogenesis ovaries (G2), displaying a peak at the ovulation-luteogenesis transition phase (G5) followed by a sharp decrease (Fig. 2a). Similarly, Mmp9 rapidly increased, peaked at G5 and decreased after G6 (Fig. 2b). Meanwhile, Mmp13 peaked at late stages of the G8-G9 protocol (Fig. 2c). Mmp14 (Fig. 2d) and Mmp19 were not modulated during the process (data not shown).
Expression profile of the Mmps’ modulators Reck, Bsg and Sparc during ovary ECM remodeling across hormone-induced estrous cycle
Figure 3 presents the mRNA expression profile of different Mmps modulators during the hormone-induced estrous cycle, the inhibitors Timps (Timp1, Timp2, and Timp3) (Fig. 3a, B, and C, respectively) and Reck (Fig. 3d) and the Mmps inducers Bsg and Sparc (Fig. 3e and f, respectively) in whole ovarian RNA extracts. Timp1 mRNA showed a clear expression peak during ovulation (G4 and G5, Fig. 3a), while Timp2 and Timp3 apparently are not significantly modulated in this model (Fig. 3b and c). Expression levels were normalized relative to the average expression of control animals treated with saline solution.
As may be observed in Fig. 1b and Table 2, the Reck protein, the only Mmp inhibitor which is anchored to the membrane by a glycosyl phosphatidyl anchor, was expressed in theca cells, ovarian stroma, and CL, but the protein expression pattern was opposite to that of Mmps, namely: no Reck protein expression was detected during early folliculogenesis and ovulation, while the highest Reck levels occurred immediately following eCG administration (G1) and during ovulation-luteogenesis transition (G5). The only exception was luteolysis (G9), during which Reck protein and the Mmps were both detected at high levels. Additionally, Reck mRNA was upregulated at the early stages of the estrous cycle during folliculogenesis, with an expression peak at G2 (Fig. 3d). Furthermore, its gene expression negatively correlated with Mmp13 collagenase expression (Additional file 1: Table S2). Conversely, the Mmp inducer Bsg displayed an mRNA expression peak during the late stage of CL maintenance (G8) (Fig. 3e).
The ovarian tissue expressed the Sparc protein in theca, granulosa and stromal cells, in oocytes and CL. The Sparc modulation pattern across the estrous cycle was similar to that of Mmps, with protein expression peaks during early folliculogenesis (G2), ovulation (G4), and luteolysis (G9), and lower expression levels between these phases (Fig. 1b, Table 2). High Sparc mRNA expression was observed at late folliculogenesis (G3) and from ovulation-luteogenesis transition to CL maintenance (G5 and G6, Fig. 3f).
In this study, we evaluated the spatio-temporal expression pattern of several Mmps and their regulating proteins in whole ovaries at several stages of the estrous cycle using immunohistochemistry to generate a panel of these players in rat ovarian dynamics. In addition, we also quantitatively assessed the mRNA expression of Mmps and their regulators at those same estrous cycle stages. Interestingly, some of the analyzed genes (Mmp9, Reck and Sparc) displayed a discrepancy between the observed mRNA and protein levels at specific phases of the induced estrous cycle. A poor correlation between mRNA and protein levels has been extensively reported highlighting the importance of an integrative approach [45,46,47] or the development of a gene specific RNA-to-protein conversion factor . Particularly, MMPs and their inhibitors have been shown to display discrepancies in their mRNA to protein correspondence in human prostate tissue and in skeletal muscle cancer in cachexia-associated matrix remodeling [49, 50] and further studies are necessary to identify the mechanisms that lead to this disparity for this particular targets.
Analysis of Mmp spatial and temporal protein distribution, through immunohistochemistry, revealed that the Mmp2 and Mmp9 gelatinases display an overall higher expression during the rat estrous cycle, when compared to the Mmp14 collagenase. Gelatinases had a widespread spatial distribution in the ovaries, being detected in the theca, stromal, and granulosa cells, oocytes and CL. Mmp14 protein expression, on the other hand, was mainly restricted to theca and stromal cells, but at late stages of the estrous cycle, during luteolysis (G9), being also detected at the CL. The Mmp14 transcript levels apparently did not vary across ovulation and CL lifespan, similarly to what has been described in the literature [24, 51]. Nevertheless, high Mmp14 protein levels were detected in CL during its regression, indicating a possible contribution of this MMP to local tissue remodeling.
Mmp2 mRNA levels displayed a 2.8-fold increase during early folliculogenesis (G2), followed by a reduction (G3) and displayed a peak during the ovulation-luteogenesis transition phase (G5). At this same stage (G5), 24 h after hCG injection, Mmp9 expression displayed an increase, peaking at luteinization (G6). Furthermore, both gelatinases displayed a reduced mRNA expression during CL maintenance (G7-G8) and regression (G9). Our results support previously published gene expression data from other groups, corroborating the increased gelatinolytic and collagenolytic activity during folliculogenesis and Mmp2 and Mmp9 expression 24 h post-hCG injection (luteogenesis, G5) in other studies [16, 23, 24, 51], indicating a role for this enzyme in the tissue remodeling associated with luteolysis.
We observed an increase in Mmp13 mRNA expression during the hormone-induced estrous cycle, which reached a peak of 150,000-fold 14 days after hCG injection (G9), suggesting that this gene is strongly correlated with luteal regression. These data are supported by others: Nothnick et al.  observed a 15-fold increase in Mmp13 mRNA expression 12 days after hCG and Liu et al. , using an adult pseudo-pregnant rat model, also showed an increase in Mmp13 expression during luteal regression .
Our data suggest that Mmp expression is tightly related to the main biological events taking place during the estrous cycle, supporting previously published data and adding evidence to the importance of proteolytic tissue-degradation for ovulation in primates and rodents . Equally important as determining Mmps expression during the estrous cycle is to evaluate the expression of their inhibitors and inducers, since the balance between all of these proteins is what determines the extent and nature of ECM remodeling.
Timp1 was the only member of the Timp family to show significant mRNA modulation in this work, presenting a sharp increase in mRNA levels during ovulation (G4), 12 h after hGC administration. Conversely, for an Mmp inhibitor, its expression positively correlates with Mmp2 gene expression. Li and Curry  reported that the expression of both Timp1 and Timp3 mRNAs was induced upon treatment with hCG in pre-ovulatory rat granulosa cells in an epidermal growth factor receptor tyrosine kinase and mitogen-activated protein kinase (MAPK) pathway-dependent fashion . This differential regulation is of extreme interest because Timps are multifunction proteins, whose Mmp-inhibition independent functions might have an important role in this model.
Takahashi and colleagues observed that the RECK MMP inhibitor has a privileged peri-cellular localization in human cells, due to its anchoring to the cell membrane by GPI . Although RECK gene is widely expressed in normal organs , no description is available to date about its expression pattern in ovaries, despite the intense MMP-directed ECM remodeling that occurs in these organs. Our analysis of Reck expression in rat ovarian dynamics, by immunohistochemistry, showed no detectable Reck protein expression in early folliculogenesis (G2) and ovulation (G4), while Mmps 2, 9 and 14 levels were elevated. Following the same pattern, Reck protein is expressed at detectable levels in theca and stromal cells at the end of folliculogenesis (G3) whereas expression of Mmps is lower. In addition, during the ovulation-luteogenesis transition phase (G5), an increase in Reck protein expression in theca and stromal cells was observed in parallel to a decrease in Mmp9 expression, but not of Mmp2, in the same cells. Also, at late stages of ovarian remodeling, Reck was expressed in the corpus luteum, during luteogenesis, CL maintenance and regression (G6-G9). Moreover, Reck mRNA reached a significant expression peak during the late stages of folliculogenesis (G3), 48 h post-eCG stimulus, opposite to the reduction observed in gelatinases mRNA levels at this point. This orchestrated modulation of Reck and Mmps expression suggests that Reck may be involved in Mmps modulation in this model. In support of this hypothesis, the literature describes that RECK controls MMPs 2, 9, and 14 activities through several mechanisms, including inhibition of their proteolytic activity and activation cascade, as well as their sequestration [36, 39, 54].
Regarding Mmp inducers, our gene expression analysis shows that Bsg mRNA levels peak at CL manteinance (G8). In agreement, other authors used a rat hormone-induced estrous cycle model to determine a Bsg mRNA expression panel, observing that this gene displayed increased expression during luteinization, which persisted until CL regression . Sparc, another Mmp inducer, was also investigated, revealing a spatio-temporal pattern modulation in rat ovaries during the estrous cycle. Spatial and temporal expression of Sparc protein was widespread, similarly to the pattern observed for Mmps. The Sparc protein profile was the similar to that of Mmps and opposite to the one displayed by Reck. According to Bagavandoss and collaborators, Sparc protein has been detected in theca and interstitial cells and its expression has been shown to increase in rat granulosa cells following hCG administration, during the final phases of ovarian development . However, in the current study, we also detected Sparc protein expression in theca, stromal, granulosa cells and oocytes in the whole hormone-induced estrous cycle, except at late folliculogenesis. It is important to note that this is the first time that Sparc protein was observed in oocytes and in the initial stages of the ovarian remodeling, namely folliculogenesis, ovulation and ovulation-luteinization transition in rat model. In addition, gene expression analyses indicate high levels of relative Sparc mRNA expression during late folliculogenesis (G3) and ovulation-luteogenesis transition (G5, G6). Taken together, these results suggest that the Sparc transcript and protein may also play a role in the estrous cycle. In this same direction, Smith et al., using an ovine model, suggested that SPARC gene is an important modulator of ovarian development and dynamics which is involved in tissue reconstruction during pregnancy .
Therefore, our results generate an expression profile panel of Mmps and their regulators, suggesting that Reck and Sparc seem to play an important role at specific steps of the process, being part of a larger picture in rat ovarian dynamics.
Here we present a global view of Mmps and their regulators during a hormone-induced estrous cycle in pre-pubertal rats (Fig. 4). Our results show that both Reck and Sparc are differentially expressed in the rat ovary model, suggesting that these proteins may play important roles in ovary ECM remodeling. For the first time, a spatio-temporal profile of Reck and Sparc was described, evidencing an opposite regulation pattern between these Mmps regulators, with Reck as a possible inhibitor and Sparc as an inducer of Mmps, during ovarian dynamics.
Equine chorionic gonadotropin
Human chorionic gonadotropin
0.1% Triton X-100 in Phosphate-buffered saline
Reversion-inducing-Cystein-rich protein with Kazal motifs
Secreted protein acid and rich in cysteine
Tissue inhibitors of metalloproteinases
Irving-Rodgers HF, Roger J, Luck MR, Rodgers RJ. Extracellular matrix of the corpus luteum. Semin Reprod Med. 2006;24:242–50.
Curry TE, Osteen KG. Cyclic changes in the matrix metalloproteinase system in the ovary and uterus. Biol Reprod. 2001;64:1285–96.
Curry TE, Osteen KG. The matrix metalloproteinase system: changes, regulation, and impact throughout the ovarian and uterine reproductive cycle. Endocr Rev. 2003;24:428–65.
McCawley LJ, Matrisian LM. Matrix metalloproteinases: they're not just for matrix anymore! Curr Opin Cell Biol. 2001;13:534–40.
Portela VM, Veiga A, Price CA. Regulation of MMP2 and MMP9 metalloproteinases by FSH and growth factors in bovine granulosa cells. Genet Mol Biol. 2009;32:516–20.
Kleiner DE, Stetler-Stevenson WG. Structural biochemistry and activation of matrix metalloproteases. Curr Opin Cell Biol. 1993;5:891–7.
Gomez DE, Alonso DF, Yoshiji H, Thorgeirsson UP. Tissue inhibitors of metalloproteinases: structure, regulation and biological functions. Eur J Cell Biol. 1997;74:111–22.
Brew K, Dinakarpandian D, Nagase H. Tissue inhibitors of metalloproteinases: evolution, structure and function. Biochim Biophys Acta. 2000;1477:267–83.
Gaddy-Kurten D, Hickey GJ, Fey GH, Gauldie J, Richards JS. Hormonal regulation and tissue-specific localization of alpha 2-macroglobulin in rat ovarian follicles and corpora lutea. Endocrinology. 1989;125:2985–95.
Smedts AM, Curry TE. Expression of basigin, an inducer of matrix metalloproteinases, in the rat ovary. Biol Reprod. 2005;73:80–7.
Huet C, Monget P, Pisselet C, Hennequet C, Locatelli A, Monniaux D. Chronology of events accompanying follicular atresia in hypophysectomized ewes. Changes in levels of steroidogenic enzymes, connexin 43, insulin-like growth factor II/mannose 6 phosphate receptor, extracellular matrix components, and matrix metalloproteinases. Biol Reprod. 1998;58:175–85.
Aston KE, Stamouli A, Thomas EJ, Vyas S, Iredale JP, Arthur MJ, Richardson MC. Effect of gonadotrophin on cell and matrix retention and expression of metalloproteinases and their inhibitor in cultured human granulosa cells modelling corpus luteum function. Mol Hum Reprod. 1996;2:26–30.
Stamouli A, O'Sullivan MJ, Frankel S, Thomas EJ, Richardson MC. Suppression of matrix metalloproteinase production by hCG in cultures of human luteinized granulosa cells as a model for gonadotrophin-induced luteal rescue. J Reprod Fertil. 1996;107:235–9.
Nothnick WB, Soloway P, Curry TE Jr. Assessment of the role of tissue inhibitor of metalloproteinase-1 (TIMP-1) during the periovulatory period in female mice lacking a functional TIMP-1 gene. Biol Reprod. 1997;56:1181–8.
Bagavandoss P. Differential distribution of gelatinases and tissue inhibitor of metalloproteinase-1 in the rat ovary. J Endocrinol. 1998;158:221–8.
Cooke RG, Nothnick WB, Komar C, Burns P, Curry TE. Collagenase and gelatinase messenger ribonucleic acid expression and activity during follicular development in the rat ovary. Biol Reprod. 1999;61:1309–16.
Robker RL, Russell DL, Espey LL, Lydon JP, O'Malley BW, Richards JS. Progesterone-regulated genes in the ovulation process: ADAMTS-1 and cathepsin L proteases. Proc Natl Acad Sci U S A. 2000;97:4689–94.
Liu K, Wahlberg P, Ny T. Coordinated and cell-specific regulation of membrane type matrix metalloproteinase 1 (MT1-MMP) and its substrate matrix metalloproteinase 2 (MMP-2) by physiological signals during follicular development and ovulation. Endocrinology. 1998;139:4735–8.
Hägglund AC, Ny A, Leonardsson G, Ny T. Regulation and localization of matrix metalloproteinases and tissue inhibitors of metalloproteinases in the mouse ovary during gonadotropin-induced ovulation. Endocrinology. 1999;140:4351–8.
Rothchild I. The regulation of the mammalian corpus luteum. Recent Prog Horm Res. 1981;37:183–298.
Smith MF, McIntush EW, Smith GW. Mechanisms associated with corpus luteum development. J Anim Sci. 1994;72:1857–72.
Niswender GD, Juengel JL, Silva PJ, Rollyson MK, McIntush EW. Mechanisms controlling the function and life span of the corpus luteum. Physiol Rev. 2000;80:1–29.
Nothnick WB, Keeble SC, Curry TE. Collagenase, gelatinase, and proteoglycanase messenger ribonucleic acid expression and activity during luteal development, maintenance, and regression in the pseudopregnant rat ovary. Biol Reprod. 1996;54:616–24.
Liu K, Olofsson JI, Wahlberg P, Ny T. Distinct expression of gelatinase a [matrix metalloproteinase (MMP)-2], collagenase-3 (MMP-13), membrane type MMP 1 (MMP-14), and tissue inhibitor of MMPs type 1 mediated by physiological signals during formation and regression of the rat corpus luteum. Endocrinology. 1999;140:5330–8.
Reich R, Daphna-Iken D, Chun SY, Popliker M, Slager R, Adelmann-Grill BC, Tsafriri A. Preovulatory changes in ovarian expression of collagenases and tissue metalloproteinase inhibitor messenger ribonucleic acid: role of eicosanoids. Endocrinology. 1991;129:1869–75.
Mann JS, Kindy MS, Edwards DR, Curry TE. Hormonal regulation of matrix metalloproteinase inhibitors in rat granulosa cells and ovaries. Endocrinology. 1991;128:1825–32.
Murray SC, Keeble SC, Muse KN, Curry TE. Regulation of granulosa cell-derived ovarian metalloproteinase inhibitor(s) by prolactin. J Reprod Fertil. 1996;107:103–8.
Kennedy JL, Muse KN, Keeble SC, Curry TE. Hormonal regulation of tissue inhibitors of metalloproteinases during follicular development in the rat ovary. Endocrine. 1996;5:299–305.
Nothnick WB, Edwards DR, Leco KJ, Curry TE. Expression and activity of ovarian tissue inhibitors of metalloproteinases during pseudopregnancy in the rat. Biol Reprod. 1995;53:684–91.
Duncan WC, McNeilly AS, Illingworth PJ. The effect of luteal "rescue" on the expression and localization of matrix metalloproteinases and their tissue inhibitors in the human corpus luteum. J Clin Endocrinol Metab. 1998;83:2470–8.
Smedts AM, Lele SM, Modesitt SC, Curry TE. Expression of an extracellular matrix metalloproteinase inducer (basigin) in the human ovary and ovarian endometriosis. Fertil Steril. 2006;86:535–42.
Bagavandoss P, Sage EH, Vernon RB. Secreted protein, acidic and rich in cysteine (SPARC) and thrombospondin in the developing follicle and corpus luteum of the rat. J Histochem Cytochem. 1998;46:1043–9.
Nuttall RK, Sampieri CL, Pennington CJ, Gill SE, Schultz GA, Edwards DR. Expression analysis of the entire MMP and TIMP gene families during mouse tissue development. FEBS Lett. 2004;563:129–34.
Podhajcer OL, Benedetti LG, Girotti MR, Prada F, Salvatierra E, Llera AS. The role of the matricellular protein SPARC in the dynamic interaction between the tumor and the host. Cancer Metastasis Rev. 2008;27:691–705.
Smith GW, Gentry PC, Long DK, Smith MF, Roberts RM. Differential gene expression within the ovine corpus luteum: identification of secreted protein acidic and rich in cysteine as a major secretory product of small but not large luteal cells. Endocrinology. 1996;137:755–62.
Takahashi C, Sheng Z, Horan TP, Kitayama H, Maki M, Hitomi K, Kitaura Y, Takai S, Sasahara RM, Horimoto A, et al. Regulation of matrix metalloproteinase-9 and inhibition of tumor invasion by the membrane-anchored glycoprotein RECK. Proc Natl Acad Sci U S A. 1998;95:13221–6.
Bornstein P, Sage EH. Matricellular proteins: extracellular modulators of cell function. Curr Opin Cell Biol. 2002;14:608–16.
Yiu GK, Chan WY, Ng SW, Chan PS, Cheung KK, Berkowitz RS, Mok SC. SPARC (secreted protein acidic and rich in cysteine) induces apoptosis in ovarian cancer cells. Am J Pathol. 2001;159:609–22.
Oh J, Takahashi R, Kondo S, Mizoguchi A, Adachi E, Sasahara RM, Nishimura S, Imamura Y, Kitayama H, Alexander DB, et al. The membrane-anchored MMP inhibitor RECK is a key regulator of extracellular matrix integrity and angiogenesis. Cell. 2001;107:789–800.
Zhang X, Healy C, Nothnick WB. Estrogen suppresses expression of the matrix metalloproteinase inhibitor reversion-inducing cysteine-rich protein with Kazal motifs (RECK) within the mouse uterus. Endocrine. 2012;42:97–106.
Zarrow MX, Quinn DL. Superovulation in the immature rat following treatment with PMS alone and inhibition of PMS-induced ovulation. J Endocrinol. 1963;26:181–8.
Jo M, Curry TE. Regulation of matrix metalloproteinase-19 messenger RNA expression in the rat ovary. Biol Reprod. 2004;71:1796–806.
Ying SY, Meyer RK. Ovulation induced by PMS and HCG in hypophysectomized immature rats. Proc Soc Exp Biol Med. 1972;139:1231–3.
Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29:e45.
Tian Q, Stepaniants SB, Mao M, Weng L, Feetham MC, Doyle MJ, Yi EC, Dai H, Thorsson V, Eng J, et al. Integrated genomic and proteomic analyses of gene expression in mammalian cells. Mol Cell Proteomics. 2004;3:960–9.
Maier T, Guell M, Serrano L. Correlation of mRNA and protein in complex biological samples. FEBS Lett. 2009;583:3966–73.
Liu Y, Beyer A, Aebersold R. On the dependency of cellular protein levels on mRNA abundance. Cell. 2016;165:535–50.
Edfors F, Danielsson F, Hallstrom BM, Kall L, Lundberg E, Ponten F, Forsstrom B, Uhlen M. Gene-specific correlation of RNA and protein levels in human cells and tissues. Mol Syst Biol. 2016;12:883.
Lichtinghagen R, Musholt PB, Lein M, Romer A, Rudolph B, Kristiansen G, Hauptmann S, Schnorr D, Loening SA, Jung K. Different mRNA and protein expression of matrix metalloproteinases 2 and 9 and tissue inhibitor of metalloproteinases 1 in benign and malignant prostate tissue. Eur Urol. 2002;42:398–406.
Devine RD, Bicer S, Reiser PJ, Velten M, Wold LE. Metalloproteinase expression is altered in cardiac and skeletal muscle in cancer cachexia. Am J Physiol Heart Circ Physiol. 2015;309:H685–91.
Jo M, Thomas LE, Wheeler SE, Curry TE. Membrane type 1-matrix metalloproteinase (MMP)-associated MMP-2 activation increases in the rat ovary in response to an ovulatory dose of human chorionic gonadotropin. Biol Reprod. 2004;70:1024–32.
Chaffin CL, Vandevoort CA. Follicle growth, ovulation, and luteal formation in primates and rodents: a comparative perspective. Exp Biol Med (Maywood). 2013;238:539–48.
Li F, Curry TE. Regulation and function of tissue inhibitor of metalloproteinase (TIMP) 1 and TIMP3 in periovulatory rat granulosa cells. Endocrinology. 2009;150:3903–12.
Welm B, Mott J, Werb Z. Developmental biology: vasculogenesis is a wreck without RECK. Curr Biol. 2002;12:R209–11.
We are grateful to Zizi de Mendonça, Debora Cristina da Costa, Ingrid da Silva Cunha, Raquel Santana and Amanda Cordeiro de Mello Pereira for their excellent technical support, to MSc. Marluce da Cunha Mantovani, MSc. Mateus Contar Adolfi, and Dr. Gustavo Gross Belchior for experimental support and advice, and to Carlos DeOcesano Pereira for help with the figures.
The following Brazilian Federal and State research agencies provided financial support for this study: São Paulo Research Foundation (FAPESP) [grant number 2013/18842–8], National Bank for Economical and Social Development (BNDES), National Council for Scientific and Technological Development (CNPq), Coordination for the Improvement of Higher Education Personnel (CAPES), National Financier of Research and Projects (FINEP), Ministry of Science, Technology and Innovation (MCTI), and Ministry of Health of Brazil/Department of Science and Technology (MS-DECIT).
Availability of data and materials
The dataset supporting the conclusions of this article is included within the article (and its Additional files).
All animal experimentation was approved by the Ethics Committee for Animal Use of the Chemistry Institute (CEUA), University of São Paulo, São Paulo, Brazil, in accordance with the National Council for the Control of Animal Experimentation (CONCEA).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Primary antibody characteristics and concentrations used in IHC assays. Table S2. Sequence and concentration of primers used in qRT-PCR assays. Table S3. Spearman correlation coefficient between the expression profiles of all analyzed genes. (DOC 90 kb)