# Design and evaluation of a unique SYBR Green real-time RT-PCR assay for quantification of five major cytokines in cattle, sheep and goats

Carinne Puech, Laurence Dedieu, Isabelle Chantal, Valérie Rodrigues

► To cite this version:

Carinne Puech, Laurence Dedieu, Isabelle Chantal, Valérie Rodrigues. Design and evaluation of a unique SYBR Green real-time RT-PCR assay for quantification of five major cytokines in cattle, sheep and goats. BMC Veterinary Research, 2015, 11, 10.1186/s12917-015-0382-0. hal-01512149

HAL Id: hal-01512149

https://hal.science/hal-01512149

Submitted on 21 Apr 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Design and evaluation of a unique SYBR Green real-time RT-PCR assay for quantification of five major cytokines in cattle, sheep and goats

Carinne Puech¹,⁴*, Laurence Dedieu², Isabelle Chantal³ and Valérie Rodrigues¹,⁴

## Abstract

**Background:** Today, when more than 60% of animal diseases are zoonotic, understanding their origin and development and identifying protective immune responses in ruminants are major challenges. Robust, efficient and cost-effective tools are preconditions to solve these challenges. Cytokines play a key role in the main mechanisms by which the immune system is balanced in response to infectious pathogens. The cytokine balance has thus become the focus of research to characterize immune response in ruminants. Currently, SYBR Green reverse transcriptase quantitative PCR (RT-qPCR) is the most widely method used to investigate cytokine gene expression in ruminants, but the conditions in which the many assays are carried out vary considerably and need to be properly evaluated. Accordingly, the quantification of gene expression by RT-qPCR requires normalization by multiple reference genes. The objective of the present study was thus to develop an RT-qPCR assay to simultaneously quantify the expression of several cytokines and reference genes in three ruminant species. In this paper, we detail each stage of the experimental protocol, check validation parameters and report assay performances, following MIQE guidelines.

**Results:** Ten novel primer sets were designed to quantify five cytokine genes (IL-4, IL-10, IL-12B, IFN-γ and TNF-α) and five reference genes (ACTB, GAPDH, H3F3A, PPIA and YWHAZ) in cattle, sheep, and goats. All the primer sets were designed to span exon-exon boundaries and use the same hybridization temperature. Each stage of the RT-qPCR method was detailed; their specificity and efficiency checked, proved and are reported here, demonstrating the reproducibility of our method, which is capable of detecting low levels of cytokine mRNA up to one copy whatever the species. Finally, we checked the stability of candidate reference gene expression, performed absolute quantification of cytokine and reference gene mRNA in whole blood samples and relative expression of cytokine mRNA in stimulated PBMC samples.

**Conclusions:** We have developed a novel RT-qPCR assay for the simultaneous relative quantification of five major cytokines in cattle, sheep and goats, and their accurate normalization by five reference genes. This accurate and easily reproducible tool can be used to investigate ruminant immune responses and is widely accessible to the veterinary research community.

**Keywords:** Cytokine expression, SYBR Green RT-qPCR, Cattle, Sheep, Goats, Immune response, Th1/Th2 response, Reference genes

<footnote>* Correspondence: carinne.puech@cirad.fr
¹INRA, UMR1309 CMAEE, Montpellier F-34398, France
²CIRAD, UMR CMAEE, Montpellier F-34398, France
Full list of author information is available at the end of the article</footnote>
## Background

Animal diseases have major economic, sanitary, health and environmental consequences worldwide. In 2013, the World Organization for Animal Health (OIE) calculated that more than 60% of animal diseases are zoonotic and that 75% of them are responsible for human emerging infections [1]. Infectious diseases are at the animal-human-ecosystem interface and are a major public health concern, particularly emerging infectious diseases [2]. The importance of ruminant livestock in sustainable agricultural systems [3] and the impact of major livestock diseases on food production is no longer the subject of debate [4]. However, developing efficient control strategies is a major challenge, and understanding disease pathogenesis and identifying protective immune responses are prerequisites, requiring robust, efficient, and cost-effective tools.

In response to pathogens, the immune system is regulated by complex mechanisms in which cytokines play a key role. They control the activation, proliferation and/or differentiation of different types of cells as well as in the secretion of antibodies and mediators [5]. Differential cytokine production and expression, in particular the Th1/Th2 balance, has thus become a widely used way to characterize disease pathogenesis in cattle [6,7], goats [8] and sheep [9,10].

The cytokine balance is routinely characterized either by quantification of secreted cytokine proteins or by the expression of cytokine messenger RNA (mRNA). A few optimized assays are currently available for the measurement of bovine and ovine cytokines secreted: ELISA kits, monoclonal antibodies and standards have been developed and marketed including interferon-gamma, (IFN-γ), tumor necrosis factor-alpha (TNF-α) [11], interleukin (IL)-10 [12], IL-12 [13] and IL-4 [14]. For the measurement of caprine cytokines secreted, most of the few specific tools available are based on cross-reactivity with other ruminants, but are not always sure [15]. Using the ELISA method for ruminants is costly and only a limited number of cytokines can be analyzed from a single sample. In addition, despite technological advances, veterinary immunologists working on ruminants lack tools to investigate host immunity [15-17]. Cost-effective tools able to quantify several cytokines in several species of ruminants would be a significant advance.

To achieve this objective, analyzing cytokine gene expression has become a widely used method to establish a cytokine profile in ruminants. Currently, reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) is the routine method used to investigate gene expression. RT-qPCR makes it possible to simultaneously characterize different genes, using small quantities of sample with high specificity, sensitivity and accuracy [18,19]. Beside the specific probe strategy, the use of SYBR Green dye is also a sensitive, robust and

reproducible method, provided that primer design and set up conditions are properly optimized. In addition, the SYBR Green method is cheaper than methods using probes, and is thus more widely and easily transferable. However, the wide range of published studies on ruminants cites many different procedures and conditions. Most of these assays are specific to only one species and some were even published without reporting all validation parameters.

The objective of our study was thus to develop an innovative and cost-effective RT-qPCR assay, to simultaneously quantify the expression of five major cytokines involved in Th1/Th2 responses, IL-4, IL-10, IL-12B, IFN-γ and TNF-α, using a single primer set and hybridization temperature for three species of ruminants, cattle, sheep and goats, in a single assay. Here we detail each stage of the experimental protocol used for the development of our SYBR Green RT-qPCR assay and offer a reliable and transferable tool. These stages range from template preparation, nucleic acid quality, reverse transcription (RT) strategy up to the PCR conditions and data reporting [20], and checking validation parameters and assay performances following the recommendations of the Minimum Information for publication of Quantitative real-time PCR Experiments (MIQE) guidelines [21]. We carefully selected five reference genes for accurate quantification of gene expression by RT-qPCR and checked the stability of their expression. The identification of internal controls to normalize gene expression is essential for relative gene expression [22], and normalization by multiple reference genes instead of one is recommended [23]. For that reason, we chose the reference genes most commonly used and/or described as stable in different types of samples. *Glyceraldehyde 3-phosphate dehydrogenase* (GAPDH) and *beta-actin* (ACTB) are the most commonly used reference genes. These genes have been described as optimal reference genes in bovine peripheral blood mononuclear cells (PBMC) [24] or whole blood [25] and thus can be used to normalize gene expression. Previous studies also showed that *Tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide* (YWAHZ) [26,27], *Peptidylprolyl Isomerase A* (PPIA) [26,28] or *H3-histone family 3A* (H3F3A) [29] are the most stable genes in bovine peripheral lymphocytes and ovine whole blood.

The SYBR Green RT-qPCR assay we developed can be used for cattle, sheep and goats to quantify the expression of five major cytokine genes (IL-4, IL-10, IL-12B, IFN-γ and TNF-α) and five candidate reference genes (GAPDH, ACTB, YWHAZ, PPIA and H3F3A). Our SYBR Green RT-qPCR assay is based on an optimized and validated protocol and is an easily reproducible and reliable tool specifically designed to investigate immune response in ruminants.
## Methods

In this paper, we use the nomenclature proposed by the Minimum Information for publication of quantitative real-time PCR experiments (MIQE) guidelines [19]. Candidate genes used for normalization are referred to as reference genes and the fractional PCR cycle used for quantification was Cq to quantification cycle.

## Animals, blood collection and ethical considerations

Samples from healthy crossbreed cattle, Suffolk sheep and Saanen goats (three animals per species) were collected in an animal housing facility (Montpellier-France). Whole blood was collected from the jugular vein either with Tempus™ Blood RNA tubes (Applied Biosystems Ltd., Warrington, UK) or heparinized BD Vacutainer® tubes (Beckton Dickinson, New Jersey, USA) according to the manufacturer's instructions. Experimental procedures for animal maintenance and blood sampling were approved by the Languedoc-Roussillon regional ethics committee (French CE-LR #36) in the Authorised Project using animals for scientific purposes #12ANI01.

## Tempus™ whole blood total RNA samples

Total RNA was extracted from whole blood of cattle, sheep and goats with Tempus™ Blood RNA System (Applied Biosystems®, Warrington, UK). After blood collection, RNA was immediately extracted using the Tempus™ Spin RNA Isolation Reagent Kit according to the manufacturer's instructions. All the samples were treated with DNase, AbsoluteRNA Wash Solution, as recommended in the kit. RNA samples were stored at -80°C until conversion into cDNA. These samples are henceforth referred to as « Tempus™ total RNA ».

