Evaluation of internal reference genes for quantitative expression analysis by real-time PCR in ovine whole blood

A novel quantitative real-time PCR with the GAPDH reference gene for peste des petits ruminants

Peste des petits ruminants (PPR) is a serious acute, highly contagious disease caused by the peste des petits ruminants virus (PPRV). This study aims to establish a qRT-PCR assay with an internal amplification control for the rapid and accurate detection of PPRV. The primers and probes for PPRV N were based on the national standard of the diagnostic techniques for PPR of China, and a pair of primers and TaqMan probes for the internal reference gene of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was designed. Optimisation of the reaction conditions, specificity, sensitivity and reproducibility tests, and clinical sample detection were conducted. The results showed that the optimal primers and probe concentrations of PPRV were 0.4 μmol/l and 0.4 μmol/l, respectively, and were 0.4 μmol/l and 0.2 μmol/l for the reference gene GAPDH, respectively. The established method has no cross-reaction with other viruses. The minimum detection limit was 6.8 copies/μl for PPRV and 190 copies/μl for GAPDH. The coefficients of variation (CV%) of PPRV and GAPDH were both lower than 2%. The results suggest that the PPRV qRT-PCR method containing internal reference genes has strong specificity, high sensitivity, and good reproducibility. The addition of internal reference genes for the sample quality control improves the accuracy of the detection.

Evaluation of reference genes for gene expression analysis by real-time quantitative PCR (qPCR) in different tissues from mice infected by Ascaris suum

Human ascariasis is the most prevalent helminth infection, affecting 445 million people worldwide. To better understand the impact of the immune system on the pathophysiology of individuals infected with Ascaris suum, mice have been used as experimental models. The RT-qPCR technique is a critical auxiliary tool of investigation used to quantify mRNA levels. However, proper normalization using reference genes is essential to ensure reliable outcomes to avoid analytical errors and false results. Despite the importance of reference genes for experimental A. suum infection studies, no specific reference genes have been identified yet. Therefore, we conducted a study to assess five potential reference genes (GAPDH, 18s, ACTB, B2M, and HPRT1) in different tissues (liver, lungs, small and large intestines) affected by A. suum larval migration in C57BL/6j mice. Tissue collection was carried out to analyze parasite burden and confirm the presence of larvae during the peak of migration in each tissue. Upon confirmation, we analyzed different genes in the tissues and found no common gene with stable expression. Our results highlight the importance of analyzing different genes and using different software programs to ensure reliable relative expression results. Based on our findings, B2M was ranked as the ideal reference gene for the liver, while 18S was the most stable gene in the lung and small intestine. ACTB, or a combination of ACTB with GAPDH, was deemed suitable as reference genes for the large intestine due to their stable expression and less variation between the control and infected groups. To further demonstrate the impact of using different reference genes, we normalized the expression of a chemokine gene (CXCL9) in all tissues. Significant differences in CXCL9 expression levels were observed between different groups in all tissues except for the large intestine. This underscores the importance of selecting appropriate reference genes to avoid overestimating target gene expression levels and encountering normalization-related issues that can lead to false results. In conclusion, our study highlights the significance of using reliable reference genes for accurate RT-qPCR analysis, especially in the context of A. suum infection studies in different tissues. Proper normalization is crucial to ensure the validity of gene expression data and avoid potential pitfalls in interpreting results.

Identification and validation of stable reference genes for expression profiling of target genes in diverse ovine tissues

Quantitative PCR (qPCR) is a widely-used technique for quantifying the expression of target genes across various tissues, as well as under different pathological and physiological conditions. One of the challenges associated with this method is the need to identify optimal reference genes (RGs) that maintain consistent expression levels under diverse experimental settings, thereby ensuring accurate biological interpretation. In this study, we conducted a thorough analysis of 18 candidate RGs (ACTB, BACH1, B2M, GAPDH, HMBS, HPRT1, PGK1, PPIA, PPIB, RPLP0, RPL19, RPS9, RPS15, RPS28, SDHA, TBP, UXT, and YWHAZ) across 10 ovine tissues (muscle, skin, kidney, liver, intestine, rumen, lung, testis, heart, and spleen) obtained from five individual sheep. We aimed to identify genes with stable expression across these tissues. A literature-based survey helped us shortlist candidate genes representing various functional classes from multiple livestock species. We employed four algorithms: geNorm, NormFinder, BestKeeper, and Delta Ct (ΔCt), to rank these genes based on their stability. A consistent trend in the rankings was observed across these different algorithms. RefFinder was then used for a comprehensive ranking, integrating the outputs from the various methods. ACTB, PPIB, BACH1, and B2M emerged as the most stable RGs, while RPS9, RPS15, and PGK1 displayed variable expression. We validated our findings through qPCR analysis of four target genes (ACTN2, CRYAB, DLK1, and TRIM54) in the skin samples from two different sheep breeds. Based on these results, we recommend ACTB, PPIB, BACH1, and B2M as reliable internal control genes for qPCR experiments involving diverse ovine tissues.

