- Open Access
A transcriptomal analysis of bovine oviductal epithelial cells collected during the follicular phase versus the luteal phase of the estrous cycle
Reproductive Biology and Endocrinology volume 13, Article number: 84 (2015)
Reproductive success depends on a functional oviduct for gamete storage, maturation, fertilization, and early embryonic development. The ovarian-derived steroids estrogen and progesterone are key regulators of oviductal function. The objective of this study was to investigate luteal and follicular phase-specific oviductal epithelial cell function by using microarray-based transcriptional profiling, to increase our understanding of mRNAs regulating epithelial cell processes, and to identify novel genes and biochemical pathways that may be found to affect fertility in the future.
Six normally cycling Angus heifers were assigned to either luteal phase (LP, n = 3) or follicular phase (FP, n = 3) treatment groups. Heifers in the LP group were killed between day 11 and 12 after estrus. Heifers in the FP group were treated with 25 mg PGF2α (Lutalyse, Pfizer, NY) at 8 pm on day 6 after estrus and killed 36 h later. Transcriptional profiling by microarray and confirmation of selected mRNAs by real-time RT-PCR analyses was performed using total RNA from epithelial cells isolated from sections of the ampulla and isthmus collected from LP and FP treatment groups. Differentially expressed genes were subjected to gene ontology classification and bioinformatic pathway analyses.
Statistical one-way ANOVA using Benjamini-hochberg multiple testing correction for false discovery rate (FDR) and pairwise comparison of epithelial cells in the ampulla of FP versus LP groups revealed 972 and 597 transcripts up- and down-regulated, respectively (P < 0.05). Within epithelial cells of the isthmus in FP versus LP groups, 946 and 817 transcripts were up- and down-regulated, respectively (P < 0.05). Up-regulated genes from both ampulla and isthmus were found to be largely involved in cholesterol biosynthesis and cell cycle pathways, while down-regulated genes were found in numerous inflammatory response pathways.
Microarray-based transcriptional profiling revealed phase of the cycle-dependent changes in the expression of mRNA within the epithelium of the oviducts’ ampulla and isthmus.
Reproductive success depends on a functional oviduct for gamete storage and maturation, fertilization, and early embryonic development. Ovarian-derived steroids are well known regulators of oviductal function. Both estrogen and the progesterone receptors are abundant in the bovine oviducts’ mucosal epithelium [1–6], yet our understanding of how the steroidal environment affects the ability of the oviduct to function remains only partially understood. Maturation of gametes and breeding will occur in an estrogen dominant environment, fertilization and early cleavage after the steroidogenic shift and later stages of cleavage and formation of the morula occur within an oviduct exposed to increasing concentrations of circulating progesterone. Increasing our understanding of the steroidal control of oviductal function is critical to the design and implementation of interventions used to manage breeding and the establishment of a pregnancy.
Functionally, the oviduct is divided into two distinct segments; the upper ampulla situated immediately below the ovarian bursa and infundibulum, and the lower isthmus which spans the interval from the ampulla to the uterus. The epithelial mucosa within these two sections consists of ciliated and secretory simple columnar epithelial cells , the ratio of which is regulated by steroid hormones [8, 9]. Steroidal regulation of processes that facilitate sperm binding , sperm release , capacitation  and hyperactivation  are all established in the literature. Epithelial cells are also an active site of biosynthesis and secretion. A 3- to 5-fold increase in the rate of oviductal secretions can be expected around the time of estrus  and amino acids including glycine, glutamate, aspartate, alanine and lysine are all found in higher concentrations in oviductal fluid than in peripheral plasma . Overall, steroid-regulated epithelial cell secretions can be considered an important mediator of the microenvironment that facilitates gamete and zygote health and early development.
While major efforts have been directed at investigating the ovary and uterus at the transcriptional level [16–18], less extensive investigation has been directed towards the oviduct. Suppressive subtractive hybridization was used in the detailed study of Bauersachs et al., , in which changes in gene expression within the oviductal epithelium were determined in heifers killed on the morning of estrus or 12 days thereafter; our report expanding on their analysis with the use of more current transcriptomal profiling technologies and the determination of spatial differences between the ampulla and the isthmus.
The objective of this descriptive study was therefore to determine global oviductal epithelial cell gene expression profiles during the follicular and luteal phases of the estrous cycle. Specifically, microarray-based transcriptional profiling was used to identify spatial and phase of the cycle-dependent changes in mRNA expression in epithelial cells isolated from the ampulla and isthmus, with the overall goals of increasing our understanding of epithelial processes and identifying novel genes that may be identified as key regulators of fertility in the future. Our results must be interpreted, however, with the knowledge that this analysis does not extract potential spatial (isthmus to ampulla and visa versa) signaling mechanisms that could affect oviductal epithelial cell gene expression profiles independent to phase of the estrous cycle. Given the size of the dataset generated by this analysis, our approach to this descriptive study is not to provide a detailed discussion of genes or processes affected by phase of the estrous cycle, but to summarize the results generated, provide our bioinformatic analyses (as Tables and Supplemental Tables) and make available our data for further analysis by others. The microarray raw data (*.cel files) collected with the GCOS software, plus the RMA-normalized and log2 transformed transcript data (Park Genomics Suite ), have been deposited into the Gene Expression Omnibus (National Center for Biotechnology Information, http://www.ncbi.nlm.nih.gov/geo) as accession number GSE63969.
Animals and tissue collection
Animal procedures involved in this study were approved by the University of Kentucky Animal Care and Use Committee. Six normally cycling Angus heifers were used and at least one spontaneous ovulation was observed in each animal prior to being included in the study. All animals were monitored for behavioral estrus (Day 0) and examined a minimum of every other day by trans-rectal ultrasonography throughout the study period, as described before [21, 22]. After ovulation was confirmed, heifers were assigned to either luteal phase (LP) or follicular phase (FP) treatment groups. Heifers in the LP group (n = 3) were killed between day 11 and 12 after estrus. The rationale for collection of LP oviducts on Day 11 or 12 post-estrus was to collect samples under a stable, high progestogenic environment, at a time when the oviduct is undergoing progesterone-dependent remodeling and repair. Heifers in the FP group (n = 3) were treated with 25 mg PGF2α (Lutalyse, Pfizer, New York, NY) at 8 pm on day 6 after estrus and killed 36 h later. This is an established protocol in which the dominant follicle of the first follicular wave of the estrous cycle is induced to differentiate into a preovulatory follicle [21–24]. The diameter of the preovulatory follicle (FP group) and the corpus luteum (LP group) prior to retrieval of the tissues was 14.6 +/− 1.0 mm and 21.5 +/− 0.8 mm, respectively. Stage of the cycle was confirmed by visual appearance of the ovaries collected at slaughter.
Heifers were killed by stunning with a captive-bolt pistol and exsanguinated in the Department of Animal and Food Sciences’ USDA approved Meat Science Lab. Immediately after exsanguination, the oviduct ipsilateral to the corpus luteum (LP group) or preovulatory follicle (FP group) was dissected free from surrounding connective tissue. Epithelial cells were isolated from small sections of the ampulla and isthmus by using a well-established technique for the collection of bovine oviductal epithelial cells [25–27]. Briefly, sections of ampulla and isthmus were gently squeezed with fine forceps under a dissecting microscope to separate epithelial cells from residual stroma. Epithelial cells were briefly centrifuged to form a pellet and then snap-frozen in liquid N2 for later extraction of RNA. Small sections of the ampulla and isthmus were also fixed for 24 h in Bouin’s fixative (Sigma-Aldrich, St. Louis, MO) then processed, sectioned and stained by the University of Kentucky Imaging Facility. Tissues were processed through an ascending series of graded ethyl alcohols, xylene and paraffin then embedded in paraffin. Embedded tissues were sectioned on a microtome at 5 μm, floated on to SuperFrost Plus slides then heat-fixed for a minimum of 2 h at 48 °C on a slide warmer. Hematoxylin and Eosin staining was performed using a ThermoShandon GLX Slide Stainer with Hematoxylin (EKI hematoxylin solution Gills III) and Eosin Y.
RNA Extraction and analysis
Epithelial cells were isolated from the ampulla and isthmus of the oviduct ispilateral to the ovary bearing the preovulatory follicle (FP group) or the corpus luteum (LP group). Total RNA was extracted from each sample of epithelial cells using TRIzol reagent (Invitrogen Corporation, Carlsbad, CA, USA) and purified with RNeasy columns (QIAGEN, Valencia, CA) according to the manufacturer’s instructions. RNA quality was analyzed by determining the RNA integrity number (RIN) using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA) at the University of Kentucky Microarray Core Facility. RNA integrity numbers were greater than 9.2 and 28S/18S rRNA absorbance ratios greater than 1.5 for all samples. RNA concentration was then determined via spectrophotometry using a NanoDrop 2000 (Thermo Fischer Scienctific-NanoDrop products, Wilmington, DE, USA). Spectrophotometry results revealed 260/280 absorbance ratios greater than 1.95 and 260/230 absorbance ratios greater than 1.5 for all samples.
The Bovine gene 1.0 ST array (GeneChip; Affymetrix, Inc., Santa Clara, CA) was used. Microarray analysis was conducted according to the manufacturer’s instructions at the University of Kentucky Microarray Core Facility, as described before [28, 29]. Briefly, RNA (3 μg/sample) was reverse transcribed to cDNA using primers containing T7 RNA polymerase, so that the resulting cDNA contained the T7 sequence. In-vitro transcription was then used for the preparation and labeling of cRNA. The biotinylated cRNA were further fragmented and used as probes to hybridize the GeneChips in the GeneChip Hybridization Oven 640, using 1 chip per RNA sample. The raw expression intensity values generated by microarray hybridization were imported into Partek Genomics Suite 6.6 (Partek Inc., St. Louis, MO). Robust Multiarray Analysis algorithm, quantile normalization, and Median Polish were applied for GeneChip background correction, log base2 transformation, conversion of expression values and probeset summarization [30, 31]. Transcripts were annotated using NetAffx annotation database for the Bovine gene 1.0 ST array and last updated in June 2014.
After data were processed for background adjustment, normalization and log2 transformed, quality of data was assessed using light intensity expression values on a per chip and per gene basis. For statistical analysis, Partek Genomics Suite 6.6 (Partek Inc.) was used to complete an F-test on least-square means to determine significance of each transcript in each comparison. Benjamini-hochberg multiple testing correction for false discovery rate (FDR) was applied and significance set to FDR adjusted P-value < 0.05. A post-hoc pairwise comparison of FP compared to LP epithelial cells from the ampulla and isthmus was completed using Fisher’s Least Significant Difference (LSD) to determine which means differed [32, 33]. Only transcripts with a fold-change value of ≥ 1.5 were included in the results. The raw data (*.cel files), plus the RMA-normalized and log2 transformed transcript data (Park Genomics Suite , have been deposited into the Gene Expression Omnibus (National Center for Biotechnology Information ) as accession number GSE63969 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE63969).
Real-time RT-PCR analysis
Real-time RT-PCR was performed using RNA samples from each tissue and phase of the estrous cycle to determine the expression of mRNA for neurotensin (NTS), binder of sperm 3 (BSP3), lactate dehydrogenase A (LDHA), cyclin D2 (CCND2), early growth response 1 (EGR1) and hydroxysteroid (17β) dehydrogenase 7 (HSD17B7) using an Eppendorf Mastercycler ep realplex 2 system (Eppendorf, Hamburg, Germany) with iQ SYBR Green Supermix (Bio-RAD, Hercules, CA), as described before . Additional validation of the Bovine gene 1.0 ST array has been reported by our laboratory . The following oligonucleotide primer pairs (5’ to 3’) were used: NTS, F: GTG TGG AAA TGT GAC AGA GCA C and R: GGT AGG CTA GAC TTT GCG GT; BSP3, F: ATT CCT GTG GTG TTC CCT CG and R: GCT CAG AGC ATC ACC TTT GC; LDHA, F: CCA ACA TGG CAG CCT TTT CC and R: ACC GCT TTC CAC TGT TCC TT; CCND2, F: CCG ACA ACT CCA TCA AGC CT and R: TGA AGT AGT GGC GCA CAG AG; EGR1, F: AGA AAG TTT GCC AGG AGC GA and R: GGA GGG ACG GAG GAG TAT GT; HSD17B7, F: ACA GCT GAA GGA CTG CTG AC and R: CCA GAC AGT GCT TCT GTT CCA; and 18S, F: CGG GGA GGT AGT GAC GAA A and R: CCG CTC CCA AGA TCC AAC TA.
Briefly, cDNA was synthesized using the SuperScript III 1st Strand Synthesis System (Invitrogen), with 0.5 μg of RNA used for each reverse transcription reaction. Real-time RT-PCR was performed with a total volume of 25 μL per reaction, with each reaction containing 5 μL of cDNA, 1 μL of a 10 μM stock of each primer (forward and reverse), 12.5 μL of 2× SYBR Green PCR Master Mix, and 5.5 μL of nuclease-free water. The typical dissociation curves of these cDNA, plus 18S as the housekeeping gene was confirmed. RT-PCR reactions were run in triplicate and gene expression was analyzed by the 2ΔΔCT method .
Gene ontology and pathway analysis
Differentially expressed transcripts were interrogated for their gene ontology classes using Partek Genomics Suite 6.6 (Partek Inc.). Partek derives biological processes, molecular functions and cellular components from geneontology.org and/or the affymetrix database. GO hierarchies leads to division of the gene list into significant classifications using Fischer’s exact test-right tailed. When the observed number of genes in a GO category is greater than expected, the GO category is enriched. Pathway analysis was completed by importing differentially expressed transcripts into QIAGEN’S Ingenuity Pathway Analysis (IPA, QIAGEN, Redwood City, CA, USA, www.qiagen.com/ingenuity). Ingenuity Pathway Analysis uses data from multiple databases to extrapolate significant pathways and Fischer’s exact test was used to determine significant pathways. Significance was set to P-value < 0.05.
Results and discussion
Real-time RT- PCR analysis of selected transcripts
The effect of stage of the estrous cycle and tissue on the expression of mRNA for NTS, BSP3, LDHA, CCND2, EGR1 and HSD17B7 was performed by real-time RT-PCR. A comparison of the results obtained by RT-PCR and microarray analysis is described in Table 1 as a validation of the microarray platform. Additional validation of this Bovine gene 1.0 ST array by real-time RT-PCR performed within our laboratory has been reported. Overall, RT-PCR revealed the same directional trends in gene expression as microarray analysis, with the magnitude of these changes typically lower after analysis using the microarray platform (ratio compression phenomena), as described by others . For example, the expression of mRNA for NTS was found to increase by 32.4- and 150-fold within epithelial cells of the ampulla, and 21.6- and 88-fold within the epithelial cells of the isthmus, in the FP versus the LP groups by microarray analysis and real-time RT-PCR, respectively.
Microarray quality control and principal component analysis
Box plots revealed mean intensity values were similar across all chips and overlapping histograms indicated the frequency of transcripts at specific intensity values for each chip were similar (Additional file 1: Figure S2 and Additional file 2: Figure S2). Quantification of signal intensity to noise revealed that spatial location and phase of the cycle, not error, accounted for the variation within the data set (mean F Ratio for attribute = 30.78, versus mean F value for error = 1.0). Consistent with this, principal component analysis (PCA), which allows for the visualization of patterns through the distribution of samples to highlight similarities and differences , revealed clear differences between LP and FP groups, as well as tissue-specific differences within the same phase of the estrous cycle (Fig. 1). Total variance (55.5 %) is the cumulative percent of variance accounted for in our datasets based upon eigenvector multivariate analysis. PC#1 (x-axis, 28 %) indicates that the largest proportion of variability is due to phase of the estrous cycle. PC#2 (y-axis, 17.8 %) indicates variability between the ampulla and isthmus. PC#3 (z-axis, 9.7 %) indicates the variability between phase of the estrous cycle and tissues.
It is well known that function of the mammalian oviduct is influenced by estradiol and progesterone, with previous studies confirming the presence of mRNA and protein for estrogen receptor alpha (ERα) and beta (ERβ) as well as the nuclear and membrane progesterone receptors (PR, PGRMC1, PGRMC2) within this organ [6, 38]. After chip normalization, a statistical one-way ANOVA and pairwise comparison (LSD test) of gene expression in epithelial cells of the ampulla and isthmus in the FP and LP groups was performed to generate a list of 2374 differentially expressed genes (DEGs). Table 2 indicates the number of DEGs in each pairwise comparison. As expected, there were a large number of DEGs identified by FP (1563 DEGs) and LP (1758 DEGs) pairwise comparisons. Among these DEGs, 947 DEGs overlapped between contrasts (Fig. 2). By pairwise comparison, 616 DEGs were exclusive to the ampulla and 811 DEGs were exclusive to the isthmus. The complete list of these transcripts is provided in Additional file 3: Table S1 (DEGs exclusive to the ampulla or isthmus) and Additional file 4: Table S2 (DEGs in common among the ampulla and isthmus).
DEGs in epithelial cells of the ampulla
The ampulla is situated immediately below the ovarian bursa and infundibulum and is the site of fertilization. Histological evaluations have shown the ampulla to have an elaborate, extensively folded epithelial layer adjacent to a thin layer of smooth muscle (Fig. 3a). Evaluation of the oviducts of cattle [14, 39, 40] confirm that this is a dynamic tissue with multiple, ongoing biological processes that include cell proliferation  and tissue remodeling . The most highly up- and down-regulated DEGs within epithelial cells of the ampulla in FP versus LP groups are provided as Table 3.
As indicated above, a detailed discussion of the role for individual genes identified by this analysis is not our objective, however examples of consistency among our profiling results and more targeted analyses performed by others is warranted. Of the up-regulated genes, the expression of mRNA for neurotensin (NTS) was found to increase by 32.4-fold within epithelial cells of the ampulla in the FP versus the LP groups. Neurotensin has multiple functions and evidence suggests that NTS plays a role in gamete and embryo transport within the oviduct . Receptors for NTS are expressed on spermatozoa, and increasing NTS administration facilitates sperm protein tyrosine phosphorylation, which is a measure of sperm capacitation . Furthermore, the acrosome reaction is promoted in capacitated spermatozoa in the presence of increasing concentrations of NTS , consistent with the increased expression of NTS observed in the follicular phase herein. In contrast to this established role of NTS, and consistent with our objective to identify novel genes that may be revealed as critical mediators of function in the future, the transcription factor early growth response 1 (EGR1) was observed to display the highest fold-change (9.1-fold) among down-regulated genes within epithelial cells of the ampulla in the FP versus the LP groups. A gonadotropin-dependent induction of EGR1 has been reported prior to ovulation in bovine follicles , however no reports to date are apparent on the function, or importance, of this transcription factor in the epithelium of the oviducts’ ampulla.
Up- and down-regulated DEGs were analyzed for enriched gene ontology classifications. Gene ontology analysis of up-regulated DEGs within epithelial cells of the ampulla in the FP versus the LP groups resulted in 117 significant biological processes, 65 significant cellular components, and 46 significant molecular functions (P < 0.05,
Additional file 5: Table S3). Cell cycle, cholesterol biosynthetic process, cell division, mitosis and protein folding were the top biological processes identified, which is not surprising considering the cell proliferation and secretory activity required to prepare the ampulla for the arrival of the gametes. With respect to cholesterol biosynthesis, steroid-dependent effects on the oviduct are documented in the literature. Cholesterol will affect the ability of spermatozoa to fertilize an oocyte , with the process of capacitation well established to require efflux of cholesterol from the plasma membrane of spermatozoa (reviewed in ). High density lipoproteins are elevated in bovine oviductal fluid during the follicular phase of the estrous cycle , however the synthesis and release of cholesterol by oviductal epithelial cells appears to be greater under a progesterone dominant environment. Esterified-cholesterol containing lipid droplets from oviductal epithelial cells were observed in greater numbers when collected from luteal-phase cows  and the concentration of cholesterol was increased in isthmic but not ampullary oviductal fluid collected from the luteal phase versus non-luteal phase animals . Further analysis of the biochemical relationship among the DEGs expressed within these top cholesterol-associated pathways may increase our understanding of the role of cholesterol during the processes of capacitation and fertilization, critical events that occur within the oviduct. Top cellular components of up-regulated DEGs within epithelial cells of the ampulla in the FP versus the LP groups included endoplasmic reticulum lumen and membrane, cytoplasm, and mitochondrion. Top molecular functions involved protein disulfide isomerase activity, FK506 binding, peptidyl-prolyl cis-trans isomerase activity, oxidoreductase activity, and dolichyl-diphosphooligosaccharide-protein glycotransferase activity.
Analysis of down-regulated DEGs within epithelial cells of the ampulla in the FP versus the LP revealed 118 significant biological processes, 18 significant cellular components, and 49 significant molecular functions (P < 0.05, Additional file 6: Table S4). The top biological processes were the innate immune response, response to nicotine, myoblast proliferation, negative regulation of MAP kinase activity, and bone morphogenetic protein (BMP) signaling pathway. With respect to immune responses, it can be postulated that increases in estradiol during the follicular phase decrease the induction of pro-inflammatory factors, with the results of our study consistent with an investigation on the effect of ovarian steroids on lipopolysaccharide (LPS)-induced responses in bovine oviductal epithelial cells in vitro . In their report, estradiol reversed the effect of LPS on pro-inflammatory gene expression and we have previously demonstrated ESR1-dependent as well as cyclic changes in the expression of the hematopoetic form of prostaglandin D synthase, a putative regulator of inflammation, in the mouse oviduct . Given that the oviductal epithelium must continually repair itself from any damage caused by exposure to gametes, seminal fluids and post-ovulatory follicular debris, steroid-dependent changes in inflammatory response mechanisms can be considered physiologically important biological processes. Lastly, top cellular components associated with the down-regulated DEGs were plasma membrane, cytoplasm, extracellular region, cytosol, and tight junction and top molecular functions were pancreatic ribonuclease activity, transmembrane signaling receptor activity, protein-L-isoaspartate (D-aspartate) O-methyltransferase activity, serine hydrolase activity, and phosphatidylinositol-3,4-bisphosphate binding.
Canonical pathway analysis of DEGs within epithelial cells of the ampulla in the FP versus the LP was then determined using QIAGEN’S Ingenuity Pathway Analysis (IPA, QIAGEN, Redwood City, www.qiagen.com/ingenuity). Top pathways up-regulated in the follicular phase largely reflected cholesterol biosynthesis (Superpathway of Cholesterol Biosynthesis, Cholesterol Biosynthesis I, Cholesterol Biosynthesis II (via 24,25-dihydrolanosterol), Cholesterol Biosynthesis III (via Desmosterol)), and Oxidative Phosphorylation (P < 0.05, Additional file 7: Table S5 and Additional file 8: Figure S3a), which was consistent with the results of the gene ontology analysis discussed above. The top pathways for down-regulated genes included Hepatic Fibrosis/Hepatic Stellate Cell Activation, Role of Pattern Recognition Receptors in Recognition of Bacteria and Viruses, Colorectal Cancer Metastasis Signaling, Ovarian Cancer Signaling, and Role of Macrophages, Fibroblasts and Endothelial Cells in Rheumatoid Arthritis (P < 0.05, Additional file 9: Table S6 and Additional file 8: Figure S3b). Again, this canonical pathway analysis is consistent with steroid-dependent regulation of inflammation and immune responses extracted by the gene ontology analysis. Worthy to note, given the increasing interest in oviductal epithelial cells as progenitors for ovarian cancer [50–52], down-regulation of the DEGs in the ovarian cancer signaling canonical pathway (MMP7, ARRB1, FZD4, FGF9, SMO, FIGF, CCND1, PDGFC, FZD7, EGFR and BCL2) could provide useful clues that can advance that important field of study. In addition, potential effects of the presence of embryos to the oviductal epithelium should not be overlooked. We observed down-regulation of 7 DEGs, including that of insulin-like growth factor binding protein 3 (IGFBP3) in epithelial cells of the ampulla in FP versus LP groups (Additional file 4: Table S2), whose expression in primary bovine oviductal epithelial cells in vitro is increased by the addition of embryos to their culture .
DEGs in epithelial cells of the isthmus
Similar to the ampulla, the mucosa of the isthmus appears well defined in histological evaluations (Fig. 3b). However, in contrast to the ampulla, the mucosa of the isthmus is located adjacent to a prominent musculature with this smooth musculature playing a key role in the movement of gametes to the site of fertilization and passage of the early embryo to the uterus [54–56]. Several of the DEGs within epithelial cells of the isthmus in the FP versus the LP were similar to those identified within the ampulla, however, there were also many DEGs unique to this spatial location (Additional file 3: Tables S1 and Additional file 4: Table S2). Among the up-regulated DEGs in the epithelium of the isthmus (Table 4), the increased expression of mRNA encoding phospholipase A2, group IVD (cytosolic) (PLA2G4D) and phospholipase A2, group IVF (PLA2G4F) is interesting, especially when considering the key role for prostaglandins within the oviduct [55, 57, 58] and the noted regulation of prostaglandin synthesis by ovarian steroids as reported by us  and others [59, 60]. Both PLA2G4D and PLA2G4F belong to the phospholipase A2 (PLA2) family, group 4. The PLA2 family of enzymes catalyzes the hydrolysis of phospholipids to liberate free fatty acids, among other molecules. Of these free fatty acids, arachidonic acid released by PLA2 enzymes acts as the precursor for the synthesis of prostaglandins [61, 62]. Consistent with our results, in ovariectomized rabbits, PLA2 activity in epithelial cells of the ampulla is reported to be increased after treatment of rabbits with estradiol . Given that PLA2G4D and PLA2G4F do not appear to be described within the oviduct, our results can be considered as new information on the local regulation of phospholipases and potentially prostaglandin secretion within this organ. Of the down-regulated DEGs expressed by epithelial cells of the isthmus, the expression of binder of sperm 3 (BSP3) is also identified as a novel transcript. BSP3 is reported to be secreted by seminal vesicles and binds to sperm , and identified within the epithelial cells of the isthmus herein. Again, the potential to modify the epithelial cell transcriptome by the presence of embryos must be acknowledged. We observed down-regulation of 10 DEGs, including that of the apoptosis regulator XIAP associated factor-1 (XAF1) in epithelial cells of the isthmus in FP versus LP groups (Additional file 4: Table S2), whose expression has been reported to be increased by the addition of embryos to cultures of primary bovine oviductal epithelial cells .
Gene ontology classification of up-regulated DEGs within epithelial cells of the isthmus in the FP versus the LP revealed 97 significant biological processes with top biological processes being protein folding, cell cycle, cell division, mitosis, and electron transport chain (P < 0.05, Additional file 10: Table S7). Similar to the ampulla, there were 57 significant cellular components affected by stage of the cycle within epithelial cells of the isthmus, including endoplasmic reticulum lumen and membrane, cytoplasm, and mitochondrion and 46 significant molecular functions. Gene ontology analysis of the down-regulated DEGs within epithelial cells of the isthmus revealed 101 significant biological processes, 26 significant cellular components, and 45 significant molecular functions (P < 0.05, Additional file 11: Table S8). Top biological processes were negative regulation of fibroblast growth factor receptor signaling pathway, calcium-independent cell-cell adhesion, inactivation of MAPK activity, brown fat cell differentiation, and SMAD protein signal transduction. Top cellular components were similar to those observed within the ampulla and included cytoplasm, plasma membrane cytosol, extracellular region, and tight junction.
Again, QIAGEN’S Ingenuity Pathway Analysis (IPA, QIAGEN, Redwood City, www.qiagen.com/ingenuity) was used to identify the canonical pathways affected by stage of the cycle within epithelial cells of the isthmus. Up-regulated DEGs within epithelial cells of the isthmus in FP versus LP were oxidative phosphorylation, mitochondrial dysfunction, superpathway of cholesterol biosynthesis, superpathway of geranylgeranyldiphosphate biosynthesis I (via mevalonate), and melavonate pathway I (P < 0.05, Additional file 12: Table S9 and Additional file 13: Figure S4A). Similar to that observed in the ampulla, cholesterol biosynthesis stands out as a key pathway with gene expression increased under an estrogen dominant environment. Consistency in response was also observed among down-regulated DEGs; top pathways for down-regulated DEGs included molecular mechanisms of cancer, basal cell carcinoma signaling, role of osteoblasts, osteoclasts and chondrocytes in rheumatoid arthritis, Role of Macrophages, Fibroblasts and Endothelial Cells in Rheumatoid Arthritis, and Corticotropin Releasing Hormone Signaling, with the complete listing provided in Additional file 14: Table S10 and Additional file 13: Figure S4B (P < 0.05).
At the genome level, estrous cycle stage-dependent effects on epithelial cell gene expression are not well defined. The current study therefore investigated changes in the expression of mRNA within the epithelium of the ampulla and isthmus of the bovine oviduct during the luteal and follicular phases of the estrous cycle. This transcriptomal profiling analysis was performed to increase our understanding of gene expression and potentially epithelial cell processes important for oviductal function and fertility, and to identify novel mRNA that may prove critical for fertility after analysis in the future.
Analysis of variance
Differentially expressed genes
False discovery rate
Bylander A, Nutu M, Wellander R, Goksor M, Billig H, Larsson DG. Rapid effects of progesterone on ciliary beat frequency in the mouse fallopian tube. Reprod Biol Endocrinol. 2010;8:48.
Wanggren K, Stavreus-Evers A, Olsson C, Andersson E, Gemzell-Danielsson K. Regulation of muscular contractions in the human Fallopian tube through prostaglandins and progestagens. Hum Reprod. 2008;23:2359–68.
Parada-Bustamante A, Orihuela PA, Rios M, Cuevas CA, Orostica ML, Velasquez LA, et al. A non-genomic signaling pathway shut down by mating changes the estradiol-induced gene expression profile in the rat oviduct. Reproduction. 2010;139:631–44.
Shao R, Nutu M, Karlsson-Lindahl L, Benrick A, Weijdegard B, Lager S, et al. Downregulation of cilia-localized Il-6R alpha by 17beta-estradiol in mouse and human fallopian tubes. Am J Physiol Cell Physiol. 2009;297:C140–51.
Valle GR, Cassali GD, Nogueira JC, Castro AC, Reis AM, Cardoso FM, et al. Nuclear estrogen and progesterone receptors in the oviduct of heifers under natural and superovulated estrous cycles. Anim Reprod Sci. 2007;101:28–37.
Saint-Dizier M, Sandra O, Ployart S, Chebrout M, Constant F. Expression of nuclear progesterone receptor and progesterone receptor membrane components 1 and 2 in the oviduct of cyclic and pregnant cows during the post-ovulation period. Reprod Biol Endocrinol. 2012;10:76.
Yaniz JL, Lopez-Gatius F, Santolaria P, Mullins KJ. Study of the functional anatomy of bovine oviductal mucosa. Anat Rec. 2000;260:268–78.
Shirley B, Reeder RL. Cyclic changes in the ampulla of the rat oviduct. J Exp Zool. 1996;276:164–73.
Reeder RL, Shirley B. Deciliation in the ampulla of the rat oviduct and effects of estrogen on the process. J Exp Zool. 1999;283:71–80.
Bureau M, Bailey JL, Sirard MA. Binding regulation of porcine spermatozoa to oviductal vesicles in vitro. J Androl. 2002;23:188–93.
Hunter RH. Sperm release from oviduct epithelial binding is controlled hormonally by peri-ovulatory graafian follicles. Mol Reprod Dev. 2008;75:167–74.
Ded L, Dostalova P, Dorosh A, Dvorakova-Hortova K, Peknicova J. Effect of estrogens on boar sperm capacitation in vitro. Reprod Biol Endocrinol. 2010;8:87.
Fujinoki M. Suppression of progesterone-enhanced hyperactivation in hamster spermatozoa by estrogen. Reproduction. 2010;140:453–64.
Carlson D, Black DL, Howe GR. Oviduct secretion in the cow. J Reprod Fertil. 1970;22:549–52.
Hugentobler SA, Diskin MG, Leese HJ, Humpherson PG, Watson T, Sreenan JM, et al. Amino acids in oviduct and uterine fluid and blood plasma during the estrous cycle in the bovine. Mol Reprod Dev. 2007;74:445–54.
Jo M, Gieske MC, Payne CE, Wheeler-Price SE, Gieske JB, Ignatius IV, et al. Development and application of a rat ovarian gene expression database. Endocrinology. 2004;145:5384–96.
Winuthayanon W, Hewitt SC, Korach KS: Uterine Epithelial Cell Estrogen Receptor Alpha-Dependent and -Independent Genomic Profiles That Underlie Estrogen Responses in Mice. Biol Reprod. 2014;91(5):110,1-10.
Minten MA, Bilby TR, Bruno RG, Allen CC, Madsen CA, Wang Z, et al. Effects of fertility on gene expression and function of the bovine endometrium. PLoS One. 2013;8, e69444.
Bauersachs S, Rehfeld S, Ulbrich SE, Mallok S, Prelle K, Wenigerkind H, et al. Monitoring gene expression changes in bovine oviduct epithelial cells during the oestrous cycle. J Mol Endocrinol. 2004;32:449–66.
Partek D. Partek documentation: turning data into discovery. St Louis: Partek Incorporated; 2009.
Bridges PJ, Fortune JE. Regulation, action and transport of prostaglandins during the periovulatory period in cattle. Mol Cell Endocrinol. 2007;263:1–9.
Bridges PJ, Komar CM, Fortune JE. Gonadotropin-induced expression of messenger ribonucleic acid for cyclooxygenase-2 and production of prostaglandins E and F2alpha in bovine preovulatory follicles are regulated by the progesterone receptor. Endocrinology. 2006;147:4713–22.
Komar CM, Berndtson AK, Evans AC, Fortune JE. Decline in circulating estradiol during the periovulatory period is correlated with decreases in estradiol and androgen, and in messenger RNA for p450 aromatase and p450 17alpha-hydroxylase, in bovine preovulatory follicles. Biol Reprod. 2001;64:1797–805.
Bridges PJ, Fortune JE. Characteristics of developing prolonged dominant follicles in cattle. Domest Anim Endocrinol. 2003;25:199–214.
Pedersen ME, Ozdas OB, Farstad W, Tverdal A, Olsaker I. Effects of bovine oviduct epithelial cells, fetal calf serum and bovine serum albumin on gene expression in single bovine embryos produced in the synthetic oviduct fluid culture system. Reprod Fertil Dev. 2005;17:751–7.
Rottmayer R, Ulbrich SE, Kolle S, Prelle K, Neumueller C, Sinowatz F, et al. A bovine oviduct epithelial cell suspension culture system suitable for studying embryo-maternal interactions: morphological and functional characterization. Reproduction. 2006;132:637–48.
Kodithuwakku SP, Miyamoto A, Wijayagunawardane MP. Spermatozoa stimulate prostaglandin synthesis and secretion in bovine oviductal epithelial cells. Reproduction. 2007;133:1087–94.
Bridges PJ, Jeoung M, Kim H, Kim JH, Lee DR, Ko C, et al. Methodology matters: IVF versus ICSI and embryonic gene expression. Reprod Biomed Online. 2011;23:234–44.
Bridges PJ, Jeoung M, Shim S, Park JY, Lee JE, Sapsford LA, et al. Hematopoetic prostaglandin D synthase: an ESR1-dependent oviductal epithelial cell synthase. Endocrinology. 2012;153:1925–35.
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, et al. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003;4:249–64.
D P. Partek documentation: turning data into discovery. St. Louis, MO: Partek Inc; 2009.
Eisenhart C. The assumptions underlying the analysis of variance. Biometrics. 1947;3:1–21.
Tamhane CA, Dunlop DD. Statistics and Data Analysis from Elementary to Intermediate. Prentice Hall. 2000;473–474.
Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–10.
Cerny KL, Garbacik S, Skees C, Burris WR, Matthews JC, Bridges PJ. Gestational form of Selenium in Free-Choice Mineral Mixes Affects Transcriptome Profiles of the Neonatal Calf Testis, Including those of Steroidogenic and Spermatogenic Pathways. Biol Trace Elem Res. 2015. doi:10.1007/s12011-015-0386-4.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods. 2001;25:402–8.
Wang Y, Barbacioru C, Hyland F, Xiao W, Hunkapiller KL, Blake J, et al. Large scale real-time PCR validation on gene expression measurements from two commercial long-oligonucleotide microarrays. BMC Genomics. 2006;7:59.
Ulbrich SE, Kettler A, Einspanier R. Expression and localization of estrogen receptor alpha, estrogen receptor beta and progesterone receptor in the bovine oviduct in vivo and in vitro. J Steroid Biochem Mol Biol. 2003;84:279–89.
Kowsar R, Hambruch N, Liu J, Shimizu T, Pfarrer C, Miyamoto A. Regulation of innate immune function in bovine oviduct epithelial cells in culture: the homeostatic role of epithelial cells in balancing Th1/Th2 response. J Reprod Dev. 2013;59:470–8.
Odau S, Gabler C, Holder C, Einspanier R. Differential expression of cyclooxygenase 1 and cyclooxygenase 2 in the bovine oviduct. J Endocrinol. 2006;191:263–74.
Gabler C, Killian GJ, Einspanier R. Differential expression of extracellular matrix components in the bovine oviduct during the oestrous cycle. Reproduction. 2001;122:121–30.
Reinecke M. Neurotensin in the human fallopian tube: immunohistochemical localization and effects of synthetic neurotensin on motor activity in vitro. Neurosci Lett. 1987;73:220–4.
Hiradate Y, Inoue H, Kobayashi N, Shirakata Y, Suzuki Y, Gotoh A, et al. Neurotensin enhances sperm capacitation and acrosome reaction in mice. Biol Reprod. 2014;91:53.
Sayasith K, Brown KA, Lussier JG, Dore M, Sirois J. Characterization of bovine early growth response factor-1 and its gonadotropin-dependent regulation in ovarian follicles prior to ovulation. J Mol Endocrinol. 2006;37:239–50.
Diaz-Fontdevila M, Bustos-Obregon E. Cholesterol and polyunsaturated acid enriched diet: effect on kinetics of the acrosome reaction in rabbit spermatozoa. Mol Reprod Dev. 1993;35:176–80.
Sheriff DS, Ali EF. Perspective on plasma membrane cholesterol efflux and spermatozoal function. J Hum Reprod Sci. 2010;3:68–75.
Ehrenwald E, Foote RH, Parks JE. Bovine oviductal fluid components and their potential role in sperm cholesterol efflux. Mol Reprod Dev. 1990;25:195–204.
Henault MA, Killian GJ. Composition and morphology of lipid droplets from oviduct epithelial cells. Anat Rec. 1993;237:466–74.
Grippo AA, Anderson SH, Chapman DA, Henault MA, Killian GJ. Cholesterol, phospholipid and phospholipase activity of ampullary and isthmic fluid from the bovine oviduct. J Reprod Fertil. 1994;102:87–93.
Reade CJ, McVey RM, Tone AA, Finlayson SJ, McAlpine JN, Fung-Kee-Fung M, et al. The fallopian tube as the origin of high grade serous ovarian cancer: review of a paradigm shift. J Obstet Gynaecol Can. 2014;36:133–40.
Sherman-Baust CA, Kuhn E, Valle BL, Shih Ie M, Kurman RJ, Wang TL, et al. A genetically engineered ovarian cancer mouse model based on fallopian tube transformation mimics human high-grade serous carcinoma development. J Pathol. 2014;233:228–37.
van der Horst PH, van der Zee M, Heijmans-Antonissen C, Jia Y, DeMayo FJ, Lydon JP, et al. A mouse model for endometrioid ovarian cancer arising from the distal oviduct. Int J Cancer. 2014;135:1028–37.
Schmaltz-Panneau B, Cordova A, Dhorne-Pollet S, Hennequet-Antier C, Uzbekova S, Martinot E, et al. Early bovine embryos regulate oviduct epithelial cell gene expression during in vitro co-culture. Anim Reprod Sci. 2014;149:103–16.
Wijayagunawardane MP, Miyamoto A. Tumor necrosis factor alpha system in the bovine oviduct: a possible mechanism for embryo transport. J Reprod Dev. 2004;50:57–62.
Wijayagunawardane MP, Miyamoto A, Taquahashi Y, Gabler C, Acosta TJ, Nishimura M, et al. In vitro regulation of local secretion and contraction of the bovine oviduct: stimulation by luteinizing hormone, endothelin-1 and prostaglandins, and inhibition by oxytocin. J Endocrinol. 2001;168:117–30.
Talo A. Myoelectrical activity and transport of unfertilized ova in the oviduct of the mouse in vitro. J Reprod Fertil. 1980;60:53–8.
Wijayagunawardane MP, Miyamoto A, Cerbito WA, Acosta TJ, Takagi M, Sato K. Local distributions of oviductal estradiol, progesterone, prostaglandins, oxytocin and endothelin-1 in the cyclic cow. Theriogenology. 1998;49:607–18.
Gauvreau D, Moisan V, Roy M, Fortier MA, Bilodeau JF. Expression of prostaglandin E synthases in the bovine oviduct. Theriogenology. 2010;73:103–11.
Szostek AZ, Siemieniuch MJ, Deptula K, Woclawek-Potocka I, Majewska M, Okuda K, et al. Ovarian steroids modulate tumor necrosis factor-alpha and nitric oxide-regulated prostaglandin secretion by cultured bovine oviductal epithelial cells. Domest Anim Endocrinol. 2011;41:14–23.
Perez Martinez S, Hermoso M, Farina M, Ribeiro ML, Rapanelli M, Espinosa M, et al. 17-beta-Estradiol upregulates COX-2 in the rat oviduct. Prostaglandins Other Lipid Mediat. 2006;80:155–64.
Balsinde J, Winstead MV, Dennis EA. Phospholipase A(2) regulation of arachidonic acid mobilization. FEBS Lett. 2002;531:2–6.
Dennis EA. Diversity of group types, regulation, and function of phospholipase A2. J Biol Chem. 1994;269:13057–60.
Morishita T, Nozaki M, Sano M, Yokoyama M, Nakamura G, Nakano H. Changes in phospholipase A2 activity of the rabbit ampullary epithelium by ovarian steroids. Prostaglandins Leukot Essent Fatty Acids. 1993;48:315–8.
Manjunath P, Sairam MR, Uma J. Purification of four gelatin-binding proteins from bovine seminal plasma by affinity chromatography. Biosci Rep. 1987;7:231–8.
This work was supported by National Institutes of Health Grant P01 HD071875, the University of Kentucky, and Kentucky Agricultural Experiment Station. The information reported in this paper (publication no. 15-07-075 ) is part of a project of the Kentucky Agricultural Experiment Station and is published with approval of the Director.
The authors declare that they have no competing interests.
K.C., L.A. and P.B. designed the study and collected the samples, K.C., E.G. and A.J.W. performed the analyses, K.C. and P.B. wrote the manuscript. All authors read and approved the final manuscript.
Supplementary Figure 1. Box plot of the log2 expression signal for each sample (microarray chip).
Supplementary Figure 2. Overlapping sample signal intensity histogram indicating the frequency of transcripts at specific signal intensity values.
Supplementary Table 1. Differentially expressed transcripts between follicular and luteal phases of the estrous cycle in epithelial cells exclusive to the ampulla and isthmus (FDR adjusted P < 0.05).
Supplementary Table 2. Differentially expressed transcripts in common between follicular and luteal phases of the estrous cycle in epithelial cells of the ampulla and isthmus (FDR adjusted P < 0.05).
Supplementary Table 3. Enriched gene ontology classifications from up-regulated DEGS within epithelial cells of the ampulla in follicular versus luteal phase groups (P < 0.05).
Supplementary Table 4. Enriched gene ontology classifications from down-regulated DEGS within epithelial cells of the ampulla in follicular versus luteal phase groups (P < 0.05).
Supplementary Table 5. Canonical pathway analysis of up-regulated DEGS within epithelial cells of the ampulla in the follicular versus luteal phase groups (P < 0.05).
Supplementary Figure 3. Top 6 Canonical pathways from up- and down-regulated differentially expressed genes within epithelial cells of the ampulla in the follicular versus luteal phase. Ingenuity Pathway Analysis software was used to determine significant pathways based on the number of significant genes expressed within the pathway using Fischer’s Exact test (P <0.05). A) Ampulla: Up-regulated pathways in the follicular versus luteal phases. B) Ampulla: Down-regulated pathways in the follicular versus luteal phases.
Supplementary Table 6. Canonical pathway analysis of down-regulated DEGS within epithelial cells of the ampulla in the follicular versus luteal phase groups (P < 0.05).
Supplementary Table 7. Enriched gene ontology classifications from up-regulated DEGS within epithelial cells of the isthmus in follicular versus luteal phase groups (P < 0.05).
Supplementary Table 8. Enriched gene ontology classifications from down-regulated DEGS within epithelial cells of the isthmus in follicular versus luteal phase groups (P < 0.05).
Supplementary Table 9. Canonical pathway analysis of up-regulated DEGS within epithelial cells of the isthmus in the follicular versus luteal phase groups (P < 0.05).
Supplementary Figure 4. Top 6 Canonical pathways from up- and down-regulated differentially expressed genes within epithelial cells of the isthmus in the follicular versus luteal phase. Ingenuity Pathway Analysis software was used to determine significant pathways based on the number of significant genes expressed within the pathway using Fischer’s Exact test (P <0.05). A) Isthmus: Up-regulated pathways in the follicular versus luteal phases. B) Isthmus: Down-regulated pathways in the follicular versus luteal phases.
Supplementary Table 10. Canonical pathways of down-regulated DEGs within epithelial cells of the isthmus in the follicular versus luteal phase groups (P < 0.05).
About this article
Cite this article
Cerny, K.L., Garrett, E., Walton, A. et al. A transcriptomal analysis of bovine oviductal epithelial cells collected during the follicular phase versus the luteal phase of the estrous cycle. Reprod Biol Endocrinol 13, 84 (2015). https://doi.org/10.1186/s12958-015-0077-1
- Oviduct epithelial cells