## Isolation, stimulation and RNA extraction of PBMC

Peripheral blood mononuclear cells (PBMC) were isolated from heparinized whole blood by density gradient centrifugation. Briefly, PBMCs were collected after centrifugation on Histopaque®-1077 (Sigma Aldrich, France) for sheep and goats and on Histopaque®-1083 (Sigma Aldrich, France) for cattle, according to the manufacturer's instructions. Cell viability was assessed by Trypan blue exclusion. PBMCs were resuspended in RPMI 1640 Medium-glutaMAX™ (Life Technologies™, USA) supplemented with 5.10⁻⁵ M βmercaptoethanol (Life Technologies™, USA), 50 µg/ml gentamycin (Life Technologies™, USA) and 10% heat inactivated fetal calf serum (FCS; Eurobio AbCys, France) and were seeded in 12-well plates at 2.10⁶ cells/ml. PBMCs were cultured in medium (unstimulated condition) or stimulated with 5 µg/ml of Concanavalin A (stimulated condition) for 36 h at 37°C and 5% CO₂. Total RNA was extracted using the RNeasy® Mini Kit (Qiagen Ltd., Crawley, UK) and treated with RNase-free DNase Set (Qiagen Ltd., Crawley, UK) for 30 min at

room temperature. RNA samples were stored at -80°C until conversion into cDNA.

## RNA quantification and quality control

The purity and quantity of RNA were assessed using a NanoDrop™ ND-1000 Spectrophotometer (Thermo Fisher Scientific, MA, USA). The A260:A280 ratio ranged from 2.1 to 2.2. RNA integrity was checked with the Agilent 2100 Bioanalyzer using the RNA 6000 Nano Assay Kit (Agilent Technologies, Inc. Santa Clara, USA). The RNA Integrity Number (RIN) ranged from 7 to 9.6.

## Reverse transcription reaction

Two hundred nanograms of total RNA from each sample were reverse transcribed with the AffinityScript QPCR cDNA synthesis kit (Agilent Technologies, Inc. Santa Clara, USA) using oligo-dT strategy, and according to the manufacturer's instructions. cDNA samples were stored as multiple aliquots at -20°C for subsequent use. No-reverse transcription (no-RT) controls (assay without reverse transcriptase) were also prepared to check for non-specific amplification.

## Cytokine and gene reference primers set design

All oligonucleotides were synthesized by Eurogentec (Seraing, Belgium).

Primers were designed by hand to span exon-exon boundaries and to fulfill the following criteria: located near the 3' end, GC% between 40 and 70%, giving an approximately 200 bp amplicon and a melting temperature between 62 and 65°C. Exon spanning primers were designed by aligning *Bos taurus* gene sequences from National Center for Biotechnology Information (NCBI) GenBank database (btau 4.6.1) with mRNA-to-genomic alignment program, Spidey (www.ncbi.nlm.nih.gov/spidey). Species homologies and exon positions of ovine and caprine primers were checked in *silico*, respectively, with *Capra hircus* (chir_1.0) and *Ovis aries* (oar_V3.1) gene sequences from NCBI GenBank. The absence of primer-dimers and secondary structures was checked in *silico* using the OligoCalculator, an online oligonucleotide sequence calculator (Sigma Aldrich Co., St. Louis, USA).

Accession numbers and the exon location of each species are listed in Table 1 and primer sequences in Table 2.

## Checking specificity with conventional PCR assays

The specificity of each primer pair was checked in preliminary conventional PCR assays with bovine, ovine and caprine cDNA synthesized from Tempus™ total RNA. All PCR reactions were performed with SureStart Taq DNA Polymerase (Agilent Technologies, Inc. Santa Clara, USA) following the manufacturer's instructions, with 10 ng of cDNA in 0.5 µM of each primer. PCR
<table><caption>Table 1 Accession number and exon locations of primers for bovine, caprine and ovine gene expression</caption><thead><tr><th rowspan="2">Target gene</th><th rowspan="2">Species</th><th rowspan="2">Accession number<sup>a</sup></th><th colspan="2">Exon position</th></tr><tr><th>Primer F</th><th>Primer R</th></tr></thead><tbody><tr><td rowspan="3">IL-4</td><td>Bos taurus</td><td>NM_173921</td><td>E2</td><td>E3-E4</td></tr><tr><td>Capra hircus</td><td>NM_001285681</td><td>E2</td><td>E3-E4</td></tr><tr><td>Ovis aries</td><td>NM_001009313</td><td>E2</td><td>E3-E4</td></tr><tr><td rowspan="3">IL-10</td><td>Bos taurus</td><td>NM_174088</td><td>E2-E3</td><td>E5</td></tr><tr><td>Capra hircus</td><td>XM_005690416</td><td>E2-E3</td><td>E4-E5</td></tr><tr><td>Ovis aries</td><td>NM_001009327</td><td>E2-E3</td><td>E4-E5</td></tr><tr><td rowspan="3">IL-12B</td><td>Bos taurus</td><td>NM_174356</td><td>E3-E4</td><td>E4-E5</td></tr><tr><td>Capra hircus</td><td>NM_001285700</td><td>E3-E4</td><td>E4-E5</td></tr><tr><td>Ovis aries</td><td>NM_001009438</td><td>E3-E4</td><td>E4-E5</td></tr><tr><td rowspan="3">INF-γ</td><td>Bos taurus</td><td>NM_174086</td><td>E3</td><td>E3-E4</td></tr><tr><td>Capra hircus</td><td>NM_001285682</td><td>E2-E3</td><td>E3</td></tr><tr><td>Ovis aries</td><td>NM_001009803</td><td>E2-E3</td><td>E3-E4</td></tr><tr><td rowspan="3">TNF-α</td><td>Bos taurus</td><td>NM_173966</td><td>E1-E2</td><td>E3-E4</td></tr><tr><td>Capra hircus</td><td>NM_001286442</td><td>E1-E2</td><td>E3-E4</td></tr><tr><td>Ovis aries</td><td>NM_001024860</td><td>E1-E2</td><td>E3-E4</td></tr><tr><td rowspan="3">GAPDH</td><td>Bos taurus</td><td>NM_001034034</td><td>E2-3</td><td>E4-5</td></tr><tr><td>Capra hircus</td><td>XM_005680968</td><td>E2-E3</td><td>E4-5</td></tr><tr><td>Ovis aries</td><td>NM_001190390</td><td>E2-E3</td><td>E4-5</td></tr><tr><td rowspan="3">H3F3A</td><td>Bos taurus</td><td>NM_00101489</td><td>E2</td><td>E2-E3</td></tr><tr><td>Capra hircus</td><td>XM_005690530</td><td>E1</td><td>E1-E2</td></tr><tr><td>Ovis aries</td><td>XM_004013633</td><td>E2</td><td>E2-E3</td></tr><tr><td rowspan="3">ACTB</td><td>Bos taurus</td><td>NM_173979</td><td>E4-E5</td><td>E5-E6</td></tr><tr><td>Capra hircus</td><td>XM_005694067</td><td>E3-E4</td><td>E4-E5</td></tr><tr><td>Ovis aries</td><td>NM_001009784</td><td>E4-E5</td><td>E5-E6</td></tr><tr><td rowspan="3">PPIA</td><td>Bos taurus</td><td>NM_178320</td><td>E4</td><td>E4-E5</td></tr><tr><td>Capra hircus</td><td>XM_005679322</td><td>E3</td><td>E3-E4</td></tr><tr><td>Ovis aries</td><td>XM_004013990</td><td>-</td><td>-</td></tr><tr><td rowspan="3">YWHAZ</td><td>Bos taurus</td><td>NM_174814</td><td>E3-E4</td><td>E4-E5</td></tr><tr><td>Capra hircus</td><td>XM_005689196</td><td>E3-E4</td><td>E4-E5</td></tr><tr><td>Ovis aries</td><td>NM_001267887</td><td>E2-E3</td><td>E3-E4</td></tr></tbody></table>

Abbreviations: F: forward; R: reverse; bp: base pair; IL, interleukin; IL12B, interleukin p40; TNF-α tumor necrosis factor-alpha; IFN-γ interferon-gamma; ACTB, beta-actin; GAPDH, Glyceraldehyde 3-phosphate dehydrogenase; H3F3A, H3-histone family 3A; PPIA, peptidylprolyl isomerase A; YWHAZ, tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide;

reactions were conducted on a SureCycler 8800 Thermal Cycler (Agilent Technologies, Santa Clara, USA). The PCR program consisted in an initial denaturation step of 10 min at 95°C, followed by 40 cycles of 30 sec at 95°C, 30 sec at 60°C, and 45 sec at 72°C, followed by a final elongation step of 10 min at 72°C. PCR products were analyzed on 2% agarose gel. Amplicon size was checked with agarose electrophoresis migration (amplicon sizes are listed in Table 2). Finally, PCR products were (i)

