Reproduction   citetrack
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS  

Reproduction (2007) 134 123-135
DOI: 10.1530/REP-07-0387
Copyright © 2007 Society for Reproduction and Fertility
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Tables
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Chen, Y.
Right arrow Articles by Roberts, R M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Chen, Y.
Right arrow Articles by Roberts, R M.

RESEARCH

A microarray analysis for genes regulated by interferon-{tau} in ovine luminal epithelial cells

Yizhen Chen1, Eric Antoniou1, Zhilin Liu1, Leonard B Hearne2 and R Michael Roberts1,3

1 Division of Animal Sciences, University of Missouri-Columbia, Missouri 65211, USA, 2 Department of Statistics, University of Missouri-Columbia, Missouri 65211, USA and 3 Department of Biochemistry, University of Missouri-Columbia, Missouri 65211, USA

Correspondence should be addressed to R M Roberts, 240b CS Bond Life Sciences Center, University of Missouri-Columbia, 1201 Rollins Street, Columbia, Missouri 65211, USA; Email: robertsrm{at}missouri.edu


    Abstract
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions
 Acknowledgements
 References
 
Interferon-{tau} (IFNT) is released by preimplantation conceptuses of ruminant species and prepares the mother for pregnancy. Although one important function is to protect the corpus luteum from the luteolytic activity of prostaglandin-F 2{alpha}, IFNT most likely regulates a range of other physiological processes in endometrium. Here, an immortalized cell line from ovine uterine luminal epithelial cells was treated with IFNT for either 8 or 24 h. RNA was subjected to cDNA microarray analysis, with RNA from untreated cells as the reference standard. Of 15 634 genes, 1274 (8%) were IFNT responsive at P<0.01 and 585 at P<0.001 to at least one treatment. Of the latter, 356 were up-regulated and 229 down-regulated. Increasing IFNT concentrations from 10 ng/ml to 10 µg/ml had minor effects, and most genes up- or down-regulated at 8 h were regulated similarly at 24 h. Although IFNT influences many genes implicated in antiviral activity and apoptosis, its action also likely regulates prostaglandin metabolism, growth factors and their receptors, apoptosis and the nuclear factor (NF)-{kappa}B cascade, extracellular matrix accretion, angiogenesis, blood coagulation, and inflammation. In particular, it increased mRNA concentrations of genes related to the vascular endothelial growth factor R2 pathway of angiogenesis and down-regulated ones associated with hypoxia. Two genes implicated in the antiluteolytic actions of IFNT (encoding cyclooxygenase-2 and the oxytocin receptor respectively) were down-regulated in response to all treatments. IFNT targets a complex range of physiological processes during the establishment of pregnancy.


    Introduction
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions
 Acknowledgements
 References
 
Interferon-{tau} (IFNT) is generally regarded as the primary signal for maternal recognition of pregnancy in the true ruminant species such as cattle and sheep and helps to maintain the structural and functional integrity of the corpus luteum (CL) at a time in a nonpregnant animal that this structure would be poised to regress (Roberts et al. 1992). By such intervention, IFNT allows progesterone production by the CL to continue and hence maintains the endometrium in a receptive state. The elongating conceptus secretes IFNT in large amounts during this period, which immediately precedes definitive attachment of trophectoderm to the uterine epithelium, i.e. ~12–21 days in sheep (Hansen et al. 1988, Farin et al. 1990) and ~14–25 days in cattle (Helmer et al. 1987, Farin et al. 1990). In its luteoprotective role, IFNT prevents or diminishes oxytocin-induced ‘pulses’ of prostaglandin F2{alpha} (PGF2{alpha}; PGF) from the uterus, most likely by suppressing the expression of oxytocin receptors (OTR) on the luminal and upper glandular epithelial cells lining the uterus (Spencer et al. 1995a). Hence, the immediate target for IFNT action is on these epithelial cells, which are exposed to very high concentrations of this protein for an extended period of time.

IFNT is a member of the Type I IFN family (Roberts et al. 1997). These cytokines have broad pleiotropic actions on their target cells (Stark et al. 1998). Type I IFNs bind to a common receptor complex consisting of two subunit polypeptides (IFNAR1 and IFNAR2; Pestka et al. 2004), both of which are highly concentrated in uterine epithelial cells relative to other cells in ovine endometrium (Rosenfeld et al. 2002). Upon binding to the receptor, Type I IFNs activate the well-studied JAK–STAT signaling pathway (Stark et al. 1998), leading to the induction of several downstream effector genes, often named interferon-stimulated genes or ISGs. In addition, Type I IFNs activate additional signal transduction pathways, including components of the MAPK pathway (Uddin et al. 2000).

Several ISGs are regulated by IFNT in the endometrium of cattle and sheep (Austin et al. 1996, Chen et al. 2006, Gray et al. 2006). In sheep, at least, increased expression of some of these genes has been noted in stromal and immune cells resident in the uterus as well as in the surface and glandular epithelium (Johnson et al. 1999a, 2000, 2001a, 2002, Choi et al. 2001). In addition, different ISGs show complex and distinct patterns of expression within the endometrium, suggesting that regulation of gene expression is complex. Most of the endometrial ISGs so far studied are ones anticipated to be up-regulated as part of the classical antiviral response. It remains unclear whether the activation of these viral-response genes has any role in preparing the uterus for pregnancy. However, at least three genes, specifically those encoding estrogen receptor-{alpha} (ESR1; Spencer et al. 1995b, Spencer & Bazer 1996), the OTR; Spencer et al. 1995a, Spencer & Bazer 1996), and cyclooxygenase-2 (PTGS2; Xiao et al. 1998, 1999, Binelli et al. 2000, Pru et al. 2001) have been implicated in the antiluteolytic action of IFNT and reported to be down-regulated following IFNT exposure. Confusingly, others have been unable to confirm this result for PTGS2, which has variously been reported to be up-regulated (Asselin et al. 1997a, 1997b, Emond et al. 2004) or unaffected by IFNT (Kim et al. 2003b). We recently showed that OTR and PTGS2 gene expression were down-regulated at day 14 of pregnancy in sheep and by treatment of ewes with IFNT and concluded that the mechanism whereby IFNT acts as an antiluteolysin still remains unclear, although a strong case can be made for the involvement of OTR (Chen et al. 2006).