Analysis of the stability of the reference genes GAPDH, SDHA and RPL-19 in sheep from a semi-arid region infected by gastrointestinal nematodes

Analyzing the stability of reference genes already described as universal is an important methodology to lead gene expression analysis because different studies have shown that the expression of universal reference genes may vary between experimental treatments. In this sense, the glyceraldehyde 3-phosphate dehydrogenase (GAPDH), Succinate dehydrogenase complex subunit A (SDHA) and Ribosomal Protein L-19 (RPL-19) reference genes (already described in other studies with sheep from different regions, breeds and infectious agents or in organisms evolutionarily close to sheep) were investigated in the abomasum, small and large intestines of resistant and susceptible crossbred sheep groups to gastrointestinal nematode infections in the Semi-arid region in Northeast of Brazil. The animals were naturally infected to determine the resistance or susceptibility status by counting eggs per gram (EPG) of feces from the gastrointestinal tract after 33 weeks of observations of infection evolution. Relative gene expression was performed by RT-qPCR methodology using Sybr green and relative gene expression stability was tested by different software programs such as REST, BestKeeper, geNorm and Normfinder. Our results showed the susceptible animals had increase in egg counts per gram of feces than resistant animals (p < 0.001), and both groups showed a mixed infection by nematodes of the genus Haemonchus, Trichostrongylus, Oesophagostomum and Trichuris. Furthermore, we show the importance of analyzing different genes in different software programs and the importance to choose ideal reference genes. In this sense, GAPDH was the most stable gene in the abomasum, whereas SDHA was the most stable in the small and large intestines. In addition, we discuss about variables which can interfere in relative expression such as breed, species, climate and tissue. However, utilizing other reference genes already described in other studies with the same and different variables should be performed.

Evaluation of the innate immune response of caprine neutrophils against Mycobacterium avium subspecies paratuberculosis in vitro

Neutrophils constitute an essential component of the innate immune response, readily killing most bacteria through phagocytosis, degranulation, and the release of neutrophil extracellular traps (NETs) among other mechanisms. These cells play an unclear role in mycobacterial infections such as Mycobacterium avium subspecies paratuberculosis (Map), the etiological agent of paratuberculosis, and its response is particularly understudied in ruminants. Herein, a wide set of techniques were adapted, or newly developed, to study the in vitro response of caprine neutrophils after Map infection. Immunofluorescence was used to demonstrate, simultaneously, chemotaxis, phagocytosis, degranulation, and NETs. The quantification of neutrophil phagocytic activity against Map at a 1:10 multiplicity of infection (MOI), through flow cytometry, showed values that varied from 4.54 to 5.63% of phagocyting neutrophils. By immunofluorescence, a 73.3 ± 14.5% of the fields showed NETs, and the mean release of DNA, attributable to NETosis, calculated through a fluorometric method, was 16.2 ± 3.5%. In addition, the RNA expression of TGF-β, TNF and IL-1β cytokines, measured through reverse transcription qPCR, was significantly higher in the two latter. Overall, neutrophil response was proportional to the number of bacteria. This work confirms that the simultaneous study of several neutrophil mechanisms, and the combination of different methodologies, are essential to reach a comprehensive understanding of neutrophil response against pathogens, demonstrates that, in vitro, caprine neutrophils display a strong innate response against Map, using their entire repertoire of effector functions, and sets the basis for further in vitro and in vivo studies on the role of neutrophils in paratuberculosis.

Screening and stability analysis of reference genes in fasting caecotrophy model in rabbits

BACKGROUND

The selection and validation of stably expressed reference genes is key for accurately quantifying the mRNA abundance of genes under different treatments. In the rabbit model of fasting caecotrophy, reports about the selection of stable reference genes are not available.

METHODS AND RESULTS

This study aims to screen suitable reference genes in different tissues (including uterus, cecum, and liver) of rabbits between control and fasting caecotrophy groups. RT-qPCR was used to analyze the expression levels of eight commonly used reference genes (including GAPDH, 18S rRNA, B2M, CYP, HPRT1, β-actin, H2afz, Ywhaz), and RefFinder (including geNorm, NormFinder, and BestKeeper) was used to analyze the expression stability of these reference genes. Our results showed that the most stable reference genes were different in different tissues and treatments. In the control and fasting caecotrophy groups, CYP, GAPDH and HPRT1 were proven to be the top stable reference genes in the uterus, cecum, and liver tissues, respectively. GAPDH and Ywhaz were proven to be the top two stable reference genes among uterus, cecum, and liver in both control and fasting caecotrophy groups.

