REGULATION OF OVULATORY GENES IN BOVINE ... - Reproduction

7 downloads 0 Views 1MB Size Report
Oct 16, 2014 - (TNFAIP6), pentraxin-3 (PTX3), epidermal growth factor like factors: amphiregulin (AREG) and. 3 epiregulin (EREG) are essential for successful ...
Page 1 of 31 Reproduction Advance Publication first posted on 16 October 2014 as Manuscript REP-14-0337

REGULATION OF OVULATORY GENES IN BOVINE GRANULOSA CELLS: LESSONS FROM siRNA SILENCING OF PTGS2 Ketan Shrestha1, Karolina Lukasik2, Anja Baufeld3, Jens Vanselow3, Uzi Moallem4 and Rina Meidan1§ 1

Department of Animal Sciences, The Robert H. Smith Faculty of Agriculture, Food and

Environment, The Hebrew University of Jerusalem, Rehovot, Israel. 2

Department of Reproductive Immunology and Pathology, Institute of Animal Reproduction

and Food Research, Polish Academy of Science, Olsztyn, Poland. 3

Department of Reproductive Biology, Leibniz Institute of Farm Animal Biology,

Dummerstorf, Germany. 4

Department of Ruminant Science, Agricultural Research Organization, Bet-Dagan, Israel.

Running title: Silencing of prostaglandin endoperoxide synthase 2. Keywords – ovulation, granulosa cells, follicle, RNA interference, endometrial cells.

Supported by a grant from Chief Scientist of the Israeli Ministry of Agriculture (RM and UM). The authors have nothing to disclose. §

Corresponding author: Prof. Rina Meidan, Department of Animal Sciences, The Robert

H. Smith Faculty of Agriculture, Food and Environment. The Hebrew University of Jerusalem, Rehovot 76100 Israel, email: [email protected], phone number: +97289489394, fax number: 97289465763.

0

Copyright © 2014 by the Society for Reproduction and Fertility.

Page 2 of 31

1

Abstract

2

Prostaglandin endoperoxide synthase 2 (PTGS2), tumor necrosis factor-alpha-induced protein-6

3

(TNFAIP6), pentraxin-3 (PTX3), epidermal growth factor like factors: amphiregulin (AREG) and

4

epiregulin (EREG) are essential for successful ovulation. We compared here the induction of

5

these ovulatory genes in bovine granulosa cells (GCs) in vivo (after LH surge) and in vitro

6

(forskolin treatment). These genes were markedly stimulated in GCs isolated from cows 21hr

7

after LH- surge. In isolated GCs, forskolin induced a distinct temporal profile for each gene.

8

Generally, there was a good agreement between the in vivo and in vitro inductions of these genes

9

except for PTX3. Lack of PTX3 induction in isolated GCs culture suggests that other follicular

10

compartments may mediate its induction by LH. Next, to study the role of PTGS2 and

11

prostaglandins in the cascade of ovulatory genes, PTGS2 was silenced with small interfering

12

RNA (siRNA). PTGS2 siRNA caused a marked and specific knockdown of PTGS2 mRNA and

13

PGE2 production (70% compared to scrambled siRNA) in bovine GCs. Importantly, PTGS2

14

silencing also reduced AREG, EREG and TNFAIP6 mRNA but not PTX3. Exogenous PGE2

15

increased AREG, EREG and TNFAIP6 mRNAs, further confirming that these genes are

16

prostanoid-dependent. A successful and specific knockdown of PTGS2 was also achieved in

17

endometrial cells (EndoCs) expressing PTGS2. Then, cholesterol-conjugated PTGS2 (chol-

18

PTGS2) siRNA that facilitates cells’ entry was examined. In EndoCs, but not in GCs, Chol-

19

PTGS2 siRNA succeeded to reduce PTGS2 and PGE2 levels even without transfection reagent.

20

PTGS2 knockdown is a promising tool to critically examine the functions of PTGS2 in the

21

reproductive tract.

22

1

Page 3 of 31

23

Introduction

24

Ovulation is triggered by LH surge. The LH receptor activates several families of

25

heterotrimeric G-proteins, but the activation of Gs and cAMP/PKA cascade is widely accepted

26

as important downstream signaling of LH action (Fan et al. 2011, Breen et al. 2013). It brings

27