<table><caption>Table 2 Primer characteristics of cytokine and reference genes for bovine, caprine and ovine gene expression</caption><thead><tr><th>Target gene</th><th>Primer sequence (5'-3')</th><th>Optimal primer concentration (nM)</th><th>Amplicon size (bp)</th></tr></thead><tbody><tr><td rowspan="2">IL-4</td><td>F-CAGCATGGAGCTGCCT</td><td>300</td><td rowspan="2">177</td></tr><tr><td>R-ACAGAACAGGTCTTGCTTGC</td><td>300</td></tr><tr><td rowspan="2">IL-10</td><td>F-CTTTAAGGGTTACCTGGGTTGC</td><td>300</td><td rowspan="2">239</td></tr><tr><td>R-CTCACTCATGGCTTTGTAGACAC</td><td>300</td></tr><tr><td rowspan="2">IL-12B</td><td>F-CAGCAGAGGCTCCTCTGAC</td><td>600</td><td rowspan="2">237</td></tr><tr><td>R-GTCTGGTTTGATGATGTCCCTG</td><td>600</td></tr><tr><td rowspan="2">INF-γ</td><td>F-CAGAGCCAAATTGTCTCCTTTC</td><td>300</td><td rowspan="2">167</td></tr><tr><td>R-ATCCACCGGAATTTGAATCAG</td><td>300</td></tr><tr><td rowspan="2">TNF-α</td><td>F-CCAGAGGGAAGAGCAGTCC</td><td>300</td><td rowspan="2">111</td></tr><tr><td>R-GGCTACAACGTGGGCTACC</td><td>300</td></tr><tr><td rowspan="2">GAPDH</td><td>F-ATCTCGCTCCTGGAAGATG</td><td>600</td><td rowspan="2">227</td></tr><tr><td>R-TCGGAGTGAACGGATTCG</td><td>300</td></tr><tr><td rowspan="2">H3F3A</td><td>F-GAGGTCTCTATACCATGGCTC</td><td>300</td><td rowspan="2">150</td></tr><tr><td>R-GTACCAGGCCTGTAACGATG</td><td>300</td></tr><tr><td rowspan="2">ACTB</td><td>F-TGGGCATGGAATCCTG</td><td>600</td><td rowspan="2">194</td></tr><tr><td>R-GGCGCGATGATCTTGAT</td><td>600</td></tr><tr><td rowspan="2">PPIA</td><td>F-TGACTTCACACGCCATAAT</td><td>300</td><td rowspan="2">180</td></tr><tr><td>R-CTTGCCATCCAACCACTC</td><td>600</td></tr><tr><td rowspan="2">YWHAZ</td><td>F-GAAAGGGATTGTGGACCAG</td><td>300</td><td rowspan="2">183</td></tr><tr><td>R-GGCTTCATCAAATGCTGTCT</td><td>300</td></tr></tbody></table>

Abbreviations: F: forward; R: reverse; bp: base pair; IL, interleukin; IL12B, interleukin p40; TNF-α tumor necrosis factor-alpha; IFN-γ interferon-gamma; ACTB, beta-actin; GAPDH, Glyceraldehyde 3-phosphate dehydrogenase; H3F3A, H3-histone family 3A; PPIA, peptidylprolyl isomerase A; YWHAZ, tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide; Amplicon sizes were determined in *Bos taurus*. Amplicon sizes in *Capra hircus* and *Ovis aries* were similar.

purified with QIAquick PCR Purification Kit (Qiagen Ltd., Crawley, UK), according to the manufacturer's instructions, (ii) quantified using a NanoDrop™ ND-1000 Spectrophotometer (Thermo Fisher Scientific, MA, USA), and (iii) sequenced (Beckman Coulter Genomics, Takeley, UK).

The absence of genomic DNA (gDNA) amplification was checked with bovine, ovine and caprine gDNA isolated from blood buffy coat, using the QIAamp DNA Mini Kit (Qiagen Ltd., Crawley, UK).

## Quantitative PCR assays

Quantitative PCR (QPCR) reactions were conducted on Mx3005P QPCR Systems™ (Agilent Technologies, Santa Clara, USA). Amplifications were performed with Brilliant SYBR® II QPCR Master kit (Agilent Technologies, Santa Clara, USA). Following the manufacturer's instructions, reactions were carried out in a final volume of 25 µl with 1.5 µl of cDNA and 1 µl of each primer at the optimized concentration. Amplifications were performed as described above for conventional PCR assays. A dissociation

<footnote>ª: GenBank accession number.</footnote>
step was included for all reactions to confirm single specific PCR product amplification and define the Tm of each amplicon. A negative control (no-template control) was included in each primer assay to check for the formation of primer-dimers.

## Linearity and efficiency of the qPCR reaction

QPCR reaction efficiency (E) was determined for each purified and quantified PCR product by performing a 10-fold serial dilution in eight points, in duplicate. Calibration curves were plotted; efficiencies and correlation coefficients were calculated by MxPro QPCR Software (Agilent Technologies, Santa Clara, USA).

## Absolute quantification and relative expression genes in Tempus™ and PBMC samples

The RT-qPCR method was checked by amplification, in duplicate, of cDNA synthesized from (i) Tempus™ total RNA from three different animals (absolute quantification), and (ii) PBMC (ConcanavalinA-stimulated and unstimulated) total RNA from three independent experiments (relative expression).

Absolute quantification results are expressed as the mean number of copies (nc) from six amplification values with standard deviation. Nc was calculated using the following equation: $Nc=(Relative amount of target (in ng) \times 10^{-9})/(DNA length(dp) \times 650) \times 6.022 \times 10^{23}$. The relative amount of target (in ng) was calculated using the calibration curve and the following equation: Relative amount of target (in ng)=10 ((number of Cq-intercept)/slope).

Relative expression was calculated with efficiency correction using the relative expression software tool (REST) [30]. Multiple reference gene normalization was chosen.

Gene expression stability of the five candidate reference genes was determined in PBMC samples and calculated using the geNorm application in Microsoft Excel [23]. Stability values (M-value) and pairwise variations (V-score) were calculated and the optimal number of reference genes required for relative expression was determined. In addition, the best stable combination of reference genes was selected with NormFinder Excel Add-In [31]. The relative expression ratio (Er) of cytokine genes in ConcanavalinA-stimulated cells (n=3) compared to unstimulated cells (n=3) was calculated using the quantification of cytokine and reference gene as number of copies. Expression variation for each cytokine gene was provided by REST [30].

## Stability of the five candidate reference genes using geNorm and NormFinder analysis

The gene expression stability of the five candidate reference genes was determined in Tempus™ and PBMC samples and calculated with two different reference gene selection algorithms. The geNorm application in

Microsoft Excel [23] and NormFinder Excel Add-In [31] calculate stability values (M-value in geNorm and p-value in NormFinder) and the reference gene with the lowest stability value is the most stable expression. NormFinder identify the best stable combination of reference genes.

As recommended by the authors of the geNorm method, an M-value <1.5 is the cut-off for suitability as reference genes for PBMC samples (heterogeneous samples) and M-value <0.5 for Tempus™ samples (homogeneous samples). In addition, using the geNorm application, we calculated the pairwise variation (V-score) to determine the optimal number of reference genes required using a V-value <0.25 as a cut-off for suitable reference genes in PBMC samples (heterogeneous samples) and V-value <0.15 in Tempus™ samples (homogeneous samples). All analyses were performed using the quantification of cytokine and reference genes as number of copies. Data were analyzed by NormFinder with groups identified (unstimulated and stimulated group).

## Reproducibility and sensitivity of the RT-qPCR method

Inter-assay variability (R) is the variation in amplification results between different runs for each sample. R was determined by amplification of cDNA synthesized from the Tempus™ total RNA, in duplicate, repeated in three independent runs. R is expressed as Cq and as the mean number of copies (nc) with standard deviation.

The limit of detection (LOD) was the highest to the lowest quantifiable number of copies performed after serial dilutions of the cDNA pool synthesized from the Tempus™ total RNA. LOD is expressed as nc.

## Results

### Primer design and control of primer specificity

Cytokine and reference gene primer sets were designed jointly for the three ruminant species, cattle, sheep and goats, according to the predefined parameters described in "Methods". The characteristics of the primer sets are listed in Table 2. All primer sets allowed specific hybridization and amplification at 60°C. Products were checked by electrophoresis migration and led to a single specific amplicon of the expected size. Also, the absence of amplification with bovine genomic DNA was confirmed.

In addition, similarities between the sequenced PCR products and the corresponding target on the mRNA reference sequences (*Bos taurus*, *Capra hircus*, *Ovis aries*) from the NCBI GenBank database were checked using the standard nucleotide BLAST program [32]. To sum up, all qPCR amplifications led to a single specific peak and confirmed the amplification of a single specific product.
## Linearity and efficiency

The optimal concentration of each primer was determined as the lowest concentration of the reverse and forward PCR primers that resulted in the lowest Cq with no formation of primer-dimer. The absence of no-specific products and primer-dimers was confirmed. Concentrations of primers are listed in Table 2. Single-peak melting curves of the PCR products defined the melting temperatures of PCR products (see Table 3 and Table 4). No amplification was detected in no-RT controls and no-template controls.

All calibration curves produced a linear standard curve and efficiency (E) ranged from 91% (for bovine TNF-α) to 102% (bovine IL-4, ovine and caprine GAPDH, caprine ACTB). All correlation coefficients were higher than 0.998 and confirmed the reaction linearity of all primer sets. PCR efficiency, the slopes and intercepts of calibration curves are listed in Table 3 (cytokine genes) and Table 4 (reference genes). Representative dissociation and calibration curves for the three species are presented in Additional file 1 (cytokine genes) and in Additional file 2 (reference genes).

## Reproducibility and sensitivity