The endometrium responds to IFNT secreted by the conceptus in a paracrine manner. The surface epithelium of the uterine lumen is probably the primary target for IFNT, although there is also evidence that this cytokine can reach cells in the underlying stroma (Johnson et al. 1999b, 2000, 2001b, Hicks et al. 2003) and even the myometrium (Ott et al. 1998, Johnson et al. 1999d, Hicks et al. 2003), suggesting that the epithelium does not act as a completely impermeable barrier to conceptus-derived proteins. Although cell lines from a number of sources have been used to explore the biological activities of IFNT, including ones from nonendometrial tissues, such as L929 (mouse fibroblast) cells, WISH (human amnion) cells, Daudi (human B-cell) cells, human 2fGH (human fibrosarcoma) cells, and U3A (STAT1-deficient 2fTGH) cells (Alexenko et al. 1997, 1999, Stewart et al. 2001b, Kim et al. 2003a), cells or tissue explants from endometrium clearly provide the best models for studying the physiological role of IFNT. In this paper, we have chosen to examine the effects on an immortalized luminal epithelial (LE) cell line initially isolated at day 5 of the estrous cycle (Johnson et al. 1999c). The LE line has been relatively well characterized and mimics at least some of the responses of the ovine uterus to IFNT in terms of ISG expression (Johnson et al. 1999c, Stewart et al. 2001a). These cells express estrogen and progesterone receptors (Johnson et al. 1999c) and both Type I IFN receptor subunits, IFNAR1 and IFNAR2 (Rosenfeld et al. 2002). They provide a convenient means for studying the effects of IFNT on gene expression in a specific ovine uterine cell type, although it would be naïve to assume that these cells can exactly recapitulate responses in vivo. Nevertheless, data from such cells can provide insight into the range of IFNT-regulated genes and pathways likely to be modulated by IFNT during pregnancy, and can be compared with existing in vivo information on gene expression changes that occur in response to the presence of a conceptus (Ishiwata et al. 2003, Bauersachs et al. 2006, Klein et al. 2006) or to exposure of the uterine tract to IFNT (Chen et al. 2006, Gray et al. 2006).

We have examined two kinds of response of the LE cells to IFNT, namely the effects of concentration and time of exposure, since both variables have been reported to be important. For example, IFNT increased PGE2 production in a dose-dependent manner (Asselin et al. 1997b, 1997c), and while low concentrations (equivalent to 5x10–8 M or less) inhibited the production of PGF in primary bovine endometrial cell cultures (Xiao et al. 1998, 1999, Asselin & Fortier 2000), high concentrations (~5x10–7 M) promoted such production (Asselin et al. 1997a, 1997c, Parent et al. 2003). In addition, IFNT reduced PTGS2 mRNA at 6 h or earlier but not at later time points (Xiao et al. 1999). To limit the number of analyses, we chose two exposure times: 8 and 24 h. The former was anticipated to provide data on ‘early’ and ‘intermediate’ genes, the second on genes whose regulation was late and possibly secondary to the initial response to IFNT.

Although gene transcript profiling for IFN-responsive genes in human cells was first performed in 1998 (Der et al. 1998), there have been few other comprehensive studies since then in any species. Genes that respond to ovine IFNT have been described in human fibroblast cells (Stewart et al. 2001b, Kim et al. 2003a). In addition, our laboratory has examined the expression of 70 genes considered likely to be IFNT responsive in endometrium of ewes that were either pregnant or had been administered IFNT (Chen et al. 2006), while, in the same species, Gray et al.(2006) used a microarray representing ~5000 ovine genes to study genes regulated during pregnancy and in response to a combination of IFNT and progesterone. In the present study, we have made the use of a 15 634 bovine unigene set to analyze IFNT responses in LE cells. There were four goals. The first was to determine whether the changes we observed earlier (Chen et al. 2006) in endometrium of ewes, including effects on genes influencing prostaglandin metabolism, would be reproduced in the epithelial cell line. The second and third goals were to determine whether raising the concentration and time of exposure to IFNT would influence gene expression profiles, recalling that the surface epithelium of the intact uterus must be exposed to unusually high amounts of the cytokine for several days. The final goal was to generate a comprehensive list of genes regulated by IFNT in ovine endometrium. Presumably, the responses of some but not necessarily all of these genes provide the basis for how IFNT mediates changes in the mother that provide for a successful establishment of pregnancy.


    Materials and Methods
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions
 Acknowledgements
 References
 
Cell culture
The LE cells were maintained in Dulbecco’s modified Eagle’s medium with F12 salts (DMEM/F12; Sigma–Aldrich Corp.), pH 7.4, containing 10% (v:v) fetal bovine serum (Gibco-BRL) and penicillin/streptomycin/amphotericin B solution (100 IU/ml, 0.1 mg/ml, 0.25 µg/ml; Gibco-BRL) as described previously (Johnson et al. 1999c).