about series of changes such as resumption of meiosis in the oocyte, reprogramming of the

28

granulosa cells (GCs) and theca cells (TCs) layer of the follicular wall. These events are

29

accompanied by the expression of new mRNAs and proteins and result in the release of a

30

fertilizable oocyte and the luteinization of follicular cells (Richards 2001, Richards 2005). In

31

preparation for oocyte expulsion, the extracellular matrix (ECM) of cumulus oocyte complex

32

(COC) undergoes remodeling (Richards 1994, Salustri et al. 1996). A number of genes

33

mediating these processes have been identified; amongst them are Prostaglandin endoperoxide

34

synthase 2 (PTGS2), tumor necrosis factor-alpha-induced protein 6 (TNFAIP6), pentraxin 3

35

(PTX3), epidermal growth factor (EGF) like factors: amphiregulin (AREG) and epiregulin

36

(EREG) (Mukhopadhyay et al. 2001, Salustri et al. 2004, Conti et al. 2006, Bridges & Fortune

37

2007).

38

LH-induced PTGS2 in the GCs layer increases prostaglandins (PGs), especially prostaglandin

39

endoperoxide 2 (PGE2) in the follicular fluid (Sirois & Richards 1992, Sirois 1994, Tsai et al.

40

1996, Liu et al. 1997). The essential role of PTGS2 is demonstrated in gene knockout

41

experiments and PTGS inhibitors studies which resulted in failed ovulation with normal

42

luteinization and corpus luteum formation (Dinchuk et al. 1995, Davis et al. 1999, Peters et al.

43

2004). Later, it was demonstrated that LH/hCG upregulate EGF-like factors in murine GCs;

44

these peptides mediated LH action in COC expansion and oocyte maturation (Park et al. 2004,

45

Conti et al. 2006, Shimada et al. 2006). During COC expansion, LH surge also induces 2

Page 4 of 31

46

TNFAIP6 which is a hyaluronan binding protein (Lee et al. 1992) critical for the stability of

47

glycosaminoglycan hyaluronan (HA) rich ECM (Ochsner et al. 2003). Study from TNFAIP6

48

knockout mice showed defective COC expansion and infertility (Fulop et al. 2003). Likewise,

49

PTX3, belonging to the long pentraxin family of inflammatory proteins, was expressed in COC

50

(Varani et al. 2002, Salustri et al. 2004). It plays protective role during the formation of HA

51

rich matrix of the COC by cross-linking HA polymers through interactions with heavy chains

52

of inter-α trypsin inhibitor (IαI) and/or TNFAIP6 (Ievoli et al. 2011). An abnormal COC

53

expansion characterized by unstable extracellular matrix was reported in PTX3 deficient mice

54

(Salustri et al. 2004). Thus, PTGS2, EGF like factors (AREG and EREG) and HA binding

55

proteins (TNFAIP6 and PTX3) are all the obligatory factors for ECM stability, COC expansion

56

and follicular rupture.

57

Most of our knowledge on ovulatory genes stems from rodents especially mice, however

58

reproductive strategies differ in rodents and large mono-ovulatory animals. Therefore, data

59

from rodents cannot always be extrapolated in mono-ovulatory mammals such as bovine (Bahr

60

& Wolf 2012). In addition, these animals are not amendable for gene deletion studies. In such

61

species, small interfering RNA (siRNA) can be used as means for gene knockdown (Fang et al.

62

2013). Therefore, we aimed here to study induction of ovulatory genes in bovine GCs under in

63

vivo and in vitro conditions. siRNA silencing of PTGS2 was used to critically examine the role

64

of PGE2 in the cascade of events in GCs that leads to ovulation. As a proof of concept we also

65

employed endometrial cells (EndoCs), these cells similarly to GCs, express PTGS2 and their

66

synthesized PGs play significant physiological roles.

67

3

Page 5 of 31

68

Materials and Methods

69

Animals and sample collection

70

Follicles were collected from Holstein-Friesian cows as described (Vanselow et al. 2010,

71

Christenson et al. 2013). Large dominant follicles before LH surge were collected after

72

carefully monitoring a growing cohort of follicles in normally cycling cows by transrectal

73

ultrasonography (Aloka SSD-500, Aloka GmbH, Meerbusch, Germany). Animals were

74

slaughtered at days 7 (n=1) and 8 (n=2) of the estrous cycle during the first follicular wave.