To evaluate the reproducibility of our qPCR assays, we performed amplification of Tempus™ samples (whole

blood samples). All genes were correctly amplified whatever the species and led to single specific product amplification at the expected Tm. Inter-assay variability was calculated using the standard deviation of the six amplification values (Cq and copy values). Variability of cytokine gene amplifications ranged between 0.0 Cq (bovine INF-γ, corresponding to 2 copies) and 0.4 Cq (bovine IL-12B, corresponding to 8 copies). Variability of reference gene amplifications ranged between 0.1 Cq (caprine YWHAZ, corresponding to 33 copies) and 0.4 Cq (bovine ACTB, corresponding to 5.1 10⁴ copies). Inter-assay variability is listed in Table 3 (cytokine genes) and Table 4 (reference genes).

Finally, we evaluated the lowest number of correctly amplified and quantified copies that led to single specific product amplification at the expected Tm (LOD) with our assay. In bovine samples, the lowest quantifiable number of copies was one copy both cytokine genes (IL-12B) and reference genes (GAPDH). In caprine samples, the lowest quantifiable number of copies was two copies for cytokine genes (IL-12B) and one copy for reference genes (GAPDH and ACTB). In ovine samples, the lowest quantifiable number of copies was one copy for cytokine genes (TNF-α) and two copies for reference genes (GAPDH). LOD are listed in Table 3 (cytokine genes) and Table 4 (reference genes).

<table><caption>Table 3 Assay performances for five cytokine genes for bovine, caprine and ovine gene expression</caption><thead><tr><th rowspan="2">Target gene</th><th rowspan="2">Species</th><th rowspan="2">Tm (°C)</th><th rowspan="2">Intercept</th><th rowspan="2">Slope</th><th rowspan="2">E (%)</th><th colspan="4">Inter-assay variability</th><th rowspan="2">LOD nc</th></tr><tr><th>Mean Cq</th><th>SD Cq</th><th>Mean nc</th><th>SD nc</th></tr></thead><tbody><tr><td rowspan="3">IL-4</td><td>cattle</td><td></td><td>2.02</td><td>-3.372</td><td>102</td><td>28.5</td><td>0.2</td><td>25</td><td>3</td><td>9</td></tr><tr><td>goat</td><td>83.9</td><td>1.37</td><td>-3.276</td><td>99</td><td>28.6</td><td>0.2</td><td>30</td><td>4</td><td>11</td></tr><tr><td>sheep</td><td></td><td>1.39</td><td>-3.283</td><td>99</td><td>28.1</td><td>0.1</td><td>38</td><td>2</td><td>5</td></tr><tr><td rowspan="3">IL-10</td><td>cattle</td><td></td><td>0.27</td><td>-3.50</td><td>95</td><td>27.6</td><td>0.2</td><td>104</td><td>16</td><td>6</td></tr><tr><td>goat</td><td>86.6</td><td>0.91</td><td>-3.33</td><td>97</td><td>28.9</td><td>0.2</td><td>30</td><td>4</td><td>7</td></tr><tr><td>sheep</td><td></td><td>1.12</td><td>-3.34</td><td>101</td><td>29.2</td><td>0.2</td><td>17</td><td>2</td><td>5</td></tr><tr><td rowspan="3">IL-12B</td><td>cattle</td><td></td><td>0.8</td><td>-3.32</td><td>97</td><td>27.3</td><td>0.4</td><td>95</td><td>8</td><td>1</td></tr><tr><td>goat</td><td>86.5</td><td>1.70</td><td>-3.28</td><td>95</td><td>35</td><td>0.6</td><td>&lt;1</td><td>-</td><td>2</td></tr><tr><td>sheep</td><td></td><td>1.59</td><td>-3.31</td><td>95</td><td>33</td><td>0.5</td><td>2</td><td>0.6</td><td>5</td></tr><tr><td rowspan="3">INF-γ</td><td>cattle</td><td></td><td>1.27</td><td>-3.449</td><td>95</td><td>27.8</td><td>0.0</td><td>79</td><td>2</td><td>5</td></tr><tr><td>goat</td><td>80.2</td><td>1.63</td><td>-3.309</td><td>95</td><td>28.3</td><td>0.1</td><td>68</td><td>5</td><td>6</td></tr><tr><td>sheep</td><td></td><td>1.02</td><td>-3.409</td><td>96</td><td>27.2</td><td>0.2</td><td>128</td><td>18</td><td>2</td></tr><tr><td rowspan="3">TNF-α</td><td>cattle</td><td></td><td>2.18</td><td>-3.417</td><td>91</td><td>25.0</td><td>0.3</td><td>837</td><td>192</td><td>2</td></tr><tr><td>goat</td><td>84.5</td><td>0.38</td><td>-3.350</td><td>99</td><td>25.1</td><td>0.1</td><td>387</td><td>35</td><td>3</td></tr><tr><td>sheep</td><td></td><td>0.04</td><td>-3.350</td><td>97</td><td>25.7</td><td>0.2</td><td>307</td><td>35</td><td>1</td></tr></tbody></table>

Abbreviations. E, reaction efficiency; Tm, melting temperature; SD: Standard deviation; LOD, limit of detection; nc: number of copies. Abbreviations for cytokine genes see Table 1.

Inter-assay variability (R) is the variation in amplification results between different runs for each sample. R was determined with amplification of cDNA synthesized from Tempus™ whole blood total RNA, in duplicate, repeated in three independent runs. R is expressed in Cq and mean number of copies (nc) with standard deviation (SD).

Limit of detection (LOD) is the lowest number of copies correctly amplified, quantified that led to single specific product amplification at the expected Tm after serial dilutions of cDNA pool synthesized from the Tempus™ total RNA. Nc was calculated using the following equation: Nc=(Relative amount of target (in ng) × 10⁻⁵)/(DNA length (dp) × 650) × 6.022 × 10²³. Relative amount of target (in ng) was calculated using the calibration curve and the following equation: Relative amount of target (in ng)=10 ((number of Cq-intercept)/slope).
<table><caption>Table 4 Assay performances for reference genes for bovine, caprine and ovine gene expression</caption><thead><tr><th rowspan="2">Target gene</th><th rowspan="2">Species</th><th rowspan="2">Tm (°C)</th><th rowspan="2">Intercept</th><th rowspan="2">Slope</th><th rowspan="2">E (%)</th><th colspan="4">Inter-assay variability</th><th rowspan="2">LOD nc</th></tr><tr><th>Mean Cq</th><th>SD Cq</th><th>Mean nc</th><th>SD nc</th></tr></thead><tbody><tr><td rowspan="3">GAPDH</td><td>cattle</td><td></td><td>2.02</td><td>-3.372</td><td>98</td><td>18.3</td><td>0.2</td><td>6.2 10<sup>4</sup></td><td>0.7 10<sup>4</sup></td><td>1</td></tr><tr><td>goat</td><td>84.5</td><td>1.37</td><td>-3.276</td><td>102</td><td>19.6</td><td>0.2</td><td>1.1 10<sup>4</sup></td><td>0.2 10<sup>4</sup></td><td>1</td></tr><tr><td>sheep</td><td></td><td>1.39</td><td>-3.283</td><td>102</td><td>24.9</td><td>0.1</td><td>285</td><td>12</td><td>2</td></tr><tr><td rowspan="3">H3F3A</td><td>cattle</td><td></td><td>0.27</td><td>-3.50</td><td>93</td><td>16.8</td><td>0.1</td><td>11.4 10<sup>4</sup></td><td>0.5 10<sup>4</sup></td><td>20</td></tr><tr><td>goat</td><td>85</td><td>0.91</td><td>-3.33</td><td>100</td><td>16.4</td><td>0.1</td><td>14.2 10<sup>4</sup></td><td>0.7 10<sup>4</sup></td><td>11</td></tr><tr><td>sheep</td><td></td><td>1.12</td><td>-3.34</td><td>99</td><td>16.5</td><td>0.1</td><td>15.4 10<sup>4</sup></td><td>1.1 10<sup>4</sup></td><td>14</td></tr><tr><td rowspan="3">ACTB</td><td>cattle</td><td></td><td>0.8</td><td>-3.32</td><td>100</td><td>15.4</td><td>0.4</td><td>19.4 10<sup>4</sup></td><td>5.1 10<sup>4</sup></td><td>8</td></tr><tr><td>goat</td><td>85.7</td><td>1.70</td><td>-3.28</td><td>102</td><td>15.7</td><td>0.2</td><td>25.6 10<sup>4</sup></td><td>3.8 10<sup>4</sup></td><td>1</td></tr><tr><td>sheep</td><td></td><td>1.59</td><td>-3.31</td><td>101</td><td>15.6</td><td>0.3</td><td>27.8 10<sup>4</sup></td><td>0.5 10<sup>4</sup></td><td>6</td></tr><tr><td rowspan="3">PPIA</td><td>cattle</td><td></td><td>1.27</td><td>-3.449</td><td>95</td><td>17.9</td><td>0.1</td><td>7.7 10<sup>4</sup></td><td>0.5 10<sup>4</sup></td><td>6</td></tr><tr><td>goat</td><td>82.3</td><td>1.63</td><td>-3.309</td><td>100</td><td>17.5</td><td>0.1</td><td>8.3 10<sup>4</sup></td><td>0.6 10<sup>4</sup></td><td>3</td></tr><tr><td>sheep</td><td></td><td>1.02</td><td>-3.409</td><td>96</td><td>17.4</td><td>0.1</td><td>8.2 10<sup>4</sup></td><td>0.5 10<sup>4</sup></td><td>8</td></tr><tr><td rowspan="3">YWHAZ</td><td>cattle</td><td></td><td>2.18</td><td>-3.417</td><td>96</td><td>18.7</td><td>0.1</td><td>7.2 10<sup>4</sup></td><td>0.7 10<sup>4</sup></td><td>9</td></tr><tr><td>goat</td><td>81.3</td><td>0.38</td><td>-3.350</td><td>99</td><td>23.8</td><td>0.1</td><td>506</td><td>33</td><td>4</td></tr><tr><td>sheep</td><td></td><td>0.04</td><td>-3.350</td><td>99</td><td>23.5</td><td>0.1</td><td>506</td><td>24</td><td>3</td></tr></tbody></table>

