Skip to main content

Differential neutrophil gene expression in early bovine pregnancy



In food production animals, especially cattle, the diagnosis of gestation is important because the timing of gestation directly affects the running of farms. Various methods have been used to detect gestation, but none of them are ideal because of problems with the timing of detection or the accuracy, simplicity, or cost of the method. A new method for detecting gestation, which involves assessing interferon-tau (IFNT)-stimulated gene expression in peripheral blood leukocytes (PBL), was recently proposed. PBL fractionation methods were used to examine whether the expression profiles of various PBL populations could be used as reliable diagnostic markers of bovine gestation.


PBL were collected on days 0 (just before artificial insemination), 7, 14, 17, 21, and 28 of gestation. The gene expression levels of the PBL were assessed with microarray analysis and/or quantitative real-time reverse transcription (q) PCR. PBL fractions were collected by flow cytometry or density gradient cell separation using Histopaque 1083 or Ficoll-Conray solutions. The expression levels of four IFNT-stimulated genes, interferon-stimulated protein 15 kDa (ISG15), myxovirus-resistance (MX) 1 and 2, and 2′-5′-oligoadenylate synthetase (OAS1), were then analyzed in each fraction through day 28 of gestation using qPCR.


Microarray analysis detected 72 and 28 genes in whole PBL that were significantly higher on days 14 and 21 of gestation, respectively, than on day 0. The upregulated genes included IFNT-stimulated genes. The expression levels of these genes increased with the progression of gestation until day 21. In flow cytometry experiments, on day 14 the expression levels of all of the genes were significantly higher in the granulocyte fraction than in the other fractions. Their expression gradually decreased through day 28 of gestation. Strong correlations were observed between the expression levels of the four genes in the granulocyte fractions obtained with flow cytometry and with density gradient separation.


The expression profiles of ISG15, MX1, MX2, and OAS1 could be a useful diagnostic biomarker of bovine gestation. Assessing the expression levels of these genes in a granulocyte fraction obtained with density gradient separation is a practical way of detecting gestation in cows within three weeks of insemination.


The early detection of gestation would help improve reproductive efficiency and bring economic benefits for the food animal industry [1, 2]. Various effective methods for gestation confirmation have been developed, such as the detection of serum and/or milk progesterone levels [35], detection of pregnancy-associated glycoproteins [6], observation of estrous behavior, and the rectal or ultrasonic detection of the conceptus/corpus luteum [7, 8].

Assessing the expression levels of interferon-tau (IFNτ)-stimulated genes (ISGs) in peripheral blood leukocytes (PBL) was recently proposed as a new approach to detecting gestation. IFNτ, which is used as a marker of pregnancy in ruminants, stimulates the expression of ISGs, such as interferon-stimulated protein 15 kDa (ISG15), interferon regulatory factors (IRF), and myxovirus-resistance (MX) proteins [915]. This new approach is generally derived from the molecular approach for analyzing the gene expression profiles of trophoblastic and endometrial tissues using microarray and quantitative real-time reverse transcription PCR (qPCR) in ruminants [14, 16, 17]. Although the PBL expression levels of some ISGs, such as ISG15, MX1, MX2, and 2′-5′-oligoadenylate synthetase (OAS1), have been found to be related to gestation [12, 13, 18], the most critical factor for gestation initiation remains to be clarified. In cows, there are reliable and useful methods (see above) for diagnosing gestation [38], and they have been applied in the field. However, diagnostic methods that produce results earlier, i.e., within about two weeks, are desired in order to shorten the open period. We hypothesized that particular PBL populations would display upregulated expression levels of genes related to pregnancy during early gestation in cows.

In the present study, we examined whether any PBL populations displayed upregulated expression levels of pregnancy-related genes during early gestation. First, the gene expression profiles of all PBLs were examined during early gestation using a custom-made bovine oligo microarray. Second, the PBLs were fractionated with flow cytometry or density gradient cell separation, and the gene expression profiles of each fraction were examined by qPCR. Finally, we searched for genes whose PBL expression levels could be used as specific markers of gestation establishment in cows.



The estrous cycles of Japanese Black x Holstein mixed cows (F1, multiparous) were synchronized via the administration of two intramuscular injections of prostaglandin F2α (0.15 mg Dalmazin, Kawasaki Pharmaceutical Co., Tokyo, Japan) every 14 days. Artificial insemination (AI) was performed on the day of estrus using frozen semen from Japanese black bulls. The day of AI was designated as day 0 of gestation, and blood was collected from the jugular vein in heparin sodium-containing vacutainers on days 0, 7, 14, 17, 21, and 28 of gestation. Blood was also collected before AI on day 0. Blood samples were used for microarray, qPCR, gradient, and flow cytometric analysis. The details of blood collection are summarized in Table 1. To analyze whole PBL, 19 cows were used. Fourteen cows were inseminated on day 0, and data were collected separately for the fertile (pregnant; n = 6) and infertile cows (not pregnant despite AI; n = 8). Blood samples were collected during the estrous cycle from the other cows as a control (not inseminated; n = 5). Only 5 out of 6 pregnant cows were used for microarray analysis. To analyze cells separated by flow cytometer or gradient, 10 cows were used. Data were collected separately for the pregnant (n = 5) and control cows (not inseminated, n = 3). Cells from these cows were used in both separation methods. The two other cows, inseminated but non-fertile, were not analyzed. Pregnancy was confirmed by transrectal ultrasonography (HS-1500 V, Honda Electronics, Aichi, Japan) at 4–5 weeks after AI. The guidelines for animal experiments outlined by the animal use and experimental committee of Iwate University (approval no.: 2007–46, 2008–59, and A201043), which approved this study, and the ethical guidelines of the National Institute of Agrobiological Sciences (approval no.: H18-036) were followed during the animal experiments.