CONCLUSIONS

Our results indicated that the combined analysis of three or more reference genes (GAPDH, HPRT1, and Ywhaz) are recommended to be used for RT-qPCR normalization in the rabbit model of fasting caecotrophy, and that GAPDH is a better choice than the other reference genes for normalizing the relative expression of target genes in different tissues of fasting caecotrophy rabbits.

Messenger RNA expression and localization of xenin in the gastrointestinal tract in sheep

The present study aimed to determine the primary sequence of ovine xenin and clarify the mRNA expression and peptide localization of xenin in the gastrointestinal tract in sheep. The colocalization of xenin and glucose-dependent insulinotropic polypeptide was also compared in the antrum and duodenum. Analysis of the nucleotide sequence of ovine xenin revealed a high degree (97.9%) of sequence homology of the sequence between sheep and cattle, and the amino acids sequence determined for ovine xenin coincided (100%) with that of other mammalian species. Real-time quantitative PCR for ovine xenin did not show regional difference in the mRNA expression ratio of xenin. In contrast to the real-time quantitative PCR results, anti-xenin positive cells were abundantly localized in the abomasal antrum (P < 0.01) and at a lesser amount in the duodenum, but no antixenin positive cells were observed in the other regions. Anti-xenin single-positive cells were in a majority in the abomasal antrum, whereas anti-xenin single-positive cells, and anti-GIP single-positive cells, and double-positive cells were even colocalized in the duodenum. These results suggest that abomasal antrum is a major source of xenin in the ovine gastrointestinal tract.

Development and evaluation of a real-time PCR panel for the detection of 20 immune markers in cattle and sheep

The establishment of a panel of immune markers is of paramount importance to understand the different transcription patterns of infectious diseases in livestock. The array of commercially available immunological assays for cattle and sheep is currently limited, due to the lack of antibodies for these species. Even though SYBR Green based real time quantitative PCR (qPCR) is the most commonly used method to study cytokine transcription in ruminants, a lack of standardization impairs its implementation in the study of different ruminant diseases. In order to obtain reliable qPCR results, several variables need to be considered: choice of reference genes for optimal normalization, variation of annealing temperature among primer sets, and assay specificity and sensitivity. In this study, we developed and validated a panel of immune markers in bovine and ovine samples using SYBR Green based qPCR in a cost-effective way with multiple primer sets optimised to amplify at a common thermal cycling temperature. Twenty primer sets were designed to quantify immune markers (IL-1b, IL-2, IL-4, IL-5, IL-6, IL-10, IL-12, IL-13, IL-15, IL-18, IL-23, TNF-α, IFN-γ, IFN-α, Ki-67, NFkB-65, TLR-3, TLR-4, TLR-8 and Rig-1) in ovine and bovine templates. For optimal normalization and selection of suitable reference genes, primer sets that measure the transcription of five reference genes were also included in the panel. The amplification efficiency, linearity and specificity was validated for all target genes. Optimal amplification conditions were achieved in both ovine and bovine samples for all gene targets, with the exception of Ki67. Relative quantification studies were performed on ovine and bovine mRNA obtained from sheep peripheral blood mononuclear cells (PBMCs) stimulated with three different treatments (PMA/Ionomycin, Concanavalin A (Con A) and pokeweed mitogen (PWM)). Pokeweed and ConA efficiently induced gene transcription of most of the targeted genes, while PMA/Ionomycin showed a weaker induction. Finally, we further assessed usability of our panel by running it on bovine monocyte derived dendritic cells (MoDCs) stimulated with different vaccines. Results confirmed the induction of a specific pro-inflammatory gene transcription pattern by rabies vaccine, which resembles the one occurring during viral infection. Altogether, we validated the efficiency and usability of an extended real-time PCR panel that gives the possibility to rapidly measure a broad spectrum of ovine and bovine immune markers by using a single set of reagents and protocol thus representing a valid and cost-effective tool for research purposes.

Characterization of MDCK cells and evaluation of their ability to respond to infectious and non-infectious stressors

