Reference genes selection for quantitative real-time PCR using RankAggreg method in different tissues of Capra hircus

Temporal Transcriptome Dynamics of Longissimus dorsi Reveals the Mechanism of the Differences in Muscle Development and IMF Deposition between Fuqing Goats and Nubian Goats

In this study, we measured the growth performance and intramuscular fat (IMF) content of the Longissimus dorsi (LD) of Fuqing goats (FQs) and Nubian goats (NBYs), which exhibit extreme phenotypic differences in terms of their production and meat quality traits. RNA-Seq analysis was performed, and transcriptome data were obtained from the LD tissue of 3-month fetuses (E3), 0-month lambs (0M), 3-month lambs (3M), and 12-month lambs (12M) to reveal the differences in the molecular mechanisms regulating the muscle development and IMF deposition between FQs and NBYs. The results showed that a higher body weight and average daily gain were observed in the NBYs at three developmental stages after birth, whereas a higher IMF content was registered in the FQs at 12M. Additionally, transcriptome profiles during the embryonic period and after birth were completely different for both FQs and NBYs. Moreover, DEGs (KIF23, CCDC69, CCNA2, MKI67, KIF11, RACGAP1, NUSAP1, SKP2, ZBTB18, NES, LOC102180034, CAPN6, TUBA1A, LOC102178700, and PEG10) significantly enriched in the cell cycle (ko04110) at E3 (FQs vs. NBYs), and DEGs (MRPS7, RPS8, RPL6, RPL4, RPS11, RPS10, RPL5, RPS6, RPL8, RPS13, RPS24, RPS15, RPL23) significantly enriched in ribosomes (ko03010) at 0M (FQs vs. NBYs) related to myogenic differentiation and fusion were identified. Meanwhile, the differences in glucose and lipid metabolism began at the E3 timepoint and continued to strengthen as growth proceeded in FQs vs. NBYs. DEGs (CD36, ADIROQR2, ACACA, ACACB, CPT1A, IGF1R, IRS2, LDH-A, PKM, HK2, PFKP, PCK1, GPI, FASN, FADS1, ELOVL6, HADHB, ACOK1, ACAA2, and ACSL4) at 3M (FQs vs. NBYs) and 12M (FQs vs. NBYs) significantly enriched in the AMPK signaling pathway (ko04152), insulin resistance (ko04931), the insulin signaling pathway (ko04910), fatty acid metabolism (ko01212), and glycolysis/gluconeogenesis (ko00010) related to IMF deposition were identified. Further, the results from this study provide the basis for future studies on the mechanisms regulating muscle development and IMF deposition in different breeds of goats, and the candidate genes identified could be used in the selection process.

Anti-inflammation-based treatment of atherosclerosis using Gliclazide-loaded biomimetic nanoghosts

In the study, a biomimetic platform for anti-inflammatory-based treatment of atherosclerotic plaque was developed. Gliclazide (GL) as an anti-inflammasome agent was encapsulated in PLGA nanoparticles (NP), which were coated by monocyte membrane using an extrusion procedure. The size and zeta potential of the nanoghost (NG) changed to 292 and - 10 nm from 189.5 to -34.1 in the core NP. In addition, the actual size of 62.5 nm with a coating layer of 5 nm was measured using TEM. The NG was also showed a sustained release profile with the drug loading content of about 4.7%. Beside to attenuated TNFα, decrease in gene expression levels of NLRP3, MyD88, NOS, IL-1β, IL-18 and caspases 1/3/8/9 in LPS-primed monocytes exposed to NG strongly indicated remarkable inflammation control. After systemic toxicity evaluation and pharmacokinetic analysis of NP and NG, intravenous NG treatment of rabbits with experimentally induced atherosclerosis revealed remarkably less plaque lesions, foam cells, lipid-laden macrophages, and pathological issues in tunica media of aorta sections. Higher expression of CD163 than CD68 in aorta of NG-treated rabbits strongly reveals higher M2/M1 macrophage polarization. The bio/hemocompatible, biomimetic and anti-inflammatory NG can be considered as a potential platform for immunotherapy of particularly atherosclerosis in the field of personalized medicine.

QuilA® adjuvanted Coxevac® sustains Th1-CD8+-type immunity and increases protection in Coxiella burnetii-challenged goats

Coxevac® is the EMA-approved veterinary vaccine for the protection of cattle and goats against Q fever, a zoonotic bacterial disease due to Coxiella burnetii. Since Coxevac® reduces bacterial shedding and clinical symptoms but does not prevent infection, novel, ready-to-use vaccine formulations are needed to increase its immunogenicity. Here, a goat vaccination-challenge model was used to evaluate the impact of the commercially available saponin-based QuilA® adjuvant on Coxevac® immunity. Upon challenge, the QuilA®-Coxevac® group showed a stronger immune response reflected in a higher magnitude of total IgG and an increase in circulating and splenic CD8⁺ T-cells compared to the Coxevac® and challenged-control groups. The QuilA®-Coxevac® group was characterized by a targeted Th1-type response (IFNγ, IP10) associated with increased transcripts of CD8⁺ and NK cells in spleens and γδ T cells in bronchial lymph nodes. Coxevac® vaccinated animals presented an intermediate expression of Th1-related genes, while the challenged-control group showed an immune response characterized by pro-inflammatory (IL1β, TNFα, IL12), Th2 (IL4 and IL13), Th17 (IL17A) and other immunoregulatory cytokines (IL6, IL10). An intriguing role was observed for γδ T cells, which were of TBX21- and SOX4-types in the QuilA®-Coxevac® and challenged control group, respectively. Overall, the addition of QuilA® resulted in a sustained Th1-type activation associated with an increased vaccine-induced bacterial clearance of 33.3% as compared to Coxevac® only. QuilA® could be proposed as a readily-applied veterinary solution to improve Coxevac® efficacy against C. burnetii infection in field settings.