Table 1 The number of samples (animals) used for each experiment

Blood cell collection for microarray analysis

After blood collection, 9 mL of whole blood was carefully layered onto 6 mL of Histopaque-1119 (Sigma, Saint Louis, MO, USA) and centrifuged at 780 × g for 45 min at room temperature. The serum fraction was removed. Five hundred microliters of the upper part of the whole blood cell pellet fraction, which included red cells, was transferred into TRIzol LS (Invitrogen, Carlsbad, CA, USA) by pipette for total RNA extraction in accordance with the manufacturer’s protocol. DNase treatment to remove genomic DNA was performed using TURBO DNA-free kits (Ambion, Austin, TX, USA).

Peripheral blood leukocyte (PBL) isolation for flow cytometry

After collection, the whole blood cell fraction was transferred to a new tube and suspended in lysis buffer (155 mM NH4Cl, 10 mM KHCO3, and 1 mM ethylenediaminetetraacetic acid) prepared at 37°C. The sample was then diluted with sorting buffer composed of Hank’s balanced salt solution (HBSS(−), Sigma) containing 2% fetal bovine serum and 10 mM 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid. The sample was inverted several times and immediately centrifuged at 1,200 × g for 10 min. The resultant pellets were resuspended in cold sorting buffer. Alternatively, some of the pellets were loaded onto Ficoll (Sigma)-Conray (Daiichi-Sankyo, Tokyo, Japan) solution (gravity: 1.072) for cell fractionation [19, 20] and RNA extraction, as described below. Approximately 2 × 105 cells were collected for the flow cytometric analysis (EPICS ALTRA MultiCOMP, Beckman Coulter, Carlsbad, CA, USA). The cells were stained with the following specific antibodies: anti-granulocyte antibody (2 μg/1 × 106 cells, MM20A, VMRD Inc., Pullman, WA, USA), anti-bovine monocyte antibody (2 μg/1 × 106 cells, BAQ151, VMRD Inc.), or anti-CD3 antibody (1 μg/1 × 106 cells, MM1A, VMRD Inc.) for 60 min on ice. The cells were then washed and resuspended in sorting buffer containing Alexa488-conjugated anti-mouse IgG (0.05 μg/1 × 106 cells, A-11001, Invitrogen). After 60 min of incubation on ice, the cells were washed, resuspended in staining buffer containing 5 μg/mL propidium iodide, and kept on ice until flow cytometry. Finally, we sorted the bovine PBL into granulocytes (G), monocytes (M), and lymphocytes (L) by flow cytometry based on the forward scatter channel (FSC) and the side scatter channel (SSC). The FSC roughly indicates a cell’s size, and the SSC shows the granularity of the cell. G display high FSC and SSC values, L display low FSC and SSC values, and M display intermediate FSC and SSC values. Based on these patterns, the PBL were sorted into four groups: all blood cells except dead cells (all), G, M, and L. Total RNA was extracted using the RNeasy Micro Kit (QIAGEN, Hilden, Germany) or the RNeasy Mini Kit (QIAGEN) with DNase treatment, in accordance with the manufacturer’s protocols, and used for the subsequent gene expression analysis.

Microarray analysis

A custom-made 15 K bovine oligo DNA microarray developed at our laboratory was used for the microarray analysis (GPL9284), which was performed according to the method of a previous report [21]. After verifying the quality of the RNA with a 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) and NanoDrop ND-1000 spectrophotometer (NanoDrop Technology Inc., Wilmington, DE, USA), we performed one-color microarray analysis. RNA integrity was confirmed; all samples had an A260/280 ratio greater than 1.8 and an RNA integrity number greater than 8.0. The oligo-microarray produced by Agilent Technologies was used in this study. Sixty-mer nucleotide probes for the customized microarray were synthesized on a glass slide. cDNA synthesis, Cy3-labeled cRNA preparation, hybridization, and the washing and scanning of the array slides were performed according to the Agilent one-color microarray-based gene expression analysis protocol. Briefly, 400 ng of total RNA from each sample were reverse-transcribed into cDNA using the Quick Amp Labeling Kit (Agilent Technologies) with an oligo dT-based primer, and then Cy3-labelled cRNA was prepared by in vitro transcription. Labeled cRNA was purified with an RNeasy Mini Kit, and the concentration and Cy3 dye incorporation (pmol Cy3/μg cRNA) were measured with a spectrophotometer. Labeled cRNA (600 ng) was fragmented and hybridized using the Gene Expression Hybridization Kit (Agilent Technologies), according to the manufacturer’s instructions. The arrays were washed using a Gene Expression Wash Pack Kit and scanned using an Agilent Microarray Scanner (Agilent Technologies). Feature Extraction ver. 9.5 (Agilent Technologies) was used for image analysis and data extraction. The microarray data from each sample were imported into GeneSpring 12 (Agilent Technologies) for further data characterization. The GEO accession numbers are as follows. Platform: GPL9284; samples: GSM1052989 to GSM1053001; series: GSE42894.