The Madin-Darby Canine Kidney (MDCK) cell line is widely used as epithelial cell model in studies ranging from viral infection to environmental pollutants, and vaccines production. However, little is known about basal expression of genes involved in innate immunity, and the ability to respond to infectious and non-infectious stressors. Therefore, the aims of our study were to evaluate the basal level of expression of pivotal genes in the innate immune response and cell cycle regulation, as well as to evaluate the ability of this cell line to respond to infectious or non-infectious stressors. As surmised in our working hypothesis, we demonstrated the constitutive expression of genes involved in the innate immune response and cell defense alike, including TLRs, Interleukins, Myd88, p65/NF-kB and p53. Moreover, we described the ability of this cell line to respond to LPS and cadmium (Cd2+) in terms of gene expression and cytokine release. These data confirm the possibility of using this cell line as a model in studies of host/pathogen interaction and response to non-infectious stressors.

A Transcriptomic Study of the Tail Fat Deposition in Two Types of Hulun Buir Sheep According to Tail Size and Sex

Simple Summary

Based on tail types, Hulun Buir sheep were divided into two lines including small and big fat-tailed, but these two lines have similar genetic background. In this study, we investigated the morphology and transcription level differences of tail fat between these two lines. The RNA-seq analyses indicated several differentially expressed genes when compared between sexes or two tail sizes. Interestingly, we also found an obvious sex difference in the fat metabolism in Hulun Buir sheep. Two different co-expression networks were only shown either in male or in female sheep. Our findings will provide theoretical background in understanding the genetic mechanism of fat deposition in sheep.

Abstract

Hulun Buir sheep of similar genetic background were divided into two lines based on tail types: Small- and big fat-tailed. To explore the molecular mechanism of fat deposition in sheep tails, we firstly evaluated the morphology and transcription level differences of tail fat between these two lines. RNA-Seq technology was used to identify differentially expressed genes (DEGs) in phenotypic extremes of tail sizes. Five comparisons were performed taking into account two factors, sex and tail type. We screened out 373 DEGs between big-tailed and small-tailed Hulun Buir sheep, and 775 and 578 DEGs between two types of tails in male and female sheep, respectively. The results showed an obvious sex difference in the fat metabolism in sheep based on gene ontology (GO), pathway, and network analyses. Intriguingly, there were two different co-expression networks only respectively shown in male and female sheep, which were insulin-related network acting on upstream pathways and PPARG-related network effect in downstream pathways. Furthermore, these two networks were linked by a classic pathway of regulating adipogenesis. This is the first study to investigate the sex differences of fat metabolism in domestic animals, and it demonstrates a new experimental way to study fat metabolism. Our findings will provide theoretical background in understanding the tail-size phenotype in sheep and can be exploited in breeding small-tailed sheep.

Whole-Genome Resequencing Reveals Loci Associated With Thoracic Vertebrae Number in Sheep

The number of vertebrae, especially thoracic vertebrae, is an important economic trait that may influence carcass length and meat production in animals. However, the genetic basis of vertebrae number in sheep is still poorly understood. To detect the candidate genes, 400 increased number of thoracic vertebrae (T14L6) and 200 normal (T13L6) Kazakh sheep were collected. We generated and sequenced 60 pools of genomic DNA (each pool prepared by mixing genomic DNA from 10 sheep with the same thoracic traits), with an average depth of coverage of 25.65×. We identified a total of 42,075,402 SNPs and 11 putatively selected genomic regions, including the VRTN gene and the HoxA gene family that regulate vertebral development. The most prominent areas of selective elimination were located in a region of chromosome 7, including VRTN, which regulates spinal development and morphology. Further investigation indicated that the expression level of the VRTN gene during fetal development was significantly higher in sheep with more thoracic vertebrae than in those with a normal number of thoracic vertebrae. A genome-wide comparison between sheep with increased and normal numbers of thoracic vertebrae showed that the VRTN gene is the major selection locus for the number of thoracic vertebrae in sheep and has the potential to be utilized in sheep breeding in the future.

Does kisspeptin participate in GABA-mediated modulation of GnRH and GnRH receptor biosynthesis in the hypothalamic-pituitary unit of follicular-phase ewes?

BACKGROUND

The inverse relationship between GnRH transcript level and GABA neurons activity has suggested that GABA at the hypothalamic level may exert a suppressive effect on subsequent steps of the GnRH biosynthesis. In the present study, we analyzed the effects of GABA type A receptor agonist (muscimol) or antagonist (bicuculline) on molecular mechanisms governing GnRH/LH secretion in follicular-phase sheep.

METHODS

ELISA technique was used to investigate the effects of muscimol and/or bicuculline on levels of post-translational products of genes encoding GnRH ligand and GnRH receptor (GnRHR) in the preoptic area (POA), anterior (AH) and ventromedial (VMH) hypothalamus, stalk/median eminence (SME), and GnRHR in the anterior pituitary (AP). Real-time PCR was chosen for determination of the effect of drugs on kisspeptin (Kiss 1) mRNA level in POA and VMH including arcuate nucleus (VMH/ARC), and on Kiss1 receptor (Kiss1r) mRNA abundance in POA-hypothalamic structures. These analyses were supplemented by RIA method for measurement of plasma LH concentration.