75

Only the largest growing, but not stagnating or regressing follicle of each animal was collected.

76

To collect dominant follicles after the LH surge, normally cycling cows were treated with

77

500µg prostaglandin F2 alpha (PGF2a) (PGF Veyx® forte, Veyx Pharma GmbH, Germany) at

78

day 8 of the cycle, to induce luteolysis. Forty-eight hours later, animals were injected with 100

79

µg of a GnRH analogue (Gonavet Veyx®,Depherelin, Veyx Pharma GmbH) to induce the LH

80

surge. The animals were slaughtered 23 hours later and the largest growing follicle of each

81

animal was isolated (n=3). GCs were collected by aspiration of the follicular fluid with an 18G

82

needle. The fluid was centrifuged (2min, 400rcf) and the sediment cells were frozen in liquid

83

nitrogen and stored at -80°C for RNA preparation. RNA isolation, reverse transcription and

84

qPCR were carried out as described previously (Christenson et al. 2013).

85

All procedures involving living animals (injections and ultrasonography examinations in cattle)

86

were performed according to the German law for animal protection (8a TierSchG i.V.m. 29

87

Tierschutz-Versuchstierverordnung). The named procedures are not subjected to specific

88

permit by the governmental authority Landesamt für Landwirtschaft, Lebensmittelsicherheit

89

und Fischerei (LALLF) Rostock. 4

Page 6 of 31

90

Isolation and culture of granulosa cells

91

Ovaries were collected at a local abattoir as previously described (Meidan et al. 1990). GCs

92

were enzymatically dispersed by using a combination of collagenase I (125 units/ml),

93

hyaluronidase III (36 units/ml) and deoxyribonuclease I (11 units/ml) in Dulbecco's modified

94

Eagle medium (DMEM)/F-12 containing 2 mM L-Glutamine and 100 µg/ml of

95

penicillin/streptomycin. Only large follicles (>10 mm in diameter) containing ≥ 4 million

96

viable cells were used. GCs were cultured overnight in medium containing 3% fetal calf serum

97

(FCS). Cells were then incubated with DMEM/F-12 containing 1% FCS (control) and with

98

PGE2 (1 µM; Cayman Chemical Co) or with forskolin (FRS; 10 µM) for various time points as

99

indicated. Unless otherwise stated all biochemical were from Sigma-Aldrich, St. Louis, MO

100

and cell culture materials were from Biological Industries, Kibbutz Beit Haemek, Israel.

101

Isolation and culture of endometrial cells (EndoCs)

102

Uterine horns from cows at late luteal phase were collected from the abattoir (Arosh et al.

103

2002). Endometrial strips were dissected and incubated in enzyme solution (18 units/ml DNAse

104

and 315 units/ml collagenase I) at 38o C for 25 minutes in shaking water bath. Cell suspension

105

was filtered with cell strainer, centrifuged and cultured in DMEM/F-12 containing 10% FCS.

106

Confluent plates were trypsinized with trypsin EDTA solution (0.25% trypsin and 0.02%

107

EDTA). Cells of passages 2-3 were utilized in this study. Based on morphology and growth in

108

culture the cells were identified as endometrial stromal cells (Cherny & Findlay 1990).

5

Page 7 of 31

109

Transfection of cells

110

For transfection experiments GCs were trypsinized with trypsin EDTA solution (0.05% trypsin

111

and 0.02% EDTA) immediately after isolation from follicles. Trypsinized cells were seeded

112

(3x105 GCs; 0.8X105 EndoCs) in 6-well plates and cultured for up to 24 hr in 3% FCS. Then

113

cells were transfected using Lipofectamine 2000 reagent (Invitrogen, Carlsbad, CA) according

114

to the manufacturer’s protocol in 1% FCS. Cells were transfected with siRNA sequence

115

targeting (50 nmol/L) naked PTGS2 siRNA or cholesterol conjugated PTGS2 siRNA or

116

scrambled

117

GUGAAAGGCUGUCCCUUUA[dT][dT],

118

UAAAGGGACAGCCUUUCAC[dT][dT] corresponding to bases 1781-1799 of the bovine

119

PTGS2 mRNA sequence (NM_174445). Sequence for cholesterol conjugated PTGS2 siRNA

120

was same as for naked PTGS2 siRNA with cholesterol conjugation in 3’end of sense strand: S,