Selection and validation of reference genes for the normalization of quantitative real-time PCR in different muscle tissues of rabbits

Background

In molecular biology studies, the selection of optimal reference genes is of vital importance for accurately quantifying gene expression. The purpose of the present study was to screen the most stable reference genes in different muscle tissues of New Zealand white rabbits and Yufeng yellow rabbits.

Methods and results

Results indicated that the most stable reference genes in the muscle tissues of New Zealand white rabbits were HPRT1, ACTB and PPIC, while HPRT1, PPIC, and RPL13A were the most stable reference genes in muscle tissues of Yufeng yellow rabbits. However, in the longissimus dorsi muscle and the abdominal wall muscle of both varieties, the most stable reference genes were HPRT1, RPL13A, and SDHA. In the quadriceps femoris muscle, the most stable reference genes were ACTB, HPRT1, and SDHA. Furthermore, the relative abundance of MYOG, MYH3 and MSTN was used to confirm the suitability and reliability of the selected most stable reference genes and the most unstable reference gene. Results revealed the same expression patterns of these myogenic genes when normalized according to the most stable genes, while normalization against the unstable reference gene altered the observed expression patterns.

Conclusions

Taken together, our results demonstrated that the most stable reference genes varied among different muscle tissues and different breeds of rabbits. However, HPRT1, PPIC and SDHA presented high stability among all examined reference genes; thus, the combined analysis of HPRT1/ PPIC/ SDHA gene provides the best reference for RT-qPCR in muscle tissues of New Zealand white rabbits and Yufeng yellow rabbits, while HPRT1 is a better choice than other reference genes when using a single reference gene to assess target gene expression. Our results provide basic data for better measuring target gene expression profiles in muscle tissues of rabbits.

Transcriptome-based selection and validation of optimal reference genes in perirenal adipose developing of goat (Capra hircus)

Brown adipose tissue (BAT) is mainly present in young mammals and is important for maintaining body temperature in neonatal mammals because of its ability to produce non-shivering thermogenesis. There is usually a large amount of BAT around the kidneys of newborn kids, but the BAT gradually "whiting" after birth. Screening and validating appropriate reference genes is a prerequisite for further studying the mechanism of goat brown adipose tissue "whiting" during the early stages. In this study, the expression stability of 17 candidate reference genes: 12 COPS8, SAP18, IGF2R, PARL, SNRNP200, ACTG1, CLTA, GANAB, GABARAP, PCBP2, CTSB, and CD151) selected based on previous transcriptome data as new candidate reference genes, 3 (PFDN5, CTNNB1, and EIF3M) recommended in previous studies, and 2 traditional reference genes (ACTB and GAPDH) was evaluated. Real-time quantitative PCR (RT-qPCR) technology was used to detect the expression level of candidate reference genes during goat BAT "whiting". Four algorithms: Normfinder, geNorm, ΔCt method, and BestKeeper, and two comprehensive algorithms: ComprFinder and RefFinder, were used to analyze the stability of each candidate reference genes. GABARAP, CLTA, GAPDH, and ACTB were identified as the most stable reference genes, while CTNNB1, CTSB, and EIF3M were the least stable. Moreover, two randomly selected target genes IDH2 and RBP4, were effectively normalized using the selected most stable reference genes. These findings collectively suggest that GABARAP, CLTA, GAPDH, and ACTB are relatively stable reference genes that can potentially be used for the development of perirenal fat in goats.

Identification of suitable reference genes for normalization of quantitative real-time PCR-based gene expression in chicken (Gallus gallus)

The recent availability of genome information greatly facilitates the fundamental research on chicken. In different organs, gene expression patterns can provide clues to understanding the biological functions. For rapid and accurate quantification of gene expression, quantitative real-time PCR (qPCR) has become one of the most widely used methods. However, the success of qPCR data normalization depends on the use of a suitable reference gene and a single reference gene is not universally suitable for all the experiments. Therefore, reference gene validation is a crucial step for different organ tissues of chicken where suitable reference genes for qPCR analysis in varieties of tissues have not been investigated exhaustively so far. In this study, we have selected 30 Gallus gallus candidate reference genes from NCBI, amplified and studied their expression profiles by qPCR in different organ tissues (breast muscle, thigh muscle, heart, liver, spleen, gizzard, and bursa) of chicken. The result showed that, for breast muscle HSP10 and RPL23, thigh muscle RPL14 and RPL13, liver ALB and HSP70, spleen ALB and GAPDH, heart CYCS and TUBA8B, gizzard RPL5 and 18S rRNA, and bursa EEF1A1 and PGK2 are most stable genes respectively. The results also showed that for different organ tissues, individual or a combination of reference genes should be selected for data normalization. In this study, we have identified and validated 30 reference genes in seven different organ tissues to provide accurate transcript normalization and quantification, which can be useful for gene expression studies in other avian species.