Abbreviations and notes, see Table 3.

## Validation of the RT-qPCR method: absolute quantification of cytokine and reference gene expression

Absolute quantification was checked by amplification, in duplicate, of cDNA synthesized from Tempus™ total RNA (whole blood samples) from three different animals of each species. All genes were correctly amplified whatever the species, and led to a single specific product amplification at the expected Tm. The mean number of copies and Cq with SD from the six amplification values are listed in Table 5.

In bovine samples, the number of cytokine gene copies detected was between 24 +/-5 copies (IL-4, corresponding to 28.7 +/-0.3 Cq) and 965 +/-199 copies (TNF-α, corresponding to 24.7 +/-0.3 Cq). In ovine and caprine samples, the number of cytokine gene copies detected was similar: between 2 +/-0.7 copies (ovine IL-12B, corresponding 33.3 +/-0.5 Cq) and 798 +/-34 copies (caprine TNF-α, corresponding to 24.0 +/-0.1 Cq).

GAPDH was the most weakly expressed reference gene with 264 +/-38 copies (ovine GAPDH, corresponding to 25.0 +/-0.2 Cq) and PPIA was the most strongly expressed reference gene with 75.6 10⁴ +/-0.8 10⁴ copies (caprine PPIA, corresponding to 17.7 +/-0.2 Cq). H3F3A and ACTB were expressed almost identically in the three species with, for example, respectively 12.5 10⁴ +/-0.7 10⁴ copies in bovine samples (corresponding to 16.7 +/-0.1 Cq) and 30.2 10⁴ +/-4.6 10⁴ copies in ovine samples (corresponding to 15.5 +/-0.2 Cq). Surprisingly, the level of GAPDH and YWHAZ was differently expressed in the three species. GAPDH was weakly expressed in ovine samples with 264 +/-38 copies (corresponding

to 25.0 +/-0.2 Cq) but strongly expressed in bovine and caprine samples with 7.6 10⁴ +/-1.2 10⁴ copies (corresponding to 18.0 +/-0.3 Cq). YWHAZ was weakly expressed in ovine and caprine samples with 482 +/-81 copies (ovine, corresponding to 23.6 +/-0.2 Cq) but strongly expressed in bovine samples with 70.9 10⁴ +/-1.9 10⁴ copies (corresponding to 18.8 +/-0.5 Cq). These differences in expression gene levels are consistent with the levels observed in unstimulated PBMC assays (see Additional file 3).

## Validation of the RT-qPCR method: relative expression of cytokine genes in stimulated PBMC samples

The relative expression of ConcanavalinA-stimulated PBMC samples compared with unstimulated PBMC samples in three independent experiments for each species was calculated using the expression level of cytokine and reference genes. The stability results are presented in Figure 1. Quantification of cytokine and reference genes in unstimulated and ConcanavalinA-stimulated PBMC samples are presented in Additional file 3.

First, the optimal number of reference genes required for relative expression was determined using the geNorm application. Based on the cut-off pairwise variation of 0.25 (V-score), in our study, the use of two reference genes was sufficient for accurate normalization whatever the species, with a V-score of between 0.12 (pairwise variation V2/3 in ovine PBMC samples) and 0.24 (pairwise variation V2/3 in caprine PBMC samples). Second, the best stable combination of two reference genes was identified by NormFinder analysis: PPIA/H3F3A genes for bovine samples
<table><caption>Table 5 Absolute quantification of cytokine and reference gene expression in whole blood samples</caption><thead><tr><th rowspan="2">Target gene</th><th rowspan="2">Species</th><th colspan="2">Cq values</th><th colspan="2">Copie number values</th></tr><tr><th>Mean Cq</th><th>SD Cq</th><th>Mean nc</th><th>SD nc</th></tr></thead><tbody><tr><td rowspan="3">IL-4</td><td>cattle</td><td>28.7</td><td>0.3</td><td>24</td><td>5</td></tr><tr><td>goat</td><td>28.5</td><td>0.2</td><td>31</td><td>3</td></tr><tr><td>sheep</td><td>28.1</td><td>0.5</td><td>39</td><td>12</td></tr><tr><td rowspan="3">IL-10</td><td>cattle</td><td>27.2</td><td>0.6</td><td>156</td><td>54</td></tr><tr><td>goat</td><td>28.2</td><td>0.5</td><td>48</td><td>17</td></tr><tr><td>sheep</td><td>28.5</td><td>0.1</td><td>28</td><td>1</td></tr><tr><td rowspan="3">IL-12B</td><td>cattle</td><td>27.9</td><td>0.6</td><td>62</td><td>21</td></tr><tr><td>goat</td><td>34.9</td><td>1.2</td><td>2</td><td>0.4</td></tr><tr><td>sheep</td><td>33.3</td><td>0.5</td><td>2</td><td>0.7</td></tr><tr><td rowspan="3">INF-γ</td><td>cattle</td><td>27.4</td><td>0.4</td><td>106</td><td>28</td></tr><tr><td>goat</td><td>27.7</td><td>0.3</td><td>103</td><td>21</td></tr><tr><td>sheep</td><td>26.3</td><td>0.3</td><td>241</td><td>52</td></tr><tr><td rowspan="3">TNF-α</td><td>cattle</td><td>24.7</td><td>0.4</td><td>965</td><td>199</td></tr><tr><td>goat</td><td>24.0</td><td>0.1</td><td>978</td><td>34</td></tr><tr><td>sheep</td><td>24.4</td><td>0.1</td><td>728</td><td>34</td></tr><tr><td rowspan="3">GAPDH</td><td>cattle</td><td>18.0</td><td>0.3</td><td>7.6 10<sup>4</sup></td><td>1.2 10<sup>4</sup></td></tr><tr><td>goat</td><td>19.1</td><td>0.1</td><td>1.6 10<sup>4</sup></td><td>0.1 10<sup>4</sup></td></tr><tr><td>sheep</td><td>25.0</td><td>0.2</td><td>264</td><td>38</td></tr><tr><td rowspan="3">H3F3A</td><td>cattle</td><td>16.7</td><td>0.1</td><td>12.5 10<sup>4</sup></td><td>0.7 10<sup>4</sup></td></tr><tr><td>goat</td><td>16.3</td><td>0.2</td><td>14.9 10<sup>4</sup></td><td>2.3 10<sup>4</sup></td></tr><tr><td>sheep</td><td>16.2</td><td>0.1</td><td>19.4 10<sup>4</sup></td><td>1.6 10<sup>4</sup></td></tr><tr><td rowspan="3">ACTB</td><td>cattle</td><td>15.1</td><td>0.3</td><td>24.5 10<sup>4</sup></td><td>4.3 10<sup>4</sup></td></tr><tr><td>goat</td><td>15.4</td><td>0.1</td><td>32.9 10<sup>4</sup></td><td>2.8 10<sup>4</sup></td></tr><tr><td>sheep</td><td>15.5</td><td>0.2</td><td>30.2 10<sup>4</sup></td><td>4.6 10<sup>4</sup></td></tr><tr><td rowspan="3">PPIA</td><td>cattle</td><td>18.4</td><td>0.2</td><td>58.0 10<sup>4</sup></td><td>0.6 10<sup>4</sup></td></tr><tr><td>goat</td><td>17.7</td><td>0.2</td><td>75.6 10<sup>4</sup></td><td>0.8 10<sup>4</sup></td></tr><tr><td>sheep</td><td>17.7</td><td>0.1</td><td>65.8 10<sup>4</sup></td><td>0.4 10<sup>4</sup></td></tr><tr><td rowspan="3">YWHAZ</td><td>cattle</td><td>18.8</td><td>0.5</td><td>70.9 10<sup>4</sup></td><td>1.9 10<sup>4</sup></td></tr><tr><td>goat</td><td>23.5</td><td>0.3</td><td>647</td><td>120</td></tr><tr><td>sheep</td><td>23.6</td><td>0.2</td><td>482</td><td>81</td></tr></tbody></table>

Abbreviations for cytokine and reference genes, see Table 1.

Absolute quantification was checked by amplification of cDNA synthesized from Tempus™ total RNA from three different animals of each species.

Absolute quantification used the standard curve method and is expressed as the number of copies (nc) and Cq with standard deviation (SD) from six amplification values.

(with a stability value of 0.024), GAPDH/H3F3A genes for caprine samples (with a stability value of 0.072) and PPIA/H3F3A genes for ovine samples (with a stability value of 0.064). Each combination of reference genes was used to calculate the relative expression of cytokine genes.

Finally, a qPCR assay enabled the relative expression of cytokine to be determined in all samples, whatever the species, even in cases of strong up-regulation of cytokine mRNA expression (i.e. caprine INF-γ) or in