Quantitative real-time RT-PCR

The expression levels of ISG15, MX1, MX2, and OAS1 were determined by qPCR according to the method described in a previous report [21]. Briefly, the DNase-treated total RNA (0.5 μg) was reverse-transcribed into cDNA using a High Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA, USA) according to the manufacturer’s instructions. The primer pairs for the SYBR Green assay were designed using the Primer Express software (Applied Biosystems). The primers used to amplify each gene are listed in Table 2. Each qPCR reaction contained cDNA template, forward and reverse primers (300 nM each), and half volume of Power SYBR Green PCR Master Mix (Applied Biosystems) according to the manufacturer’s instruction. A non-template control was included. PCR and the resulting relative increase in reporter fluorescent dye emission were monitored in real time using an ABI7300 real-time PCR system (Applied Biosystems). The thermal cycling conditions included initial sample incubation at 50°C for 2 min and 95°C for 10 min, followed by 40 cycles of 95°C for 15 s and 60°C for 1 min. The relative differences in the initial amounts of each cDNA species were determined by comparing their threshold cycle (Ct) values. To quantify the mRNA copy number, standard curves for each gene were generated by serial dilution of the plasmid containing the corresponding cDNA. The dissociation curve for detecting the SYBR Green-based objective amplicon was confirmed because SYBR Green also detects double-stranded DNA, including primer dimers, contaminating DNA, and PCR products from misannealed primers. Contaminating DNA or primer dimers appear as a peak separate from the desired amplicon peak. The mRNA copy number of each gene was determined at each time point. In addition, the ratio of the mRNA expression level of each gene to that of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was also calculated to adjust for variations in the qPCR reaction. We analyzed other genes, such as ribosomal protein L27 and beta-actin, as other internal makers for qPCR, but there was no significant difference. All values are presented as mean ± standard error of the mean (SEM).

Table 2 Primer sequences for quantitative RT-PCR

Granulocyte collection by gradient centrifugation

To collect the G fraction without using flow cytometry, whole blood samples (10 mL) were diluted with an equal volume of HBSS(−) and loaded onto Histopaque 1083 solution (Sigma) or Ficoll-Conray solution in a 50-mL conical tube (BD Falcon), then centrifuged at 1,000 × g for 30 min. After centrifugation, the upper (plasma) and middle (mononucleate cells) fractions were removed, then added 2 mL HBSS(−) and added 0.7% NH4Cl up to 50 mL. After a 5 min incubation on ice, the sample was centrifuged at 1,800 × g for 5 min, and the supernatant was removed. One milliliter of 0.7% NH4Cl was added to the tube, and the sample was incubated for 5 min on ice. The tube was centrifuged at 1,800 × g for 5 min. The resultant pellet was transferred to a new microtube (Eppendorf) for RNA extraction using TRIzol (Invitrogen) and gene analysis. In this experiment, blood samples were collected from ten cows: five pregnant cows (four subjected to blood sample collection throughout the experiment, one subjected to blood sample collection only on day 14) from days 0 to 28 of gestation, two non-pregnant cows (data were not included in the analysis), and three cyclic cows as a control.

Statistical analyses

Microarray data were analyzed statistically with Student’s t-test and summarized using GeneSpring 12 (Agilent Technologies). All qPCR data and extracted microarray data were analyzed with one-way ANOVA followed by the Turkey-Kramer or Dunnett’s multiple comparison test using JMP 7 software (SAS Institute Inc., Cary, NC, USA). The correlations between the data in microarray and qPCR were analyzed with the Spearman’s rank correlation analysis (correlation coefficient: rs) using JMP software because the data were non-parametric. The correlations between the flow cytometric collection and gradient collection were determined by the Pearson product–moment correlation coefficient using JMP software, and the correlation coefficients were calculated as r-values.


Global gene expression profiles of PBL

Gene expression in PBL during early gestation was assessed using an oligo-microarray. On days 14 and 21, the expression of 72 and 28 genes, respectively, was significantly increased, compared with expression on day 0 (see Additional file 1: Table S1 and Additional file 2: Table S2). Some genes that respond to IFNτ (ISG15, and MX2) were only significantly upregulated on day 21.

Microarray analysis showed that in the whole PBL population the expression levels of ISG15, MX1, MX2, and OAS1 increased as gestation progressed (Figure 1). qPCR confirmed the expression of the four ISGs in whole PBL; the copy numbers of ISG15 and MX2 were significantly higher on day 21 of gestation than on day 0 of gestation, and MX1 and OAS1 were increased but not significant (Figure 1B). However, when examined the expression levels of some of the genes in Additional file 1: Tables S1 and Additional file 2: Table S2 (see Additional files 1 and 2) except abovementioned ISGs more closely, we found that there were no common significant genes (over 2 fold changes) between days 14 and 21 of gestation except three geens; FERM domain containing 4A (NM_001192267), low density lipoprotein receptor-related protein 1 (XR_082763) and EST (EE896098), by microarray. The microarray gene expression data and qPCR data showed a strong correlation for ISG15 (rs = 0.709, p = 0.0067), MX1 (rs = 0.571, p = 0.0413), MX2 (rs = 0.692, p = 0.0087), and OAS1 (rs = 0.923, p < 0.001).