Screening and validation of reference genes for qRT-PCR of bovine skeletal muscle-derived satellite cells

The accuracy of sixteen commonly used internal reference genes was assessed in skeletal muscle-derived satellite cells of Qinchuan cattle at different stages of proliferation and induction of differentiation to determine the most suitable ones. Quantitative real-time PCR and three commonly used algorithmic programs, GeNorm, NormFinder and BestKeeper, were used to evaluate the stability of expression of the candidate internal reference genes (GAPDH, ACTB, PPIA, LRP10, HPRT1, YWHAZ, B2M, TBP, EIF3K , RPS9, UXT, 18S rRNA, RPLP0, MARVELD, EMD and RPS15A) in skeletal muscle-derived satellite cells at 0, 12, 24, 36 and 48 h of growth and after differentiation for 0, 2, 4, 6 and 8 days. The expression of two satellite cell marker genes, CCNA2 and MYF5, was used for validation analysis. The results of the software analyses showed that GAPDH and RPS15A were the most stable reference gene combinations during in vitro proliferation of bovine skeletal muscle-derived satellite cells, RPS15A and RPS9 were the most stable reference gene combinations during in vitro induction of differentiation of the cells, and PPIA was the least stable reference gene during proliferation and differentiation and was not recommended. This study lays the foundation for the selection of reference genes for qRT-PCR during the proliferation and induction of differentiation of bovine skeletal muscle-derived satellite cells.

Genome-Wide Identification of Reference Genes for Reverse-Transcription Quantitative PCR in Goat Rumen

Simple Summary

The rumen plays an essential role as a digestive organ and serves as the primary site of energy substrate absorption for the productive ruminants. Understanding gene expression profiles is necessary to explore the intrinsic regulatory mechanisms of rumen development in goats. The selection of suitable reference genes (RGs) was the primary assay before the real-time quantitative PCR (RT-qPCR). We identified sixteen genome-wide candidate RGs for normalization of gene expression assessments in goat rumen tissues. We demonstrate that the RGs selected (RPS4X and RPS6) were more stably expressed than the commonly used HKGs (ACTB and GAPDH) in goat rumen tissues, suggesting that the ribosomal protein gene family may be another source for the RG pool.

Abstract

As the largest chamber of the ruminant stomach, the rumen not only serves as the principal absorptive surface and nutrient transport pathway from the lumen into the animal, but also plays an important short-chain fatty acid (SCFA) metabolic role in addition to protective functions. Accurate characterization of the gene expression profiles of genes of interest is essential to the exploration of the intrinsic regulatory mechanisms of rumen development in goats. Thus, the selection of suitable reference genes (RGs) is an important prerequisite for real-time quantitative PCR (RT-qPCR). In the present study, 16 candidate RGs were identified from our previous transcriptome sequencing of caprine rumen tissues. The quantitative expressions of the candidate RGs were measured using the RT-qPCR method, and the expression stability of the RGs was assessed using the geNorm, NormFinder, and BestKeeper programs. GeNorm analysis showed that the M values were less than 0.5 for all the RGs except GAPT4, indicating that they were stably expressed in the rumen tissues throughout development. RPS4X and RPS6 were the two most stable RGs. Furthermore, the expressions of two randomly selected target genes (IGF1 and TOP2A), normalized by the selected most stable RGs (RPS4X and RPS6), were consistent with the results of RNA sequencing, while the use of GAPDH and ACTB as RGs resulted in altered profiles. Overall, RPS4X and RPS6 showed the highest expression stability and the lowest coefficients of variation, and could be used as the optimal reference combination for quantifying gene expression in rumen tissues via RT-qPCR analysis.

Evaluation of candidate reference genes for quantitative real-time PCR analysis in a male rat model of dietary iron deficiency

BACKGROUND

Quantitative real-time polymerase chain reaction (qPCR) is a reliable and efficient method for quantitation of gene expression. Due to the increased use of qPCR in examining nutrient-gene interactions, it is important to examine, develop, and utilize standardized approaches for data analyses and interpretation. A common method used to normalize expression data involves the use of reference genes (RG) to determine relative mRNA abundance. When calculating the relative abundance, the selection of RG can influence experimental results and has the potential to skew data interpretation. Although common RG may be used for normalization, often little consideration is given to the suitability of RG selection for an experimental condition or between various tissue or cell types. In the current study, we examined the stability of gene expression using BestKeeper, comparative delta quantitation cycle, NormFinder, and RefFinder in a variety of tissues obtained from iron-deficient and pair-fed iron-replete rats to determine the optimal selection among ten candidate RG.

RESULTS

Our results suggest that several commonly used RG (e.g., Actb and Gapdh) exhibit less stability compared to other candidate RG (e.g., Rpl19 and Rps29) in both iron-deficient and iron-replete pair-fed conditions. For all evaluated RG, Tfrc expression significantly increased in iron-deficient animal livers compared to the iron-replete pair-fed controls; however, the relative induction varied nearly 4-fold between the most suitable (Rpl19) and least suitable (Gapdh) RG.