cases of decreased expression of cytokine mRNA (i.e. bovine IL-10). Cytokine expression results are presented in the whisker-box plots provided by REST [30] (see Figure 2).

## Stability of the five reference genes and comparison between geNorm and NormFinder analyses and samples

The expression of the five selected reference genes (GAPDH, ACTB, YWHAZ, PPIA and H3F3A) was determined in all samples, in both PBMC and Tempus™ (whole blood) samples, using two statistical approaches, geNorm application [23] and NormFinder software [31]. Figure 1 shows the ranked lists of genes and pairwise variations for the PBMC samples and Figure 3, for the whole blood samples.

In the PBMC samples, *M*-values were lower than 1.5 and between 0.455 (bovine ACTB) and 0.999 (ovine H3F3A), and *V*-values were lower than 0.25 and between 0.12 (caprine V2/3) and 0.24 (ovine V2/3). In whole blood samples, *M*-values were lower than 0.5 and between 0.152 (caprine PPIA) and 0.374 (ovine H3F3A), and *V*-values were lower than 0.15 and between 0.039 (caprine V2/3) and 0.062 (ovine V2/3).

In bovine samples, the most stable reference gene in both whole blood and PBMC samples was ACTB, whatever the analysis. Also, the ranked list of reference genes, provided by both software, was the same for PBMC and for whole blood samples.

In ovine samples, stability results were very similar. For PBMC samples, the ranked lists of reference genes, produced by geNorm and NormFinder were the same and the most stable reference gene was ACTB. In whole blood samples, the three most stable genes were GAPDH, ACTB and PPIA, with very similar stability values (0.209 < *M*-value < 0.220).

In caprine samples, the ranked lists of reference genes for samples and analyses differed. In PBMC samples, ACTB and YWHAZ were the two most stable reference genes according to both geNorm and NormFinder analyses. In the whole blood samples, PPIA was the most stable reference gene and the three most stable genes were the same, with very close stability values (0.152 < *M*-value < 0.176).

## Discussion

We developed a novel SYBR Green RT-qPCR assay for the simultaneous relative quantification of five major cytokines (IL-4, IL-10, IL-12B, IFN-γ and TNF-α) in cattle, sheep and goats, and accurate normalization with five candidate reference genes (GAPDH, ACTB, H3F3A, PPIA and YWHAZ). The main advantage of this original work is the design of an optimal primer set for three dominant species with the same hybridization temperature. To our knowledge, this is the first SYBR Green assay
Stability values

<table><thead><tr><th rowspan="2">gene name</th><th colspan="3">NormFinder</th></tr><tr><th>M-value</th><th>gene name</th><th>p-value</th></tr></thead><tbody><tr><td>ACTB</td><td>0.455</td><td>ACTB</td><td>0.128</td></tr><tr><td>YWHAZ</td><td>0.463</td><td>YWHAZ</td><td>0.154</td></tr><tr><td>GAPDH</td><td>0.518</td><td>GAPDH</td><td>0.222</td></tr><tr><td>PPIA</td><td>0.706</td><td>PPIA</td><td>0.378</td></tr><tr><td>H3F3A</td><td>0.721</td><td>H3F3A</td><td>0.387</td></tr></tbody></table>

Pairwise Variations

A

<figure><img src="image_2.png" /></figure>

Stability values

<table><thead><tr><th rowspan="2">gene name</th><th colspan="3">NormFinder</th></tr><tr><th>M-value</th><th>gene name</th><th>p-value</th></tr></thead><tbody><tr><td>ACTB</td><td>0.646</td><td>ACTB</td><td>0.094</td></tr><tr><td>YWHAZ</td><td>0.747</td><td>YWHAZ</td><td>0.335</td></tr><tr><td>PPIA</td><td>0.807</td><td>PPIA</td><td>0.371</td></tr><tr><td>GAPDH</td><td>0.967</td><td>GAPDH</td><td>0.544</td></tr><tr><td>H3F3A</td><td>0.999</td><td>H3F3A</td><td>0.570</td></tr></tbody></table>

Pairwise Variations

B

<figure><img src="image_3.png" /></figure>

Stability values

<table><thead><tr><th rowspan="2">gene name</th><th colspan="3">NormFinder</th></tr><tr><th>M-value</th><th>gene name</th><th>p-value</th></tr></thead><tbody><tr><td>ACTB</td><td>0.559</td><td>ACTB</td><td>0.202</td></tr><tr><td>YWHAZ</td><td>0.572</td><td>YWHAZ</td><td>0.265</td></tr><tr><td>H3F3A</td><td>0.688</td><td>GAPDH</td><td>0.343</td></tr><tr><td>GAPDH</td><td>0.695</td><td>H3F3A</td><td>0.368</td></tr><tr><td>PPIA</td><td>0.814</td><td>PPIA</td><td>0.433</td></tr></tbody></table>

Pairwise Variations

C

<figure><img src="image_4.png" /></figure>

**Figure 1** Stability of five candidate reference genes in PBMC samples using geNorm and NormFinder analysis. Abbreviations. PBMC, peripheral blood mononuclear cells. Abbreviations for reference genes, see Table 1. PBMCs were cultured with medium (unstimulated condition) or stimulated with 5 µg/ml of Concanavalin A (stimulated condition) for 36 h. Stability of the five reference genes were calculated in bovine (A), caprine (B) and ovine (C) PBMC samples using geNorm application and NormFinder software. Stability values were calculated as M-value in geNorm and p-value in NormFinder. The reference genes were presented in ranked list of the most stable gene to the least stable gene. The pairwise variation (V-score) was performed with geNorm application to determine the optimal number of required reference genes for relative expression in bovine (A), caprine (B) and ovine (C) PBMC samples.

