Heifer synchronization and ovariectomy
A cohort of cyclic heifers (presence of corpus luteum > 16 mm and pre-ovualtory follicle > 12 mm) had their estrous cycle synchronized. The heifers received an intravaginal progesterone-releasing device (1.38 g progesterone; Eazi-Breed CIDR, Zoetis, Parsippany-Troy Hills, NJ, USA) and an injection of GnRH agonist (100 µg Factrel®, Zoetis, Parsippany-Troy Hills, NJ, USA) intramuscularly. The intravaginal device was removed 5 days later, and heifers were given PGF2α analog (25 mg dinoprost tromethamine, Lutalyse, Zoetis, Parsippany-Troy Hills, NJ, USA) intramuscularly at the time of the intravaginal device removal and again 24 h later. Heifers were examined daily by transrectal ultrasonography to assess response to the synchronization program. Heifers that initiate follicle development and ovulate after the first GnRH treatment and had luteolysis after PGF2α were eligible to have an ovariectomy. The procedure to remove the ovaries was performed at 48 h after the second PGF2α injection when the dominant follicles reached at least 12 mm in diameter. Ovariectomy was performed via colpotomy in the standing position under sedation with 15 mg intramuscular xylazine (AnaSed® Injection, Lloyd, Inc., Shenandoah, IA, USA) and caudal epidural anesthesia with 5 mL of 2% (w/v) lidocaine HCl (VetOne®, Boise, Idaho, USA) performed. Heifers were treated pre-operatively with ceftiofur crystalline free acid (6.6 mg/kg; Excede®, Zoetis, Parsippany-Troy Hills, NJ, USA) injected subcutaneously in the base of the ear to prevent infection. Flunixin meglumine (2.2 mg/kg; Norbrook® Inc., Overland Park, KS, USA) was administered intravenously daily for two to control inflammation and pain. An incision was made in the dorsolateral aspect of the vaginal fornix, and the peritoneum was manually punctured after blunt dissection through the adventitia. The mesovarium was manually compressed, and the ovary containing the dominant follicle was removed using a chain écraseur. Ovaries were placed in ice-cold phosphate-buffered saline solution containing 2% antibiotic-antimycotic mixture (25 µg/mL amphotericin B, 10,000 units/mL penicillin, 10,000 µg/mL streptomycin; Gibco, Gaithersburg, MD, USA) for transport to the laboratory.
Isolation of follicular wall cells and treatment allocation
The pre-ovulatory follicle was identified and dissected away from each ovary for use in the tissue culture system (Fig. 1). Follicular fluid was aspirated to facilitate further dissection and frozen at − 80 °C. The follicles were dissected into quarters, and the theca interna with adherent granulosa cells was peeled from the theca externa and surrounding stromal tissue. The remaining follicle wall preparations (theca interna and granulosa cells) were cut into 26 pieces (average weight: 5.3 ± 0.7 mg), 24 of which were transferred to a costar 24-well plate (1 piece/well; Cambridge, MA, USA) for tissue culture, as previously described (21, 22). The remainder of the tissue was flash-frozen in liquid nitrogen and maintained at − 80 °C until mRNA extraction.
The dissected 24 follicle wall pieces containing theca and granulosa cells were randomly allocated to receive a culture medium that was either supplemented with purified bovine NGF (100 ng/mL, n = 12) or left untreated (control, n = 12). The single plate was incubated at 37 °C in a humidified incubator gassed with 5% CO2:95% air for 72 h. The NGF used was purified from bovine seminal plasma, as described previously (13). The follicle wall pieces were cultured in 0.5 mL of medium consisting of Eagle’s MEM (Invitrogen, Carlsbad, CA, USA) supplemented with 1% L-glutamine (Gibco, Gaithersburg, MD, USA), 1% nonessential amino acids (Sigma-Aldrich, St. Louis, MO, USA), 1% penicillin-streptomycin (Sigma-Aldrich), 1% ITS (10 ng/ml insulin, 5.5 ng/ml transferrin, 5.5 ng/ml selenium, Sigma-Aldrich), 10% fetal bovine serum (FBS, Atlanta Biologicals, Lawrenceville, GA), 40 ng/mL cortisol (Sigma-Aldrich), 4 ng/mL human recombinant LH (Dr. A. F. Parlow, National Hormone and Peptide Program, Harbor- UCLA Medical Center, Torrance, CA, USA), and 4 ng/mL human recombinant FSH (Dr. A. F. Parlow, National Hormone and Peptide Program).
Hormone assays
The aliquots (0.5 mL) of culture medium at 3, 6, 12, 24, 48, and 72 h of culture were collected and preserved at − 20 °C for subsequent steroid assays. Progesterone, testosterone, and estradiol-17β secretions by granulosa and theca cells into the culture media were assessed using immunoassays (Immulite 2000 XPi platform; Siemens Medical Solutions, Malvern, PA, USA, Inc.). Total hormone production for each well was calculated by multiplying the measured concentration by the volume of media (0.5 mL) and then dividing by tissue weight (mg). Intra-assay coefficient of variations were 4.0% (testosterone), 2.4% (progesterone), and 3.1% (estradiol-17β). Inter-assay coefficient of variations were 12% (testosterone), 19% (progesterone), and 15% (estradiol-17β). The progesterone assay had a detection range of 0.2 to 40 ng/mL and a sensitivity of 0.1 ng/mL. The testosterone assay had a detection range of 20 to 1600 ng/mL and a sensitivity of 15 ng/dL. The estradiol-17β assay had a detection range of 20-2000 pg/mL and a sensitivity of 15 pg/mL.
Quantitative real-time PCR analyses
At the completion of the 72-h culture period, follicle wall pieces were weighed and flash-frozen and kept at -80 °C until RNA extraction. Follicular tissue mRNA expression was determined for LH/ choriogonadotropin receptor (LHCGR), FSH receptor (FSHR), PGE synthase (PGES), vascular endothelial growth factor A isoform 121 (VEGFA121), fibroblast growth factor 2 (FGF2), estrogen receptor 1 (ESR1), steroidogenic acute regulatory protein (STAR), cytochrome P450 family 11 subfamily A member 1 (CYP11A1), cytochrome P450 family 17 subfamily A member 1 (CYP17A1), cytochrome P450 family 19 subfamily A member 1 (CYP19A1), hydroxyl-delta-5-steroid dehydrogenase 3-beta (HSD3B), and hydroxysteroid 17-beta dehydrogenase (HSD17B). Primers were designed for the constitutively expressed mRNAs, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal protein L 15 (RPL15), and ribosomal protein L 19 (RPL19), with the expression value of each gene normalized to the mean values of these genes (Table 1). Relative expression values were obtained by determining the PCR amplification efficiency (E = 2) to the power of the delta-delta threshold cycle (ΔΔCt) obtained from the ΔCt least square mean differences of pairwise comparisons between initial and cultured tissue (23).
Table 1
List of genes and primers used for quantitative real-time PCR
Target genes and abbreviations | NCBI Sequence | Primer | Primer Sequence |
Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) | NM_001034034 | Forward Reverse | 5’-GGCGCCAAGAGGGTCAT-3’ 5’-ACGCCCATCACAAACATGG-3’ |
Ribosomal protein L 15 (RPL15) | AY786141 | Forward Reverse | 5’-TGGAGAGTATTGCGCCTTCTC-3’ 5’-CACAAGTTCCACCACACTATTGG-3’ |
Ribosomal protein L 19 (RPL19) | NM_001040515 | Forward Reverse | 5’-CAGACGATACCTGAATCTAAGAAGA-3’ 5’-TGAGAATCCGCTTGTTTTTGAA-3’ |
Follicular stimulating hormone receptor (FSHR) | NM_174061 | Forward Reverse | 5′-CGACTCTGTCACTGCTCTAACGG-3′ 5′-CGTCAATTCCTTTGGCATAGGTGG-3′ |
Luteinizing hormone/ choriogonadotropin receptor (LHCGR) | NM_174381 | Forward Reverse | 5’-CAGTCCCCCGCTTTCTCAT-3’ 5’-GTAGAGCCCCATGCAGAAGTCT-3’ |
Steroidogenic acute regulatory protein (STAR) | XR_083945 | Forward Reverse | 5’-GGATTAACCAGGTTCGGCG-3’ 5’-CTCTCCTTCTTCCAGCCCTC-3’ |
Cytochrome P450 family 11 subfamily A member 1 (CYP11A1) | NM_176644 | Forward Reverse | 5’-GCCACATCGAGAACTTCCAGAAG-3’ 5’-CTGGTGTGGAACATCTTGTAGACG-3’ |
Hydroxyl-delta-5-steroid dehydrogenase 3-beta (HSD3B) | NM_174343 | Forward Reverse | 5′-TGTTGGTGGAGGAGAAGGATCTG-3′ 5′-TGGGTACCTTTCACATTGACGTTC-3′ |
Hydroxysteroid 17-beta dehydrogenase (HSD17B) | NM_001102365 | Forward Reverse | 5′-TTGTGCGAGAGTCTGGCGATTCT-3′ 5′-AGGAATCGCTCGGTGGTGAAGTA-3′ |
Cytochrome P450 family 17 subfamily A member 1 (CYP17A1) | NM_174304 | Forward Reverse | 5’-TGTGGCCCCTACGCTGAT-3’ 5’-CGCCAATGCTGGAGTCAAT-3’ |
Cytochrome P450 family 19 subfamily A member 1 (CYP19A1) | NM_174305 | Forward Reverse | 5′-GTCCGAAGTTGTGCCTATTGCCAGC-3′ 5′- CCTCCAGCCTGTCCAGATGCTTGG − 3′ |
Estrogen receptor 1 (ESR1) | NM_001001443 | Forward Reverse | 5’-AGGGAAGCTCCTATTTGCTCC-3’ 5’-CGGTGGATGTGGTCCTTCTCT-3’ |
Fibroblast growth factor 2 (FGF2) | NM_174056 | Forward Reverse | 5’-GAACGGGGGCTTCTTCCT-3’ 5’-CCCAGTTCGTTTCAGTGCC-3’ |
Prostaglandin E synthase (PGES) | NM_174443 | Forward Reverse | 5’-AGGACGCTCAGAGACATGGA-3’ 5’-TTCGGTCCGAGGAAAGAGTA-3’ |
Vascular endothelial growth factor A isoform 121 (VEGFA121) | NM_174216 | Forward Reverse | 5’- CCGTCCCATTGAGACCCTG-3’ 5’-CGGCTTGTCACAATTTTTCTTGTC-3’ |
Follicle tissue lysis and RNA extraction were conducted according to the manufacturer's recommendations (PureLink RNA Mini Kit, Invitrogen, Carlsbad, CA, USA). Isolated RNA was evaluated for concentration and purity using a NanoDrop One Spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). A maximum of 2 µg of mRNA was used to synthesize complementary DNA using a commercial kit (High-capacity cDNA Reverse Transcription Kit, Applied Biosystems, Foster City, CA, USA) supplemented with RNase inhibitor (RNase Inhibitor, human placenta, New England BioLabs, Ipswich, MA, USA). Complementary DNA was used for quantitative real-time reverse transcription PCR using a 7500 Real-Time PCR Detection System (Applied Biosciences) with PowerUp™ SYBR™ Green Master Mix (2X; Applied Biosciences). All assays were carried out in triplicate for each target mRNA. The amplification conditions were as follows: 50 °C for 2 min, 95 °C for 2 min, and 40 cycles at 95 °C for 15 s and 60 °C for 60 s.
Statistical analyses
Data are presented as percentage mean ± SEM. All statistical analyses were performed using R Version 3.4.3 (https://www.r-project.org/). Normality was confirmed using a Shapiro-Wilk test of the residuals. Non-normal data were transformed using Tukey’s Ladder of Powers. If transformation did not result in a normalized population, a Kruskal-Wallis rank-sum test was performed. Analysis of variance was applied to parametric data using a general linear mixed model with repeated measures applied for hormone data. The covariance structure that resulted in the smallest Bayesian information criterion was selected from the mixed models. Significance was declared at P ≤ 0.05.