RESULTS

The study demonstrated that muscimol and bicuculline significantly decreased or increased GnRH biosynthesis in all analyzed structures, respectively, and led to analogous changes in plasma LH concentration. Similar muscimol- and bicuculline-related alterations were observed in levels of GnRHR. However, the expression of Kiss 1 and Kiss1r mRNAs in selected POA-hypothalamic areas of either muscimol- and bicuculline-treated animals remained unaltered.

CONCLUSIONS

Our data suggest that GABAergic neurotransmission is involved in the regulatory pathways of GnRH/GnRHR biosynthesis and then GnRH/LH release in follicular-phase sheep conceivably via indirect mechanisms that exclude involvement of Kiss 1 neurons.

Selection of reference genes for expression analysis in mouse models of acute alcoholic liver injury

Investigations of hepatic gene expression are crucial for determining the molecular factors involved in acute alcoholic liver injury. The results of liver molecular investigations may reveal etiologically important genomic alterations. Therefore, it is necessary to normalize gene expression data to identify stable genes, which may be used as a reference under different experimental conditions. The aim of the present study was to apply reverse transcription‑quantitative polymerase chain reaction analysis and use analysis software to investigate the expression stability of candidate reference genes in hepatic tissues from mice with acute alcoholic liver injury. The acute alcoholic liver injury models were established by the intragastric administration of alcohol (5 mg/kg) in Imprinting Control Region mice. Total RNA was isolated from the mouse livers, following which the expression levels of seven reference genes, β-actin, glyceraldehyde 3-phosphate dehydrogenase (Gadph), glucuronidase β, hypoxanthine phosphoribosyltransferase 1 (Hprt1), 18S ribosomal RNA, TATA binding protein and β‑2 microglobulin, were examined, and gene expression stability was assessed using the geNorm, NormFinder and BestKeeper tools. The geNorm analysis revealed that the gene with the lowest variability was Hprt1. Hprt1 and Gapdh were validated as the optimal reference gene pair in all samples from all groups. The NormFinder and BestKeeper results showed that Hprt1 was the most stable gene in all samples. Alcohol induces endoplasmic reticulum (ER) stress, causing changes in the expression levels of ER stress‑associated genes. The stability of Hprt1 was verified by the expression analysis of ER stress‑associated genes, and gene expression levels in the ethanol groups were upregulated, with a significant difference in expression, compared with those in the control group. Therefore, Hprt1 was selected as the most stable gene, and Hprt1 and Gapdh were determined to be the optimum gene pair in mouse models of acute alcoholic liver injury. The reliability of the Hprt1 gene was confirmed by expression analysis of ER stress‑associated genes.

The influence of dopaminergic system inhibition on biosynthesis of gonadotrophin-releasing hormone (GnRH) and GnRH receptor in anoestrous sheep; hierarchical role of kisspeptin and RFamide-related peptide-3 (RFRP-3)

This study aimed to explain how prolonged inhibition of central dopaminergic activity affects the cellular processes governing gonadotrophin-releasing hormone (GnRH) and LH secretion in anoestrous sheep. For this purpose, the study included two experimental approaches: first, we investigated the effect of infusion of sulpiride, a dopaminergic D2 receptor antagonist (D2R), on GnRH and GnRH receptor (GnRHR) biosynthesis in the hypothalamus and on GnRHR in the anterior pituitary using an immunoassay. This analysis was supplemented by analysis of plasma LH levels by radioimmunoassay. Second, we used real-time polymerase chain reaction to analyse the influence of sulpiride on the levels of kisspeptin (Kiss1) mRNA in the preoptic area and ventromedial hypothalamus including arcuate nucleus (VMH/ARC), and RFamide-related peptide-3 (RFRP-3) mRNA in the paraventricular nucleus (PVN) and dorsomedial hypothalamic nucleus. Sulpiride significantly increased plasma LH concentration and the levels of GnRH and GnRHR in the hypothalamic–pituitary unit. The abolition of dopaminergic activity resulted in a significant increase in transcript level of Kiss1 in VMH/ARC and a decrease of RFRP-3 in PVN. The study demonstrates that dopaminergic neurotransmission through D2R is involved in the regulatory pathways of GnRH and GnRHR biosynthesis in the hypothalamic–pituitary unit of anoestrous sheep, conceivably via mechanisms in which Kiss1 and RFRP-3 participate.