CONCLUSION

These results indicate the selection and use of RG should be empirically determined and RG selection may vary across experimental conditions and biological tissues.

Selection of reference genes for quantitative real-time PCR normalization in European quail tissues

Coturniculture has been standing out as an industrial poultry activity in several countries around the world because of the several adaptive advantages of quails. Research that considers the analysis of gene expression can enhance this activity. This study aimed to analyze the stability of reference genes (RGs) in different tissues of quails (both males and females) for the recommendation of use in gene expression studies by the quantitative reverse transcription-polymerase chain reaction (RT-qPCR). The expression stability of ten RGs (ACTA1, ACTB, B2M, GAPDH, HMBS, SDHA, HPRT1, MRPS27, MRPS30, and RPL5) was analyzed in four tissues (breast muscle, abdominal fat, liver, and intestine), and assessed using the statistical tools geNorm, NormFinder, comparative ΔC_q method, and BestKeeper. The HPRT1 gene was the most stable in all quail tissues tested, followed by MRPS27 and MRPS30 in breast muscle, B2M and RPL5 in abdominal fat, HMBS and B2M in the liver, and RPL5 and HMBS in the intestine. These results may help studies using RT-qPCR assays to assess quail tissues from both sexes because they provide data on the most stable genes, which should be tested as candidate RGs for other experimental conditions.

Understanding seasonal weight loss tolerance in dairy goats: a transcriptomics approach

Background

Seasonal weight loss (SWL) is a very important limitation to the production of ruminants in the Mediterranean and Tropical regions. In these areas, long dry seasons lead to poor pastures with low nutritional value. During the dry season, ruminants, particularly those raised in extensive production systems, lose around 30% of their body weight. Seasonal weight loss has important consequences on animal productive performance and health. In this study, RNA sequencing was used to characterize feed restriction effects in dairy goat of 2 breeds with different SWL tolerance: Majorera (tolerant) and Palmera (susceptible). Nine Majorera and ten Palmera goats were randomly distributed in a control and a restricted group: Majorera Control (adequately fed; MC; n = 4), Palmera Control (adequately fed; PC; n = 6), Majorera Restricted (feed restricted; ME; n = 5) and Palmera Restricted (feed restricted; PE; n = 4). On day 22 of the trial, mammary gland biopsies were collected for transcriptomics analysis.

Results

From these samples, 24,260 unique transcripts were identified. From those, 82 transcripts were differentially expressed between MC and ME, 99 between PC and PE, twelve between both control groups and twenty-nine between both restricted groups.

Conclusions

Feed restriction affected several biochemical pathways in both breeds such as: carbohydrate and lipid transport; intracellular trafficking, RNA processing and signal transduction. This research also highlights the importance or involvement of the genes in tolerance (ENPP1, S-LZ, MT2A and GPNB) and susceptibility (GPD1, CTPS1, ELOVL6 and NR4A1) to SWL with respectively higher expression in the Majorera restriced group and the Palmera restricted group in comparison to the control groups. In addition, results from the study may be extrapolated to other dairy ruminant species.

Using RNA-Seq to Identify Reference Genes of the Transition from Brown to White Adipose Tissue in Goats

Simple Summary

Brown adipose tissue (BAT) plays important roles in unique non-shivering thermogenesis. It is necessary to select reference genes during the transition process from brown (BAT) to white adipose tissue (WAT) for Quantitative PCR (qPCR) analysis. In this study, CTNNB, PFDN5 and EIF3M, selected from RNA sequencing data, were the most suitable reference genes. The present study provides a detailed analysis of the expression stability of reference genes for the study of gene expression profiling during the transition process from BAT to WAT.

Abstract

Brown adipose tissues have unique non-shivering thermogenesis functions, can be found in newborn ruminate animals, and then are gradually replaced by white adipose tissues in adulthood. For the purpose of exploring the intrinsic mechanism underlying the conversion process from brown (BAT) to white adipose tissue (WAT), it is necessary to utilize Quantitative PCR (qPCR) to study gene expression profiling. In this study, we identified reference genes that were consistently expressed during the transformation from goat BAT to WAT using RNA-seq data. Then, twelve genes were evaluated as candidate reference genes for qPCR in goat perirenal adipose tissue using three tools (geNorm, Normfinder, and BestKeeper). In addition, the selected reference genes were used to normalize the gene expression of PGC-1α and GPAT4. It was found that traditional reference genes, such as GAPDH, RPLP0, HPRT1, and PPIA were not suitable for target gene normalization. In contrast, CTNNB, PFDN5, and EIF3M, selected from RNA sequencing data, showed the least variation and were recommended as the best reference genes during the transformation from BAT to WAT.

Influence of heat stress on reference genes stability in heart and liver of two chickens genotypes

Introduction