Figure 1
figure 1

Microarray (A) and qPCR (B) analysis of IFNτ-stimulated gene expression in PBL. (A) Blood was collected from cows on days 0 (n = 5), 14 (n = 3), and 21 (n = 5) after artificial insemination (AI). Three and five cows samples were anlyzed on day 14 and 21 of gestation, respectively. Each sample was applied to an oligo microarray. Data are shown as mean ± SEM. (B) Open circles represent estrous cycle cows (control, n = 5). Crosses and solid circles represent infertile cows (n = 8) and pregnant cows (n = 6), respectively, which were subjected to AI on day 0 and diagnosed via ultrasonography. Data are shown as mean ± SEM. The copy number of each mRNA is indicated to confirm the minimum amount of mRNA produced. Values significantly different from the value on day 0 of gestation (before AI) are shown with an asterisk. P values less than 0.05 were considered significant. d: day of gestation and/or estrous cycle.

Fractionation of peripheral blood leukocytes and the gene expression levels of each fraction

We hypothesized that the expression levels of genes involved in gestation would vary between different types of PBL. Therefore, the PBL were fractionated using the Ficoll-Conray method, and the fractions were subjected to flow cytometry. The specificity of each fraction was confirmed with flow cytometric analysis using antibodies specific for different PBL cell types (see Additional file 3: Figure S1). Thereafter, cells were identified as G, M, and L according to their FSC and SSC scatter plot profiles without using antibodies (Figure 2) were analyzed in each fraction until day 28 of gestation using qPCR (Figure 3). On day 14, all four genes displayed significantly higher expression levels (p < 0.05) in the G fraction than in the other fractions, including the non-separated blood cells (all, Figure 3). After day 14, the expression of IGS15, MX1, MX2, and OAS1 in the G fraction gradually decreased through day 28 of gestation.

Figure 2
figure 2

Characteristics of cells in the Ficoll-Conray fractions. After obtaining cell fractions with the Histopaque 1119 and Ficoll-Conray gradient methods, the characteristics of the cells and their proportions were determined by flow cytometry using a range of cell collection gates (SSC and FSC). H: whole PBL obtained with Histopaque 1119; FL: lower fraction obtained with Ficoll-Conray; FM: middle fraction obtained with Ficoll-Conray. L: lymphocytes; M: monocytes; G: granulocytes. The numbers indicate the frequency of each cell type. Representative figures are shown.

Figure 3
figure 3

mRNA profiles of each PBL fraction. (A) Profiles of ISG15 (a to d), MX1 (e to h), MX2 (i to l), and OAS1 (m to p) mRNA levels in the four PBL fractions: all blood cells except dead cells (all), granulocytes (G), monocytes (M), and lymphocytes (L) in early pregnant (closed circles; n = 4) and cyclic (open circles; n = 3) cows. Values were adjusted for GAPDH mRNA expression. They represent the relative fold change compared with the value on day 0. The error bars represent the SEM. Values significantly different from the value on day 0 of gestation are shown with an asterisk. P values less than 0.05 were considered significant. (B) ISG15 (a), MX1 (b), MX2 (c), and OAS1 (d) mRNA levels in the PBL fractions on day 14 of pregnancy (n = 5). A: all fractions; G: G fraction; M: M fraction; L: L fraction. Values represent the mean values adjusted for GAPDH mRNA expression, and the error bars represent the SEM. Values significantly different from the value in all fractions are shown with an asterisk (p < 0.05).

Simplification of the granulocyte collection method

When we separated the whole blood cell population into three fractions (upper, middle, and lower fractions) using Ficoll-Conray or Histopaque solutions, we found about 90% of the G in the lower fraction. The accuracy of our G sorting method was confirmed by flow cytometry using anti-bovine G antibody and morphological analysis. The gene expression levels of the four target genes in the G fractions obtained with flow cytometry and with Ficoll-Conray separation displayed high correlation coefficients (see Additional file 4: Figure S2); the correlation coefficients (r) for ISG15, MX1, MX2, and OAS1 were 0.892 (p = 0.0070), 0.883 (p = 0.0084), 0.893 (p = 0.0068), and 0.822 (p = 0.0232), respectively.


ISGs were detected in PBL, and their expression pattern correlated with the onset of gestation. Our results confirm those from previous studies [1215]. The present study showed G fraction was a specific expression PBL cell for ISGs even though other PBL cells expressed ISGs. ISG expression in G may be a suitable diagnostic indicator of early gestation in cows. Methods for the early diagnosis of pregnancy in cattle have been sought because such methods will improve the reproductive and running efficiency of the cattle industry [1, 2]. Although many technologies can diagnose gestation, including rectal palpation, ultrasonography, and milk progesterone tests, these methods are generally applicable three weeks after fertilization [38]. The present study describes a new method for determining pregnancy within three weeks, specifically 14–21 days after AI. This method assesses the expression of ISGs in the G fraction of PBL using gradient separation and qPCR. The technique is more specific and practical than the traditional methods of gestation diagnosis mentioned above.