A tissue-based approach to selection of reference genes for quantitative real-time PCR in a sheep osteoporosis model

BACKGROUND

In order to better understand the multifactorial nature of osteoporosis, animal models are utilized and compared to healthy controls. Female sheep are well established as a model for osteoporosis induced by ovariectomy, calcium and vitamin D low diet, application of steroids, or a combination of these treatments. Transcriptional studies can be performed by applying quantitative real time PCR (RT-qPCR). RT-qPCR estimates mRNA-levels of target genes in relation to reference genes. A chosen set of reference genes should not show variation under experimental conditions. Currently, no standard reference genes are accepted for all tissue types and experimental conditions. Studies examining reference genes for sheep are rare and only one study described stable reference in mandibular bone. However, this type of bone differs from trabecular bone where most osteoporotic fractures occur. The present study aimed at identifying a set of reference genes for relative quantification of transcriptional activity of ovine spine bone and ovine in vitro differentiated mesenchymal stromal cells (MSC) for reliable comparability.

METHODS

Twelve candidate reference genes belonging to different functional classes were selected and their expression was measured from cultured ovMSCs (n = 18) and ovine bone samples (n = 16), respectively. RefFinder was used to rank the candidate genes.

RESULTS

We identified B2M, GAPDH, RPL19 and YWHAZ as the best combination of reference genes for normalization of RT-qPCR results for transcriptional analyses of these ovine samples.

CONCLUSION

This study demonstrates the importance of applying a set of reference genes for RT-qPCR analysis in sheep. Based on our data we recommend using four identified reference genes for relative quantification of gene expression studies in ovine bone or for in vitro experiments with osteogenically differentiated ovine MSCs.

Genome-wide analysis of circular RNAs in prenatal and postnatal pituitary glands of sheep

Circular RNAs (circRNAs) are a class of animal non-coding RNAs and play an impor-tant role in animal growth and development. However, the expression and function of circRNAs in the pituitary gland of sheep are unclear. Transcriptome profiling of circRNAs in the pituitary gland of sheep may enable us to understand their biological functions. In the present study, we identified 10,226 circRNAs from RNA-seq data in the pituitary gland of prenatal and postnatal sheep. Reverse transcription PCR and DNA sequencing analysis confirmed the presence of several circRNAs. Real-time RT-PCR analysis showed that sheep circRNAs are resistant to RNase R digestion and are expressed in prenatal and postnatal pituitary glands. GO and KEGG enrichment analysis showed that host genes of differentially expressed circRNAs are involved in the regulation of hormone secretion as well as in several pathways related to these processes. We determined that numerous circRNAs interact with pituitary-specific miRNAs that are involved in the biologic functions of the pituitary gland. Moreover, several circRNAs contain at least one IRES element and open reading frame, indicating their potential to encode proteins. Our study provides comprehensive expression profiles of circRNAs in the pituitary gland, thereby offering a valuable resource for circRNA biology in sheep.

Genome-wide analysis of circular RNAs in prenatal and postnatal muscle of sheep

Circular RNAs (circRNAs), a type of non-coding RNA with circular structure, were generated by back splicing and widely expressed in animals and plants. Recent studies have shown that circRNAs extensively participate in cell proliferation, cell differentiation, cell autophagy and other biological processes. However, the role and expression of circRNAs in the development and growth of muscle have not been studied in sheep. In our study, we first used RNA-seq to study the circRNAs in prenatal and postnatal longissimus dorsi muscle of sheep. A total of 6113 circRNAs were detected from the RNA-seq data. Several circRNAs were identified using reverse transcription PCR, DNA sequencing and RNase R digestion experiments. The expression levels of several circRNAs in prenatal and postnatal muscle were confirmed by Real-Time RT-PCR. The gene ontology (GO) and KEGG enrichment analysis of the host gene of the circRNAs showed that these circRNAs were mainly involved in the growth and development of muscle related signaling pathways. These circRNAs might sponge microRNAs (miRNAs) in predicted circRNA-miRNA-mRNA networks. The circRNAs expression profiles in muscle provided an important reference for the study of circRNAs in sheep.

Assessment of reference genes for reliable analysis of gene transcription by RT-qPCR in ovine leukocytes