121

GUGAAAGGCUGUCCCUUUA[dA][dT][CholTEG],

122

UAAAGGGACAGCCUUUCAC[dA][dT]. Scrambled siRNA sequence - negative control was

123

S,

124

ACGUGACACGUUCGGAGAA[dT][dT]. A day after transfection PTGS2 was induced in

125

GCs by FRS (10 µM) and in EndoCs by phorbol 12-myristate 13-acetate (PMA; 100ng/ml),

126

then cells were harvested and total RNA was extracted from cells 48hr after transfection.

127

For treatment of GCs and EndoCs in absence of transfection reagent, cells were incubated with

128

cholesterol-conjugated PTGS2 siRNA for 24hr and total RNA was extracted from cells after

129

48hr.

siRNA.

Sequence

of

naked

PTGS2

UUCUCCGAACGUGUCACGU[dT][dT],

siRNA

was

sense

antisense

(S), (AS),

AS,

and

AS,

6

Page 8 of 31

130

RNA extraction and real-time PCR

131

Total RNA was isolated from cells using TriFast reagent (Peqlab Biotechnologie GmbH,

132

Erlangen, Germany) according to the manufacturer's instructions. 1 µg of total RNA was

133

reverse transcribed using M-MuLV Reverse Transcriptase (200 units/uL), M-MuLV RT Buffer

134

(New England Biolabs, Ipswich, MA), random primer (100 nM), oligo-dT (100 µM), and

135

dNTPs mix (100 mM) (Bioline Reagents Limited, London). Real-time PCRs were performed

136

using the Mx3000P system (Stratagene, Garden Grove, CA), using Platinum SYBR Green

137

(SuperMix, Invitrogen, Carlsbad, CA), as previously described (Klipper et al. 2009). Gene

138

expression (PTGS2, TNFAIP6, PTX3, AREG and EREG) was analyzed by quantitative real-

139

time PCR. The glyceraldehyde 3-phosphate dehydrogenase (GAPDH) gene was used as the

140

housekeeping gene. The threshold cycle number (Ct) was used to quantify the relative

141

abundance of the gene; arbitrary units were calculated as 2-∆Ct=2- (Ct target gene - Ct housekeeping gene).

142

The primer sequences used were as follows: GAPDH (NM_001034034) forward, 5’-

143

GTCTTCACTACCATGGAGAAGG -3’, reverse, 5’- TCATGGATGACCTTGGCCAG -3’;

144

PTGS2 (NM_174445) forward, 5’- CAGCGGTGCAGCAAATCCTTG -3’,

145

CTGTGTTGGGAGTGGGTTTCA

146

CTATAGCTGCTTTCGTCTCTGC -3’, reverse, 5’- CGTTCTTCAGCGACACCTTCA -3’;

147

EREG (XM_002688367) forward, 5’- CTGCTGCTCGTCCTGGTTTTC -3’, reverse, 5’-

148

GCTGTGCAGTTATCTCCCGAC

149

ATGGCTTGAACAAGCAGCAGG -3’, reverse, 5’- GCCATCCACCCAGCAGCACA -3’;

150

PTX3 (NM_001076259) forward, 5’- AGCCTCTTGCCTCGTCCCCTC -3’, reverse, 5’-

151

TCTGAGTTCTCCGCCGACACT -3’.

-3’;

-3’;

AREG

(NM_001099092)

TNFAIP6(NM_001007813)

reverse, 5’-

forward,

forward,

5’-

5’-

152 7

Page 9 of 31

153

Prostaglandin E2 analysis

154

Media from the cell cultures were collected on the day of RNA isolation. Levels of PGE2 were

155

measured by PGE2 enzyme immunoassay kit- monoclonal (Cayman Chemical Co) according to

156

manufacturer’s instruction. The standard curve ranged from 6.4 pg/ml to 1000 pg/ml. Cross

157

reactivity of the assay with PGE2 was 100% and 0.01 % for PGF-2a.

158

Statistical analyses

159

Data are presented as means ± SEM; experiments were repeated at least three times.

160

Expressions of specific mRNA transcripts were normalized relative to the abundance of

161

GAPDH mRNA. Data were analyzed by Student t-test. Differences were considered significant

162

at P < 0.05. Asterisks represent significant differences from their respective controls. *P