Microarray analysis is a very useful method for assessing gene expression profiles, and it has become a popular method for examining physiological events in food production animals [16, 17, 2224]. However, recent studies have shown that while microarray analysis is a useful method for determining the gene expression profiles of cattle reproductive tissues, it might be impractical for the routine diagnosis of gestation [1214]. Although we previously described some preliminary results on ISG expression obtained using PBL, the expression of ISGs varied widely [14]. Green et al. [13] carefully examined whether it was possible to use the gene expression profile of leukocytes for gestation diagnosis and concluded that assessing the expression levels of ISGs is a reliable method for pregnancy detection, especially in heifers. In the present study, only multiparous cows were used, but ISG expression appeared to be a suitable indicator. Another group confirmed that the expression of ISGs in peripheral blood mononuclear cells reflects embryonic mortality during early gestation [18]. These findings generally agree with those obtained in the present study, in which ISG expression was more reliable in the G fraction. In addition, microarray analyses in the study by Green et al. [13] and in the present study detected various genes related to gestation events, but the microarray expression profiles obtained in the two studies differed. Differences between strains and in the parous status of the cows may account for the discrepancies.

The PBL gene expression profile is considered a potential source of markers that reflect the physiological changes that occur in early gestation because PBL are involved in various immunological responses. Recent transcriptome analysis studies have shown that PBL are an excellent source of biological markers of physiological conditions in humans and other animals [2527]. For example, gene expression profiles display extensive, rapid, and complex changes that reflect the changes in the levels of active humoral substances in the circulation, such as those caused by infection, inflammation, cytokines, and hormones [28, 29]. Thus, PBL gene expression profiles might be a good source of biomarkers of changes of physiological significance [15].

We analyzed PBL to identify the cell subset that is most sensitive to gestation signals from the conceptus in mother cows. The most important issue is the development of an accurate and simple method for collecting target blood cells. Consistent with previous studies, microarray data suggest that various genes could be used as biomarkers of pregnancy, including ISGs [13, 2124, 3032]. IFNτ stimulates various genes in the bovine endometrium, and some of them are expressed in PBL [1215]. The endometrial response to the presence of a fertilized embryo or the conceptus is surprisingly rapid, and previous studies have suggested that such stimuli affect endometrial cells and blood cells. In the present study, we selected ISG15, MX1, MX2, and OAS1 as potential markers of gestation. The expression levels of these genes increased from days 14 to 21 after AI, but declined by day 28. The changes in the expression levels of these genes might be related to changes in the concentration of IFNτ during gestation [33]. ISG15, MX1, MX2, and OAS1 displayed significantly increased expression levels in the G fraction, especially from days 14 to 21 of gestation, and displayed less significant increases in their expression levels in the whole blood cell population.

The reason and mechanism for the specific stimulation of G by IFNτ are unclear because G have short lives (less than a week) in the peripheral blood. In addition, the increases in the expression levels of ISG15, MX1, MX2, and OAS1 started between days 7 to 14 of gestation, when the bovine conceptus has not begun producing IFNτ [33]. A low level of IFNτ may be sufficient to enhance these genes. Interestingly, IFNτ is a 23-kDa protein that is only produced by the conceptus, which secretes it into the uterus. How then does IFNτ reach blood cells? There is clear evidence that IFNτ passes into the blood via the uterine vein [10, 34], though a detailed study is needed examine whether this large molecule is able to pass from the uterine cavity to the peripheral blood stream. Another question is whether the expression of ISGs is induced by IFNτ or viral infection. ISGs respond to signals from the conceptus and from viruses, such as the bovine viral diarrhea virus (BVDV), which causes infectious diarrhea and mucosal disease in calves [35, 36]. It is difficult to eliminate the effects of BVDV infection on the expression of ISGs following the start of conception. However, G have a short life in the blood stream, perhaps less than a week, and thus they may quickly and specifically respond to gestational signals from the conceptus. In fact, IFNτ signals through the Janus family of tyrosine kinases (JAK) and the signal transducer and activator of transcription (STAT) pathway [37]. The phosphorylation of STAT1 is a critical step in the cell signaling response to type I interferons, including IFNτ. The phosphorylation intensity of STAT1 is different in various peripheral blood white cells; the intensity in L is higher than that in G [38]. These results suggest a possible way to establish whether the signal cascade derives from BVDV because the expression intensities of ISGs were specifically higher in G than in L in the present study. This may mean that the responsiveness and sensitivity of PBL to the ISG signaling pathway depends on the cell type. Therefore, ISG expression in the G fraction may be a possible indicator for early gestation diagnosis in cattle, but it does not provide any direct evidence about how IFNτ stimulates the expression of ISGs in circulating blood cells, specifically G. Further examinations are required to clarify the sensitivity of G to viral infections and gestation signals, mainly IFNτ. We also detected more than 20 new genes other than ISGs that increased after the onset of gestation (Additional file 1: Table S1 and Additional file 2: Table S2). They may be new indicator genes. Further studies are necessary to find new makers.

Flow cytometric analysis showed that the main IFNτ-responsive blood cell during early gestation in bovines was G. This is the first evidence that G specifically respond to gestation signals, including IFNτ. However, although cell fractionation with a flow cytometer is an accurate method, it is not practical for everyday use. A simpler and more convenient method is required for routine use. Thus, we used two density gradient solutions, Ficoll-Conray and Histopaque, as a cell collection method. Both solutions worked well. The accuracy of these methods was confirmed by morphological examinations using Giemsa staining and flow cytometry. Gradient separation is an accurate and simple way to isolate PBL, especially the G fraction.