Real-time polymerase chain reaction (RT-qPCR) is an important tool for analyzing gene expression. However, before analyzing the expression of target genes, it is crucial to normalize the reference genes, in order to find the most stable gene to be used as an endogenous control. A gene that remains stable in all samples under different treatments is considered a suitable normalizer. In this sense, we aimed to identify stable reference genes for normalization of target genes in the heart and liver tissues from two genetically divergent groups of chickens (Cobb 500^® commercial line and Peloco backyard chickens) under comfort and acute heat stress environmental conditions. Eight reference genes (ACTB, HPRT1, RPL5, EEF1, MRPS27, MRPS30, TFRC and LDHA) were analyzed for expression stability. The samples were obtained from 24 chickens, 12 from the backyard Peloco and 12 from the Cobb 500® line, exposed to two environmental conditions (comfort and heat stress). Comfort temperature was 23°C and heat stress temperature was 39.5°C for one hour. Subsequently, the animals were euthanized, and heart and liver tissue fragments were collected for RNA extraction and amplification. To determine the stability rate of gene expression, three different statistical algorithms were applied: BestKeeper, geNorm and NormFinder, and to obtain an aggregated stability list, the RankAgregg package of R software was used.

Results

The most stable genes using BestKeeper tool, including the two factors (genetic group and environmental condition), were LDHA, RPL5 and MRPS27 for heart tissue, and TFRC, RPL5 and EEF1 for liver tissue. Applying geNorm algorithm, the best reference genes were RPL5, EEF1 and MRPS30 for heart tissue and LDHA, EEF1 and RPL5 for liver. Using the NormFinder algorithm, the best normalizer genes were EEF1, RPL5 and LDHA in heart, and EEF1, RPL5 and ACTB in liver tissue. In the overall ranking obtained by RankAggreg package, considering the three algorithms, the RPL5, EEF1 and LDHA genes were the most stable for heart tissue, whereas RPL5, EEF1 and ACTB were the most stable for liver tissue.

Conclusion

According to the RankAggreg tool classification based on the three different algorithms (BestKeeper, geNorm and NormFinder), the most stable genes were RPL5, EEF1 and LDHA for heart tissue and RPL5, EEF1 and ACTB for liver tissue of chickens subjected to comfort and acute heat stress environmental conditions. However, the best reference genes may vary depending on the experimental conditions of each study, such as different breeds, environmental stressors, and tissues analyzed. Therefore, the need to perform priori studies to assay the best reference genes at the outset of each study is emphasized.

Selection of suitable reference genes for gene expression studies in myxosporean (Myxozoa, Cnidaria) parasites

Myxozoans (Cnidaria: Myxozoa) are an extremely diversified group of endoparasites some of which are causative agents of serious diseases in fish. New methods involving gene expression studies have emerged over the last years to better understand and control myxozoan diseases. Quantitative RT-PCR is the most extensively used approach for gene expression studies. However, the accuracy of the results depends on the normalization of the data to reference genes. We studied the expression of eight commonly used reference genes, adenosylhomocysteinase (AHC1), beta actin (ACTB), eukaryotic translation elongation factor 2 (EF2), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), hypoxanthine-guanine phosphoribosyltransferase 1 (HPRT1), DNA-directed RNA polymerase II (RPB2), 18S ribosomal RNA (18S), 28S ribosomal RNA (28S) across different developmental stages of three myxozoan species, Sphaerospora molnari, Myxobolus cerebralis and Ceratonova shasta, representing the three major myxozoan linages from the largest class Myxosporea. The stable reference genes were identified using four algorithms: geNorm, NormFinder, Bestkeeper and ΔCq method. Additionally, we analyzed transcriptomic data from S. molnari proliferative and spore-forming stages to compare the relative amount of expressed transcripts with the most stable reference genes suggested by RT-qPCR. Our results revealed that GAPDH and EF2 are the most uniformly expressed genes across the different developmental stages of the studied myxozoan species.

Selection and validation of reference genes for quantitative Real-Time PCR in Arabis alpina

Arabis alpina is a perennial arctic-alpine plant and an upcoming model organism for genetics and molecular biology for the Brassicaceae family. One essential method for most molecular approaches is the analysis of gene expression by reverse-transcription quantitative Real-Time PCR (RT-qPCR). For the normalisation of expression data in RT-qPCR experiments, it is essential to use reliable reference genes that are not affected under a wide range of conditions. In this study we establish a set of 15 A. alpina reference genes that were tested under different conditions including cold, drought, heat, salt and gibberellic acid treatments. Data analyses with geNORM, BestKeeper and NormFinder revealed the most stable reference genes for the tested conditions: RAN3, HCF and PSB33 are most suitable for cold treatments; UBQ10 and TUA5 for drought; RAN3, PSB33 and EIF4a for heat; CAC, TUA5, ACTIN 2 and PSB33 for salt and PSB33 and TUA5 for gibberellic acid treatments. CAC and ACTIN 2 showed the least variation over all tested samples. In addition, we show that two reference genes are sufficient to normalize RT-qPCR data under our treatment conditions. In future studies, these reference genes can be used for an adequate normalisation and thus help to generate high quality RT-qPCR data in A. alpina.

Selection and validation of suitable reference genes for qPCR gene expression analysis in goats and sheep under Peste des petits ruminants virus (PPRV), lineage IV infection