that enables quantification of bovine, ovine and caprine cytokines with a single set of reagents using the same protocol. Due to the lack of tools to investigate host immunity by quantification of secreted cytokines, analysis of cytokine gene expression has become a widely used method to establish a cytokine profile in ruminants. Many infectious diseases (for example Johne's disease) affect many animal species and threaten livestock especially in the case of multi-species grazing. Our assay could thus greatly improve research by scientists involved in ruminant diseases, particularly to characterize the immune response to pathogens affecting cattle, sheep, and goats.

RT-qPCR is one of the most widely used methods to investigate cytokine gene expression in ruminants thanks to its high specificity, sensitivity and accuracy. It is also the most suitable for expression analysis of miscellaneous

cytokines from a single sample and a small quantity of template material. However, to ensure accurate reproducible quantitative data, strict standard operating procedures have to be followed and many reviews have demonstrated the importance of the RT-qPCR workflow [21,33,34].

For that reason, we developed the novel SYBR Green RT-qPCR assay presented here taking particular care to maximize the assay, detailing each step of the experimental protocol, from primer design, sample preparation or reverse transcription to the choice and normalization by reference genes, and reporting on the performance of the assays.

We first focused our attention on the criteria needed for primer design, particularly primer location straddling exons, and the choice of a single hybridization
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

**Figure 2** Cytokine relative expression in bovine, caprine and ovine stimulated PBMC samples. Abbreviations. PBMC, peripheral blood mononuclear cells. Abbreviations for reference genes, see Table 1. PBMCs were cultured with medium (unstimulated condition) or stimulated with 5 µg/ml of Concanavalin A (stimulated condition) for 36 h. The relative expression ratio (Er) of cytokine genes in Concanavalin A-stimulated cells compared to unstimulated cells was calculated (amplification of three independent experiments in duplicate) in bovine (A), caprine (B) and ovine (C) samples using the relative expression software tool (REST) with the most stable combination of two reference genes identified with NormFinder software. REST Software uses randomization technique and error bars represents distribution of permutated expression data. Black boxes represent up-regulated cytokine expression and white boxes down-regulated cytokine expression.
Stability values

<table><thead><tr><th colspan="2">geNorm</th><th colspan="2">NormFinder</th></tr><tr><th>gene name</th><th><i>M</i>-value</th><th>gene name</th><th><i>p</i>-value</th></tr></thead><tbody><tr><td>ACTB</td><td>0.202</td><td>ACTB</td><td>0.043</td></tr><tr><td>GAPDH</td><td>0.209</td><td>GAPDH</td><td>0.043</td></tr><tr><td>H3F3A</td><td>0.251</td><td>H3F3A</td><td>0.130</td></tr><tr><td>PPIA</td><td>0.270</td><td>PPIA</td><td>0.161</td></tr><tr><td>YWHAZ</td><td>0.350</td><td>YWHAZ</td><td>0.234</td></tr></tbody></table>

Pairwise Variations

A

<figure><img src="image_8.png" /></figure>

Stability values

<table><thead><tr><th colspan="2">geNorm</th><th colspan="2">NormFinder</th></tr><tr><th>gene name</th><th><i>M</i>-value</th><th>gene name</th><th><i>p</i>-value</th></tr></thead><tbody><tr><td>PPIA</td><td>0.152</td><td>PPIA</td><td>0.021</td></tr><tr><td>H3F3A</td><td>0.169</td><td>ACTB</td><td>0.021</td></tr><tr><td>ACTB</td><td>0.176</td><td>H3F3A</td><td>0.062</td></tr><tr><td>YWHAZ</td><td>0.213</td><td>YWHAZ</td><td>0.129</td></tr><tr><td>GAPDH</td><td>0.350</td><td>GAPDH</td><td>0.237</td></tr></tbody></table>

Pairwise Variations

B

<figure><img src="image_9.png" /></figure>

Stability values

<table><thead><tr><th colspan="2">geNorm</th><th colspan="2">NormFinder</th></tr><tr><th>gene name</th><th><i>M</i>-value</th><th>gene name</th><th><i>p</i>-value</th></tr></thead><tbody><tr><td>GAPDH</td><td>0.209</td><td>PPIA</td><td>0.061</td></tr><tr><td>ACTB</td><td>0.219</td><td>GAPDH</td><td>0.079</td></tr><tr><td>PPIA</td><td>0.220</td><td>ACTB</td><td>0.101</td></tr><tr><td>YWHAZ</td><td>0.298</td><td>YWHAZ</td><td>0.164</td></tr><tr><td>H3F3A</td><td>0.374</td><td>H3F3A</td><td>0.2445</td></tr></tbody></table>

Pairwise Variations

C

<figure><img src="image_10.png" /><figcaption>**Figure 3** Stability of five candidate reference genes in whole blood samples using geNorm and NormFinder analysis. Abbreviations for reference genes, see Table 1. Whole blood (or Tempus™) samples were performed from three different animals of each species. Stability of the five reference genes were calculated in bovine (A), caprine (B) and ovine (C) whole blood samples using geNorm application and NormFinder software. Stability values were calculated as *M*-value in geNorm and *p*-value in NormFinder. The reference genes were presented in ranked list of the most stable gene to the least stable gene. The pairwise variation (V-score) was performed with geNorm application to determine the optimal number of required reference genes for relative expression in bovine (A), caprine (B) and ovine (C) whole blood samples.</figcaption></figure>

temperature. S. Taylor [33] underlines the importance of primer design and the choice of target sequence to ensure specific efficient amplification of the products. Indeed, a design that spans exon-exon boundaries of at least one of the two primers, in addition to the DNAse treatment, prevents amplification of the target from contaminating gDNA. In our study, all primers were carried out fulfilling our criteria-with the exception of five primers that could not be designed to span exon-exon boundaries-to achieve a GC% of between 40% and 70% and a melting temperature of between 62°C and 65°C. This was the case for forward primers of H3F3A, PPIA, IL-4, and INF-γ genes, and for the reverse primer of the IL-10 gene. In these cases, the

second primer was then designed to span exon-exon boundaries. Thus, absence of amplification in no-RT controls and with genomic DNA, reaction linearity higher than 0.998 and high PCR efficiencies of all primer sets helped design sound primers and ensured that contaminating gDNA amplification did not occur.

Using the same hybridization temperature made it possible to analyze many samples in the same run, to compare the level of expression of a particular gene between different samples and to perform simultaneous normalization. With the common primer design, the same hybridization temperature had the additional advantage of enabling different samples of three ruminant species to be analyzed in the same run. This reduced the time and cost
of the analysis, and also improved reproducibility, as demonstrated by our low run-to-run variation.

We then carefully selected five reference genes, the most widely used (GAPDH and ACTB) and/or those reported to be the most stable, for accurate normalization in both different types of samples and in the three target species. Accurate normalization is essential to eliminate sample-to-sample variations and to obtain reliable results, particularly in relative gene expression studies [21] for which multiple normalization genes with three to five reference genes could be required [23].

Before selecting suitable reference genes for normalization, the stability of the candidate reference genes (GAPDH, ACTB, H3F3A, PPIA and YWHAZ) has to be confirmed and the optimal number of reference genes required for an accurate relative expression has to be determined. In fact, there is no consensus on which method should be used to examine reference gene expression stability. We chose to evaluate the gene expression stability of the five reference genes using geNorm application [23] and NormFinder software [31] as two complementary statistical approaches. We showed that the most stable genes were similar or very close with the geNorm application and NormFinder software. In whole blood and PBMC assays, stability analyses allowed confirmation of the stability of the five candidate reference genes, which two were sufficient for an accurate normalization of our assays. While many studies recommend multiple normalization genes with three to five reference genes [33], we performed new REST analyses and compared the results of cytokine expression using the best combination genes selected by NormFinder software and the two or three most stable genes selected by geNorm application. Our results showed that gene expression was similar whatever the combination of genes selected (see results in Additional file 4), suggesting that a panel of five reference genes could be sufficient and underlining the importance of selecting adequate reference genes and checking their stability. In addition, differences and similarities in the expression stability of the reference genes were compiled. In our ovine PBMC assay, ACTB, YWHAZ and PPIA were the three most stable genes. These results were confirmed in our whole blood assay, where ACTB and PPIA were two of the three most stable reference genes. Several studies reported YWHAZ to be the most stable reference gene in ovine samples [27,35-37]. In addition to the above mentioned results, in our study, ACTB was also the most stable reference genes in bovine PBMC and whole blood samples, followed by YWHAZ and GAPDH. These results are partially in agree with Sheridan's results using whole blood samples and with Spalenza and Robinson's results using PBMC samples [24-26]. Our results emphasize the importance of following strict standard

operating procedures and maximizing the assay to prevent technical variations. Finally, our results demonstrate the advantages of our RT-qPCR workflow, especially in the case of limited sample availability and funds.

In addition, we revealed an original difference in the expression level of GAPDH and YWHAZ genes among the three species, in both unstimulated PBMC and in whole blood samples. GAPDH and YWHAZ are involved in two different biological processes and have sophisticated functions in organisms in which their expression profile might fluctuate. GAPDH is an important enzyme for energy metabolism while YWHAZ plays a role in protein folding. The absence of any correlation between them in the variations we observed could exclude a GAPDH/YWHAZ co-regulation.

Our results emphasize the need to select different reference genes for accurate normalization, i.e. the most stable gene, or better yet, the most stable gene combination.

We demonstrated the sensitivity of the method and the correct cytokine gene expression quantification. In the case of low level mRNA, up-regulation, or a decrease in expression genes, an efficient RT-qPCR must be able to assess the cytokine profile regardless of the status of the animal or of the experimental conditions. Our assay ensured the detection of low levels of cytokine mRNA up to one or two copies. This was the case, for example, for IL-12B expression. We observed a low level of IL-12B expression in bovine and ovine unstimulated PBMC samples, as previously reported by Weiss in bovine macrophages and by Budhia in ovine PBMC [38]. Also, despite the low levels of IL-12B expression, 3 +/−1 copies and 8 +/−2 copies were obtained in bovine and ovine un-stimulated PBMC samples respectively, up-regulation was determined relative to bovine and ovine stimulated PBMC samples by REST.

Relative expression results showed that IL-10 is up-regulated in ovine stimulated PBMC and not or only weakly down-regulated in bovine and caprine stimulated PBMC. We are unable to explain this difference in expression genes between species but we can compare with IL-12 expression, which was more strongly up-regulated in bovine and caprine stimulated PBMC.

Our SYBR Green RT-qPCR assay and the method we used to develop it, from primer design, through sample preparation and including validation, could be used to develop new cytokine primer sets and to increase the panel of cytokines and reference genes. These reference genes will also be useful for any gene expression studies whatever the species of ruminant.

## Conclusions

We developed a novel SYBR Green RT-qPCR assay for the simultaneous relative quantification of five major cytokines (IL-4, IL-10, IL-12B, IFN-γ and TNF-α) in
cattle, sheep and goats, and their accurate normalization with five candidate reference genes. This novel SYBR Green RT-qPCR assay, which we validated in whole blood and PBMC samples, carries out accurate and reproducible quantitative data and ensures the detection of low levels of cytokine mRNA. This original assay is a robust, efficient, cost-effective, widely accessible tool to identify immunological markers for infections in three species of ruminants, cattle, sheep and goats.

## Additional files

**Additional file 1: Representative dissociation and calibration curves of qPCR for reference genes.** Dissociation curves were plotted by MxPro QPCR Software. Single-peak melting curves defined the melting temperatures of PCR products and confirmed the amplification of a single specific product. QPCR reaction efficiency was determined for each purified and quantified PCR product by performing a 10-fold serial dilution in eight points, in duplicate. Calibration curves were plotted by MxPro QPCR Software.

**Additional file 2: Representative dissociation and calibration curves of qPCR for cytokine genes.** Dissociation curves were plotted by MxPro QPCR Software. Single-peak melting curves defined the melting temperatures of PCR products and confirmed the amplification of a single specific product. QPCR reaction efficiency was determined for each purified and quantified PCR product by performing a 10-fold serial dilution in eight points, in duplicate. Calibration curves were plotted by MxPro QPCR Software.

**Additional file 3: Quantification of cytokine and reference gene expression in unstimulated/stimulated PBMC assays.** Abbreviations for cytokine and reference genes, see Table 1. Quantification of cytokine and reference gene expression was checked by amplification of cDNA synthesized from Concanavalin A-stimulated and unstimulated PBMC total RNA from three independent experiments for each species. Quantification of gene expression used the standard curve method and is expressed as the number of copies (nc) and Cq with standard deviation (SD) from three amplification values.

**Additional file 4: Cytokine gene expression using three combinations of reference genes selected by NormFinder and geNorm analyses.** Abbreviations for cytokine and reference genes, see Table 1. The relative expression ratio of cytokine genes in Concanavalin A-stimulated cells compared to unstimulated cells was calculated (amplification of three independent experiments in duplicate) in bovine (A), caprine (B) and ovine (C) samples using the relative expression software tool (REST) using three different combinations of reference genes selected by geNorm application an NormFinder software.

## Abbreviations

RT-qPCR: Reverse transcriptase quantitative polymerase chain reaction;
IL: Interleukin; IL-12B: Interleukin p40; TNF-α: Tumor necrosis factor-alpha;
IFN-γ: Interferon-gamma; ACTB: Beta-actin; GAPDH: Glyceraldehyde 3-phosphate dehydrogenase; H3F3A: H3-histone family 3A; PPIA: Peptidylprolyl Isomerase A;
YWHAZ: Tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide; MIQE: Minimum Information for publication of Quantitative Real-time PCR Experiments; E: Reaction efficiency; PBMC: Peripheral blood mononuclear cells; REST: Relative expression software tool; Er: Relative expression ratio.

## Competing interests

The authors declare that they have no competing interests.

## Authors' contributions

CP conceived the study, designed the experiments, performed all the experiments, analyzed the data and drafted the manuscript. LD participated in writing the manuscript. IC participated in developing the proposal. VR conceived the study and corrected the manuscript. All the authors have read and approved the final manuscript.

## Acknowledgements

The authors would like to thank Virginie Dupuy for helpful discussions and her critical review of this paper. We thank also Noël Richard and Thierry Maccota for their technical assistance in animal care and sampling.

## Author details

¹INRA, UMR1309 CMAEE, Montpellier F-34398, France. ²CIRAD, DGD-RS-Dist, Montpellier F-34398, France. ³CIRAD, UMR Interyrp, Montpellier F-34398, France. ⁴CIRAD, UMR CMAEE, Montpellier F-34398, France.

Received: 2 October 2014 Accepted: 27 February 2015

Published online: 17 March 2015

## References

1. OIE. 2012. Animal health: a multifaceted challenge. [http://www.oie.int/fileadmin/Home/eng/Media_Center/docs/pdf/Key_Documents/ANIMAL-HEALTH-EN-FINAL.pdf]

2. Lefrançois T, Pineau T. Public health and livestock: emerging diseases in food animals. Animal Frontiers. 2014;4:4–6.

3. Hackmann TJ, Spain JN. Invited review: ruminant ecology and evolution: perspectives useful to ruminant livestock research and production. J Dairy Sci. 2010;93:1320–34.

4. Grace D, Mutua F, Ochungo P, Kruska R, Jones K, Brierley L, et al. Mapping of Poverty and Likely Zoonoses Hotspots. 2012.

5. Mosmann TR, Coffman RL. TH1 and TH2 cells: different patterns of lymphokine secretion lead to different functional properties. Annu Rev Immunol. 1989;7:145–73.

6. Boddu-Jasmine HC, Witchell J, Vordermeier M, Wangoo A, Goyal M. Cytokine mRNA expression in cattle infected with different dosages of Mycobacterium bovis. Tuberculosis. 2008;88:610–5.

7. Coussens PM, Verman N, Coussens MA, Elfman MD, McNulty AM. Cytokine gene expression in peripheral blood mononuclear cells and tissues of cattle infected with mycobacterium avium subsp. paratuberculosis: evidence for an inherent proinflammatory gene expression pattern. Infect Immun. 2004;72:1409–22.

8. Singh PK, Singh SV, Saxena VK, Singh MK, Singh AV, Sohal JS. Expression profiles of different cytokine genes in peripheral blood mononuclear cells of goats infected experimentally with native strain of Mycobacterium avium subsp. paratuberculosis. Anim Biotechnol. 2013;24:187–97.

9. Smeed JA, Watkins CA, Rhind SM, Hopkins J. Differential cytokine gene expression profiles in the three pathological forms of sheep paratuberculosis. BMC Vet Res. 2007;3:18.

10. Channappanavar R, Singh KP, Singh R, Umeshappa CS, Ingale SL, Pandey AB. Enhanced proinflammatory cytokine activity during experimental bluetongue virus-1 infection in Indian native sheep. Vet Immunol Immunopathol. 2012;145:485–92.

11. Kwong LS, Thom M, Sopp P, Rocchi M, Wattegedera S, Entrican G, et al. Production and characterization of two monoclonal antibodies to bovine tumour necrosis factor alpha (TNF-alpha) and their cross-reactivity with ovine TNF-alpha. Vet Immunol Immunopathol. 2010;135:320–4.

12. Kwong LS, Hope JC, Thom ML, Sopp P, Duggan S, Bembridge GP, et al. Development of an ELISA for bovine IL-10. Vet Immunol Immunopathol. 2002;85:213–23.

13. Hope JC, Kwong LS, Entrican G, Wattegedera S, Vordermeier HM, Sopp P, et al. Development of detection methods for ruminant interleukin (IL)-12. J Immunol Methods. 2002;266:117–26.

14. Hope JC, Kwong LS, Thom M, Sopp P, Mwangi W, Brown WC, et al. Development of detection methods for ruminant interleukin (IL)-4. J Immunol Methods. 2005;301:114–23.

15. Hope JC, Sopp P, Wattegedera S, Entrican G. Tools and reagents for caprine immunology. Small Rumin Res. 2012;103:23–7.

16. Entrican G, Lunney JK, Rutten VP, Baldwin CL. A current perspective on availability of tools, resources and networks for veterinary immunology. Vet Immunol Immunopathol. 2009;128:24–9.

17. Entrican G. New technologies for studying immune regulation in ruminants. Vet Immunol Immunopathol. 2002;87:485–90.

18. Giulietti A, Overbergh L, Valcx D, Decallonne B, Bouillon R, Mathieu C. An overview of real-time quantitative PCR: applications to quantify cytokine gene expression. Methods. 2001;25:386–401.

19. Bustin SA. Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. J Mol Endocrinol. 2000;25:169–93.
20. Bustin SA, Nolan T. Pitfalls of quantitative real-time reverse-transcription polymerase chain reaction. J Biomol Tech. 2004;15:155-66.

21. Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55:611-22.

22. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2=ΔΔCT method. Methods. 2001;25:402-8.

23. Vandesompele J, Preter KD, Pattyn F, Poppe B, Roy NV, Paepe AD, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3: research0034.

24. Robinson TL, Sutherland IA, Sutherland J. Validation of candidate bovine reference genes for use with real-time PCR. Vet Immunol Immunopathol. 2007;115:160-5.

25. Sheridan MP, Browne JA, MacHugh DE, Costello E, Gormley E. Impact of delayed processing of bovine peripheral blood on differential gene expression. Vet Immunol Immunopathol. 2012;145:199-205.

26. Spalenza V, Girolami F, Bevilacqua C, Riondato F, Rasero R, Nebbia C, et al. Identification of internal control genes for quantitative expression analysis by real-time PCR in bovine peripheral lymphocytes. Vet J. 2011;189:278-83.

27. Peletto S, Bertuzzi S, Campanella C, Modesto P, Maniaci MG, Bellino C, et al. Evaluation of internal reference genes for quantitative expression analysis by real-time PCR in ovine whole blood. Int J Mol Sci. 2011;12:7732-47.

28. Taraktsoglou M, Szalabska U, Magee DA, Browne JA, Sweeney T, Gormley E, et al. Transcriptional profiling of immune genes in bovine monocyte-derived macrophages exposed to bacterial antigens. Vet Immunol Immunopathol. 2011;140:130-9.

29. Meade KG, Gormley E, Doyle MB, Fitzsimons T, O'Farrelly C, Costello E, et al. Innate gene repression associated with Mycobacterium bovis infection in cattle: toward a gene signature of disease. BMC Genomics. 2007;8:400.

30. Pfaffl MW, Horgan GW, Dempfle L. Relative expression software tool (REST®) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 2002;30:e36.

31. Andersen CL, Jensen JL, Ørntoft TF. Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004;64:5245-50.

32. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403-10.

33. Taylor S, Wakem M, Dijkman G, Alsarraj M, Nguyen M. A practical approach to RT-qPCR-publishing data that conform to the MIQE guidelines. Methods. 2010;50:S1-5.

34. Derveaux S, Vandesompele J, Hellemans J. How to do successful gene expression analysis using real-time PCR. Methods. 2010;50:227-30 [The Ongoing Evolution of qPCR].

35. Jarczak J, Kaba J, Bagnicka E. The validation of housekeeping genes as a reference in quantitative Real Time PCR analysis: application in the milk somatic cells and frozen whole blood of goats infected with caprine arthritis encephalitis virus. Gene. 2014;549:280-5.

36. Vorachek WR, Hugejiletu, Bobe G, Hall JA. Reference gene selection for quantitative PCR studies in sheep neutrophils. Int J Mol Sci. 2013;14:11484-95.

37. Zang R, Bai J, Xu H, Zhang L, Yang J, Yang L, et al. Selection of suitable reference genes for real-time quantitative PCR studies in Lanzhou Fat-tailed Sheep (Ovis aries). Asian J Animal Vet Adv. 2011;6:789-804.

38. Budhia S, Haring LF, McConnell I, Blacklaws BA. Quantitation of ovine cytokine mRNA by real-time RT-PCR. J Immunol Methods. 2006;309:160-72.

Submit your next manuscript to BioMed Central
and take full advantage of:

• Convenient online submission

• Thorough peer review

• No space constraints or color figure charges

• Immediate publication on acceptance

• Inclusion in PubMed, CAS, Scopus and Google Scholar

• Research which is freely available for redistribution

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central