Flow cytometric separation of PBL suggested that G are most sensitive to IFNτ stimulation, and that the expression levels of ISGs in the G fraction might be good biomarkers for diagnosing pregnancy in bovines very early during gestation. The density gradient cell method for the collection of PBL is a useful way of identifying G. It represents a simpler and more practical method for diagnosing bovine pregnancy in cows within three weeks of insemination.



Bovine viral diarrhea virus






Interferon-stimulated gene








Quantitative real-time reverse transcription PCR


2′-5′-oligoadenylate synthetase


Peripheral blood leukocytes


Signal transducer and activator of transcription.


  1. Oltenacu PA, Ferguson JD, Lednor AJ: Economic evaluation of pregnancy diagnosis in dairy cattle: a decision analysis approach. J Dairy Sci. 1990, 73: 2826-2831. 10.3168/jds.S0022-0302(90)78970-9.

    Article  CAS  PubMed  Google Scholar 

  2. de Vries A, Bartolome J, Broaddus B: What is early pregnancy diagnosis worth?. Proceedings of the 2nd Florida Dairy Road Show: 7–11 February 2005. Edited by: Vries A. 2005, Florida, 31-41.

    Google Scholar 

  3. Nakao T, Sugihashi A, Kawata K, Saga N, Tsunoda N: Milk progesterone levels in cows with normal or prolonged estrous cycles, referenced to an early pregnancy diagnosis. Jpn J Vet Sci. 1983, 45: 495-499. 10.1292/jvms1939.45.495.

    Article  CAS  Google Scholar 

  4. Morton JM, Wynn PC: Assessing ovulation detection performance in commercial dairy herds using progesterone concentrations from limited numbers of strategically collected milk samples. J Dairy Sci. 2010, 93: 3019-3030. 10.3168/jds.2009-2938.

    Article  CAS  PubMed  Google Scholar 

  5. Okumu LA, Forde N, Fahey AG, Fitzpatrick E, Roche JF, Crowe MA, Lonergan P: The effect of elevated progesterone and pregnancy status on mRNA expression and localisation of progesterone and oestrogen receptors in the bovine uterus. Reproduction. 2010, 140: 143-153. 10.1530/REP-10-0113.

    Article  CAS  PubMed  Google Scholar 

  6. Zoli AP, Guilbault LA, Delahaut P, Ortiz WB, Beckers JF: Radioimmunoassay of a bovine pregnancy-associated glycoprotein in serum: its application for pregnancy diagnosis. Biol Reprod. 1992, 46: 83-92. 10.1095/biolreprod46.1.83.

    Article  CAS  PubMed  Google Scholar 

  7. Fricke PM: Scanning the future–ultrasonography as a reproductive management tool for dairy cattle. J Dairy Sci. 2002, 85: 1918-1926. 10.3168/jds.S0022-0302(02)74268-9.

    Article  CAS  PubMed  Google Scholar 

  8. Romano JE, Thompson JA, Kraemer DC, Westhusin ME, Forrest DW, Tomaszweski MA: Early pregnancy diagnosis by palpation per rectum: influence on embryo/fetal viability in dairy cattle. Theriogenology. 2007, 67: 486-493. 10.1016/j.theriogenology.2006.08.011.

    Article  PubMed  Google Scholar 

  9. Rempel LA, Francis BR, Austin KJ, Hansen TR: Isolation and sequence of an interferon-tau-inducible, pregnancy- and bovine interferon-stimulated gene product 15 (ISG15)-specific, bovine ubiquitin-activating E1-like (UBE1L) enzyme. Biol Reprod. 2005, 72: 365-372. 10.1095/biolreprod.104.033027.

    Article  CAS  PubMed  Google Scholar 

  10. Bott RC, Ashley RL, Henkes LE, Antoniazzi AQ, Bruemmer JE, Niswender GD, Bazer FW, Spencer TE, Smirnova NP, Anthony RV, et al: Uterine vein infusion of interferon tau (IFNT) extends luteal life span in ewes. Biol Reprod. 2010, 82: 725-735. 10.1095/biolreprod.109.079467.

    Article  CAS  PubMed  Google Scholar 

  11. Yang L, Wang XL, Wan PC, Zhang LY, Wu Y, Tang DW, Zeng SM: Up-regulation of expression of interferon-stimulated gene 15 in the bovine corpus luteum during early pregnancy. J Dairy Sci. 2010, 93: 1000-1011. 10.3168/jds.2009-2529.

    Article  CAS  PubMed  Google Scholar 

  12. Han H, Austin KJ, Rempel LA, Hansen TR: Low blood ISG15 mRNA and progesterone levels are predictive of non-pregnant dairy cows. J Endocrinol. 2006, 191: 505-512. 10.1677/joe.1.07015.

    Article  CAS  PubMed  Google Scholar 

  13. Green JC, Okamura CS, Poock SE, Lucy MC: Measurement of interferon-tau (IFN-tau) stimulated gene expression in blood leukocytes for pregnancy diagnosis within 18-20d after insemination in dairy cattle. Anim Reprod Sci. 2010, 121: 24-33.

    Article  CAS  PubMed  Google Scholar 

  14. Hashizume K: Diagnosis of gestation using genes expression profiles in bovine peripheral white blood cells. Livestock Technol (Chikusan Gijutu, in Japanese). 2010, 5: 6-10.

    Google Scholar 

  15. Ott TL, Gifford CA: Effects of early conceptus signals on circulating immune cells: lessons from domestic ruminants. Am J Reprod Immunol. 2010, 64: 245-254. 10.1111/j.1600-0897.2010.00912.x.

    Article  CAS  PubMed  Google Scholar 

  16. Ishiwata H, Katsuma S, Kizaki K, Patel OV, Nakano H, Takahashi T, Imai K, Hirasawa A, Shiojima S, Ikawa H, et al: Characterization of gene expression profiles in early bovine pregnancy using a custom cDNA microarray. Mol Reprod Dev. 2003, 65: 9-18. 10.1002/mrd.10292.

    Article  CAS  PubMed  Google Scholar 

  17. Ushizawa K, Herath CB, Kaneyama K, Shiojima S, Hirasawa A, Takahashi T, Imai K, Ochiai K, Tokunaga T, Tsunoda Y, Tsujimoto G, Hashizume K: cDNA microarray analysis of bovine embryo gene expression profiles during the pre-implantation period. Reprod Biol Endocrinol. 2004, 2: 77-10.1186/1477-7827-2-77.

    Article  PubMed Central  PubMed  Google Scholar 

  18. Matsuyama S, Kojima T, Kato S, Kimura K: Relationship between quantity of IFNT estimated by IFN-stimulated gene expression in peripheral blood mononuclear cells and bovine embryonic mortality after AI or ET. Reprod Biol Endocrinol. 2012, 10: 21-10.1186/1477-7827-10-21.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Okeda T, Ono J, Takaki R, Todo S: Simple mrthod for the collection of pancreatic islets by the use of Ficoll-Conray gradient. Endocrinologia Japonica. 1979, 26: 495-499. 10.1507/endocrj1954.26.495.

    Article  CAS  PubMed  Google Scholar 

  20. Takaki R, Ono J: Preparation and culture of isolated islets and dissociated islet cells from rodent pancreas. J Tissue Culture Method. 1985, 9: 61-66. 10.1007/BF01797776.

    Article  Google Scholar 

  21. Sugawara K, Kizaki K, Herath CB, Hasegawa Y, Hashizume K: Transforming growth factor beta family expression at the bovine feto-maternal interface. Reprod Biol Endocrinol. 2010, 8: 120-10.1186/1477-7827-8-120.

    Article  PubMed Central  PubMed  Google Scholar 

  22. Klein C, Bauersachs S, Ulbrich SE, Einspanier R, Meyer HH, Schmidt SE, Reichenbach HD, Vermehren M, Sinowatz F, Blum H, et al: Monozygotic twin model reveals novel embryo-induced transcriptome changes of bovine endometrium in the preattachment period. Biol Reprod. 2006, 74: 253-264. 10.1095/biolreprod.105.046748.

    Article  CAS  PubMed  Google Scholar 

  23. Forde N, Carter F, Spencer TE, Bazer FW, Sandra O, Mansouri-Attia N, Okumu LA, McGettigan PA, Mehta JP, McBride R, et al: Conceptus-induced changes in the endometrial transcriptome: how soon does the cow know she is pregnant?. Biol Reprod. 2011, 85: 144-156. 10.1095/biolreprod.110.090019.

    Article  CAS  PubMed  Google Scholar 

  24. Bauersachs S, Mitko K, Ulbrich SE, Blum H, Wolf E: Transcriptome studies of bovine endometrium reveal molecular profiles characteristic for specific stages of estrous cycle and early pregnancy. Exp Clin Endocrinol Diabetes. 2008, 116: 371-384. 10.1055/s-2008-1076714.

    Article  CAS  PubMed  Google Scholar 

  25. Huck B, Steck T, Habersack M, Dietl J, Kämmerer U: Pregnancy associated hormones modulate the cytokine production but not the phenotype of PBMC-derived human dendritic cells. Eur J Obstet Gyn Reprod Biol. 2005, 122: 85-94. 10.1016/j.ejogrb.2005.02.017.

    Article  CAS  Google Scholar 

  26. Ledderose C, Heyn J, Limbeck E, Kreth S: Selection of reliable reference genes for quantitative real-time PCR in human T cells and neutrophils. BMC Res Notes. 2001, 4: 427-

    Article  Google Scholar 

  27. Galvão KN, Felippe MJ, Brittin SB, Sper R, Fraga M, Galvão JS, Caixeta L, Guard CL, Ricci A, Gilbert RO: Evaluation of expression by blood monocytes of lactating Holstein cows with or without postpartum uterine disease. Theriogenology. 2012, 77: 356-372. 10.1016/j.theriogenology.2011.08.008.

    Article  PubMed  Google Scholar 

  28. Burton JL, Madsen SA, Chang LC, Weber PS, Buckham KR, van Dorp R, Hickey MC, Earley B: Gene expression signatures in neutrophils exposed to glucocorticoids: a new paradigm to help explain “neutrophil dysfunction” in parturient dairy cows. Vet Immunol Immunopath. 2005, 105: 197-219. 10.1016/j.vetimm.2005.02.012.

    Article  CAS  Google Scholar 

  29. Sukumaran B, Carlyon JA, Cai JL, Berliner N, Fikrig E: Early transcriptional response of human neutrophils to Anaplasma phagocytophilum infection. Infect Immun. 2005, 73: 8089-8099. 10.1128/IAI.73.12.8089-8099.2005.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Austin KJ, Carr AL, Pru JK, Hearne CE, George EL, Belden EL, Hansen TR: Localization of ISG15 and conjugated proteins in bovine endometrium using immunohistochemistry and electron microscopy. Endocrinology. 2004, 145: 967-975.

    Article  CAS  PubMed  Google Scholar 

  31. Song G, Bazer FW, Spencer TE: Pregnancy and interferonτ regulate RSAD2 and IFIH1 expression in the ovine uterus. Reproduction. 2007, 133: 285-295. 10.1530/REP-06-0092.

    Article  CAS  PubMed  Google Scholar 

  32. Bauersachs S, Ulbrich SE, Zakhartchenko V, Minten M, Reichenbach M, Reichenbach HD, Blum H, Spencer TE, Wolf E: The endometrium responds differently to cloned versus fertilized embryos. Proc Natl Acad Sci U S A. 2009, 106: 5681-5686. 10.1073/pnas.0811841106.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Roberts MR, Chen Y, Ezashi T, Walker AM: Interferons and the maternal-conceptus dialog in mammals. Semin Cell Dev Biol. 2008, 19: 170-177. 10.1016/j.semcdb.2007.10.007.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Oliveira JF, Henkes LE, Ashley RL, Purcell SH, Smirnova NP, Veeramachaneni DN, Anthony RV, Hansen TR: Expression of interferon (IFN)-stimulated genes in extrauterine tissues during early pregnancy in sheep is the consequence of endocrine IFN-tau release from the uterine vein. Endocrinology. 2008, 149: 1252-1259.

    Article  CAS  PubMed  Google Scholar 

  35. Hansen TR, Smirnova NP, Van Campen H, Shoemaker ML, Ptitsyn AA, Bielefeldt-Ohmann H: Maternal and fetal response to fetal persistent infection with bovine viral diarrhea virus. Am J Reprod Immunol. 2010, 64: 295-306. 10.1111/j.1600-0897.2010.00904.x.

    Article  CAS  PubMed  Google Scholar 

  36. Weiner CM, Smirnova NP, Webb BT, Van Campen H, Hansen TR: Interferon stimulated genes, CXCR4 and immune cell responses in peripheral blood mononuclear cells infected with bovine viral diarrhea virus. Res Vet Sci. 2012, 93: 1081-1088. 10.1016/j.rvsc.2012.01.011.

    Article  CAS  PubMed  Google Scholar 

  37. Samuel CE: Antiviral actions of interferons. Clin Microbiol Rev. 2001, 14: 778-809. 10.1128/CMR.14.4.778-809.2001.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Tochizawa S, Ohmoto Y, Mori T: A novel modification of a flow cytometric assay of phosphorylated STAT1 in whole blood lymphocytes for rapid detection of interferon-alpha signal in vivo. J Immunol Methods. 2006, 313: 29-37. 10.1016/j.jim.2006.03.015.

    Article  CAS  PubMed  Google Scholar 