Identification of suitable candidate reference genes is an important prerequisite for validating the gene expression data obtained from downstream analysis of RNA sequencing using quantitative real time PCR (qRT-PCR). Though existence of a universal reference gene is myth, commonly used reference genes can be assessed for expression stability to confer their suitability to be used as candidate reference genes in gene expression studies. In this study, we evaluated the expression stability of ten most commonly used reference genes (GAPDH, ACTB, HSP90, HMBS, 18S rRNA, B2M, POLR2A, HPRT1, ACAC, YWHAZ) in fourteen different Peste des petits ruminants virus (PPRV) infected tissues of goats and sheep. RefFinder and RankAggreg software were used to deduce comprehensive ranking of reference genes. Our results suggested HMBS and B2M in goats and HMBS and HPRT1 in sheep can be used as suitable endogenous controls in gene expression studies of PPRV infection irrespective of tissues and condition as a whole, thus eliminating the use of tissue specific/ condition specific endogenous controls. We report for the first time suitable reference genes for gene expression studies in PPRV infected tissues. The reference genes determined here can be useful for future studies on gene expression in sheep and goat infected with PPRV, thus saving extra efforts and time of repeating the reference gene determination and validation.

Fibroblast growth factor 21 induces lipolysis more efficiently than it suppresses lipogenesis in goat adipocytes

Fibroblast growth factor 21 (FGF21) potentially regulates glucose and lipid metabolism in energy homeostasis. We investigated dynamic changes in goat adipocytes treated with 75 nM FGF21 for 24, 36 and 48 h. Compared to controls, FGF21-treated adipocytes displayed smaller lipid droplets and altered levels of the mRNA transcripts encoding several lipolysis genes. The genes with elevated mRNA levels included: ATGL, HSL, CPT-1, and UCP1, and this was observed mainly at 24 and 36 h (P < 0.05). Some gene expression was attenuated including lipogenesis genes, such as SREBP1, PPARγ, C/EBPα, and ACC. This attenuation was observed mainly at 24 h (P < 0.05). Among the genes that were significantly induced or inhibited, ATGL, PGC1α, and C/EBPα were observed a significant effect at 48 h (P < 0.05). In addition, FGF21 treatment greatly increased number of mitochondria and the expression of genes implicated in mitochondrial biogenesis, such as PGC1α, NRF1, and TFAM. These results suggest that FGF21 treatment induced lipolysis more effectively than it suppressed lipogenesis in goat adipocytes, and that mitochondrial biogenesis plays an important role in these cells.

Validation of Common Housekeeping Genes as Reference for qPCR Gene Expression Analysis During iPS Reprogramming Process

Induced pluripotent stem cell (iPS) reprogramming allows to turn a differentiated somatic cell into a pluripotent cell. This process is accompanied by many changes in fundamental cell properties, such as energy production, cell-to-cell interactions, cytoskeletal organization, and others. Real-time quantitative polymerase chain reaction (RT-qPCR) can be used as a quantitative method of gene expression analysis to investigate iPS reprogramming but it requires a validation of reference genes for the accurate assessment of target genes’ expression. Currently, studies evaluating the performance of reference genes during iPS reprogramming are lacking. In this study we analysed the stability of 12 housekeeping genes during 20 days of iPS reprogramming of murine cells based on statistical analyses of RT-qPCR data using five different statistical algorithms. This study reports strong variations in housekeeping gene stability during the reprogramming process. Most stable genes were Atp5f1, Pgk1 and Gapdh, while the least stable genes were Rps18, Hprt, Tbp and Actb. The results were validated by a proof-of-point qPCR experiment with pluripotent markers Nanog, Rex1 and Oct4 normalized to the best and the worst reference gene identified by the analyses. Overall, this study and its implications are particularly relevant to investigations on the cell-state and pluripotency in iPS reprogramming.

Influence of heat stress, sex and genetic groups on reference genes stability in muscle tissue of chicken

Quantitative RT-PCR is an important technique for assessing gene expression. However, a proper normalization of reference genes prior to expression analyses of target genes is necessary. The best normalizer is that gene which remains stable in all samples from different treatments. The aim of this study was to identify stable reference genes for normalization of target genes in muscle tissue from three genetically divergent chickens groups (Peloco, Cobb 500^® and Caneluda) under environmental (heat stress and comfort) and sex influence. Expressions of ten reference genes were tested for stability in breast muscular tissue (Pectoralis major muscle). Samples were obtained from 36 males and females of two backyard breeds (Caneluda and Peloco) and one commercial line (Cobb 500^®) under two environments. The heat stress and comfort temperature were 39 and 23°C, respectively. Animals were housed in the Animal Science Department at Universidade Estadual do Sudoeste da Bahia. We analyzed the expression data by four statistical tools (SLqPCR, NormFinder, Bestkeeper and Comparative CT). According to these tools, genes stability varied according to sex, genetic group and environment, however, some genes remained stable in all analyzes. There was no difference between the most stable genes for sex effect, being MRPS27 more stable for both males and females. In general, MRPS27 was the most stable gene. Within the three genetic groups, the most stable genes were RPL5, HMBS and EEF1 to Cobb 500^®, Peloco and Caneluda, respectively. Within the environment, the most stable gene under comfort and heat stress conditions was HMBS and MRPS27, respectively. BestKeeper and Comparative Ct were less correlated (28%) and SLqPCR and NormFinder were the most correlated (98%). MRPS27, RPL5 and MRPS30 genes were considered stable according the overall ranking and can be used as normalizer of relative expression of target genes in muscle tissue of chickens under heat stress.

