Acessibilidade / Reportar erro

Reference gene selection for quantitative PCR in liver, skeletal muscle, and jejunum of Bos indicus cattle

ABSTRACT

The objective of the present study was to evaluate the stability of candidate reference genes and select the genes that can be used for normalizing real-time polymerase chain reaction (PCR) in the liver, skeletal muscle, and jejunum tissues of Nellore or Nellore × Angus steers fed different diets. Fourteen purebred and 14 crossbred steers were used, in which half of the animals of each genetic group received a diet containing whole shelled corn (WSC) and the other half whole shelled corn and sugarcane bagasse (WSCB). Stability was analyzed by the RefFinder program. To validate the selection of candidate reference genes, the expression of target genes was evaluated using the different groups of reference genes. The most stable genes were 18S, ACTB, and CASC3 for skeletal muscle; HMBS, ACTB, and 18S for the liver; and GAPDH, ACTB, and CASC3 for the jejunum, regardless of breed and diet provided. Possible errors caused in data analyses were clarified comparing the more and less stable genes as reference for normalization of the target genes FASN, ACOX, SCD1, MGAM, and SLC2A1. The use of the more stable and less stable sets of reference genes may lead to different conclusions in respect to the expression profile of the target studied gene. The selection of more suitable reference genes for each experiment is of utmost importance to ensure the reliability of gene expression studies so that they can be applied in practice.

bovine; housekeeping genes; normalization; RT-qPCR; tissues

1. Introduction

Elucidation of various molecular mechanisms, previously unknown, occurred after discovery of the polymerase chain reaction (PCR) method in the scientific community, which has become prominent in numerous research studies involving animal production (Oliveira et al., 2014Oliveira, D. M.; Chalfun-Junior, A.; Chizzotti, M. L.; Barreto, H. G.; Coelho, T. C.; Paiva, L. V.; Coelho, C. P.; Teixeira, P. D.; Schoonmaker, J. P. and Ladeira, M. M. 2014. Expression of genes involved in lipid metabolism in the muscle of beef cattle fed soybean or rumen-protected fat, with or without monensin supplementation. Journal of Animal Science 92:5426-5436. https://doi.org/10.2527/jas.2014-7855
https://doi.org/10.2527/jas.2014-7855...
; Teixeira et al., 2017Teixeira, P. D.; Oliveira, D. M.; Chizzotti, M. L.; Chalfun-Junior, A.; Coelho, T. C.; Gionbelli, M. P.; Paiva, L. V.; Carvalho, J. R. R. and Ladeira, M. M. 2017. Subspecies and diet affect the expression of genes involved in lipid metabolism and chemical composition of muscle in beef cattle. Meat Science 133:110-118. https://doi.org/10.1016/j.meatsci.2017.06.009
https://doi.org/10.1016/j.meatsci.2017.0...
; Lopes et al., 2020Lopes, R. C.; Sampaio, C. B.; Trece, A. S.; Teixeira, P. D.; Gionbelli, T. R. S.; Santos, L. R.; Costa, T. C.; Duarte, M. S. and Gionbelli, M. P. 2020. Impacts of protein supplementation during late gestation of beef cows on maternal skeletal muscle and liver tissues metabolism. Animal 14:1867-1875. https://doi.org/10.1017/S1751731120000336
https://doi.org/10.1017/S175173112000033...
).

The relative quantification of reverse transcription PCR (RT-qPCR) is an accurate and sensitive method to quantify mRNA levels of target genes. However, data processing can affect interpretation of results (Lisowski et al., 2008Lisowski, P.; Pierzchała, M.; Gościk, J.; Pareek, C. S. and Zwierzchowski, L. 2008. Evaluation of reference genes for studies of gene expression in the bovine liver, kidney, pituitary, and thyroid. Journal of Applied Genetics 49:367-372. https://doi.org/10.1007/bf03195635
https://doi.org/10.1007/bf03195635...
), being necessary a normalization to correct inevitable experimental variations. The most common method of data normalization is the reference genes, which generally have a relatively stable expression pattern in tissues, since they are required for basic cellular processes (Chapman and Waldenström, 2015Chapman, J. R. and Waldenström, J. 2015. With reference to reference genes: a systematic review of endogenous controls in gene expression studies. PLoS One 10:e0141853. https://doi.org/10.1371/journal.pone.0141853
https://doi.org/10.1371/journal.pone.014...
). However, it is becoming increasingly apparent that the stability of some usual reference genes can change according to species, tissue types, and experimental conditions (Lisowski et al., 2008Lisowski, P.; Pierzchała, M.; Gościk, J.; Pareek, C. S. and Zwierzchowski, L. 2008. Evaluation of reference genes for studies of gene expression in the bovine liver, kidney, pituitary, and thyroid. Journal of Applied Genetics 49:367-372. https://doi.org/10.1007/bf03195635
https://doi.org/10.1007/bf03195635...
; Kaur et al., 2018Kaur, R.; Sodhi, M.; Sharma, A.; Sharma, V. L.; Verma, P.; Swami, S. K.; Kumari, P. and Mukesh, M. 2018. Selection of suitable reference genes for normalization of quantitative RT-PCR (RT-qPCR) expression data across twelve tissues of riverine buffaloes (Bubalus bubalis). PLoS One 13:e0191558. https://doi.org/10.1371/journal.pone.0191558
https://doi.org/10.1371/journal.pone.019...
; Lozano-Villegas et al., 2021Lozano-Villegas, K. J.; Rodríguez-Hernández, R.; Herrera-Sánchez, M. P.; Uribe-García, H. F.; Naranjo-Gómez, J. S.; Otero-Arroyo, R. J. and Rondón-Barragán, I. S. 2021. Identification of reference genes for expression studies in the whole-blood from three cattle breeds under two states of livestock weather safety. Animals 11:3073. https://doi.org/10.3390/ani11113073
https://doi.org/10.3390/ani11113073...
).

Therefore, selecting reference genes with less variation has become a requirement for RT-qPCR analysis, and an important factor for obtaining reliable results. If gene expression is normalized with a reference gene that fluctuate under experimental conditions, relevant changes in expression of the target genes can easily be missed or overemphasized (Chapman and Waldenström, 2015Chapman, J. R. and Waldenström, J. 2015. With reference to reference genes: a systematic review of endogenous controls in gene expression studies. PLoS One 10:e0141853. https://doi.org/10.1371/journal.pone.0141853
https://doi.org/10.1371/journal.pone.014...
).

Various studies have already described selection of reference genes associated with different experimental conditions in bovines and have selected reference genes related to the mammary gland, fat, liver (Saremi et al., 2012Saremi, B.; Sauerwein, H.; Dänicke, S. and Mielenz, M. 2012. Technical note: Identification of reference genes for gene expression studies in different bovine tissues focusing on different fat depots. Journal of Dairy Science 95:3131-3138. https://doi.org/10.3168/jds.2011-4803
https://doi.org/10.3168/jds.2011-4803...
; Bonnet et al., 2013Bonnet, M.; Bernard, L.; Bes, S. and Leroux, C. 2013. Selection of reference genes for quantitative real-time PCR normalisation in adipose tissue, muscle, liver and mammary gland from ruminants. Animal 7:1344-1353. https://doi.org/10.1017/S1751731113000475
https://doi.org/10.1017/S175173111300047...
), skeletal muscle (Mberema and Sparagano, 2017Mberema, C. H. H. and Sparagano, O. A. E. 2017. Expression stability of reference genes in the skeletal muscles of beef cattle. African Journal of Biotechnology 16:261-267. https://doi.org/10.5897/AJB2016.15140
https://doi.org/10.5897/AJB2016.15140...
), somatic cells (Muhaghegh-Dolatabady et al., 2017Muhaghegh-Dolatabady, M.; Hossainy-Dolatabady, H.; Heidari Arjlo, E. and Mahmoudi, R. 2017. Validation of reference genes for real time PCR normalization in milk somatic cells of Holstein dairy cattle. Iranian Journal of Applied Animal Science 7:229-234.), blood (Lozano-Villegas et al., 2021Lozano-Villegas, K. J.; Rodríguez-Hernández, R.; Herrera-Sánchez, M. P.; Uribe-García, H. F.; Naranjo-Gómez, J. S.; Otero-Arroyo, R. J. and Rondón-Barragán, I. S. 2021. Identification of reference genes for expression studies in the whole-blood from three cattle breeds under two states of livestock weather safety. Animals 11:3073. https://doi.org/10.3390/ani11113073
https://doi.org/10.3390/ani11113073...
), and kidney, pituitary gland, and thyroid (Lisowski et al., 2008Lisowski, P.; Pierzchała, M.; Gościk, J.; Pareek, C. S. and Zwierzchowski, L. 2008. Evaluation of reference genes for studies of gene expression in the bovine liver, kidney, pituitary, and thyroid. Journal of Applied Genetics 49:367-372. https://doi.org/10.1007/bf03195635
https://doi.org/10.1007/bf03195635...
). However, there are no studies regarding the reference genes more suitable to muscle, liver, and jejunum in Zebu cattle under nutritional stress.

The hypothesis of this study is that the stability of reference genes differs among genetic groups and diets, and that this may result in errors in analysis of gene expression. Thus, this study aimed to evaluate the stability of candidate reference genes and select the genes that can be used for normalizing RT-PCR in the liver, skeletal muscle, and jejunum tissues of Nellore or Nellore × Angus steers fed different diets.

2. Material and Methods

2.1. Animals, diets, and sample collection

Research on animals was conducted according to ethical standards and approved by the local Ethics and Animal Welfare Committee (protocol 056/15).

Fourteen Nellore purebred steers and 14 Angus × Nellore crossbred steers with mean age of 25 months and body weight of 353±25.3 kg were housed in individual pens into a completely randomized design in a 2 × 2 factorial arrangement (two diets and two breeds). Half of the animals of each genetic group received a diet containing whole shelled corn (WSC) and the other half received whole shelled corn and sugarcane bagasse (WSCB) (Table 1).

Table 1
Percentage of ingredients and chemical composition of experimental diets

The animals were slaughtered after the experimental period of 96 days preceded by a 20-day adaptation period. Muscle, liver, and jejunum samples were collected soon after slaughter, labeled, and placed in 5-mL cryogenic tubes, frozen, and transported in liquid nitrogen until being stored in an ultra-freezer (−80 °C).

2.2. Design and efficiency of the primers

The primers of the reference genes were designed by means of sequences registered and published in the GenBank public database, at the National Center for Biotechnology Information (NCBI) platform. Characterization of the genes and the Open Reading Frames (ORF) of the selected sequences were obtained by the ORFinder tool of the NCBI. The primers were designed using the Oligo Perfect Designer software based on the sequence accessed in the GenBank and analyzed using OligoAnalizer 3.1.

Primers were then commercially synthesized (Invitrogen, Carlsbad, CA, USA), and the coefficient of regression (R2) and its efficiencies (E%) (Table 2) were determined from a pool of samples diluted at ratios of 1:5, 1:25, 1:125, 1:625, and 1:3125.

Table 2
Description of the candidate reference genes and of the target genes for bovines

2.3. RNA extraction and cDNA synthesis

RNA extraction was performed using two methods, in accordance with the tissue studied to avoid degradation and also to get better quality and quantity. RNA was extracted from muscle samples (80-100 mg) using the reagent Trizol (Invitrogen, Gaithersburg, MD, USA) according to manufacturer’s instructions and was homogenized (IKA® T18 ULTRA-TURRAX® Basic Homogenizer, Wilmington, NC, USA) together with 1 mL of Trizol reagent. Adaptations were made in relation to the amounts of chloroform (300 µL) and isopropanol (400 µL), and the step of chloroform addition was performed twice. Samples were treated with DNase DNA free (Ambion, Austin, TX, USA) to eliminate possible contaminations.