Download references


The authors sincerely thank Drs. Kei Imai, Yutaka Hashiyada, and Mr. Hideo Matsuda at the National Livestock Breeding Center, Fukushima, Japan for collecting blood samples. This research was supported by a grant (Hogakenkyu 19658101) from the Ministry of Education, Culture, Sport, Science, and Technology of Japan.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Kazuyoshi Hashizume.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contribution

KK collected samples, conducted microarray, planned the study, and drafted the manuscript. AS-K collected samples and conducted microarray, flow cytometric analysis, and qRT-PCR. TF carried out flow cytometric analysis. TT and MH collected samples and prepared animals. KH designed the experiments, prepared samples, and wrote the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Table S1. Genes whose expression increased significantly between D0 and D14 of gestation. Fold change from D0 to D14 assessed by microarray analysis (n = 3). (DOCX 33 KB)


Additional file 2: Table S2. Genes whose expression increased significantly between D0 and D21 of gestation. Fold change from D0 to D21 assessed by microarray analysis (n = 5). (DOCX 30 KB)


Additional file 3: Figure S1. Separation of blood cells with specific antibodies using flow cytometry. H: whole blood white cells; FL: lower fraction obtained with Ficoll-Conray separation; FM: middle fraction obtained with Ficoll-Conray separation; NS: no stain; IgG: collected with IgG; G: collected with anti-granulocyte antibody; M7: collected with anti-monocyte antibody; CD3: collected with anti-CD3 antibody. The numbers indicate the percentage of cells in each antibody-derived fraction. They are the mean of four independent examinations (n = 4). Representative figures are shown. (TIFF 3 MB)


Additional file 4: Figure S2. Correlations between gene expression profiles in granulocyte fractions obtained by flow cytometry and Ficoll-Conray separation. Blood samples were divided and simultaneously analyzed using the two different methods. G: granulocyte, FL: flow cytometry. Data were normalized to GAPDH expression (n = 7). (TIFF 3 MB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Kizaki, K., Shichijo-Kizaki, A., Furusawa, T. et al. Differential neutrophil gene expression in early bovine pregnancy. Reprod Biol Endocrinol 11, 6 (2013).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: