- Open Access
The effect of progesterone replacement on gene expression in the corpus luteum during induced regression and late luteal phase in the bonnet monkey (Macaca radiata)
Reproductive Biology and Endocrinologyvolume 9, Article number: 20 (2011)
In higher primates, although LH/CG play a critical role in the control of corpus luteum (CL) function, the direct effects of progesterone (P4) in the maintenance of CL structure and function are unclear. Several experiments were conducted in the bonnet monkey to examine direct effects of P4 on gene expression changes in the CL, during induced luteolysis and the late luteal phase of natural cycles.
To identify differentially expressed genes encoding PR, PR binding factors, cofactors and PR downstream signaling target genes, the genome-wide analysis data generated in CL of monkeys after LH/P4 depletion and LH replacement were mined and validated by real-time RT-PCR analysis. Initially, expression of these P4 related genes were determined in CL during different stages of luteal phase. The recently reported model system of induced luteolysis, yet capable of responsive to tropic support, afforded an ideal situation to examine direct effects of P4 on structure and function of CL. For this purpose, P4 was infused via ALZET pumps into monkeys 24 h after LH/P4 depletion to maintain mid luteal phase circulating P4 concentration (P4 replacement). In another experiment, exogenous P4 was supplemented during late luteal phase to mimic early pregnancy.
Based on the published microarray data, 45 genes were identified to be commonly regulated by LH and P4. From these 19 genes belonging to PR signaling were selected to determine their expression in LH/P4 depletion and P4 replacement experiments. These 19 genes when analyzed revealed 8 genes to be directly responsive to P4, whereas the other genes to be regulated by both LH and P4. Progesterone supplementation for 24 h during the late luteal phase also showed changes in expression of 17 out of 19 genes examined.
These results taken together suggest that P4 regulates, directly or indirectly, expression of a number of genes involved in the CL structure and function.
In mammals, the secretion of progesterone (P4) by corpus luteum (CL) is absolutely essential for establishment and, in some species, maintenance of pregnancy. In higher primates, LH and chorionic gonadotropin (CG) have been suggested to be the principal trophic factors responsible for P4 secretion in the CL . Whether P4 plays a role in the maintenance of structure and function of CL has not been fully elucidated in higher primates. Rothchild postulated that P4 is the primary stimulus of its own secretion and that intraluteal P4, among other effects such as control of structural integrity and steroidogenic capacity, is responsible for regulation of production of luteolysin, the prostaglandin (PG) F2α, within the CL[2, 3]. More recent studies have provided several lines of evidence, some of them with mechanistic insights, in support of the direct effects of P4 on CL. Expression of progesterone receptor (PR) isoforms in CL have been reported in several mammalian species [4–6]. Several studies have suggested that by way of its proliferative and anti-apoptotic actions, P4 functions as survival factor of CL in human , rat [8–10], and cattle [11, 12]. Like other steroid nuclear receptors, PR utilizes a plethora of cofactors termed coactivators or corepressors to regulate gene expression . The PR cofactors identified to date include coactivators like (SRC1-3, CBP/p300 and NCOA1-3)  and corepressors like NCOR1-2 involved in modulation of PR activity in vivo [14, 15].
We have recently standardized a GnRH R antagonist-induced luteolysis model system in the monkey in which induced luteolysis could be reversed by exogenous LH administration . Employing this model system, microarray analysis of differentially expressed genes in CL tissue during induced luteolysis and LH replacement following induced luteolysis has been determined . The GnRH R antagonist-induced regressed CL with ablated LH action, yet capable of responding to LH replacement, affords an ideal situation to examine direct effects of P4 on CL structure and function. Extensive tissue remodeling with breakdown and renewal of extracellular matrix (ECM) that occurs during spontaneous luteolysis requires participation of matrix metalloproteinases (MMPs) and their tissue inhibitors, TIMPs [17–19]. It has been reported that decrease in P4 levels during late luteal phase of the non fertile cycle is associated with changes in expression of ECM regulators . It remains to be determined whether P4 regulates expression of tissue proteinases, especially following conception. The purpose of this study was to examine effects of P4 action on gene expression changes and function of CL. Experiments were carried out to determine expression of genes in CL tissue that encode different elements of PR complex and few downstream targets of PR activation throughout the luteal phase, after LH/P4 depletion, after P4 replacement and after P4 supplementation during the late luteal phase.
Oligonucleotide primers were synthesized by Sigma-Genosys, Bangalore, India. DyNAzyme™ II DNA polymerase (F-501L) was purchased from Finnzymes, Espoo, Finland. Moloney murine leukemia virus (MMuLV) reverse transcriptase (RT) and 100 bp DNA ladder were obtained from MBI Fermentas GmbH (St. Leon-Rot, Germany). Power SYBR® Green PCR master mix was obtained from Applied Biosystems, Foster City, CA, USA. GnRH R antagonist [Cetrorelix®; (CET)] was a kind gift from Asta Medica, Frankfurt, Germany. ALZET® Osmotic pump Model 2ML1 (infusion rate 10 μl/h) was obtained from Alza Corporation, Palo Alto, CA. Antibodies specific to phospho-p38 (9211), phospho-p42/44 MAPK (9101), p38 MAPK (9212), p42/44 MAPK (9102), pMKK3/6 (9231), p38 MAPK assay kit (9820), MMP-9 (G657) and Phototype-HRP Western detection system with horseradish peroxidase-linked anti-rabbit IgG (7071) were purchased from Cell Signaling Technology, Inc. Danvers, MA. Antibodies specific to PR (sc-538), NCOA1 (sc-8995), NCOA2 (also designated as GRIP1; sc-8996), NCOA3 (sc-25742) were procured from Santa Cruz Biotechnology, Santa Cruz, CA. Crystalline P4 (P0130) and all other reagents were purchased from Sigma Aldrich Corp. St. Louis, MO or sourced locally.
Animal protocols, blood samples and CL collection
Experimental protocols involving monkeys in this study were approved by the Institutional Animal Ethics Committee of the Indian Institute of Science, Bangalore. Adult female bonnet monkeys (Macaca radiata) weighing 3.3-5.1 kg were utilized for the study. The general care and housing of monkeys have been described elsewhere . In this study, one day after occurrence of peak E2 surge was designated as day 1 of the luteal phase, and CL was collected on designated days of the luteal phase and/or after administration of different treatments (see below). To retrieve CL from experimental monkeys, ovaries were accessed by performing laparotomy under aseptic conditions on ketamine hydrochloride (15 mg/kg BW) and/or pentobarbital sodium (8-12 mg/kg BW) anesthetized monkeys. Under sterile conditions, the excised CL was transferred to a petri dish containing filter paper, wiped dry, weighed, cut into 4-5 pieces, placed in individual sterile cryovials, snap-frozen in liquid nitrogen and stored at -70°C until analysis.
Experiment 1: Examination of differentially expressed genes during induced luteolysis and LH replacement experiments
In order to identify genes associated with LH and P4 receptor signaling, the microarray data of induced luteolysis (LH/P4 depletion) and rescue of CL function by exogenous LH replacement studies deposited (#GSE7827 and #GSE8371) by us previously were examined. The main focus was on identifying changes in expression of genes associated with PR and genes downstream of PR signaling. For this purpose, various cofactors and corepressors (NCOA1-3, NCOR1-2, FKBP4 and FKBP5) of PR were selected for analysis. Expression of MMP-2 and 9, TIMP-3, ADAMTS-1, PGRMC1, BMP5, CDH2, WNT7A, HOXA1, IHH genes regarded as downstream target of gonadotropin and/or PR signaling in CL and in traditional P4 target tissues such as uterus were also examined [19, 22–26]. The expanded forms of gene symbols used in this study are provided in Table 1. Real time RT-PCR analysis was carried out on RNA samples isolated from CL tissues collected during luteolysis [VEH (n = 5 animals) and LH/P4 depletion (n = 5 animals)] and LH rescue (LH replacement, n = 3 animals each for PBS and LH treatments) of CL function experiments as reported previously .
Experiment 2: Expression of various components of PR signaling complex and PR target genes throughout the luteal phase
To study expression of PR, various components of PR signaling complex and some of the downstream target genes of PR signaling, corpora lutea (n = 3 animals/stage of luteal phase) were collected from monkeys at early (day 5), mid (day 8), and late (day 14) stage of the luteal phase of the menstrual cycle as reported previously . Also, CL (n = 3 animals) was collected from monkeys on day 1 of menstrual cycle (d1M), a time point when luteolytic events are manifested. Blood samples were collected from monkeys on the day of CL collection for determination of P4 concentration.
Experiment 3: Effects of P4 replacement on expression of PR, coactivators, corepressors and P4target genes during induced luteolysis
We have recently reported that replacement of LH post LH/P4 depletion leads to brisk and sustained increase in P4 concentration suggesting rescue of CL function following reestablishment of LH levels . In the present experiments, we determined the direct effects of P4 on CL during absence of luteotrophic stimulus. The different components of the PR signaling complex and expression of few genes considered as target of PR activation were examined. In experiment 3.1, with a view to mimic mid luteal phase circulating P4 concentration, exogenous P4 was administered through implantation of P4 filled ALZET pumps to monkeys depleted of LH/P4 (P4 replacement model). For this purpose, pilot experiments were carried out in adult female monkeys to examine the feasibility of providing P4 as continuous infusion by way of implantation of P4 filled ALZET pump. Keeping in mind the limitation on use of organic solvents compatible with the pump, initially 16 mg of P4 was dissolved in 300 μl of ethanol and the solution was made up to 2 ml by propylene glycol and the entire solution was transferred to 2 ML1 ALZET pump. It was determined that three, 32.5 mg of P4 filled pumps that provided infusion of 487.5 μg/30 μl/h, were required to be implanted in anesthetized monkeys to establish circulating P4 concentration higher or in the range of mid luteal phase concentration. In order to determine the time course of P4 secretion immediately before and at different time intervals after injection of CET (LH/P4 depletion), monkeys on day 7 of the luteal phase (n = 3 animals) were administered CET (150 μg/kg BW, s.c.) and blood samples were collected twice daily until onset of menses. In experiment 3.2, three groups of monkeys (n = 3-4 animals/group) were administered 5.25% glucose (vehicle for CET treatment; VEH; group 1), CET [150 μg/kg BW; (LH/P4 depletion); CET; group 2] and CET (150 μg/kg BW) followed 24 h later with implantation of P4 filled three ALZET pumps designed to infuse 487.5 μg of P4/30 μl/h for 24 h [CET+P4; (P4 replacement); group 3]. Blood samples were collected immediately before and at different time intervals throughout the experiment. In group 3, ALZET pumps were removed 24 h later, and further blood samples were collected to monitor circulating P4 levels. Corpora lutea were retrieved from monkeys of all three groups (VEH, CET and CET+P4).
Experiment 4: Effects of P4supplementation on CL function during late luteal phase
An experiment was conducted to gain insight into direct effects of increased P4 during rescue of CL function that occurs during late luteal phase of the fertile menstrual cycle. On day 13 of the luteal phase of non-mated females, three P4 filled ALZET pumps (97.5 mg of P4) were implanted for 24 h, CL (n = 3 animals) was harvested and the pumps were removed 24 h later. Blood samples were collected immediately before and at different time intervals during implantation and after retrieval of CL as well as removal of implants for determining the P4 secretion pattern. For purposes of comparison, CL (n = 3 animals) was harvested from untreated monkeys (untreated control) on day 14 of the luteal phase of the menstrual cycle.
Total RNA was isolated from CL tissues obtained from different experiments using TRI reagent according to manufacturer's instructions. RNA samples were analyzed using NanoDrop ND-1000 UV-VIS spectrophotometer and samples with A260/A280 values ~1.8 -1.9 were selected for further analyses.
Real-time RT-PCR analysis
Real time RT-PCR analysis was carried out as described previously . Total RNA (1 μg) extracted from CL of monkeys from various experiments were reverse transcribed using MMuLV RT in Eppendorf mastercycler, epgradient PCR machine. Briefly, In a 20 μl total reaction, 1 μg of total RNA along with 1 μl of Oligo dT primer were incubated at 65°C for 10 min and then snap chilled on ice for 5 min. To this RNA-Oligo dT hybrid, a cocktail containing 4 μl of 5X RT buffer [250 mM Tris HCl (pH 8.3 at 25°C), 250 mM KCl, 20 mM MgCl2 and 10 mM DTT], 10 mM dNTPs and DEPC treated water were added followed by addition of 200 units (1 μl) of Revert Aid™MMuLV RT enzyme. The reaction mixture was incubated at 42°C for 60 min. The resulting cDNAs were used as template for real time RT-PCR analysis using gene specific primers. The primer pairs were designed using Rhesus macaque sequences submitted at NCBI and ENSEMBL using Primer Express software v2.0 (Applied Biosystems, Foster City, CA, USA) spanning the exon-exon junction. The diluted cDNA samples equivalent to 10 ng of total RNA were subjected to analysis in Applied Biosystems 7500 Fast Real Time PCR system with SDS v 1.4 program employing Power SYBR green 2X PCR master mix. The 10 μl real time RT-PCR mixture contained cDNA equivalent to 10 ng of total RNA, 5 μl of PCR master mix and 5 μM each of forward and reverse gene specific primers. The initial denaturation was carried out at 95°C for 10 min, with further 40 cycles of denaturation (95°C for 30 sec), annealing (specific annealing temperature for 30 sec) and extension (72°C for 30 sec) with a final extension of 5 min at 72°C. PCR reactions were carried out in duplicates in 96 well plates. For each gene a no template control (NTC) was included and dissociation/melting curves were generated to determine the specificity of primers. Real time RT-PCR data for each gene was normalized using L19 expression level as internal control within each cDNA sample. The fold change in expression of the genes was determined using the ΔΔCt method, which calculates the fold change using the formula: Fold change = 2-ΔΔCt, where Ct= Threshold cycle i.e. the cycle number at which the relative fluorescence of test samples increases above the background fluorescence and ΔΔCt= [Ct gene of interest (treated sample)-Ct of L19 (treated sample)] - [Ct gene of interest (control sample)-Ct of L19 (control sample)]. The list of genes and details of the primers employed in the real time RT-PCR analysis along with the annealing temperature and expected product size are provided in Table 2.
E2 and P4 concentrations in serum were determined by specific RIA as reported previously .
Immunoblot analyses and in vitro P38MAPK (P38) assays of CL tissue lysates were carried out as per published procedures. Assays were carried out with equal amount of protein lysates from the luteal tissue of various treatments as reported previously from the laboratory .
Data were expressed as mean ± SEM. Statistical evaluation of mean differences of serum P4 concentrations and real-time RT-PCR expression among different experimental groups were analyzed by one-way ANOVA, followed by the Newman-Keuls multiple comparison tests (PRISM GraphPad, version 4.0; GraphPad Software Inc., San Diego, CA) and Student's t- test to compare between two groups. A P value of <0.05 was considered statistically significant. The linear regression analysis on microarray and real time RT-PCR data was done as reported earlier .
Identification and validation of differentially expressed genes in LH/P4depletion and LH replacement experiments
The experimental details of CL collection from LH/P4 depletion (CET-induced luteolysis) and LH replacement (CET+LH) models and the microarray data comparing these models deposited in NCBI's Gene Expression Ominibus (#GSE7827 and #GSE8371) have been described previously . The published microarray data was mined for P4← responsive genes, and 45 differentially expressed genes were identified based on criteria of >1.5 fold change in both LH/P4 depletion and LH replacement experiments. We narrowed down to the 45 genes by comparing our previous data on LH replacement studies with P4 -responsive genes reported in literature for other tissues [28, 29]. From the 45 genes identified, 19 genes belonging to steroidogenesis and P4 target/regulation related genes (PR, PR binding proteins, cofactors, corepressors and some of the genes considered as downstream targets of PR activation) considered necessary for CL structure and function were selected for further studies. It should be noted that of the 19 genes, expression of only 12 genes was examined during different stages of the luteal phase in experiment 2 and expression of all 19 genes was determined in experiments 3 and 4.
Figure 1A shows microarray and real time RT-PCR fold changes in expression of genes associated with PR signaling/regulation in CL tissue from monkeys of LH/P4← depletion and LH replacement models. Validation of microarray data from LH/P4← depletion model by real time RT-PCR analyses provided a good correlation in expression changes of most genes examined except for PR, NCOR1, and HOXA1 which showed higher expression changes in real time RT-PCR analysis, while expression of WNT7A and AREG genes showed inverse correlation by real time RT-PCR analysis (Figure 1A). The linear regression analysis on log2-transformed data of microarray and real time RT-PCR analyses suggested a good correlation (R = 0.5013, P = 0.0484). A number of genes that were examined to be differentially expressed after LH/P4← depletion were found to be reversed following LH replacement (Figure 1B). Expression of P4← responsive genes like PR, NCOA3, TIMP3, BMP5 and CDH2 was down regulated after LH replacement, whereas expression of genes like NCOR1, NCOR2 and MMP-9 was up regulated after LH replacement (Figure 1B). The microarray data and real time RT-PCR data for LH replacement experiments showed good correlation except for the observation of profound down regulation in BMP5 expression by microarray analysis than real time RT-PCR analysis (R = 0.7028, P = 0.0312).
Expression of PR cofactors and its downstream target genes throughout the luteal phase
Real time RT-PCR analysis results of PR, NCOA1-3, NCOR1-2, FKBP5, TIMP 3, MMP-2, ADAMTS-1, IHH and HOXA1 in CL collected from different stages of luteal phase are represented in Figure 2A. In order to gain information on the fold expression changes during different stages, expression of each gene at early (E) stage was set as 1 and expression at other stages of luteal phase was represented as fold change in relation to the early stage. Corresponding to the waxing and waning levels of circulating P4 throughout different stages of the luteal phase, the expression of PR mRNA was not statistically different (P > 0.05) across all the stages of luteal phase (Figure 2A). The protein expression for PR isoform B was low at early stage and its expression was high at all the other stages examined, with higher expression at mid luteal phase compared to early stage (1 vs 1.82 ± 0.2 fold; early vs mid stage, respectively; P < 0.05) (Figure 2B). However, the protein expression for PR isoform A was low at early and mid stage, and became higher at later stages of the luteal phase with late luteal stage being significantly higher (1 vs 1.45 ± 0.06 fold, early vs late stage, respectively; P < 0.05). The fold change in expression of NCOA1 mRNA and protein was low at early and mid stage (Figure 2A), and was significantly higher (P < 0.05) in CL collected from monkeys on day 1 of menses (1 vs 3.62 ± 0.04 fold, early stage vs d1 M, respectively; P < 0.05). Expression of NCOA2 and NCOR2 mRNA was higher (P < 0.05) in the late stage CL compared to their expression at early, mid and d1 M stages (Figure 2A), and the protein expression for NCOA2 was high in mid stage CL. Expression changes of FKBP5 and NCOA3 mRNA did not vary significantly during different stages of luteal phase (P > 0.05; Figure 2A). The protein expression of NCOA3 was high in mid stage CL compared to early stage CL (1 vs 1.2 ± 0.05 fold, early vs mid stage, respectively; P < 0.05, Figure 2B). The fold change in expression of NCOR1 showed a progressive increase with expression being higher in CL of late stage and d1 M compared to other stages (Figure 2A). The expression of proteinase family genes such as TIMP3 and MMP-2 were higher in CL at late luteal stage compared to early and mid luteal stages. Expression of ADAMTS-1, IHH and HOXA1 in CL tissues collected at different stages of luteal phase did not vary significantly (P > 0.05; Figure 2A).
Effects of P4replacement on expression of PR, coactivators and corepressors during induced luteolysis
Circulating P4 concentrations were 3.8 ± 0.3, 1.5 ± 0.07, 0.6 ± 0.1, 0.56 ± 0.01, 0.5 ± 0.04, 0.3 ± 0.05 and 0.6 ± 0.1 ng/ml, respectively at 0, 12, 24, 36, 48, 60 and 72 h post LH/P4 depletion (CET treatment; Figure 3A). All monkeys showed menses 96-120 h post LH/P4 depletion. A brief outline of the experiment involving exogenous P4 replacement is presented in Figure 3B. The P4 concentration in monkeys combined from all three groups was 3.5 ± 0.6 ng/ml prior to VEH (n = 3) or CET (n = 7) treatment (Figure 3C), and were maintained higher (4.1 ± 0.3 and 4.98 ± 0.3 ng/ml) at 12 and 24 h after administration of VEH, whereas P4 concentrations declined significantly 12 and 24 h post LH/P4 depletion (1.5 ± 0.07 and 0.6 ± 0.1 ng/ml at 12 and 24 h respectively; P < 0.05, Figure 3C). Replacement of P4 by way of implantation of P4 containing ALZET pumps at 24 h post CET treatment resulted in significant (P < 0.05) increase in P4 concentrations within 6 h post treatment and remained high throughout the 24 h duration of implantation. Upon removal of P4 containing ALZET pumps and lutectomy, P4 concentrations declined to reach 1.38 ± 0.6 and 0.76 ± 0.17 ng/ml at 12 and 24 h, respectively (Figure 3C).
To analyze gene expression changes post P4 replacement, 19 genes were selected for real time RT-PCR analyses and the results were presented in Figure 4A. These genes were selected for analysis based on literature data from studies on classic P4 target organs, P4 responsive cell lines as well as from our microarray data of the LH/P4 depletion (CET induced luteolysis) experiments [16, 28, 29]. Initially the fold change in expression of PR mRNA was low in VEH treated monkeys, increased to 15 fold in LH/P4 depletion model, and then decreased significantly (P < 0.05) after P4 replacement (Figure 4A) suggesting PR to be a P4 target gene. Similar to pattern of PR mRNA expression, the protein expression of PR isoform B was low initially in CL of VEH treated monkeys and increased post LH/P4 depletion (1 vs 2.72 ± 0.04 fold VEH vs CET, respectively; P < 0.05; Figure 4B), while reduced post P4 replacement. However, the protein expression of PR isoform A was low in CL of VEH treated monkeys and remained low thereafter in both LH/P4 depletion and P4 replacement conditions (Figure 4B). The fold change in expression of PGRMC1 and FKBP4 mRNA was significantly lower (P < 0.05) after LH/P4 depletion and no significant effect was observed after P4 replacement (Figure 4A). Expression of NCOA1 mRNA decreased marginally after LH/P4 depletion while P4 replacement increased its expression significantly by 3-4 fold (P < 0.05; Figure 4A) suggesting NCOA1 mRNA to be regulated by P4 directly. However, protein expression for NCOA1 increased significantly after LH/P4 depletion (1 vs 2.38 ± 0.04 fold, VEH vs CET treatment; P < 0.05) (Figure 4B), but no significant change was observed following P4 replacement (P > 0.05). Messenger levels of NCOA2 and FKBP5 in CL tissue did not significantly change after LH/P4 depletion or P4 replacement (P > 0.05; Figure 4A). On the other hand, the protein expression for NCOA2 was significantly (P < 0.05) higher after LH/P4 depletion (1 vs 1.87 ± 0.05 fold, VEH vs CET treatment; P < 0.05). However, P4 replacement resulted in significant decrease in protein levels (1 vs 0.38 ± 0.10 fold, VEH vs CET+P4 treatment; P < 0.05) compared to CL obtained from VEH or CET treated monkeys indicating indirect/direct regulation of NCOA2 by P4 at the protein level. Messenger levels of NCOA3 in CL tissue did not change significantly after LH/P4 depletion or P4 replacement (P > 0.05; Figure 4A). However, the protein expression for NCOA3 was low in CL from VEH treated monkeys and increased significantly after LH/P4 depletion (1 vs 3.10 ± 0.50 fold, VEH vs CET treatment; P < 0.05) while decreased after P4 replacement (Figure 4B) suggesting a role for P4 in the regulation of NCOA3 protein levels. The fold change in expression of NCOR1 mRNA increased significantly 5-8 fold after LH/P4 depletion and became lower following P4 replacement (P < 0.05; Figure 4A). Similarly expression of NCOR2 mRNA increased (but not statistically significant) after LH/P4 depletion, but had no effect following P4 replacement (Figure 4A). Expression of WNT7A mRNA increased significantly (P < 0. 05; Figure 4A) after LH/P4 depletion, and the expression decreased significantly (P < 0.05) after P4 replacement demonstrating regulation of WNT7A mRNA by P4. Expression of AREG mRNA decreased significantly (P < 0.05) after LH/P4 depletion and remained low post P4 replacement (Figure 4A). The fold change in expression of BMP5 increased 3 fold after LH/P4 depletion, but the expression decreased significantly (P < 0.05) upon P4 replacement (Figure 4A) suggesting BMP5 gene to be a P4 responsive gene. Though the fold change in expression of cadherin (CDH2) tended to be higher (but not statistically significant) after LH/P4 depletion showed no significant change after P4 replacement (Figure 4A). Expression of HOXA1 and IHH mRNA were increased significantly (P < 0.05) after LH/P4 depletion, while expression of HOXA1 mRNA decreased after P4 replacement, expression of IHH mRNA did not change significantly after P4 replacement (Figure 4A).
The fold expression change of ADAMTS-1 mRNA after LH/P4 depletion was lower compared to VEH treatment (P < 0.05) and the expression increased significantly (P < 0.05) following P4 replacement (Figure 4A). Fold change in expression of MMP-9 mRNA increased post LH/P4 depletion, and decreased significantly post P4 replacement (P < 0.05; Figure 4A) indicating indirect/direct regulation of both ADAMTS-1 and MMP-9 mRNAs by P4 at the transcription level. Further, the fold change in expression of other genes like MMP-2 and TIMP3 increased significantly following LH/P4 depletion, but P4 replacement had no effect (P < 0.05; Figure 4A).
Analysis of luteal tissue lysates for changes in P38 activation indicated that P38 activity was lower in CL after LH/P4 depletion, however, P4 replacement increased P38 activation significantly (P < 0.05; Figure 5). In vitro MAPK assays revealed lower levels of pATF-2, a substrate for activated P38 kinase, in CL after LH/P4 depletion, whereas, the effect of LH/P4 depletion on pATF-2 levels was prevented by P4 replacement (Figure 5). The kinetics of P38 activation was correlated with activation of its upstream kinase, pMKK3/6 protein levels. LH/P4 depletion decreased the pMKK3/6 protein expression significantly (P < 0.05), while P4 replacement increased its levels significantly (P < 0.05; Figure 5) indicating other components of MAP kinase signaling to be active post P4 replacement. Expression of other MAP kinase like ERK was also determined. Phosphorylated ERK levels did not change significantly in CL tissues both after LH/P4 depletion and P4 replacement (P > 0.05; Figure 5).
Effects of P4supplementation on gene expression changes during the late luteal phase
Circulating P4 concentrations in monkeys before and after implantation of P4 containing ALZET pumps and in untreated control monkeys are shown in Figure 6A. Implantation of P4 pumps resulted in significant (P < 0.05) increase in concentrations of P4 i.e., 4.4 ± 0.64, 4.2 ± 0.37, 4.4 ± 0.35 ng/ml at 6, 12, 24 h post P4 supplementation, respectively and the corresponding P4 concentrations from untreated control monkeys were 0.70 ± 0.1, 0.6 ± 0.1 and 0.49 ± 0.03 ng/ml at 6, 12 and 24 h, respectively (Figure 6A). Real time RT-PCR analysis showed >2 fold up regulation of PR, NCOA1, NCOA3, NCOR2, FKBP4, FKBP5, CDH2, HOXA1 and AREG expression in P4 supplemented monkeys, while expression of BMP5, MMP-2, PGRMC1, and ADAMTS-1 was significantly down regulated (P < 0.05; Figure 6B). Immunoblot analysis for PR (both isoforms A&B) protein showed no significant change after P4 supplementation, whereas, immunoblot analysis for NCOA1 and NCOA2 proteins showed significantly lower (P < 0.05) expression after P4 supplementation (NCOA1, 1 vs 0.25 ± 0.2 fold, for NCOA2; 1 vs 0.62 ± 0.2 fold; Figure 6C). Immunoblot analysis for protein expression of NCOA3 and MMP-9 showed lower levels after P4 supplementation (not statistically significant P > 0.05; Figure 6C).
Together with earlier reports for presence of PR in CL of many other species, the results of expression of PR, various binding proteins, cofactors and P4 regulation of expression of few genes in the present study further confirm P4 actions within the CL. The importance of P4 in the regulation of CL structure and function has been demonstrated in rats. It was shown that P4 administration during postpartum period had an anti apoptotic action in the CL . Although the rat CL does not express nuclear PR, it was suggested that P4 mediates its action by binding to its membrane receptors . In recent years, the membrane receptor-mediated actions of P4 appear to have emerged as an important mechanism for activation of P4 signal transduction pathway. It was reported that rat luteal cells express membrane and progestin binding proteins, progestin membrane receptor (PMR) α, PMR β, PMR γ, membrane component 1 (PRMC1/PGMRC1) and Rda 288 . In a recent study, it was demonstrated that PGMRC1-dependent mechanism appears to promote human granulose/luteal cell survival  confirming the novel membrane-bound progestin receptors with defined actions. In the present study, expression of PGRMC1 was demonstrated in the CL and its expression became lower after LH/P4 depletion. This finding was similar to the earlier observation reported previously in the rhesus macaque . Further restoration of PGRMC1 expression was observed after LH replacement in LH/P4 ablated rhesus macaque . In contrast, in the present study PGRMC1 expression was not restored following LH replacement in LH/P4 depleted monkeys. The difference in the findings between the two studies could be attributed to differences in the treatment protocols employed and/or differences in the response to LH replacement observed between the two macaques.
Rothchild  aptly described CL as the most unique endocrine gland in the body, since there appears to be large interspecies variations in the mechanisms of regulation of its function. In higher primates, it is well established that the sole trophic stimulus for P4 secretion is the pulsatile secretion of pituitary LH. It is possible that P4 may not have an important role in its own secretion, but it might have an important role in maintaining structural integrity of the CL capable of responding to the luteotrophic stimulus. Hutchison and Zeleznik  demonstrated the rescue of CL function in the hypothalamus lesioned monkeys after reestablishment of pulsatile GnRH treatment indicating the resilience of CL tissue to the deprivation of luteotrophic support for long periods. In our studies, it was confirmed that CL was responsive to luteotrophic support for up to 48 h following inhibition of pituitary LH secretion . We utilized LH/P4 depleted and LH replacement model systems that did not require extreme surgical procedures, yet allow examining the CL function during LH/P4 depletion period as well as hormone replacement period. During induced luteolysis, even though LH replacement had profound effects on expression of several genes , but this could be the result of both direct as well as indirect effects of LH. In the present study, examination of direct effects of P4 on CL tissue revealed changes in expression of PR and NCOR1, but expression of many other genes involved in PR activation were not significantly affected. Since, the duration of P4 replacement lasted only for 24 h, for the effect to be apparent a longer duration of P4 replacement might be required. Alternatively, perhaps the effects of P4 on expression of genes may be observed only when high intraluteal P4 levels as seen in functional CL is achieved. In the present study, only circulating P4 levels were mimicked but more studies are required to achieve the high intra luteal P4 levels and determine its effect on the luteal expression of genes.
Tissue proteinases and genes associated with tissue remodelling namely, TIMPs, MMPs and ADAMTS-1 have been examined in the CL tissues of higher primates. In the present study, the findings of higher expression of TIMP3 and MMP-9 was observed during latter part of the luteal phase which is in accordance with similar findings reported for CL of the rhesus macaque . It was previously reported that ADAMTS-1 expression was high in early CL but declined thereafter at other stages of the luteal phase in the rhesus macaque . In the present study, however, no discernible pattern of expression of ADAMTS-1 was observed during the luteal phase. It should be pointed out that in the rhesus macaque the highest expression was seen only at day 2 of the luteal phase, while in the present study the day 2 CL was not examined. Treatment with GnRH R antagonists [i.e. LH/P4 depletion models] in both the macaques resulted in decreased expression of ADAMTS-1 at mid luteal phase. Interestingly, LH replacement, but not steroid replacement, prevented the decreased ADAMTS-1 expression seen after GnRH R antagonist treatment in the rhesus macaque, but in the present study both LH (observed only in real time RT-PCR analysis) and P4 replacement treatments restored or increased its expression. The reason for different findings on ADAMTS-1 expression after P4 (in the present study) or progestin (rhesus monkey; ) treatments is difficult to explain and perhaps, related to treatment protocol employed in both these studies. Also, the unexpected decrease in ADAMTS-1 expression seen after P4 supplementation at late luteal phase observed in the present study is difficult to explain. Additional experiments are necessary to address the regulation of ADAMTS-1 expression in the CL by LH and P4.
Stouffer and Young, 2004  reported a differential expression of MMP-9 with regulation of its expression to be at different levels involving both transcriptional and translational mechanisms in the CL of the rhesus macaque. In the present study, expression of MMP-9 in LH/P4 depletion and LH replacement studies remained high, but P4 replacement decreased its expression, suggesting MMP-9 expression appear to be regulated by P4. In accordance with our data, regulation of MMP-9 expression by P4 had been reported by other investigations in studies involving human endometrial explants and rabbit cervix [33, 34]. It should be pointed out that MMP-9 expression during spontaneous luteolysis [data not shown] and P4 supplementation during late luteal phases did not show significant change indicating that regulation of MMP-9 expression may be dependent on circulating LH as well as P4 concentrations. The decrease in expression of MMP-2 mRNA post P4 replacement during late luteal phase was similar to the reported decrease in MMP-2 expression in the human CL after treatment with hCG to mimic early pregnancy . It is possible that LH and P4 may synergistically act to regulate MMP-2 mRNA expression during the late luteal phase following P4 replacement as both are present in this model system, however, the molecular mechanisms for this synergistic action needs to be clarified. Decreased expression of pro-MMP-9, pro-MMP-2 and MMP-2 in endometrial cancer cell lines after in vitro administration of medroxy progesterone acetate, an synthetic progesterone preparation,  suggests that ECM remodelling could well be controlled by P4 as well as by LH in the CL tissue.
Expression of WNT7A, AREG, BMP5, HOXA1 and IHH, which are regarded as P4 target genes in uterus and mammary gland were examined in the CL after LH/P4 depletion and P4 replacement. WNT7A coordinates a variety of cell and developmental pathways and reported to guide hormonal responses during postnatal uterine growth . In the present study, the observation that P4 decreased WNT7A expression was similar to finding reported in the uterus . The observation of down regulated expression of BMP5 upon P4 replacement in the present study is consistent with the findings in uterus. Both WNT7A and BMP5 pathways appear to be interlinked since BMPs are shown to be induced by WNT signaling in the neurons . BMPs have been shown to stimulate steroidogenesis and proliferation in porcine theca cells , but the role of P4 regulated WNT and BMP genes on the structure and function of CL remains to be determined. The observation of low levels of phosphorylated P38 in CL post LH/P4 depletion, and its reversal after P4 replacement are consistent with similar findings reported in uterus and mammary gland after P4 treatment in mice . By activating the MAP kinase pathway, P4 appears to activate a unique biologically important molecule to control CL structure and function. It should be pointed out that changes in dynamics of MAP kinase activities have been shown to regulate CL function . PR is known to be phosphorylated by various kinases including MAP kinases in response to various physiological stimuli . The activation of MAP kinases have been shown to regulate expression and function of coactivators like NCOA1  and NCOA3 . Ligand-mediated down regulation and extensive loss of nuclear steroid receptor protein is well recognized and during this process, different phosphorylation pathways may affect the proteasome-mediated degradation of nuclear receptors including PR . It is reported that NCOA3 is phosphorylated by P38 and the phosphorylated NCOA3 is subjected to degradation by proteasome . Whether activation of P38 by P4 would have a role in regulating similar processes in the CL tissue remains to be addressed, however, the results in the present study suggest operation of similar mechanism. Gene expression analysis of 19 genes belonging to PR signaling during LH/P4 depletion and P4 replacement experiments revealed expression of 8 genes [(PR, NCOA1, NCOR1, WNT7A, BMP5, HOXA1, ADAMTS-1 and MMP-9) and protein expression of 3 genes] to be directly regulated by P4. The expression of other 11 genes appear to be regulated by both LH and P4 in the CL (see Table 3). In order to further extend the observations of direct actions of P4 and to further delineate role of P4 in the regulation of CL function, circulating P4 concentrations were increased for 24 h to mimic early pregnancy levels during the late luteal phase of the non fertile cycle. Expression of many of the genes examined was observed to be differentially expressed indicating that P4 may have an important role in the regulation of CL function during early pregnancy. Also, it should be pointed out that the general housekeeping activity that include anti-apoptotic actions, structural integrity of luteal cells and perhaps regulation of expression of few genes associated with overall maintenance of CL structure may be attributed to direct or indirect actions of P4.
In summary, experiments were performed to examine direct effects of P4← on the regulation of expression of genes in the monkey CL. The first experiment comprised of analysis of previously published microarray data to identify differentially expressed genes considered as target of P4 action and validation of some of these genes by real time RT-PCR analysis. In the second experiment, expression of many of the genes identified in the first experiment was examined throughout the luteal phase. In the third experiment, the direct effects of P4← on expression of many of these genes was examined following P4← replacement in monkeys depleted for circulating endogenous LH and P4← by way of inhibition of pituitary gonadotropin secretion. In the fourth experiment, the effects of P4← supplementation on expression of genes during the late luteal phase with declining endogenous P4← levels were examined. This experiment was done to determine whether increased P4← seen during early pregnancy regulate expression of these genes. The results from these studies indicated that P4← appears to regulate expression of many of the genes in the CL.
Experiments were carried out to assess the direct effects of P4 on CL function in monkeys. Expression of some of the genes involved in PR signaling and genes considered as targets of LH and P4 was analyzed in LH/P4 depletion and P4 replacement model as well as P4 supplementation model during the late luteal phase, and the results indicated that expression of number of genes appeared to be regulated directly or indirectly by P4. These results suggest that replacement of P4 during LH/P4 depletion (induced luteolysis model) is suitable for assessing the effects of P4 on CL function. Further, these studies suggest that CL could serve as a target tissue for examining genomic and non genomic actions of P4.
Stouffer RL, Xu F, Duffy DM: Molecular control of ovulation and luteinization in the primate follicle. Front Biosci. 2007, 12: 297-307. 10.2741/2065.
Rothchild I: The regulation of the mammalian corpus luteum. Recent Prog Horm Res. 1981, 37: 183-298.
Rothchild I: The corpus luteum revisited: are the paradoxical effects of RU486 a clue to how progesterone stimulates its own secretion?. Biol Reprod. 1996, 55 (1): 1-4. 10.1095/biolreprod55.1.1.
Stouffer RL, Duffy DM: Receptors for sex steroids in the primate corpus luteum New insight into gonadotropin and steroid action. Trends Endocrinol Metab. 1995, 6 (3): 83-89. 10.1016/1043-2760(94)00214-O.
Stouffer RL: Progesterone as a mediator of gonadotrophin action in the corpus luteum: beyond steroidogenesis. Hum Reprod Update. 2003, 9 (2): 99-117. 10.1093/humupd/dmg016.
Bishop CV, Hennebold JD, Stouffer RL: The effects of luteinizing hormone ablation/replacement versus steroid ablation/replacement on gene expression in the primate corpus luteum. Mol Hum Reprod. 2009, 15 (3): 181-193. 10.1093/molehr/gap005.
Engmann L, Losel R, Wehling M, Peluso JJ: Progesterone regulation of human granulosa/luteal cell viability by an RU486-independent mechanism. J Clin Endocrinol Metab. 2006, 91 (12): 4962-4968. 10.1210/jc.2006-1128.
Goyeneche AA, Calvo V, Gibori G, Telleria CM: Androstenedione interferes in luteal regression by inhibiting apoptosis and stimulating progesterone production. Biol Reprod. 2002, 66 (5): 1540-1547. 10.1095/biolreprod66.5.1540.
Goyeneche AA, Deis RP, Gibori G, Telleria CM: Progesterone promotes survival of the rat corpus luteum in the absence of cognate receptors. Biol Reprod. 2003, 68 (1): 151-158. 10.1095/biolreprod.102.007898.
Telleria CM, Goyeneche AA, Cavicchia JC, Stati AO, Deis RP: Apoptosis induced by antigestagen RU486 in rat corpus luteum of pregnancy. Endocrine. 2001, 15 (2): 147-155. 10.1385/ENDO:15:2:147.
Okuda K, Korzekwa A, Shibaya M, Murakami S, Nishimura R, Tsubouchi M, Woclawek-Potocka I, Skarzynski DJ: Progesterone is a suppressor of apoptosis in bovine luteal cells. Biol Reprod. 2004, 71 (6): 2065-2071. 10.1095/biolreprod.104.028076.
Rueda BR, Hendry IR, Hendry IW, Stormshak F, Slayden OD, Davis JS: Decreased progesterone levels and progesterone receptor antagonists promote apoptotic cell death in bovine luteal cells. Biol Reprod. 2000, 62 (2): 269-276. 10.1095/biolreprod62.2.269.
McKenna NJ, O'Malley BW: Combinatorial control of gene expression by nuclear receptors and coregulators. Cell. 2002, 108 (4): 465-474. 10.1016/S0092-8674(02)00641-4.
Wagner BL, Norris JD, Knotts TA, Weigel NL, McDonnell DP: The nuclear corepressors NCoR and SMRT are key regulators of both ligand- and 8-bromo-cyclic AMP-dependent transcriptional activity of the human progesterone receptor. Mol Cell Biol. 1998, 18 (3): 1369-1378.
Jackson TA, Richer JK, Bain DL, Takimoto GS, Tung L, Horwitz KB: The partial agonist activity of antagonist-occupied steroid receptors is controlled by a novel hinge domain-binding coactivator L7/SPA and the corepressors N-CoR or SMRT. Mol Endocrinol. 1997, 11 (6): 693-705. 10.1210/me.11.6.693.
Priyanka S, Jayaram P, Sridaran R, Medhamurthy R: Genome-wide gene expression analysis reveals a dynamic interplay between luteotropic and luteolytic factors in the regulation of corpus luteum function in the bonnet monkey (Macaca radiata). Endocrinology. 2009, 150 (3): 1473-1484. 10.1210/en.2008-0840.
Smith MF, McIntush EW, Ricke WA, Kojima FN, Smith GW: Regulation of ovarian extracellular matrix remodelling by metalloproteinases and their tissue inhibitors: effects on follicular development, ovulation and luteal function. J Reprod Fertil Suppl. 1999, 54: 367-381.
Young KA, Hennebold JD, Stouffer RL: Dynamic expression of mRNAs and proteins for matrix metalloproteinases and their tissue inhibitors in the primate corpus luteum during the menstrual cycle. Mol Hum Reprod. 2002, 8 (9): 833-840. 10.1093/molehr/8.9.833.
Young KA, Stouffer RL: Gonadotropin and steroid regulation of matrix metalloproteinases and their endogenous tissue inhibitors in the developed corpus luteum of the rhesus monkey during the menstrual cycle. Biol Reprod. 2004, 70 (1): 244-252. 10.1095/biolreprod.103.022053.
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 (7): 2470-2478. 10.1210/jc.83.7.2470.
Suresh PS, Medhamurthy R: Dynamics of circulating concentrations of gonadotropins and ovarian hormones throughout the menstrual cycle in the bonnet monkey: role of inhibin A in the regulation of follicle-stimulating hormone secretion. Am J Primatol. 2009, 71 (10): 817-824. 10.1002/ajp.20709.
Young KA, Tumlinson B, Stouffer RL: ADAMTS-1/METH-1 and TIMP-3 expression in the primate corpus luteum: divergent patterns and stage-dependent regulation during the natural menstrual cycle. Mol Hum Reprod. 2004, 10 (8): 559-565. 10.1093/molehr/gah079.
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 USA. 2000, 97 (9): 4689-4694. 10.1073/pnas.080073497.
Peluso JJ, Liu X, Saunders MM, Claffey KP, Phoenix K: Regulation of ovarian cancer cell viability and sensitivity to cisplatin by progesterone receptor membrane component-1. J Clin Endocrinol Metab. 2008, 93 (5): 1592-1599. 10.1210/jc.2007-2771.
Jeong JW, Lee KY, Han SJ, Aronow BJ, Lydon JP, O'Malley BW, DeMayo FJ: The p160 steroid receptor coactivator 2, SRC-2, regulates murine endometrial function and regulates progesterone-independent and -dependent gene expression. Endocrinology. 2007, 148 (9): 4238-4250. 10.1210/en.2007-0122.
Khatua A, Wang X, Ding T, Zhang Q, Reese J, DeMayo FJ, Paria BC: Indian hedgehog, but not histidine decarboxylase or amphiregulin, is a progesterone-regulated uterine gene in hamsters. Endocrinology. 2006, 147 (9): 4079-4092. 10.1210/en.2006-0231.
Yadav VK, Medhamurthy R: Dynamic changes in mitogen-activated protein kinase (MAPK) activities in the corpus luteum of the bonnet monkey (Macaca radiata) during development, induced luteolysis, and simulated early pregnancy: a role for p38 MAPK in the regulation of luteal function. Endocrinology. 2006, 147 (4): 2018-2027. 10.1210/en.2005-1372.
Ace CI, Okulicz WC: Microarray profiling of progesterone-regulated endometrial genes during the rhesus monkey secretory phase. Reprod Biol Endocrinol. 2004, 2: 54-10.1186/1477-7827-2-54.
Richer JK, Jacobsen BM, Manning NG, Abel MG, Wolf DM, Horwitz KB: Differential gene regulation by the two progesterone receptor isoforms in human breast cancer cells. J Biol Chem. 2002, 277 (7): 5209-5218. 10.1074/jbc.M110090200.
Cai Z, Stocco C: Expression and regulation of progestin membrane receptors in the rat corpus luteum. Endocrinology. 2005, 146 (12): 5522-5532. 10.1210/en.2005-0759.
Peluso JJ, Liu X, Gawkowska A, Johnston-MacAnanny E: Progesterone activates a progesterone receptor membrane component 1-dependent mechanism that promotes human granulosa/luteal cell survival but not progesterone secretion. J Clin Endocrinol Metab. 2009, 94 (7): 2644-2649. 10.1210/jc.2009-0147.
Hutchison JS, Zeleznik AJ: The corpus luteum of the primate menstrual cycle is capable of recovering from a transient withdrawal of pituitary gonadotropin support. Endocrinology. 1985, 117 (3): 1043-1049. 10.1210/endo-117-3-1043.
Imada K, Ito A, Sato T, Namiki M, Nagase H, Mori Y: Hormonal regulation of matrix metalloproteinase 9/gelatinase B gene expression in rabbit uterine cervical fibroblasts. Biol Reprod. 1997, 56 (3): 575-580. 10.1095/biolreprod56.3.575.
Marbaix E, Donnez J, Courtoy PJ, Eeckhout Y: Progesterone regulates the activity of collagenase and related gelatinases A and B in human endometrial explants. Proc Natl Acad Sci USA. 1992, 89 (24): 11789-11793. 10.1073/pnas.89.24.11789.
Di Nezza LA, Jobling T, Salamonsen LA: Progestin suppresses matrix metalloproteinase production in endometrial cancer. Gynecol Oncol. 2003, 89 (2): 325-333. 10.1016/S0090-8258(03)00089-1.
Carta L, Sassoon D: Wnt7a is a suppressor of cell death in the female reproductive tract and is required for postnatal and estrogen-mediated growth. Biol Reprod. 2004, 71 (2): 444-454. 10.1095/biolreprod.103.026534.
Kasai M, Satoh K, Akiyama T: Wnt signaling regulates the sequential onset of neurogenesis and gliogenesis via induction of BMPs. Genes Cells. 2005, 10 (8): 777-783. 10.1111/j.1365-2443.2005.00876.x.
Brankin V, Quinn RL, Webb R, Hunter MG: BMP-2 and -6 modulate porcine theca cell function alone and co-cultured with granulosa cells. Domest Anim Endocrinol. 2005, 29 (4): 593-604. 10.1016/j.domaniend.2005.04.001.
Han SJ, Tsai SY, Tsai MJ, O'Malley BW: Distinct temporal and spatial activities of RU486 on progesterone receptor function in reproductive organs of ovariectomized mice. Endocrinology. 2007, 148 (5): 2471-2486. 10.1210/en.2006-1561.
Lange CA, Shen T, Horwitz KB: Phosphorylation of human progesterone receptors at serine-294 by mitogen-activated protein kinase signals their degradation by the 26 S proteasome. Proc Natl Acad Sci USA. 2000, 97 (3): 1032-1037. 10.1073/pnas.97.3.1032.
Wu RC, Qin J, Yi P, Wong J, Tsai SY, Tsai MJ, O'Malley BW: Selective phosphorylations of the SRC-3/AIB1 coactivator integrate genomic reponses to multiple cellular signaling pathways. Mol Cell. 2004, 15 (6): 937-949. 10.1016/j.molcel.2004.08.019.
Gianni M, Parrella E, Raska I, Gaillard E, Nigro EA, Gaudon C, Garattini E, Rochette-Egly C: P38MAPK-dependent phosphorylation and degradation of SRC-3/AIB1 and RARalpha-mediated transcription. EMBO J. 2006, 25 (4): 739-751. 10.1038/sj.emboj.7600981.
We are grateful to Dr. Basavanagouda and staff of Primate Research Laboratory for assistance with surgeries and Mr. Kunal BS for help rendered in the preparation of MS. Financial support from Department of Biotechnology, Government of India, to conduct these studies is gratefully acknowledged.
The authors declare that they have no competing interests.
PS and RM participated in designing, conducting experiments, analysis of results and preparation of manuscript. KCJ participated in the preparation of manuscript. All authors read and approved the final manuscript.