The liver and jejunum RNA were extracted using the Promega SV RNA Isolation kit (Promega, Madison, WI, USA) according to manufacturers. Tissue samples were weighed (130 mg) and homogenized (IKA® T18 ULTRA-TURRAX® Basic Homogenizer, Wilmington, NC, USA with 1 mL of RNA Lysis Buffer (RLA).

Samples were then quantified in a nanospectrophotometer (DeNovix DS-11, Wilmington, DE, USA) at A260 nm for the purpose of determining the amount (ng/μL) and purity (260/280 and 260/230) of the extracted RNA. After these procedures, the RNA samples were stored at −80 ℃.

Two methods were used for cDNA synthesis: the High-Capacity cDNA Reverse Transcription kit (Applied Biosystems, Foster City, CA, USA) for muscle samples, and the GoScript™ reverse transcription kit (Promega, Madison, WI, USA) for liver and jejunum samples, according to manufacturers’ specifications.

2.4. Quantitative gene expression analysis

For quantitative gene expression analysis by RT-qPCR, the Mastercycler® ep Realplex (Eppendorf) model was used with the SYBR Green detection system (Applied Biosystems, Foster City, CA, USA) and the cDNA obtained from RNA extracted from the liver, muscle, and jejunum of bovines. The thermal reaction conditions were 2 min at 50 ℃ and 10 min at 95 ℃, followed by 40 cycles of 15 s at 95 ℃ and 1 min at 60 ℃, ending with the melting curve. For each reaction, 1 μL of cDNA at the concentration of 1:5, 0.3 μL of each primer (1.2 µM), and 5.0 μL of Master Mix SYBR Green were used, for a final volume of 10.0 μL/sample on a MicroAmp Optical 96-well reaction plate (Applied Biosystems, Foster City, CA, USA).

The entire RT-qPCR experiment for each gene under study was conducted from cDNAs obtained from seven biological replications, with three technical replicates for each one. The Cq was determined by the number of cycles in which the fluorescence generated within a reaction crosses the baseline (threshold).

2.5. Data and statistical analysis

2.5.1. Evaluation of reference genes

The data were calibrated by the GenEX program (MultiD, Gothenburg, Sweden) for correcting the values according to the efficiency found for each gene evaluated. The SigmaPlot 12.0 software was employed for construction of boxplot diagrams to illustrate the levels and variations in expression of candidate reference genes. Data were then analyzed using the RefFinder online program, and thus the more stable genes for each condition were determined. Pairwise evaluation was carried out by the GeNorm algorithm to determine the minimum number of genes necessary; we needed to add one gene as a reference when the standard deviation (SD) was greater than 0.15.

2.5.2. Validation of expression of reference genes

For validation of the reference genes, a completely randomized design in a 2 × 2 factorial arrangement was used, with four treatments [Nellore steers fed a whole shelled corn diet (NWSC); Nellore steers fed a whole shelled corn and sugarcane bagasse diet (NWSCB); Nellore × Angus steers fed a whole shelled corn diet (NAWSC); and Nellore × Angus steers fed a whole shelled corn and sugarcane bagasse diet (NAWSCB)] and seven replications per treatment. Relative expression of target genes were calculated according to the method described by Pfaffl (2001)Pfaffl, M. W. 2001. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Research 29:e45. https://doi.org/10.1093/nar/29.9.e45
https://doi.org/10.1093/nar/29.9.e45...
, which is based on Cq values that are corrected for the amplification efficiency of each primer pair.

Gene expressions were analyzed using the PROC MIXED procedure of SAS (Statistical Analysis System, version 9.2), through the model:

γ i j = μ + α i + β j + α β i j + e i j (1)

in which γij is the response variable measured, μ is the overall mean, αi is the i-th level of the fixed effect of diet (two levels), βj is the j-th level of the fixed effect of breed (two levels), αβij is the interaction of breed (βj) and diet (αi), and eij is the random error associated with γij.

A Shapiro-Wilk test was used to analyze data normality, when gene expression data did not exhibit normal distribution, the data were transformed using the RANK procedure of SAS. After that, analysis of variance (ANOVA) and Tukey’s test were used on the validation data of reference genes to determine statistical differences among the means of the different treatments.

3. Results

Levels of gene expression were determined from values of the Quantification cycle (Cq) of the seven reference genes evaluated in the different tissues of the animals, regardless of their genetic group and diets. Simultaneous evaluation of the effect of breeds and diets on muscle, liver, and jejunum tissue (Figure 1) showed that the lowest mean Cq values were 11.9, 11.9, and 13.7 and the highest were 25.7, 27, and 27.4, respectively.

Figure 1
Variation in the expression of candidate reference genes in the breeds in different tissues.

The ranking of candidate reference genes was obtained by the RefFinder program from the geometric mean of the results from the software programs GeNorm, BestKeeper, NormFinder, and Delta Cq. The minimum number of genes for tissues, breeds, and the combination between variables were determined by pairwise analysis by the GeNorm algorithm (Vandesompele et al., 2002Vandesompele, J.; De Preter, K.; Pattyn, F.; Poppe, B.; Van Roy, N.; De Paepe, A. and Speleman, F. 2002. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biology 3. https://doi.org/10.1186/gb-2002-3-7-research0034
https://doi.org/10.1186/gb-2002-3-7-rese...
). Vandesompele et al. (2002)Vandesompele, J.; De Preter, K.; Pattyn, F.; Poppe, B.; Van Roy, N.; De Paepe, A. and Speleman, F. 2002. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biology 3. https://doi.org/10.1186/gb-2002-3-7-research0034
https://doi.org/10.1186/gb-2002-3-7-rese...
proposed 0.15 as the cutoff value for the paired variation. After inclusion of the second reference gene, a substantial reduction in the V value (0.15) was observed in GeNorm, showing the need for at least two and at most four reference genes for normalization of the study for the different conditions of this study (Tables 3-5). For the muscle of pure and crossbred animals fed WSC four reference genes were needed (Table 3); on the other hand, for liver (Table 4) and jejunum (Table 5), the minimum number of genes was three.

Table 3
Reference genes and minimum number of genes recommended for normalization in real time PCR studies in pure- (Nellore) and crossbred (Nellore × Angus) animals in skeletal muscle
Table 4
Reference genes and minimum number of genes recommended for normalization in real time PCR studies in pure- (Nellore) and crossbred (Nellore × Angus) animals in the liver
Table 5
Reference genes and minimum number of genes recommended for normalization in real time PCR studies in pure- (Nellore) and crossbred (Nellore × Angus) animals in the jejunum

The results obtained showed that the ranking (more stable to the less stable gene) for the muscle according to the RefFinder program were Ribosomal Protein 18S (18S), Beta-actin (ACTB), Cancer susceptibility candidate 3 (CASC3), Eukaryotic translation elongation factor 1 alpha 2 (EEF1A2), Glyceraldehyde-3-phosphate (GAPDH), Hydroxymethylbilane synthase mRNA (HMBS), and Ubiquitin C (UBC) (Table 6), being 18S, ACTB, CASC3, and EEF1A2 (Table 3) the most stable genes. The use of CASC3 in purebred animals and of EEF1A2 was expendable in comparing the different genetic groups and diets simultaneously (Table 3). For the liver, the ranking of the reference genes (RefFinder method) was HMBS, ACTB, 18S, CASC3, GAPDH, UBC, and EEF1A2 (Table 7) and the most stable genes were HMBS, ACTB, 18S, and GAPDH (Table 4). In the crossbred animals receiving the two diets or only WSC, GAPDH was more stable than the 18S gene that was selected in purebred animals. Finally, in the jejunum, the most stable reference genes were GAPDH, ACTB, CASC3, and 18S (Table 5), considering that the ranking of the reference genes were GAPDH, ACTB, CASC3, UBC, HMBS, 18S, and EEF1A2 (Table 8). In general, there was a small effect of diet on the reference genes in the tissues.

Table 6
Ranking of reference genes in the muscle of pure- and crossbred animals, regardless of the diet, by the different evaluation methods

Table 7
Ranking of reference genes in the liver of pure- and crossbred animals, regardless of the diet, by the different evaluation methods

Table 8
Ranking of reference genes in the jejunum of pure- and crossbred animals, regardless of the diet, by the different evaluation methods

4. Discussion

A reference gene should have expression near the target gene; however, that is not always possible. Thus, the general directives recommend that the Cq values be from 15 to 30 (Wan et al., 2010Wan, H.; Zhao, Z.; Qian, C.; Sui, Y.; Malik, A. A. and Chen, J. 2010. Selection of appropriate reference genes for gene expression studies by quantitative real-time polymerase chain reaction in cucumber. Analytical biochemistry 399:257-261. https://doi.org/10.1016/j.ab.2009.12.008
https://doi.org/10.1016/j.ab.2009.12.008...
). In spite of the wide variation found in the Cq of the genes tested (Figure 1), the values found in this study are at the established interval (15 < Cq < 30), except the 18S gene. Although it had very low Cq, the 18S gene was considered stable for muscle and liver according to the RefFinder program.

In addition, other studies have also found a wide variation in Cq in evaluating reference genes in plants or animals under different experimental conditions (Ling et al., 2014Ling, H.; Wu, Q.; Guo, J.; Xu, L. and Que, Y. 2014. Comprehensive selection of reference genes for gene expression normalization in sugarcane by real time quantitative RT-PCR. PLoS One 9:e97469. https://doi.org/10.1371/journal.pone.0097469
https://doi.org/10.1371/journal.pone.009...
; Olias et al., 2014Olias, P.; Adam, I.; Meyer, A.; Scharff, C. and Gruber, A. D. 2014. Reference genes for quantitative gene expression studies in multiple avian species. PLoS One 9:e99678. https://doi.org/10.1371/journal.pone.0099678
https://doi.org/10.1371/journal.pone.009...
; Robledo et al., 2014Robledo, D.; Hernández-Urcera, J.; Cal, R. M.; Pardo, B. G.; Sánchez, L.; Martínez, P. and Viñas, A. 2014. Analysis of qPCR reference gene stability determination methods and a practical approach for efficiency calculation on a turbot (Scophthalmus maximus) gonad dataset. BMC Genomics 15:648. https://doi.org/10.1186/1471-2164-15-648
https://doi.org/10.1186/1471-2164-15-648...
; Gentile et al., 2016Gentile, A.-M.; Lhamyani, S.; Coín-Aragüez, L.; Oliva-Olivera, W.; Zayed, H.; Vega-Rioja, A.; Monteseirin, J.; Romero-Zerbo, S.-Y.; Tinahones, F.-J.; Bermúdez-Silva, F.-J. and El Bekay, R. 2016. RPL13A and EEF1A1 are suitable reference genes for qPCR during adipocyte differentiation of vascular stromal cells from patients with different BMI and HOMA-IR. PLoS One 11:e0157002. https://doi.org/10.1371/journal.pone.0157002
https://doi.org/10.1371/journal.pone.015...
; Die et al., 2017Die, J. V.; Baldwin, R. L.; Rowland, L. J.; Li, R.; Oh, S.; Li, C.; Connor, E. E. and Ranilla, M.-J. 2017. Selection of internal reference genes for normalization of reverse transcription quantitative polymerase chain reaction (RT-qPCR) analysis in the rumen epithelium. PLoS One 12:e0172674. https://doi.org/10.1371/journal.pone.0172674
https://doi.org/10.1371/journal.pone.017...
; Fernandes-Brum et al., 2017Fernandes-Brum, C. N.; Garcia, B. O.; Moreira, R. O.; Ságio, S. A.; Barreto, H. G.; Lima, A. A.; Freitas, N. C.; Lima, R. R.; Carvalho, C. H. S. and Chalfun-Júnior, A. 2017. A panel of the most suitable reference genes for RT-qPCR expression studies of coffee: screening their stability under different conditions. Tree Genetics & Genomes 13:131. https://doi.org/10.1007/s11295-017-1213-1
https://doi.org/10.1007/s11295-017-1213-...
).

Some studies have already described the selection of reference genes associated with different tissues and conditions in bovines (Table 9) and showed that the choice of genes is highly dependent on experimental conditions and that most of the studies were carried out with milk cattle and taurine cattle. Saremi et al. (2012)Saremi, B.; Sauerwein, H.; Dänicke, S. and Mielenz, M. 2012. Technical note: Identification of reference genes for gene expression studies in different bovine tissues focusing on different fat depots. Journal of Dairy Science 95:3131-3138. https://doi.org/10.3168/jds.2011-4803
https://doi.org/10.3168/jds.2011-4803...
found different reference genes recommended for the liver and muscle tissues of bovines fed conjugated linoleic acid (CLA). Lisowski et al. (2008)Lisowski, P.; Pierzchała, M.; Gościk, J.; Pareek, C. S. and Zwierzchowski, L. 2008. Evaluation of reference genes for studies of gene expression in the bovine liver, kidney, pituitary, and thyroid. Journal of Applied Genetics 49:367-372. https://doi.org/10.1007/bf03195635
https://doi.org/10.1007/bf03195635...
evaluated different cattle breeds in different developmental stages and likewise found that two or more reference genes are necessary and that they differ according to the tissue/organ evaluated (liver, kidney, pituitary gland, or thyroid gland).

Table 9
Studies carried out on validation of references genes in bovines

Studies carried out prior to publication of the Minimum Information for Publication of Quantitative Real-Time PCR Experiments (MIQE guidelines) used only one reference gene (Hamalainen et al., 2001Hamalainen, H. K.; Tubman, J. C.; Vikman, S.; Kyrölä, T.; Ylikoski, E.; Warrington, J. A. and Lahesmaa, R. 2001. Identification and validation of endogenous reference genes for expression profiling of T helper cell differentiation by quantitative real-time RT-PCR. Analytical Biochemistry 299:63-70. https://doi.org/10.1006/abio.2001.5369
https://doi.org/10.1006/abio.2001.5369...
; Deindl et al., 2002Deindl, E.; Boengler, K.; van Royen, N. and Schaper, W. 2002. Differential expression of GAPDH and β-actin in growing collateral arteries. Molecular and Cellular Biochemistry 236:139-146. https://doi.org/10.1023/A:1016166127465
https://doi.org/10.1023/A:1016166127465...
; Glare et al., 2002Glare, E. M.; Divjak, M.; Bailey, M. J. and Walters, E. H. 2002. β-Actin and GAPDH housekeeping gene expression in asthmatic airways is variable and not suitable for normalising mRNA levels. Thorax 57:765-770. https://doi.org/10.1136/thorax.57.9.765
https://doi.org/10.1136/thorax.57.9.765...
) to normalize RT-qPCR data. However, after its publication, the use of two or more reference genes has been recommended to ensure stable expression under all experimental conditions (Bustin et al., 2009Bustin, S. A.; Benes, V.; Garson, J. A.; Hellemans, J.; Huggett, J.; Kubista, M.; Mueller, R.; Nolan, T.; Pfaffl, M. W.; Shipley, G. L.; Vandesompele, J. and Wittwer, C. T. 2009. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clinical Chemistry 55:611-622. https://doi.org/10.1373/clinchem.2008.112797
https://doi.org/10.1373/clinchem.2008.11...
). In this study, as a consequence of the data analysis by the GeNorm software, the need for at least two, and up to four, reference genes was shown for normalization of the study for the different breeds and diets studied and for the skeletal muscle, liver, and jejunum tissues. Likewise, Pérez et al. (2008)Pérez, R.; Tupac-Yupanqui, I. and Dunner, S. 2008. Evaluation of suitable reference genes for gene expression studies in bovine muscular tissue. BMC Molecular Biology 9:79. https://doi.org/10.1186/1471-2199-9-79
https://doi.org/10.1186/1471-2199-9-79...
and Saremi et al. (2012)Saremi, B.; Sauerwein, H.; Dänicke, S. and Mielenz, M. 2012. Technical note: Identification of reference genes for gene expression studies in different bovine tissues focusing on different fat depots. Journal of Dairy Science 95:3131-3138. https://doi.org/10.3168/jds.2011-4803
https://doi.org/10.3168/jds.2011-4803...
showed the need for three genes for muscle; and Lisowski et al. (2008)Lisowski, P.; Pierzchała, M.; Gościk, J.; Pareek, C. S. and Zwierzchowski, L. 2008. Evaluation of reference genes for studies of gene expression in the bovine liver, kidney, pituitary, and thyroid. Journal of Applied Genetics 49:367-372. https://doi.org/10.1007/bf03195635
https://doi.org/10.1007/bf03195635...
and Saremi et al. (2012)Saremi, B.; Sauerwein, H.; Dänicke, S. and Mielenz, M. 2012. Technical note: Identification of reference genes for gene expression studies in different bovine tissues focusing on different fat depots. Journal of Dairy Science 95:3131-3138. https://doi.org/10.3168/jds.2011-4803
https://doi.org/10.3168/jds.2011-4803...
showed the need for two to three reference genes for the liver. However, we did not find studies that investigated the reference genes more suitable for the jejunum of bovines, which is important for studies of gene expression of membrane transport proteins.

A review of literature regarding reference genes in gene expression studies showed that since the publication of the MIQE, the most common choice of reference genes were ACTB (used in 38% of the studies) and GAPDH (37%); the 18S gene (12%) was a little less used in vertebrates (Chapman and Waldenström, 2015Chapman, J. R. and Waldenström, J. 2015. With reference to reference genes: a systematic review of endogenous controls in gene expression studies. PLoS One 10:e0141853. https://doi.org/10.1371/journal.pone.0141853
https://doi.org/10.1371/journal.pone.014...
). For the tissues evaluated, at least one of the three genes cited above was stable, making them strong candidates for reference genes in studies with RT-qPCR in zebu cattle and their crosses. The ACTB was among the more stable reference genes under all the conditions in this study (Tables 3-5). One of the reasons for this is that ACTIN is the most abundant protein in eukaryotic cells (Nakajima-Iijima et al., 1985Nakajima-Iijima, S.; Hamada, H.; Reddy, P. and Kakunaga, T. 1985. Molecular structure of the human cytoplasmic beta-actin gene: interspecies homology of sequences in the introns. Proceedings of the National Academy of Sciences of the United States of America 82:6133-6137. https://doi.org/10.1073/pnas.82.18.6133
https://doi.org/10.1073/pnas.82.18.6133...
). There are at least six ACTIN isoforms in vertebrates: four types in muscle (skeletal, cardiac, vascular smooth, and stomach smooth actins) and two non-muscular types (β-cytoplasmic and γ-cytoplasmic actins) (Vandekerckhove and Weber, 1979Vandekerckhove, J. and Weber, K. 1979. The complete amino acid sequence of actins from bovine aorta, bovine heart, bovine fast skeletal muscle, and rabbit, slow skeletal muscle: a protein-chemical analysis of muscle actin differentiation. Differentiation 14:123-133. https://doi.org/10.1111/j.1432-0436.1979.tb01021.x
https://doi.org/10.1111/j.1432-0436.1979...
), and they are expressed in all types of tissues, as shown in the results of this study, making this gene an optimal normalizer (Tables 3-5).

The EEF1A2 gene was considered the least stable reference gene under most of the conditions studied in the liver and jejunum, and UBC was the least stable in muscle (Tables 6-8). Both are genes expressed under diverse conditions; however, in this study, there was greater variation in relation to the tissues studied. The EEF1A gene has two isoforms, EEF1A1 and EEF1A2, that codify the second most abundant protein (after ACTIN) and catalyzes bonding of tRNA to the ribosomal site. The isoform EEF1A2 is present especially in the brain, heart, and muscle (Lee and Surh, 2009Lee, M. H. and Surh, Y. J. 2009. eEF1A2 as a putative oncogene. Annals of the New York Academy of Sciences 1171:87-93. https://doi.org/10.1111/j.1749-6632.2009.04909.x
https://doi.org/10.1111/j.1749-6632.2009...
). The expression of this gene in muscle was greater and more stable, with a mean Cq value of 17.6 (Figure 1a), compared with liver (Figure 1b) and jejunum (Figure 1c) tissues, with mean values of 27 and 26.2, respectively. In the latter two tissues, this gene was selected as one of the three least stable (Tables 7 and 8), whereas for muscle, the gene was stable in the purebred and crossbred animals, regardless of the diet fed (Table 6).

In muscle, the UBC gene was highly variable in the breeds and diets studied, and thus, it is not considered a good normalizer for this tissue (Figure 1a). A possible explanation for this is that the gene codifies a protein that participates in selective degradation of short-lived proteins in eukaryotic cells and is involved in immune response and cell development and death (Hershko and Ciechanover, 1998Hershko, A. and Ciechanover, A. 1998. The Unibiquitin System. Annual Review of Biochemistry 67:425-479. https://doi.org/10.1146/annurev.biochem.67.1.425
https://doi.org/10.1146/annurev.biochem....
).

Based on the results that show differences in the stability of candidate reference genes, statistical analyses were performed using the more stable and less stable genes found in this study to reveal possible errors in analysis of gene expression of target genes.

Fatty Acid Synthase (FASN) and Acyl-CoA Oxidase 1 (ACOX) are genes involved in synthesis and oxidation of saturated fatty acids, and the protein coded by ACOX is the first enzyme of beta-oxidation of fatty acids in the peroxisomes (Schulz, 1996Schulz, H. 1996. Chapter 3 - Oxidation of fatty acids. p.75-99. In: New comprehensive biochemistry. Vance, D. E. and Vance, J. E., eds. Elsevier. https://doi.org/10.1016/S0167-7306(08)60510-4
https://doi.org/10.1016/S0167-7306(08)60...
). The gene FASN, in its turn, codifies the enzyme responsible for the final steps of de novo synthesis of long-chain fatty acids from acetyl-CoA, malonyl-CoA, and NADPH (Ladeira et al., 2016Ladeira, M. M.; Schoonmaker, J. P.; Gionbelli, M. P.; Dias, J. C. O.; Gionbelli, T. R. S.; Carvalho, J. R. R. and Teixeira, P. D. 2016. Nutrigenomics and beef quality: A review about lipogenesis. International Journal of Molecular Sciences 17:918. https://doi.org/10.3390/ijms17060918
https://doi.org/10.3390/ijms17060918...
). When the more stable genes were used for normalization of the study (Figure 2), both FASN and ACOX target genes in the muscle tissue were more expressed in WSC diets, and FASN was also more expressed in crossbred animals. However, upon normalizing the data with the less stable genes, statistical differences were not found for the diets studied, and only tendencies were found for the breeds in both genes. Therefore, these results show that the wrong choice of reference genes would lead to misinterpretation and distortion of results.

Figure 2
Expression level of target genes FASN and ACOX in muscle, normalized by the more stable and less stable reference genes.

The target gene Stearoyl-CoA Desaturase (SCD1) was evaluated in the liver (Figure 3), and Maltase-Glucoamylase 2 (MGAM) and Solute Carrier Family 2 Member 1 (SLC2A1) in the jejunum (Figure 4). The SCD1 gene is present throughout the organism and is highly expressed in the liver (Ntambi and Miyazaki, 2004Ntambi, J. M. and Miyazaki, M. 2004. Regulation of stearoyl-CoA desaturases and role in metabolism. Progress in Lipid Research 43:91-104. https://doi.org/10.1016/S0163-7827(03)00039-0
https://doi.org/10.1016/S0163-7827(03)00...
). It is responsible for codifying an enzyme responsible for desaturation of long-chain fatty acids, thus acting in regulation of membrane fluidity (Ntambi et al., 1988Ntambi, J. M.; Buhrow, S. A.; Kaestner, K. H.; Christy, R. J.; Sibley, E.; Kelly, T. J. and Lane, M. D. 1988. Differentiation-induced gene expression in 3T3-L1 preadipocytes. Characterization of a differentially expressed gene encoding stearoyl-CoA desaturase. Journal of Biological Chemistry 263:17291-17300. https://doi.org/10.1016/S0021-9258(19)77834-X
https://doi.org/10.1016/S0021-9258(19)77...
; Ntambi and Miyazaki, 2004Ntambi, J. M. and Miyazaki, M. 2004. Regulation of stearoyl-CoA desaturases and role in metabolism. Progress in Lipid Research 43:91-104. https://doi.org/10.1016/S0163-7827(03)00039-0
https://doi.org/10.1016/S0163-7827(03)00...
). Evaluation of expression of the SCD1 gene in the liver showed that upon using the more stable genes, statistical effects of the diets and breeds would not be found. In contrast, with the use of the less stable genes, the inclusion of sugarcane bagasse in the diet would increase expression of the gene in the tissue (Figure 3). In the jejunum, the MGAM gene is expressed in the brush border of the small intestine membrane and codifies an enzyme responsible for hydrolyzing maltose into glucose (Mochizuki et al., 2010Mochizuki, K.; Honma, K.; Shimada, M. and Goda, T. 2010. The regulation of jejunal induction of the maltase–glucoamylase gene by a high‐starch/low‐fat diet in mice. Molecular Nutrition & Food Research 54:1445-1451. https://doi.org/10.1002/mnfr.200900467
https://doi.org/10.1002/mnfr.200900467...
). Statistical differences were not found in expression of the MGAM gene in the breeds and diets evaluated when the less stable reference genes were used. However, upon using the more stable ones, there was a tendency that the diets containing sugarcane bagasse would result in greater expression of the gene in the jejunum (Figure 4). Finally, the SLC2A1, which is expressed in the enterocytes of the small intestine and codifies the sodium-glucose linked transporter (SGLT1), has its expression similar for diets and breeds, regardless of selection of more stable or less stable genes (Figure 4). In this case, there was statistical difference for diets in both groups of reference genes used as normalizers in this study. However, the relative expression with the use of the less stable genes was much greater, which can also lead to poor interpretation of data by authors.

Figure 3
Expression level of target gene SCD1 in the liver, normalized by the more stable and less stable reference genes.

Figure 4
Expression level of target genes MGAM and SLC2A1 in the jejunum, normalized by the more stable and less stable reference genes.

It can be inferred that the choice of reference genes highly affects the results of analyses of gene expression that can occur due to differences brought about under determined metabolic conditions of determined tissues. Thus, there is the possibility of results exhibiting wide variation in gene expression. In the present study, there were few variations in the expression and stability of reference genes among the genetic groups and the diets studied within a specific tissue. However, analysis of the expression of the reference genes among the tissues showed that expressions of GAPDH and EEF1A2 were greater and less stable in muscle compared with the liver and jejunum. Therefore, these genes would likely not be considered ideal for normalization of studies comparing these three tissues.

5. Conclusions

The more stable genes for zebu animals and their crosses were 18S, ACTB, and CASC3 for muscle tissue; HMBS, ACTB, and 18S for liver tissue; and GAPDH, ACTB, and CASC3 for jejunum tissue, regardless of the breed and diet. Comparison of the more stable and less stable genes as a reference for normalization of target genes shows that the wrong choice of these genes will lead to misleading interpretation of the results.

Acknowledgments

The authors would like to thank Conselho Nacional de Desenvolvimento Científico e Tecnológico (grant number: 457779/2014-4), Fundação de Amparo à Pesquisa do Estado de Minas Gerais (Belo Horizonte, Minas Gerais, Brazil) (grant number: CVZ–PPM: 00441-14), and Cargill (Itapira, São Paulo, Brazil) for funding this study.

References

  • Bionaz, M. and Loor, J. J. 2007. Identification of reference genes for quantitative real-time PCR in the bovine mammary gland during the lactation cycle. Physiological Genomics 29:312-319. https://doi.org/10.1152/physiolgenomics.00223.2006
    » https://doi.org/10.1152/physiolgenomics.00223.2006
  • Bonnet, M.; Bernard, L.; Bes, S. and Leroux, C. 2013. Selection of reference genes for quantitative real-time PCR normalisation in adipose tissue, muscle, liver and mammary gland from ruminants. Animal 7:1344-1353. https://doi.org/10.1017/S1751731113000475
    » https://doi.org/10.1017/S1751731113000475
  • Bustin, S. A.; Benes, V.; Garson, J. A.; Hellemans, J.; Huggett, J.; Kubista, M.; Mueller, R.; Nolan, T.; Pfaffl, M. W.; Shipley, G. L.; Vandesompele, J. and Wittwer, C. T. 2009. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clinical Chemistry 55:611-622. https://doi.org/10.1373/clinchem.2008.112797
    » https://doi.org/10.1373/clinchem.2008.112797
  • Chapman, J. R. and Waldenström, J. 2015. With reference to reference genes: a systematic review of endogenous controls in gene expression studies. PLoS One 10:e0141853. https://doi.org/10.1371/journal.pone.0141853
    » https://doi.org/10.1371/journal.pone.0141853
  • Deindl, E.; Boengler, K.; van Royen, N. and Schaper, W. 2002. Differential expression of GAPDH and β-actin in growing collateral arteries. Molecular and Cellular Biochemistry 236:139-146. https://doi.org/10.1023/A:1016166127465
    » https://doi.org/10.1023/A:1016166127465
  • Die, J. V.; Baldwin, R. L.; Rowland, L. J.; Li, R.; Oh, S.; Li, C.; Connor, E. E. and Ranilla, M.-J. 2017. Selection of internal reference genes for normalization of reverse transcription quantitative polymerase chain reaction (RT-qPCR) analysis in the rumen epithelium. PLoS One 12:e0172674. https://doi.org/10.1371/journal.pone.0172674
    » https://doi.org/10.1371/journal.pone.0172674
  • Fernandes-Brum, C. N.; Garcia, B. O.; Moreira, R. O.; Ságio, S. A.; Barreto, H. G.; Lima, A. A.; Freitas, N. C.; Lima, R. R.; Carvalho, C. H. S. and Chalfun-Júnior, A. 2017. A panel of the most suitable reference genes for RT-qPCR expression studies of coffee: screening their stability under different conditions. Tree Genetics & Genomes 13:131. https://doi.org/10.1007/s11295-017-1213-1
    » https://doi.org/10.1007/s11295-017-1213-1
  • Gentile, A.-M.; Lhamyani, S.; Coín-Aragüez, L.; Oliva-Olivera, W.; Zayed, H.; Vega-Rioja, A.; Monteseirin, J.; Romero-Zerbo, S.-Y.; Tinahones, F.-J.; Bermúdez-Silva, F.-J. and El Bekay, R. 2016. RPL13A and EEF1A1 are suitable reference genes for qPCR during adipocyte differentiation of vascular stromal cells from patients with different BMI and HOMA-IR. PLoS One 11:e0157002. https://doi.org/10.1371/journal.pone.0157002
    » https://doi.org/10.1371/journal.pone.0157002
  • Glare, E. M.; Divjak, M.; Bailey, M. J. and Walters, E. H. 2002. β-Actin and GAPDH housekeeping gene expression in asthmatic airways is variable and not suitable for normalising mRNA levels. Thorax 57:765-770. https://doi.org/10.1136/thorax.57.9.765
    » https://doi.org/10.1136/thorax.57.9.765
  • Hamalainen, H. K.; Tubman, J. C.; Vikman, S.; Kyrölä, T.; Ylikoski, E.; Warrington, J. A. and Lahesmaa, R. 2001. Identification and validation of endogenous reference genes for expression profiling of T helper cell differentiation by quantitative real-time RT-PCR. Analytical Biochemistry 299:63-70. https://doi.org/10.1006/abio.2001.5369
    » https://doi.org/10.1006/abio.2001.5369
  • Hershko, A. and Ciechanover, A. 1998. The Unibiquitin System. Annual Review of Biochemistry 67:425-479. https://doi.org/10.1146/annurev.biochem.67.1.425
    » https://doi.org/10.1146/annurev.biochem.67.1.425
  • Kaur, R.; Sodhi, M.; Sharma, A.; Sharma, V. L.; Verma, P.; Swami, S. K.; Kumari, P. and Mukesh, M. 2018. Selection of suitable reference genes for normalization of quantitative RT-PCR (RT-qPCR) expression data across twelve tissues of riverine buffaloes (Bubalus bubalis). PLoS One 13:e0191558. https://doi.org/10.1371/journal.pone.0191558
    » https://doi.org/10.1371/journal.pone.0191558
  • Kishore, A.; Sodhi, M.; Khate, K.; Kapila, N.; Kumari, P. and Mukesh, M. 2013. Selection of stable reference genes in heat stressed peripheral blood mononuclear cells of tropically adapted Indian cattle and buffaloes. Molecular and Cellular Probes 27:140-144. https://doi.org/10.1016/j.mcp.2013.02.003
    » https://doi.org/10.1016/j.mcp.2013.02.003
  • Ladeira, M. M.; Schoonmaker, J. P.; Gionbelli, M. P.; Dias, J. C. O.; Gionbelli, T. R. S.; Carvalho, J. R. R. and Teixeira, P. D. 2016. Nutrigenomics and beef quality: A review about lipogenesis. International Journal of Molecular Sciences 17:918. https://doi.org/10.3390/ijms17060918
    » https://doi.org/10.3390/ijms17060918
  • Lee, M. H. and Surh, Y. J. 2009. eEF1A2 as a putative oncogene. Annals of the New York Academy of Sciences 1171:87-93. https://doi.org/10.1111/j.1749-6632.2009.04909.x
    » https://doi.org/10.1111/j.1749-6632.2009.04909.x
  • Ling, H.; Wu, Q.; Guo, J.; Xu, L. and Que, Y. 2014. Comprehensive selection of reference genes for gene expression normalization in sugarcane by real time quantitative RT-PCR. PLoS One 9:e97469. https://doi.org/10.1371/journal.pone.0097469
    » https://doi.org/10.1371/journal.pone.0097469
  • Lisowski, P.; Pierzchała, M.; Gościk, J.; Pareek, C. S. and Zwierzchowski, L. 2008. Evaluation of reference genes for studies of gene expression in the bovine liver, kidney, pituitary, and thyroid. Journal of Applied Genetics 49:367-372. https://doi.org/10.1007/bf03195635
    » https://doi.org/10.1007/bf03195635
  • Lopes, R. C.; Sampaio, C. B.; Trece, A. S.; Teixeira, P. D.; Gionbelli, T. R. S.; Santos, L. R.; Costa, T. C.; Duarte, M. S. and Gionbelli, M. P. 2020. Impacts of protein supplementation during late gestation of beef cows on maternal skeletal muscle and liver tissues metabolism. Animal 14:1867-1875. https://doi.org/10.1017/S1751731120000336
    » https://doi.org/10.1017/S1751731120000336
  • Lozano-Villegas, K. J.; Rodríguez-Hernández, R.; Herrera-Sánchez, M. P.; Uribe-García, H. F.; Naranjo-Gómez, J. S.; Otero-Arroyo, R. J. and Rondón-Barragán, I. S. 2021. Identification of reference genes for expression studies in the whole-blood from three cattle breeds under two states of livestock weather safety. Animals 11:3073. https://doi.org/10.3390/ani11113073
    » https://doi.org/10.3390/ani11113073
  • Mberema, C. H. H. and Sparagano, O. A. E. 2017. Expression stability of reference genes in the skeletal muscles of beef cattle. African Journal of Biotechnology 16:261-267. https://doi.org/10.5897/AJB2016.15140
    » https://doi.org/10.5897/AJB2016.15140
  • Mochizuki, K.; Honma, K.; Shimada, M. and Goda, T. 2010. The regulation of jejunal induction of the maltase–glucoamylase gene by a high‐starch/low‐fat diet in mice. Molecular Nutrition & Food Research 54:1445-1451. https://doi.org/10.1002/mnfr.200900467
    » https://doi.org/10.1002/mnfr.200900467
  • Muhaghegh-Dolatabady, M.; Hossainy-Dolatabady, H.; Heidari Arjlo, E. and Mahmoudi, R. 2017. Validation of reference genes for real time PCR normalization in milk somatic cells of Holstein dairy cattle. Iranian Journal of Applied Animal Science 7:229-234.
  • Nakajima-Iijima, S.; Hamada, H.; Reddy, P. and Kakunaga, T. 1985. Molecular structure of the human cytoplasmic beta-actin gene: interspecies homology of sequences in the introns. Proceedings of the National Academy of Sciences of the United States of America 82:6133-6137. https://doi.org/10.1073/pnas.82.18.6133
    » https://doi.org/10.1073/pnas.82.18.6133
  • Ntambi, J. M.; Buhrow, S. A.; Kaestner, K. H.; Christy, R. J.; Sibley, E.; Kelly, T. J. and Lane, M. D. 1988. Differentiation-induced gene expression in 3T3-L1 preadipocytes. Characterization of a differentially expressed gene encoding stearoyl-CoA desaturase. Journal of Biological Chemistry 263:17291-17300. https://doi.org/10.1016/S0021-9258(19)77834-X
    » https://doi.org/10.1016/S0021-9258(19)77834-X
  • Ntambi, J. M. and Miyazaki, M. 2004. Regulation of stearoyl-CoA desaturases and role in metabolism. Progress in Lipid Research 43:91-104. https://doi.org/10.1016/S0163-7827(03)00039-0
    » https://doi.org/10.1016/S0163-7827(03)00039-0
  • Olias, P.; Adam, I.; Meyer, A.; Scharff, C. and Gruber, A. D. 2014. Reference genes for quantitative gene expression studies in multiple avian species. PLoS One 9:e99678. https://doi.org/10.1371/journal.pone.0099678
    » https://doi.org/10.1371/journal.pone.0099678
  • Oliveira, D. M.; Chalfun-Junior, A.; Chizzotti, M. L.; Barreto, H. G.; Coelho, T. C.; Paiva, L. V.; Coelho, C. P.; Teixeira, P. D.; Schoonmaker, J. P. and Ladeira, M. M. 2014. Expression of genes involved in lipid metabolism in the muscle of beef cattle fed soybean or rumen-protected fat, with or without monensin supplementation. Journal of Animal Science 92:5426-5436. https://doi.org/10.2527/jas.2014-7855
    » https://doi.org/10.2527/jas.2014-7855
  • Pérez, R.; Tupac-Yupanqui, I. and Dunner, S. 2008. Evaluation of suitable reference genes for gene expression studies in bovine muscular tissue. BMC Molecular Biology 9:79. https://doi.org/10.1186/1471-2199-9-79
    » https://doi.org/10.1186/1471-2199-9-79
  • Pfaffl, M. W. 2001. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Research 29:e45. https://doi.org/10.1093/nar/29.9.e45
    » https://doi.org/10.1093/nar/29.9.e45
  • Robledo, D.; Hernández-Urcera, J.; Cal, R. M.; Pardo, B. G.; Sánchez, L.; Martínez, P. and Viñas, A. 2014. Analysis of qPCR reference gene stability determination methods and a practical approach for efficiency calculation on a turbot (Scophthalmus maximus) gonad dataset. BMC Genomics 15:648. https://doi.org/10.1186/1471-2164-15-648
    » https://doi.org/10.1186/1471-2164-15-648
  • Saremi, B.; Sauerwein, H.; Dänicke, S. and Mielenz, M. 2012. Technical note: Identification of reference genes for gene expression studies in different bovine tissues focusing on different fat depots. Journal of Dairy Science 95:3131-3138. https://doi.org/10.3168/jds.2011-4803
    » https://doi.org/10.3168/jds.2011-4803
  • Schulz, H. 1996. Chapter 3 - Oxidation of fatty acids. p.75-99. In: New comprehensive biochemistry. Vance, D. E. and Vance, J. E., eds. Elsevier. https://doi.org/10.1016/S0167-7306(08)60510-4
    » https://doi.org/10.1016/S0167-7306(08)60510-4
  • Teixeira, P. D.; Oliveira, D. M.; Chizzotti, M. L.; Chalfun-Junior, A.; Coelho, T. C.; Gionbelli, M. P.; Paiva, L. V.; Carvalho, J. R. R. and Ladeira, M. M. 2017. Subspecies and diet affect the expression of genes involved in lipid metabolism and chemical composition of muscle in beef cattle. Meat Science 133:110-118. https://doi.org/10.1016/j.meatsci.2017.06.009
    » https://doi.org/10.1016/j.meatsci.2017.06.009
  • Vandekerckhove, J. and Weber, K. 1979. The complete amino acid sequence of actins from bovine aorta, bovine heart, bovine fast skeletal muscle, and rabbit, slow skeletal muscle: a protein-chemical analysis of muscle actin differentiation. Differentiation 14:123-133. https://doi.org/10.1111/j.1432-0436.1979.tb01021.x
    » https://doi.org/10.1111/j.1432-0436.1979.tb01021.x
  • Vandesompele, J.; De Preter, K.; Pattyn, F.; Poppe, B.; Van Roy, N.; De Paepe, A. and Speleman, F. 2002. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biology 3. https://doi.org/10.1186/gb-2002-3-7-research0034
    » https://doi.org/10.1186/gb-2002-3-7-research0034
  • Varshney, N.; Mohanty, A. K.; Kumar, S.; Kaushik, J. K.; Dang, A. K.; Mukesh, M.; Mishra, B. P.; Kataria, R.; Kimothi, S. P.; Mukhopadhyay, T. K.; Malakar, D.; Prakash, B. S.; Grover, S. and Batish, V. K. 2012. Selection of suitable reference genes for quantitative gene expression studies in milk somatic cells of lactating cows (Bos indicus). Journal of Dairy Science 95:2935-2945. https://doi.org/10.3168/jds.2011-4442
    » https://doi.org/10.3168/jds.2011-4442
  • Wan, H.; Zhao, Z.; Qian, C.; Sui, Y.; Malik, A. A. and Chen, J. 2010. Selection of appropriate reference genes for gene expression studies by quantitative real-time polymerase chain reaction in cucumber. Analytical biochemistry 399:257-261. https://doi.org/10.1016/j.ab.2009.12.008
    » https://doi.org/10.1016/j.ab.2009.12.008

Publication Dates

  • Publication in this collection
    18 July 2022
  • Date of issue
    2022

History

  • Received
    22 July 2021
  • Accepted
    6 Apr 2022
Sociedade Brasileira de Zootecnia Universidade Federal de Viçosa / Departamento de Zootecnia, 36570-900 Viçosa MG Brazil, Tel.: +55 31 3612-4602, +55 31 3612-4612 - Viçosa - MG - Brazil
E-mail: rbz@sbz.org.br