Validation of reference genes for cryopreservation studies with the gorgonian coral endosymbiont Symbiodinium

Quantification by real-time RT-PCR requires a stable internal reference known as a housekeeping gene (HKG) for normalising the mRNA levels of target genes. The present study identified and validated stably expressed HKGs in post-thaw Symbiodinium clade G. Six potential HKGs, namely, pcna, gapdh, 18S rRNA, hsp90, rbcl, and ps1, were analysed using three different algorithms, namely, GeNorm, NormFinder, and BestKeeper. The GeNorm algorithm ranked the candidate genes as follows in the order of decreasing stability: pcna and gapdh > ps1 > 18S rRNA > hsp90 > rbcl. Results obtained using the NormFinder algorithm also showed that pcna was the most stable HKG and ps1 was the second most stable HKG. We found that the candidate HKGs examined in this study showed variable stability with respect to the three algorithms. These results indicated that both pcna and ps1 were suitable for normalising target gene expression determined by performing real-time RT-PCR in cryopreservation studies on Symbiodinium clade G. The results of the present study would help future studies to elucidate the effect of cryopreservation on gene expression in dinoflagellates.

The 60S ribosomal protein L13 is the most preferable reference gene to investigate gene expression in selected organs from turkeys and chickens, in context of different infection models

Evaluation of reference genes for expression studies in chickens and turkeys is very much limited and unavailable for various infectious models. In this study, eight candidate reference genes HMBS, HPRT1, TBP, VIM, TFRC, RPLP0, RPL13 and RPS7 were evaluated by five different algorithms (GeNorm, NormFinder, BestKeeper©, delta CT, RefFinder) to assess their stability. In order to analyze a broad variation of tissues, spleen, liver, caecum and caecal tonsil of different aged specific pathogen free (SPF) layer chickens and commercial turkeys, uninfected or infected with the extracellular pathogen Histomonas meleagridis, were included. For tissue samples from SPF chickens RPL13 and TBP were found to be the most stable reference genes. Further testing of RPL13 and TBP in the same organs of uninfected and infected SPF broiler chickens with the intracellular pathogen fowl aviadenovirus confirmed this finding. In tissue samples from turkeys, a stable expression of RPL13 and TFRC genes was noticed. Overall, the determined reference genes should be considered whenever gene expression studies in spleen, liver, caecum and caecal tonsil of chickens and turkeys are performed.

Comparative Transcriptome Analysis Reveals that a Ubiquitin-Mediated Proteolysis Pathway Is Important for Primary and Secondary Hair Follicle Development in Cashmere Goats

Background

The fleece of cashmere goats contains two distinct populations of fibers, a short and fine non-medullated insulating cashmere fiber and a long and coarse medullated guard hair. The former is produced by secondary follicles (SFs) and the later by primary follicles (PFs). Evidence suggests that the induction of PFs and SFs may require different signaling pathways. The regulation of BMP2/4 signaling by noggin and Edar signaling via Downless genes are essential for the induction of SFs and PFs, respectively. However, these differently expressed genes of the signaling pathway cannot directly distinguish between the PFs and SFs.

Results

In this study, we selected RNA samples from 11 PFs and 7 SFs that included 145,525 exons. The pathway analysis of 4512 differentially expressed exons revealed that the most statistically significant metabolic pathway was related to the ubiquitin–mediated proteolysis pathway (UMPP) (P<3.32x 10⁻⁷). In addition, the 51 exons of the UMPP that were differentially expressed between the different types of hair follicle (HFs) were compared by cluster analysis. This resulted in the PFs and SFs being divided into two classes. The expression level of two selected exons was analyzed by qRT-PCR, and the results indicated that the expression patterns were consistent with the deep sequencing results obtained by RNA-Seq.

Conclusions

Based on the comparative transcriptome analysis of 18 HFs from cashmere goats, a large number of differentially expressed exons were identified using a high-throughput sequencing approach. This study suggests that UMPP activation is a prominent signaling pathway for distinguishing the PFs and SFs of cashmere goats. It is also a meaningful contribution to the theoretical basis of the biological study of the HFs of cashmere goats and other mammals.

Identification of Reference Genes for Quantitative Gene Expression Studies in a Non-Model Tree Pistachio (Pistacia vera L.)