With the availability of genetic sequencing data, quantitative reverse transcription PCR (RT-qPCR) is increasingly being used for the quantification of gene transcription across species. Too often there is little regard to the selection of reference genes and the impact that a poor choice has on data interpretation. Indeed, RT-qPCR provides a snapshot of relative gene transcription at a given time-point, and hence is highly dependent on the stability of the transcription of the reference gene(s). Using ovine efferent lymph cells and peripheral blood mono-nuclear cells (PBMCs), the two most frequently used leukocytes in immunological studies, we have compared the stability of transcription of the most commonly used ovine reference genes: YWHAZ, RPL-13A, PGK1, B2M, GAPDH, HPRT, SDHA and ACTB. Using established algorithms for reference gene normalization "geNorm" and "Norm Finder", PGK1, GAPDH and YWHAZ were deemed the most stably transcribed genes for efferent leukocytes and PGK1, YWHAZ and SDHA were optimal in PBMCs. These genes should therefore be considered for accurate and reproducible RT-qPCR data analysis of gene transcription in sheep.

Evaluation of reference gene suitability for quantitative expression analysis by quantitative polymerase chain reaction in the mandibular condyle of sheep

Reference genes are commonly used as a reliable approach to normalize the results of quantitative polymerase chain reaction (qPCR), and to reduce errors in the relative quantification of gene expression. Suitable reference genes belonging to numerous functional classes have been identified for various types of species and tissue. However, little is currently known regarding the most suitable reference genes for bone, specifically for the sheep mandibular condyle. Sheep are important for the study of human bone diseases, particularly for temporomandibular diseases. The present study aimed to identify a set of reference genes suitable for the normalization of qPCR data from the mandibular condyle of sheep. A total of 12 reference genes belonging to various functional classes were selected, and the expression stability of the reference genes was determined in both the normal and fractured area of the sheep mandibular condyle. RefFinder, which integrates the following currently available computational algorithms: geNorm, NormFinder, BestKeeper, and the comparative ΔCt method, was used to compare and rank the candidate reference genes. The results obtained from the four methods demonstrated a similar trend: RPL19, ACTB, and PGK1 were the most stably expressed reference genes in the sheep mandibular condyle. As determined by RefFinder comprehensive analysis, the results of the present study suggested that RPL19 is the most suitable reference gene for studies associated with the sheep mandibular condyle. In addition, ACTB and PGK1 may be considered suitable alternatives.

Selection of reference genes for gene expression studies related to intramuscular fat deposition in Capra hircus skeletal muscle

The identification of suitable reference genes is critical for obtaining reliable results from gene expression studies using quantitative real-time PCR (qPCR) because the expression of reference genes may vary considerably under different experimental conditions. In most cases, however, commonly used reference genes are employed in data normalization without proper validation, which may lead to incorrect data interpretation. Here, we aim to select a set of optimal reference genes for the accurate normalization of gene expression associated with intramuscular fat (IMF) deposition during development. In the present study, eight reference genes (PPIB, HMBS, RPLP0, B2M, YWHAZ, 18S, GAPDH and ACTB) were evaluated by three different algorithms (geNorm, NormFinder and BestKeeper) in two types of muscle tissues (longissimus dorsi muscle and biceps femoris muscle) across different developmental stages. All three algorithms gave similar results. PPIB and HMBS were identified as the most stable reference genes, while the commonly used reference genes 18S and GAPDH were the most variably expressed, with expression varying dramatically across different developmental stages. Furthermore, to reveal the crucial role of appropriate reference genes in obtaining a reliable result, analysis of PPARG expression was performed by normalization to the most and the least stable reference genes. The relative expression levels of PPARG normalized to the most stable reference genes greatly differed from those normalized to the least stable one. Therefore, evaluation of reference genes must be performed for a given experimental condition before the reference genes are used. PPIB and HMBS are the optimal reference genes for analysis of gene expression associated with IMF deposition in skeletal muscle during development.

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

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12917-015-0382-0) contains supplementary material, which is available to authorized users.

Identification of stably expressed reference genes for RT-qPCR data normalization in defined localizations of cyclic bovine ovaries

Ovaries are highly complex organs displaying morphological, molecular and functional differences between their cortical zona parenchymatosa and medullary zona vasculosa, and also between the different cyclic luteal stages. Objective of the present study was to validate expression stability of twelve putative reference genes (RGs) in bovine ovaries, considering the intrinsic heterogeneity of bovine ovarian tissue with regard to different luteal stages and intra-ovarian localizations. The focus was on identifying RGs, which are suitable to normalize RT-qPCR results of ovaries collected from clinical healthy cattle, irrespective of localization and the hormonal stage. Expression profiles of twelve potential reference genes (GAPDH, ACTB, YWHAZ, HPRT1, SDHA, UBA52, POLR2C, RPS9, ACTG2, H3F3B, RPS18 and RPL19) were analysed. Evaluation of gene expression differences was performed using genorm, normfinder, and bestkeeper software. The most stably expressed genes according to genorm, normfinder and bestkeeper approaches contained the candidates H3F3B, RPS9, YWHAZ, RPS18, POLR2C and UBA52. Of this group, the genes YWHAZ, H3F3B and RPS9 could be recommended as best-suited RGs for normalization purposes on healthy bovine ovaries irrespective of the luteal stage or intra-ovarian localization.