Recombinant ovIFNT
Recombinant protein was produced as described elsewhere (Ealy et al. 1998a, 1998b). Antiviral assays to assess potency of the preparation were performed on MDBK cell (ATCC#CCL22, American Type Culture Collection; Ealy et al. 1998b). The antiviral activity of this preparation of IFNT was 2x108 international units/mg.

Treatment of cells with IFNT
For the microarray experiments, three stocks of LE cells, frozen on different occasions, were used. Each experiment, which was performed thrice, employed three replicate groups of cells representing the three different cell stocks. Both experiments had an identical experimental design. In each, cells were allowed to reach ~80% confluence in 75 cm tissue culture flasks. The medium was then replaced with 5.0 ml fresh DMEM/F12 containing either recombinant ovine IFNT at two different concentrations (10 ng/ml; ~5x10–10 M and 10 µg/ml; ~5x10–7 M) or with the IFNT vehicle (1xPBS) as a control (Parent et al. 2003). The lowest concentration of IFNT was close to the KD for the ovine receptor, while the highest was supersaturating, but probably closer to IFNT concentrations in the uterine lumen at the time of maximal production by the conceptus. Cells were exposed to IFNT for either 8 or 24 h before isolation of RNA. There were, therefore, 36 samples (9 per treatment) of RNA prepared (see below). Cultures of pooled, untreated (control) cells that were harvested at the 8-h time point provided the reference RNA for both the microarray and real-time PCR analyses (Table 1Go).


View this table:
[in this window]
[in a new window]

 
Table 1 Real-time PCR analysis of mRNA expression in ovine luminal epithelial cells treated with interferon-{tau} (IFNT).
 
Microarray analysis
The DNA array slides were made by spotting cDNA PCR products in 3xSSC solutions onto the poly-L-lysine coated glass slides. Each array contained 17 692 cDNA probes from an available cDNA library of 16 156 probes (Smith et al. 2001) and a bovine ovary EST library with 1536 probes (Takasuga et al. 2001). In addition, the new library represented 15 634 unique sequences.

RNA was isolated from individual plates of cells using RNA STAT-60 (Tel-Test, Friendswood, TX, USA) according to manufacturer’s recommendations. The quantity of RNA for each extraction was determined by measuring absorbance at 260 nm; purity was calculated from the ratios of absorbance at 260 and 280 nm; and quality was evaluated by agarose gel electrophoresis. Extracted RNA was stored at –80 °C. For the microarray analysis, the triplicate samples for each of the four experimental treatments were pooled (RNA ratio 1:1:1) to provide four pools of RNA. Since the experiment was performed thrice, there were 12 RNA sample pools in total, which were each analyzed through a comparison with a single sample of control RNA from untreated cells. Thus, each RNA sample represents pooled material from three different isolates of LE cells, thereby minimizing any confounding effect of variation among different batches of cells.

The microarray analysis utilized also-called Reference Design approach in which we compared RNA from the various treatment groups with the control RNA, following the microarray protocol from TIGR (The Institute for Genomic Research; Hegde et al. 2000). RNA (15 µg) was reverse transcribed using SuperScript III reverse transcriptase (Invitrogen) to provide incorporation of the reactive amine derivative of dUTP, 5-(3-aminoallyl)-2'-deoxyuridine 5'-triphosphate (Sigma) into the cDNA. After removing free nucleotides, the aminoallyl-labeled samples were coupled to dye. The cDNAs from treated and control were labeled separately with either cyanine 3 (Cy3) or Cy5 dyes by standard procedures (Amersham Biosciences). The resulting labeled probes were purified, mixed together, and then redissolved in 25 µl hybridization solution (50% formamide, 5xSSC, 0.1% SDS, 10 µg ovine Cot-1 DNA, and 10 µg poly-A in water). The samples were denatured and then applied to the microarrays. After incubation at 42 °C for 16 h, the arrays were sequentially washed with 1xSSC/0.2% SDS, 0.1xSSC/0.2% SDS, and 0.1xSSC solutions. The arrays were then dried and scanned using a Genepix 4000B scanner (Axon Instruments, Union City, CA, USA) and analyzed by the Genepix Pro software. As a control, the dye labeling was reversed in one of the hybridizations (dye-swap control).

Data analysis
The intensity value of each spot for Cy3 and Cy5 fluorescence (media foreground intensity minus media background intensity) was analyzed by R/MAANOVA (MicroArray ANOVA; Wu et al. 2002), which runs in the statistical software program ‘R’ (Beckers et al. 1988) from www.r-project.org. Briefly, data were first normalized by means of the regional Lowess algorithm to account for systematic intensity dependence and spatial variation in microarray log ratios, which is to correct the bias caused largely by differences of dye incorporation (Berger et al. 2004). Then, a two-stage ANOVA model was applied to the data. In the first stage, the model was


Formula

where Y is the log intensity reading for a particular gene on a certain array, e.g. array i, dye j of RNA variety k; µ is the overall mean; Ai is the effect of array i (i=1,...,12), Dj is the effect of dye j (j=1,2); and AD is the effect of array and dye interactions. This is the cross-gene stage of the calculation, which removes the overall effects of array, dye effects, and the effects of array and dye interactions.

The residual of the first stage is used as the input for the second stage, which is the gene-specific model.


Formula

Here, G is the average effect of the gene, AG is the array-by-gene variation, and DG is the dye-by-gene variation. VG is the effects of the treatment on the expression of the gene of interest. The error {varepsilon} is the residual. A permutation F-test (Cui & Churchill 2003) was applied by MAANOVA, and the adjusted P values were calculated for every gene to indicate if expression was changed by treatment. Finally, a two-dimensional hierarchical clustering analysis was performed using GeneCluster (www.broad.mit.edu/cancer/software/genecluster2/gc2.html) and TreeView software (Eisen et al. 1998; http://rana.lbl.gov/EisenSoftware.htm; Fig. 1Go). We combined use of the DAVID Bioinformatics Resources Data Base (2006; http://david.abcc.ncifcrf.gov/; Dennis et al. 2003, Brunner et al. 2004) with manual examination of regulated genes to provide functional annotation clustering.


Figure 1
View larger version (25K):
[in this window]
[in a new window]

 
Figure 1 Two-way hierarchical cluster analysis of gene expression in control- and IFNT-treated LE cells. Here we show only the 585 genes that are differentially regulated (P<0.001). The gene set depicts two major groups (up- or down-regulated), which can be further divided into five subgroups (subgroups 1–5). The data clustered according to treatment time (8 and 24 h) and dose of IFNT used (L, low; H, high).

 
Real-time PCR analysis
The relative expression levels of five selected ovine genes were determined by real-time PCR in part to confirm the reproducibility of the microarray and also to provide additional quantitative data for genes of particular interest. Primers specific for the genes for ISG-15/17, originally known as ubiquitin cross-reactive protein (ISG15), PTGS2, OTR, insulin-like growth factor-2 (IGF2), and hypoxia-inducible factor 1 subunit {alpha} (HIF1A) were designed with Primer Express software (Applied Biosystems, Foster City, CA, USA). The primer pairs are listed in Supplementary Table 1 which can be viewed online at www.reproduction-online.org/supplemental/.

Real-time PCR was performed and analyzed by using an ABI prism 7500 sequence detection system (Applied Biosystems). The protocol for real-time RT-PCR analysis has been described earlier (Chen et al. 2006). Relative expression in comparison with ribosomal protein L19 (RPL19) was calculated by the comparative CT method (threshold cycle). The relative quantitation of candidate gene expression in each treatment groups was determined by using the comparative CT method (2{Delta}{Delta}CT) as described in user bulletin 2 ABI Prism 7500 Sequence Detection System.

Each of the five transcripts was amplified in triplicate in a single PCR run. The entire experiment was then repeated with RNA from two additional experiments from the same treatment to provide a total of three replicates. Real-time PCR quantification of gene expression level in each treatment was the mean of three real-time RT-PCR experiments. Differences in relative mRNA expression between experimental groups were assessed by one-way ANOVA, followed by pairwise comparison and by using a least significant difference test (Brunner et al. 2004). All experimental data are shown as the mean±S.E.M. Values were considered significantly different at P<0.05.


    Results
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions
 Acknowledgements
 References
 
Identification of IFN-{tau}responsive genes
Using the linear mixed model described in the Material and Methods, we identified 585 genes that had altered expression level at an adjusted P value <0.001 following IFNT treatments (Supplementary Table 2 which can be viewed online at www.reproduction-online.org/supplemental/). At least another ~700 genes (a total of 1274 or 8% of all genes) were IFNT responsive if the threshold of significance was raised to an adjusted P value <0.01 (data available upon request). Among the 585 genes in Supplementary Table 2, 567 have annotations that are based on sequence similarity with known genes. We applied a two-way (genes and treatments) hierarchical clustering method (Eisen et al. 1998) on the 585 genes regulated by IFNT relative to the expression of the same genes in untreated cells that had been isolated 8 h after the addition of IFNT to the treatment groups (Fig. 1Go). The treatments clustered into two main branches, which corresponded to the 8 and 24 h treatments. The genes regulated by the low (10 ng/ml) and high (10 µg/ml) treatment groups formed two groups within these main branches of the tree, i.e. according to the duration of the treatments. Two major groups of genes can be easily recognized, one that consists of genes that were up-regulated in one or more of the treatments and a second that consists of down-regulated genes. A few genes were up-regulated at one time point and down-regulated at the other. In general, a concentration of 10 ng/ml IFNT was about as effective as 10 µg/ml. Based on visual inspection of the expression patterns, the genes illustrated in Fig. 1Go could be divided arbitrarily into five subgroups, although there were distinct expression patterns visible within each grouping.

Members of subgroup 1 (Supplementary Table 3 which can be viewed online at www.reproduction-online.org/supplemental/), which contains 68 genes, were generally expressed more strongly at 8 h than at 24 h (Fig. 1Go), although there were several exceptions. The genes in subgroup 1 encode products with a broad range of biochemical roles. Among them are the genes for transcription factors associated with downstream action of Type 1 IFN, including RNF31 (interferon-stimulated transcription factor 3-{gamma} 48 kDa) and interferon regulatory factor 7 (IRF7), three guanylate-binding proteins (GBP4, GBP5, and GEM), and a few genes encoding inflammatory cytokines, e.g. TNFSF10 and IL6ST. The extent of up-regulation of the subgroup 1 genes by IFNT at both the low and high concentrations was relatively modest when compared with many of the genes in subgroup 2 discussed below.

Subgroup 2 comprises 252 genes (Supplementary Table 4 which can be viewed online at www.reproduction-online.org/supplemental/) whose expression was, in general, sustained over 24 h (Fig. 1Go). Subgroup 2 contained the majority of the genes that were up-regulated fivefold or more by IFNT. Many of these highly up-regulated genes have long been recognized as responsive to Type I IFN, and many are included in the original list of genes identified by microarrays in human cells (Der et al. 1998, Kim et al. 2003a). They include Mx1, ISG15, interferon-induced protein 6-16 (G1P3), interferon-induced protein with tetratricopeptide repeats (IFIT2), 2'-5'-oligoadenylate synthetase (OAS1), IFN-induced protein-35 (IFI35), multiple (MHC) genes from classes I and II, and ß-2-microglobulin (PDLIM3). Many of the proteins encoded by these subgroup 2 genes contribute to the antiviral response of cells targeted by Type 1 IFN and are ‘classical’ ISGs (Stark et al. 1998).

Subgroup 3 consists of 31 genes that were modestly, but significantly up-regulated by IFNT at 8 h but provided inconsistent temporal patterns of expression thereafter. The genes did not appear to fall into any particular functional class (Supplementary Table 5 which can be viewed online at www.reproduction-online.org/supplemental/).

Subgroup 4 is comprised 87 genes that were down-regulated by all IFNT treatments (Fig. 1Go; Supplementary Table 6 which can be viewed online at www.reproduction-online.org/supplemental/). IFN responses were, in general, sustained over the 24-h exposure to IFNT. This cohort contains several genes of potential interest in endometrial responses to pregnancy, including ones encoding metalloproteinases, the OTR, prostaglandin synthase, and IGF-2, which are discussed later.

The final subgroup contains 123 genes that were generally more strongly down-regulated at 24 h than at 8 h (Supplementary Table 7 which can be viewed online at www.reproduction-online.org/supplemental/). Again the responses to the high IFNT concentration were only marginally different from those at the low concentration.

Validation of gene expression by real-time PCR
To validate our results from the multiple microarray analyses, we performed real-time PCR experiments on five selected genes. ISG15 was chosen as a positive control for validating the experimental procedure because it has been previously shown to be up-regulated by IFNT in ovine LE cells (Johnson et al. 1999c). We selected PTGS2, IGF2, and hypoxia-inducible factor-1{alpha} subunit (HIF1A), because we had previously demonstrated these genes to be down-regulated in endometrium from ewes that were either pregnant or, alternatively, treated with IFNT (Chen et al. 2006). OTR was selected because of its importance in controlling endometrial responses to oxytocin. In general, the results from real-time PCR experiments showed the same trends as the microarray experiments (Table 2Go), although the magnitude of fold-changes was sometimes quite different. One possible explanation is that the real-time PCR experiments were performed on un-pooled samples of RNA obtained from individual plates of cells, while the microarrays employed pooled RNA from replicate samples within an experiment.


View this table:
[in this window]
[in a new window]

 
Table 2 Regulated genes in prostaglandin metabolic pathway.
 

    Discussion
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions
 Acknowledgements
 References
 
These experiments, employing microarrays based on bovine EST sequences, allowed us to collect data on over 15 000 genes, the majority of which were at least partially annotated by comparison with human data. Our study identified a large number of genes that have altered expression patterns when LE cells were treated with IFNT. We did not employ a cutoff based on a particular value of fold-change, e.g. twofold, since there is no reason to assume that values of less than this arbitrary value are of minor significance. Also in absence of a detailed time-course, it is not possible to infer whether a low fold-change at 8 h, for example, does not reflect a more major change at an earlier or later time point. We chose to concentrate on the ‘top’ 585 genes (of which 567 were at least partially annotated) where the changes were highly significant (adjusted P<0.001). Of these genes, 356 were up-regulated, while the reminder showed a reduction in expression. Although some of these genes were late responders, most were fully up- or down-regulated within 8 h. In addition, the lowest concentration of IFNT, which was deliberately chosen to be close to the KD for the ovine receptor (Knickerbocker & Niswender 1989, Li & Roberts 1994), was generally about as effective as the highest concentration (~100xKD) in causing changes in gene expression, although there were some exceptions (Supplementary Table 2). This result is consistent with the assumption that Type I IFNs are effective at low receptor occupancy (Alexenko et al. 1997).

Other experiments employing microarrays to study the effects of IFNA and IFNB on human cells also indicated that large numbers of genes are responsive to Type I IFN (Der et al. 1998). Kim et al. (2003a) identified 101 IFNT-responsive genes in human U3A (STAT1-deficient 2fTGH) cells. Most recently, Gray et al.(2006) identified 180 genes that were regulated by IFNT independently of progesterone in day 14 ovine endometrium. Yan et al.(2004) in a proteomics analysis performed on human liver cells reported that as many as 1364 polypeptides were regulated by IFNA. In addition, these data illustrate the many downstream effects of Type I IFN on their targets. As in all previous studies, the majority of the IFN-responsive genes shown in Supplementary Table 2 were difficult to categorize, possibly reflecting the broad, pleiotropic effects of Type I IFN on their target cells. Many of these same IFNT-responsive genes, particularly ones implicated in an antiviral response, have been observed to be regulated during pregnancy in cattle and sheep, although it is clear that early maternal responses to the presence of a conceptus are likely to be more than simply a result of exposure of LE cells to IFNT. However, a number of interesting patterns of gene expression change that might relate to the effects of IFNT during pregnancy could be discerned.

Although our results agree closely with an in vivo IFNT study performed on a limited number of selected genes (~70 vs 15 000 in the present study; Chen et al. 2006), they are difficult to compare with another, recent, in vivo study that utilized microarray technology (Gray et al. 2006). The latter utilized ~5000 ovine ESTs and a relatively complex experimental design, in which IFNT effects were distinguished from those of progesterone by exposing endometrium of some ewes to an antiprogestin. Gray et al.(2006) described 180 genes regulated by IFNT, but many of these were not represented in our gene list (Supplementary Table 2) and vice versa. Some of the more interesting differences are remarked upon later in the discussion that follows.

IFNT regulates the prostaglandin synthesis pathway
The production and release of uterine PGF is crucial to controlling the regression of the CL during early pregnancy in cattle and sheep (McCracken et al. 1999). Conversely, PGE2 is luteotrophic (Magness et al. 1981, Thibodeaux et al. 1992). As shown in Fig. 2Go, the production of these prostaglandins is controlled by several enzymes, including PTGS2, PGE synthase (PGES), PGF synthase (AKR1C1), and PG 15-dehydrogenase (PTGES), and their metabolism by several other gene products, including 15-hydroxy prostaglandin dehydrogenase (PGDH) and the prostaglandin transporter (SLC21A2). Of these, PTGS2, which is responsible for the conversion of arachidonic acid into PGH2, is considered to catalyze the rate-limiting step in the biosynthesis of prostaglandins in endometrial tissue (Simmons et al. 2004). In addition, an increase in the PGE2:PGF ratio in favor of PGE2 occurs during pregnancy in sheep and cattle (Silvia et al. 1984, Payne & Lamming 1994). The regulation of PTGS2 expression by IFNT has been exceedingly controversial, with some reporting an increase and others a decrease in expression of either the enzyme or the mRNA (see Introduction). In the present study, it is clear that PTGS2 and PGF synthase (AKR1C1) gene expression are down-regulated in LE cells in response to IFNT, while transcripts for PGE2 synthase (PGES) increase in concentration. We obtained essentially identical results when comparing endometrium from day 14 pregnant and nonpregnant ewes and after ewes at day 14 of their estrous cycles had been treated with IFNT (Chen et al. 2006). These data contrast with other data (Gray et al. 2006) in which PTGS2 gene expression was up-regulated in response to IFNT treatment. We have no explanation for these differences in outcome. Changes in expression of other genes with a role in prostaglandin metabolism and action were also evident. Transcripts for PGDH, the PGE receptor (PTGER4), and the prostaglandin transporter (SLC21A2), for example, increased. Overall, our data are generally consistent with the view that IFNT treatment favors PGE production over PGF (Fig. 2Go).


Figure 2
View larger version (17K):
[in this window]
[in a new window]

 
Figure 2 Prostaglandin metabolism in mammalian cells. The diagram is based on one downloaded from http://www.genome.jp/kegg/. Up or down arrows indicate whether the gene is positively or negatively regulated by IFNT in the LE cells.

 
None of the genes in Fig. 2Go thought to be involved in PGF metabolism appeared to be differentially responsive as IFNT concentrations were raised from 10 ng/ml to 10 µg/ml. Thus, we were unable to provide any supporting evidence for the view (Asselin & Fortier 2000) that prostaglandin metabolism, and specifically the production of PGF, might differ according to the concentrations of IFNT employed, with low concentrations favoring lowered PGF production and high concentrations favoring an increase. On the other hand, transcript profiling can be misleading since it provides only information relating to mRNA concentrations and nothing about how well these messages are translated or the concentration of relevant proteins in the cells.

IFNT up-regulate components of signal transduction pathways
Type I IFN signaling has been studied extensively. Since all Type I IFN share the same receptor complex, IFNT is likely to activate some or all of the signal pathways employed by IFNA and IFNB, including the JAK–STAT, and MAPK pathways (de Veer et al. 2001). Our data indicate that IFNT, as expected, up-regulated genes encoding key components of these two signaling pathways, STAT1 and RNF31 in the JAK–STAT pathway and MAP3K14 (Supplementary Tables 2 and 3) from the MAPK pathway. Components of the two pathways have also been identified as IFNA-responsive genes in human cell lines (Der et al. 1998, de Veer et al. 2001). At least 16 transcription factors, presumably downstream of IFNT, increased in expression, and at least three were down-regulated (Supplementary Table 8 which can be viewed online at www.reproduction-online.org/supplemental/). With the exception of STAT1, RNF31, IRF7, and IRF9, however, the relationship of these transcription factors to the action of Type 1 IFNs has not been explored.

IFNT promotes increased expression of inducers of apoptosis and intracellular proteolysis
A DAVID annotation analysis indicated that 24 genes (P<0.00 006) in Supplementary Table 1 were associated with control of apoptosis, with the majority up-regulated. We have listed some of these genes in Table 3Go. They include BAD (Bcl-2 antagonist of cell death), (Yang et al. 1995), BID (BH3-interacting domain death agonist), which induces caspases and apoptosis and counters the protective effect of Bcl-2 (Wang et al. 1996), and BNIP3L (Bcl2/adenovirus E1B 19 kDa protein-interacting protein 3-like), a well-characterized inducer of apoptosis (Opferman & Korsmeyer 2003). The genes for caspase 8 (CASP8) and several proteosome components were also up-regulated. The ability of Type I IFNs to cause apoptosis (Chawla-Sarkar et al. 2003) and induce genes involved in apoptosis (Der et al. 1998, Leaman et al. 2003) is well established in human target cells and cultured epithelial cells from bovine endometrium (Wang et al. 2003) although we have not examined whether cell death resulted in the LE cells used here. Similarly, a link between IFN action and induction of proteosomal activity has been described by others (Li & Hassel 2001, Liu et al. 2003). These changes in gene expression may reflect a potential ability of conceptus IFNT to cause remodeling of the uterine epithelium and underlying stromal tissue prior to firm attachment of the trophoblast to the uterine wall.


View this table:
[in this window]
[in a new window]

 
Table 3 Regulated genes related to apoptosis.
 
Interestingly, at least seven genes (P<0.005) (FADD, TNFSF10, SCOTIN, TRIM38, BST2, CXXC5, EEFID) associated with control of the NF-{kappa}B cascade were up-regulated by IFNT under all four experimental treatments (Supplementary Table 1). Somewhat similar changes on genes playing a role in protein degradation and the NF-{kappa}B cascade have been noted when comparing endometrium of pregnant and nonpregnant cattle (Bauersachs et al. 2006).

IFNT regulates the expression of angiogenesis factors
The conceptus as it begins to form a close association with the uterine wall recruits a localized blood supply in the underlying maternal endometrium, which is accompanied by local edema in the regions where the trophoblast makes contact (Spencer et al. 2004). The formation of new blood vessels involves the coordinated actions of proteolytic enzymes, protease inhibitors, extracellular matrix components, and cell adhesion molecules. IFNA and IFNB are generally considered to be angiostatic, in part by down-regulating the production of FGF2 (Fidler 2000). On the other hand, IFNA has also been reported to increase expression of several angiogenic factors (de Veer et al. 2001). The response to IFNT in LE cells was also mixed (Table 4Go). Some genes were down-regulated, including FGF2, angiopoietin-1 (ANGPT1), VE-cadherin (CDH5), and HIF1A and EPAS1 (hypoxia-inducible factor-2{alpha}), which direct responses to low O2 conditions, while other genes implicated in blood vessel formation were up-regulated, e.g. angiopoietin-2 and -4 (ANG2 and ANG4), and vascular endothelial growth factors B and C (VEGFB and VEGFC). These ostensibly conflicting data are not easy to interpret and suggest complex effects of IFNT on angiogenesis in the endometrium. One possibility is that IFNT activates certain angiogenic pathways, such as the one that involves VEGF–VEGFR2 system, while at the same time it down-regulates pathways normally initiated by hypoxia. Such a complementary system for promoting angiogenesis has recently been reported in tumors (Casanovas et al. 2005). In the pregnant uterus, we propose that LE and possibly the upper, contiguous glandular epithelial cells, which are responsible for transducing signals received from the trophoblast to the underlying tissue, are stimulated to produce increased amounts of angiogenic factors such as VEGFs and therefore may direct the growth of new vessels from the existing capillary network bed.


View this table:
[in this window]
[in a new window]

 
Table 4 Regulated genes related to angiogenesis.
 
In addition to its effects on genes associated with angiogenesis, IFNT up-regulated seven genes linked to blood coagulation (CD93, F2RL2, PLAU, TAUI, ANXA5, PLSCR1). The implications of such a link are unclear.

IFNT effects on the IGF system
IGF2 is believed to play a major role in controlling nutrient balance between the mother and the conceptus (Reik et al. 2003). Its gene is paternally imprinted, i.e. expressed in a mono-allelic manner from the paternal chromosome. Interestingly, IGF2 was down-regulated in LE cells in response to IFNT (Table 5Go), a result consistent with that reported by us for in vivo effects (Chen et al. 2006). By contrast, IGFBP2, IGFBP3, and IGFBP4 became up-regulated in the LE cells. One possibility is that IFNT from the conceptus functions to repress or deregulate the local maternal IGF2 activity, possibly to gain control of maternal nutrient supply and angiogenesis.


View this table:
[in this window]
[in a new window]

 
Table 5 Regulated genes in the insulin-like growth factor system.
 
IFNT regulates expression of MMPs, TIMPs, and extracellular matrix components
Proteins of the matrix metalloproteinase (MMP) family are involved in the degradation of extracellular matrix during many normal as well as pathological processes (Fata et al. 2000) and are also able to promote cell growth indirectly through interactions with cytokines and growth factors (Curry & Osteen 2003). Tissue inhibitors of metalloproteinases (TIMPs) counteract and therefore modulate the activities of the MMPs and may also act as autocrine/paracrine factors in cellular proliferation and differentiation. IFNT down-regulated genes for several MMPs, up-regulated TIMP1, TIMP2, TIMP3, and TIMP4, and increased transcript concentrations for a number of genes encoding extracellular matrix components (Table 6Go). Others have demonstrated that IFNT inhibits production of MMPs in primary cultures of bovine endometrial cells (Salamonsen et al. 1994, Hashizume et al. 2003). In addition, these observations suggest that IFNT promotes extracellular matrix accretion rather than breakdown in the ovine uterus during the preimplantation period. By contrast, matrix degradation appears to play a role in the decidualization process during implantation in mouse (Das et al. 1997). Decidualization does not occur in species, such as the sheep, with a minimally invasive trophoblast. Conceivably, the action of IFNT is to prevent decidualization rather than to induce it.


View this table:
[in this window]
[in a new window]

 
Table 6 Genes related to extracellular matrix.
 
Other genes implicated in the luteoprotective action of IFNT
As noted in our earlier in vivo work (Chen et al. 2006), IFNT can rapidly down-regulate OTR, which is consistent with the work of others over the last decade (Mirando et al. 1993, Lamming et al. 1995, Spencer et al. 1995a, 1995b). On the other hand, we failed to obtain a sufficiently strong hybridization signal for the mRNA of either ESR1 or ESR2 to determine whether these genes were significantly regulated in our system by IFNT. One well-accepted model for the luteoprotective action of IFNT is that the cytokine causes reduction of OTR number on the uterine epithelium indirectly as an outcome of preventing the normal up-regulation of ESR1 (Spencer et al. 1995a, 1996), whereas in the nonpregnant cycle the rise in OTR allows the hormone, oxytocin to mediate PGF release (Robinson et al. 1999). Clearly, our results are consistent with one part of this model, namely the down-regulation of OTR in response to IFNT, but could not address the second.


    Conclusions
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions
 Acknowledgements
 References
 
The present study provides a comprehensive, descriptive profiling of gene expression in ovine LE cells in responsive to IFNT. Out of more than 15 000 genes examined, we identified over 1200 genes that were regulated (P<0.01). Of these, we concentrated on the genes whose transcripts changed most significantly (P<0.001). Many of the genes identified in the present study as IFN responsive have been noted earlier by others in human cells in response to IFNA and IFNB (Der et al. 1998, Kim et al. 2003a, Yan et al. 2004). However, our list is more comprehensive and has also expanded our understanding of the likely biological activities of IFNT on a physiologically relevant cell type. Our results suggest that IFNT, like other Type I IFN, up-regulates many genes that undoubtedly play a role in the antiviral activities of this class of cytokine. However, it is also clear that IFNT action on the uterus likely involves effects on prostaglandin metabolism, apoptosis, extracellular matrix accretion, angiogenesis and blood coagulation, inflammation, and growth factor production and activity. The recognition that such processes are likely to be modulated as the conceptus begins to attach to the uterine wall provides new insights into how IFNT might act during early pregnancy and the likely complexity of its role in maternal recognition of pregnancy in ruminants.


    Acknowledgements
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions
 Acknowledgements
 References
 
We thank Dr Tom Spencer, Texas A & M University, for supplying the LE cells and Norma McCormack for editing the manuscript and preparing the Figures and Tables. The authors declare that there is no conflict of interest that would prejudice the impartiality of this scientific work. Research was supported by a grant from NIH Grant HD 21896 to RMR.


    Footnotes
 
Received 21 December 2006
First decision 2 February 2007
Accepted 2 April 2007


    References
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions
 Acknowledgements
 References
 

Alexenko AP, Leaman DW, Li J & Roberts RM 1997 The antiproliferative and antiviral activities of IFN-{tau} variants in human cells. Journal of Interferon and Cytokine Research 17 769–779.[ISI][Medline]

Alexenko AP, Ealy AD & Roberts RM 1999 The cross-species antiviral activities of different IFN-{tau} subtypes on bovine, murine, and human cells: contradictory evidence for therapeutic potential. Journal of Interferon and Cytokine Research 19 1335–1341.[CrossRef][ISI][Medline]

Asselin E & Fortier MA 2000 Detection and regulation of the messenger for a putative bovine endometrial 9-keto-prostaglandin E(2) reductase: effect of oxytocin and interferon-{tau}. Biology of Reproduction 62 125–131.[Abstract/Free Full Text]

Asselin E, Drolet P & Fortier MA 1997a Cellular mechanisms involved during oxytocin-induced prostaglandin F2{alpha} production in endometrial epithelial cells in vitro: role of cyclooxygenase-2. Endocrinology 138 4798–4805.[Abstract/Free Full Text]

Asselin E, Lacroix D & Fortier MA 1997b IFN-{tau} increases PGE2 production and COX-2 gene expression in the bovine endometrium in vitro. Molecular and Cellular Endocrinology 132 117–126.[CrossRef][ISI][Medline]

Asselin E, Bazer FW & Fortier MA 1997c Recombinant ovine and bovine interferons {tau} regulate prostaglandin production and oxytocin response in cultured bovine endometrial cells. Biology of Reproduction 56 402–408.[Abstract]

Austin KJ, Ward SK, Teixeira MG, Dean VC, Moore DW & Hansen TR 1996 Ubiquitin cross-reactive protein is released by the bovine uterus in response to interferon during early pregnancy. Biology of Reproduction 54 600–606.[Abstract]

Bauersachs S, Ulbrich SE, Gross K, Schmidt SE, Meyer HH, Wenigerkind H, Vermehren M, Sinowatz F, Blum H & Wolf E 2006 Embryo-induced transcriptome changes in bovine endometrium reveal species-specific and common molecular markers of uterine receptivity. Reproduction 132 319–331.[Abstract/Free Full Text]

Beckers RA, Chambers JM & Wilks AR 1988 The New S Language, New York: Chapman & Hall.

Berger J, Hautaniemi S, Jarvinen A-K, Edgren H, Mitra S & Astola J 2004 Optimized LOWESS normalization parameter selection for DNA microarray data. BMC Bioinformatics 5 194.[CrossRef][Medline]

Binelli M, Guzeloglu A, Badinga L, Arnold DR, Sirois J, Hansen TR & Thatcher WW 2000 Interferon-{tau} modulates phorbol ester-induced production of prostaglandin and expression of cyclooxygenase-2 and phospholipase-A(2) from bovine endometrial cells. Biology of Reproduction 63 417–424.[Abstract/Free Full Text]

Brunner A, Yakovlev I & Strauss S 2004 Validating internal controls for quantitative plant gene expression studies. BMC Plant Biology 4 14.[CrossRef][Medline]

Casanovas O, Hicklin DJ, Bergers G & Hanahan D 2005 Drug resistance by evasion of antiangiogenic targeting of VEGF signaling in late-stage pancreatic islet tumors. Cancer Cell 8 299–309.[CrossRef][ISI][Medline]

Chawla-Sarkar M, Lindner DJ, Liu YF, Williams BR, Sen GC, Silverman RH & Borden EC 2003 Apoptosis and interferons: role of interferon-stimulated genes as mediators of apoptosis. Apoptosis 8 237–249.[CrossRef][ISI][Medline]

Chen Y, Green JA, Antoniou E, Ealy AD, Mathialagan N, Walker AM, Avalle MP, Rosenfeld CS, Hearne LB & Roberts RM 2006 Effect of interferon-{tau} administration on endometrium of nonpregnant ewes: a comparison with pregnant ewes. Endocrinology 147 2127–2137.[Abstract/Free Full Text]

Choi Y, Johnson GA, Burghardt RC, Berghman LR, Joyce MM, Taylor KM, Stewart MD, Bazer FW & Spencer TE 2001 Interferon regulatory factor-two restricts expression of interferon-stimulated genes to the endometrial stroma and glandular epithelium of the ovine uterus. Biology of Reproduction 65 1038–1049.[Abstract/Free Full Text]

Cui X & Churchill GA 2003 Statistical tests for differential expression in cDNA microarray experiments. Genome Biology 4 210.[CrossRef][Medline]

Curry TE Jr & Osteen KG 2003 The matrix metalloproteinase system: changes, regulation, and impact throughout the ovarian and uterine reproductive cycle. Endocrine Reviews 24 428–465.[Abstract/Free Full Text]

Das SK, Yano S, Wang J, Edwards DR, Nagase H & Dey SK 1997 Expression of matrix metalloproteinases and tissue inhibitors of metalloproteinases in the mouse uterus during the peri-implantation period. Developmental Genetics 21 44–54.[CrossRef][ISI][Medline]

Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC & Lempicki RA 2003 DAVID: database for annotation, visualization, and integrated discovery. Genome Biology 4 P3.[CrossRef][Medline]

Der SD, Zhou A, Williams BR & Silverman RH 1998 Identification of genes differentially regulated by interferon {alpha}, ß, or {gamma} using oligonucleotide arrays. PNAS 95 15623–15628.[Abstract/Free Full Text]

Ealy AD, Alexenko AP, Keisler DH & Roberts RM 1998a Loss of the signature six carboxyl amino acid tail from ovine interferon-{tau} does not affect biological activity. Biology of Reproduction 58 1463–1468.[Abstract/Free Full Text]

Ealy AD, Green JA, Alexenko AP, Keisler DH & Roberts RM 1998b Different ovine interferon-{tau} genes are not expressed identically and their protein products display different activities. Biology of Reproduction 58 566–573.[Abstract/Free Full Text]

Eisen MB, Spellman PT, Brown PO & Botstein D 1998 Cluster analysis and display of genome-wide expression patterns. PNAS 95 14863–14868.[Abstract/Free Full Text]

Emond V, MacLaren LA, Kimmins S, Arosh JA, Fortier MA & Lambert RD 2004 Expression of cyclooxygenase-2 and granulocyte-macrophage colony-stimulating factor in the endometrial epithelium of the cow is up-regulated during early pregnancy and in response to intrauterine infusions of interferon-{tau}. Biology of Reproduction 70 54–64.[Abstract/Free Full Text]

Farin CE, Imakawa K, Hansen TR, McDonnell JJ, Murphy CN, Farin PW & Roberts RM 1990 Expression of trophoblastic interferon genes in sheep and cattle. Biology of Reproduction 43 210–218.[Abstract]

Fata JE, Ho AT, Leco KJ, Moorehead RA & Khokha R 2000 Cellular turnover and extracellular matrix remodeling in female reproductive tissues: functions of metalloproteinases and their inhibitors. Cellular and Molecular Life Sciences 57 77–95.[CrossRef][ISI][Medline]

Fidler IJ 2000 Angiogenesis and cancer metastasis. Cancer Journal 6 S134–S141.[ISI][Medline]

Gray CA, Abbey CA, Beremand PD, Choi Y, Farmer JL, Adelson DL, Thomas TL, Bazer FW & Spencer TE 2006 Identification of endometrial genes regulated by early pregnancy, progesterone, and interferon {tau} in the ovine uterus. Biology of Reproduction 74 383–394.[Abstract/Free Full Text]

Hansen TR, Imakawa K, Polites HG, Marotti KR, Anthony RV & Roberts RM 1988 Interferon RNA of embryonic origin is expressed transiently during early pregnancy in the ewe. Journal of Biological Chemistry 263 12801–12804.[Abstract/Free Full Text]

Hashizume K, Takahashi T, Shimizu M, Todoroki J, Shimada A, Hirata M, Sato T & Ito A 2003 Matrix-metalloproteinases-2 and -9 production in bovine endometrial cell culture. Journal of Reproduction and Development 49 45–53.[CrossRef][ISI]

Hegde P, Qi R, Abernathy K, Gay C, Dharap S, Gaspard R, Hughes JE, Snesrud E, Lee N & Quackenbush J 2000 A concise guide to cDNA microarray analysis. Biotechniques 29 548–562.[ISI][Medline]

Helmer SD, Hansen PJ, Anthony RV, Thatcher WW, Bazer FW & Roberts RM 1987 Identification of bovine trophoblast protein-1, a secretory protein immunologically related to ovine trophoblast protein-1. Journal of Reproduction and Fertility 79 83–91.[Abstract]

Hicks BA, Etter SJ, Carnahan KG, Joyce MM, Assiri AA, Carling SJ, Kodali K, Johnson GA, Hansen TR, Mirando MA et al. 2003 Expression of the uterine Mx protein in cyclic and pregnant cows, gilts, and mares. Journal of Animal Science 81 1552–1561.[Abstract/Free Full Text]

Ishiwata H, Katsuma S, Kizaki K, Patel OV, Nakano H, Takahashi T, Imai K, Hirasawa A, Shiojima S, Ikawa H et al. 2003 Characterization of gene expression profiles in early bovine pregnancy using a custom cDNA microarray. Molecular Reproduction and Development 65 9–18.[CrossRef][ISI][Medline]

Johnson G, Austin K, Collins A, Murdoch W & Hansen T 1999a Endometrial ISG17 mRNA and a related mRNA are induced by interferon-{tau} and localized to glandular epithelial and stromal cells from pregnant cows. Endocrine 10 243–252.[ISI][Medline]

Johnson GA, Austin KJ, Collins AM, Murdoch WJ & Hansen TR 1999b Endometrial ISG17 mRNA and a related mRNA are induced by interferon-{tau} and localized to glandular epithelial and stromal cells from pregnant cows. Endocrine 10 243–252.[ISI][Medline]

Johnson GA, Burghardt RC, Newton GR, Bazer FW & Spencer TE 1999c Development and characterization of immortalized ovine endometrial cell lines. Biology of Reproduction 61 1324–1330.[Abstract/Free Full Text]

Johnson GA, Spencer TE, Hansen TR, Austin KJ, Burghardt RC & Bazer FW 1999d Expression of the interferon {tau} inducible ubiquitin cross-reactive protein in the ovine uterus. Biology of Reproduction 61 312–318.[Abstract/Free Full Text]

Johnson GA, Spencer TE, Burghardt RC, Joyce MM & Bazer FW 2000 Interferon-{tau} and progesterone regulate ubiquitin cross-reactive protein expression in the ovine uterus. Biology of Reproduction 62 622–627.[Abstract/Free Full Text]

Johnson G, Stewart M, Gray C, Choi Y, Burghardt R, Yu-Lee L, Bazer F & Spencer T 2001a Effects of the estrous cycle, pregnancy, and interferon {tau} on 2', 5'-oligoadenylate synthetase expression in the ovine uterus. Biology of Reproduction 64 1392–1399.[Abstract/Free Full Text]

Johnson GA, Stewart MD, Gray CA, Choi Y, Burghardt RC, Yu-Lee LY, Bazer FW & Spencer TE 2001b Effects of the estrous cycle, pregnancy, and interferon {tau} on 2',5'-oligoadenylate synthetase expression in the ovine uterus. Biology of Reproduction 64 1392–1399.[Abstract/Free Full Text]

Johnson GA, Joyce MM, Yankey SJ, Hansen TR & Ott TL 2002 The interferon stimulated genes (ISG) 17 and Mx have different temporal and spatial expression in the ovine uterus suggesting more complex regulation of the Mx gene. Journal of Endocrinology 174 R7–R11.[Abstract]

Kim S, Choi Y, Bazer FW & Spencer TE 2003a Identification of genes in the ovine endometrium regulated by interferon {tau} independent of signal transducer and activator of transcription 1. Endocrinology 144 5203–5214.[Abstract/Free Full Text]

Kim S, Choi Y, Spencer TE & Bazer FW 2003b Effects of the estrous cycle, pregnancy and interferon {tau} on expression of cyclooxygenase two (COX-2) in ovine endometrium. Reproductive Biology and Endocrinology 1 58.[CrossRef]

Klein C, Bauersachs S, Ulbrich SE, Einspanier R, Meyer HH, Schmidt SE, Reichenbach HD, Vermehren M, Sinowatz F, Blum H et al. 2006 Monozygotic twin model reveals novel embryo-induced transcriptome changes of bovine endometrium in the preattachment period. Biology of Reproduction 74 253–264.[Abstract/Free Full Text]

Knickerbocker JJ & Niswender GD 1989 Characterization of endometrial receptors for ovine trophoblast protein- 1 during the estrous cycle and early pregnancy in sheep. Biology of Reproduction 40 361–369.[Abstract]

Lamming GE, Wathes DC, Flint AP, Payne JH, Stevenson KR & Vallet JL 1995 Local action of trophoblast interferons in suppression of the development of oxytocin and oestradiol receptors in ovine endometrium. Journal of Reproduction and Fertility 105 165–175.[Abstract]

Leaman DW, Chawla-Sarkar M, Jacobs B, Vyas K, Sun Y, Ozdemir A, Yi T, Williams BR & Borden EC 2003 Novel growth and death related interferon-stimulated genes (ISGs) in melanoma: greater potency of IFN-ß compared with IFN-{alpha}2. Journal of Interferon and Cytokine Research 23 745–756.[CrossRef][ISI][Medline]

Li XL & Hassel BA 2001 Involvement of proteasomes in gene induction by interferon and double-stranded RNA. Cytokine 14 247–252.[CrossRef][ISI][Medline]

Li J & Roberts RM 1994 Interferon-{tau} and interferon-{alpha} interact with the same receptors in bovine endometrium. Use of a readily iodinatable form of recombinant interferon-{tau} for binding studies. Journal of Biological Chemistry 269 13544–13550.[Abstract/Free Full Text]

Liu M, Li XL & Hassel BA 2003 Proteasomes modulate conjugation to the ubiquitin-like protein, ISG15. Journal of Biological Chemistry 278 1594–1602.[Abstract/Free Full Text]

Magness RR, Huie JM, Hoyer GL, Huecksteadt TP, Reynolds LP, Seperich GJ, Whysong G & Weems CW 1981 Effect of chronic ipsilateral or contralateral intrauterine infusion of prostaglandin E2 (PGE2) on luteal function of unilaterally ovariectomized ewes. Prostaglandins and Medicine 6 389–401.[CrossRef][ISI][Medline]

McCracken JA, Custer EE & Lamsa JC 1999 Luteolysis: a neuroendocrine-mediated event. Physiological Reviews 79 263–323.[Abstract/Free Full Text]

Mirando MA, Harney JP, Zhou Y, Ogle TF, Ott TL, Moffatt RJ & Bazer FW 1993 Changes in progesterone and oestrogen receptor mRNA and protein and oxytocin receptors in endometrium of ewes after intrauterine injection of ovine trophoblast interferon. Journal of Molecular Endocrinology 10 185–192.[Abstract]

Opferman JT & Korsmeyer SJ 2003 Apoptosis in the development and maintenance of the immune system. Nature Immunology 4 410–415.[CrossRef][ISI][Medline]

Ott TL, Yin J, Wiley AA, Kim HT, Gerami-Naini B, Spencer TE, Bartol FF, Burghardt RC & Bazer FW 1998 Effects of the estrous cycle and early pregnancy on uterine expression of Mx protein in sheep (Ovis aries). Biology of Reproduction 59 784–794.[Abstract/Free Full Text]

Parent J, Villeneuve C, Alexenko AP, Ealy AD & Fortier MA 2003 Influence of different isoforms of recombinant trophoblastic interferons on prostaglandin production in cultured bovine endometrial cells. Biology of Reproduction 68 1035–1043.[Abstract/Free Full Text]

Payne JH & Lamming GE 1994 The direct influence of the embryo on uterine PGF2{alpha} and PGE2 production in sheep. Journal of Reproduction and Fertility 101 737–741.[Abstract]

Pestka S, Krause CD & Walter MR 2004 Interferons, interferon-like cytokines, and their receptors. Immunological Reviews 202 8–32.[CrossRef][ISI][Medline]

Pru JK, Rueda BR, Austin KJ, Thatcher WW, Guzeloglu A & Hansen TR 2001 Interferon-{tau} suppresses prostaglandin F2{alpha} secretion independently of the mitogen-activated protein kinase and nuclear factor {kappa}B pathways. Biology of Reproduction 64 965–973.[Abstract/Free Full Text]

Reik W, Constancia M, Fowden A, Anderson N, Dean W, Ferguson-Smith A, Tycko B & Sibley C 20