The tree species, Pistacia vera (P. vera) is an important commercial product that is salt-tolerant and long-lived, with a possible lifespan of over one thousand years. Gene expression analysis is an efficient method to explore the possible regulatory mechanisms underlying these characteristics. Therefore, having the most suitable set of reference genes is required for transcript level normalization under different conditions in P. vera. In the present study, we selected eight widely used reference genes, ACT, EF1α, α-TUB, β-TUB, GAPDH, CYP2, UBQ10, and 18S rRNA. Using qRT-PCR their expression was assessed in 54 different samples of three cultivars of P. vera. The samples were collected from different organs under various abiotic treatments (cold, drought, and salt) across three time points. Several statistical programs (geNorm, NormFinder, and BestKeeper) were applied to estimate the expression stability of candidate reference genes. Results obtained from the statistical analysis were then exposed to Rank aggregation package to generate a consensus gene rank. Based on our results, EF1α was found to be the superior reference gene in all samples under all abiotic treatments. In addition to EF1α, ACT and β-TUB were the second best reference genes for gene expression analysis in leaf and root. We recommended β-TUB as the second most stable gene for samples under the cold and drought treatments, while ACT holds the same position in samples analyzed under salt treatment. This report will benefit future research on the expression profiling of P. vera and other members of the Anacardiaceae family.

Selection of Reliable Reference Genes for Gene Expression Analysis under Abiotic Stresses in the Desert Biomass Willow, Salix psammophila

Salix psammophila is a desert shrub willow that has extraordinary adaptation to abiotic stresses and plays an important role in maintaining local ecosystems. Moreover, S. psammophila is regarded as a promising biomass feedstock because of its high biomass yields and short rotation coppice cycle. However, few suitable reference genes (RGs) for quantitative real-time polymerase chain reaction (qRT-PCR) constrain the study on normalization of gene expression in S. psammophila until now. Here, we investigated the expression stabilities of 14 candidate RGs across tissue types and under four abiotic stress treatments, including heat, cold, salt, and drought treatments. After calculation of PCR efficiencies, three different software, NormFinder, geNorm, and BestKeeper were employed to analyze systematically the qRT-PCR data, and the outputs were merged by RankAggreg software. The optimal RGs selected for gene expression analysis were EF1α (Elongation factor-1 alpha) and OTU (OTU-like cysteine protease family protein) for different tissue types, UBC (Ubiquitin-conjugating enzyme E2) and LTA4H (Leukotriene A-4 hydrolase homolog) for heat treatment, HIS (Histone superfamily protein H3) and ARF2 (ADP-ribosylation factor 2) for cold treatment, OTU and ACT7 (Actin 7) for salt treatment, UBC and LTA4H for drought treatment. The expression of UBC, ARF2, and VHAC (V-type proton ATPase subunit C) varied the least across tissue types and under abiotic stresses. Furthermore, the relative genes expression profiles of one tissue-specific gene WOX1a (WUSCHEL-related homeobox 1a), and four stress-inducible genes, including Hsf-A2 (Heat shock transcription factors A2), CBF3 (C-repeat binding factor 3), HKT1 (High-Affinity K⁺ Transporter 1), and GST (Glutathione S-transferase), were conducted to confirm the validity of the RGs in this study. These results provided an important RGs application guideline for gene expression characterization in S. psammophila.

Transcription analysis of genes involved in lipid metabolism reveals the role of chromium in reducing body fat in animal models

Chromium was proposed to be an essential trace element over 50 years ago and has been accepted as an essential element for over 30 years. The recent studies indicated that the addition of supra nutritional amounts of chromium to the diet can only be considered as having pharmacological effects. However, the precise mechanism through which chromium acts on lipid, carbohydrate, protein and nucleic acid metabolism are relatively poor studied. To uncover, at least partially, the role of chromium in lipid metabolism, in this study, we evaluated the expression status of eight important genes, involved in fat biosynthesis and lipid metabolism, in four different tissue types (liver, subcutaneous fat, visceral fat, and longissimus muscle) in domestic goat kids feeding on three different chromium levels. The quantitative real-time PCR (RT-PCR) was established for expression analyses with HSP90 gene was used as reference gene. The results showed that supplementation of goats with 1.5mg/day chromium significantly decreases the expression of the ACC1, DGAT1, FABP4, FAS, HSL, LEP genes, but does not affect the expression of the LPL and SCD1 genes in all studied tissues. This study highlights, for the first time, the role of supra nutritional levels of chromium in lipid biosynthesis and metabolism. These findings are of especial importance for improving meat quality in domestic animals.

Identification of Suitable Reference Genes for Real Time Quantitative Polymerase Chain Reaction Assays on Pectoralis major Muscle in Chicken (Gallus gallus )

Thirteen reference genes were investigated to determine their stability to be used as a housekeeping in gene expression studies in skeletal muscle of chickens. Five different algorithms were used for ranking of reference genes and results suggested that individual rankings of the genes differed among them. The stability of the expression of reference genes were validated using samples obtained from the Pectoralis major muscle in chicken. Samples were obtained from chickens in different development periods post hatch and under different nutritional diets. For gene expression calculation the ΔΔCt approach was applied to compare relative expression of pairs of genes within each of 52 samples when normalized to mitochondrially encoded cytochrome c oxidase II (MT-CO2) target gene. Our findings showed that hydroxymethylbilane synthase (HMBS) and hypoxanthine phosphoribosyl transferase 1 (HPRT1) are the most stable reference genes while transferrin receptor (TFRC) and beta-2-microglobulin (B2M) ranked as the least stable genes in the Pectoralis major muscle of chickens. Moreover, our results revealed that HMBS and HPRT1 gene expression did not change due to dietary variations and thus it is recommended for accurate normalization of RT-qPCR data in chicken Pectoralis major muscle.

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.