Evaluation of internal reference genes for quantitative expression analysis by real-time reverse transcription-PCR in somatic cells from goat milk

Reverse transcription (RT) quantitative real-time PCR (qPCR) is the most accurate and easy-to-perform technique to measure the expression level of a selected gene of interest by quantifying mRNA transcripts. The use of reference genes is commonly accepted as the most reliable approach to normalize RT-qPCR data and reduce possible errors generated in the quantification of gene expression. The optimal number and choice of reference genes are experimentally validated for specific tissues or cell types and experimental designs. To date, data on qPCR normalization in goats are scarce and the most suitable reference genes in this species have been identified for only a limited number of tissues. The aim of this study was to determine an optimal combination of stably expressed reference genes in caprine milk somatic cells (MSC) from healthy and infected mammary glands. For the purpose, we performed RT-qPCR for 10 commonly used reference genes from various functional classes and then determined their expression level in MSC from goats intramammary challenged with Staphylococcus aureus and in MSC from healthy controls, with a view to select genes whose stability would be unaffected under infection conditions. The geNorm and NormFinder algorithms were used for validating the reference genes. Furthermore, to demonstrate the importance of normalization of gene expression with appropriate reference genes, we tested the effect of using a combination of the least stable genes for expression analysis evaluation. On the basis of our evaluation, we recommend the use of a panel of reference genes that should include G6PD, YWHAZ, and ACTB for caprine MSC gene expression profiling. The expression of the 2 genes of interest, pentraxin-related protein (PTX3) and secreted phosphoprotein 1 (SPP1), was evaluated by RT-qPCR in all samples collected pre- and postinfection, and the recommended reference genes were used to normalize the data. Our study provides a validated panel of optimal reference genes for the identification of genes differentially expressed by qRT-PCR in caprine MSC. Moreover, we provided a set of intron-spanning primer sequences that could be suitable for gene expression experiments using SYBR Green chemistry on other caprine tissues and cells.

Reference gene selection for quantitative PCR studies in sheep neutrophils

Reference genes are essential for studying mRNA expression with quantitative PCR (qPCR). We investigated 11 potential neutrophil reference genes (RPL19, GAPDH, ACTB, B2M, HPRT, G6PD, TFRC, PGK1, YWHAZ, SDHA and GYPC) for sheep under disease conditions of foot rot (FR) and with or without Se supplementation. Initial screening was based on gene expression level (<28 Cq cycles) and variability (SD < 1.5 Cq cycles) and excluded TFRC, GYPC and HPRT from further analysis. Expression stability of the remaining genes was evaluated using four software programs: geNorm, NormFinder, BestKeeper and the comparative delta Cq method. The neutrophil reference genes, G6PD, YWHAZ, GAPDH, RPL19 and SDHA, consistently ranked among the top five most stable genes under these experimental conditions. The SDHA gene expression was not stable in FR-diseased sheep receiving Se treatment and, thus, cannot be recommended as a reference gene. The commonly used genes, PGK1, ACTB and B2M, were not reliable reference genes, underscoring the need to validate neutrophil reference genes under different experimental conditions. Multiple references genes rather than a single gene may provide more robust and reliable results. The best pair of reference genes was SDHA/G6PD in healthy sheep and GADPH/YWHAZ in FR-diseased sheep.

Selection of suitable reference genes for normalization of quantitative real-time PCR in cartilage tissue injury and repair in rabbits

When studying the altered expression of genes associated with cartilage regeneration by quantitative real-time RT-PCR (RT-qPCR), reference genes with highly stable expression during different stages of chondrocyte developmental are necessary to normalize gene expression accurately. Until now, no reports evaluating expression changes of commonly used reference genes in rabbit articular cartilage have been published. In this study, defects were made in rabbit articular cartilage, with or without insulin-like growth factor 1 (IGF-1) treatment, to create different chondrocyte living environments. The stability and intensity of the expressions of the candidate reference genes glyceraldehyde-3-phosphate dehydrogenase (GAPDH), 18S Ribosomal RNA (18S rRNA), cyclophilin (CYP), hypoxanthine phosphoribosyl transferase (HPRT1), and beta-2-microglobulin (B2M) were evaluated. The data were analyzed by geNorm and NormFinder. B2M and 18S rRNA were identified to be suitable reference genes for rabbit cartilage tissues.