With Reference to Reference Genes: A Systematic Review of Endogenous Controls in Gene Expression Studies

RPLP0/TBP are the most stable reference genes for human dental pulp stem cells under osteogenic differentiation

BACKGROUND

Validation of the reference gene (RG) stability during experimental analyses is essential for correct quantitative real-time polymerase chain reaction (RT-qPCR) data normalisation. Commonly, in an unreliable way, several studies use genes involved in essential cellular functions [glyceraldehyde-3-phosphate dehydrogenase (GAPDH), 18S rRNA, and β-actin] without paying attention to whether they are suitable for such experimental conditions or the reason for choosing such genes. Furthermore, such studies use only one gene when Minimum Information for Publication of Quantitative Real-Time PCR Experiments guidelines recommend two or more genes. It impacts the credibility of these studies and causes distortions in the gene expression findings. For tissue engineering, the accuracy of gene expression drives the best experimental or therapeutical approaches.

AIM

To verify the most stable RG during osteogenic differentiation of human dental pulp stem cells (DPSCs) by RT-qPCR.

METHODS

We cultivated DPSCs under two conditions: Undifferentiated and osteogenic differentiation, both for 35 d. We evaluated the gene expression of 10 candidates for RGs [ribosomal protein, large, P0 (RPLP0), TATA-binding protein (TBP), GAPDH, actin beta (ACTB), tubulin (TUB), aminolevulinic acid synthase 1 (ALAS1), tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta (YWHAZ), eukaryotic translational elongation factor 1 alpha (EF1a), succinate dehydrogenase complex, subunit A, flavoprotein (SDHA), and beta-2-microglobulin (B2M)] every 7 d (1, 7, 14, 21, 28, and 35 d) by RT-qPCR. The data were analysed by the four main algorithms, ΔCt method, geNorm, NormFinder, and BestKeeper and ranked by the RefFinder method. We subdivided the samples into eight subgroups.

RESULTS

All of the data sets from clonogenic and osteogenic samples were analysed using the RefFinder algorithm. The final ranking showed RPLP0/TBP as the two most stable RGs and TUB/B2M as the two least stable RGs. Either the ΔCt method or NormFinder analysis showed TBP/RPLP0 as the two most stable genes. However, geNorm analysis showed RPLP0/EF1α in the first place. These algorithms’ two least stable RGs were B2M/GAPDH. For BestKeeper, ALAS1 was ranked as the most stable RG, and SDHA as the least stable RG. The pair RPLP0/TBP was detected in most subgroups as the most stable RGs, following the RefFinfer ranking.

CONCLUSION

For the first time, we show that RPLP0/TBP are the most stable RGs, whereas TUB/B2M are unstable RGs for long-term osteogenic differentiation of human DPSCs in traditional monolayers.

RPLP0/TBP are the most stable reference genes for human dental pulp stem cells under osteogenic differentiation

BACKGROUND

Validation of the reference gene (RG) stability during experimental analyses is essential for correct quantitative real-time polymerase chain reaction (RT-qPCR) data normalisation. Commonly, in an unreliable way, several studies use genes involved in essential cellular functions [glyceraldehyde-3-phosphate dehydrogenase (GAPDH), 18S rRNA, and β-actin] without paying attention to whether they are suitable for such experimental conditions or the reason for choosing such genes. Furthermore, such studies use only one gene when Minimum Information for Publication of Quantitative Real-Time PCR Experiments guidelines recommend two or more genes. It impacts the credibility of these studies and causes distortions in the gene expression findings. For tissue engineering, the accuracy of gene expression drives the best experimental or therapeutical approaches.

AIM

To verify the most stable RG during osteogenic differentiation of human dental pulp stem cells (DPSCs) by RT-qPCR.

METHODS

We cultivated DPSCs under two conditions: Undifferentiated and osteogenic differentiation, both for 35 d. We evaluated the gene expression of 10 candidates for RGs [ribosomal protein, large, P0 (RPLP0), TATA-binding protein (TBP), GAPDH, actin beta (ACTB), tubulin (TUB), aminolevulinic acid synthase 1 (ALAS1), tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta (YWHAZ), eukaryotic translational elongation factor 1 alpha (EF1a), succinate dehydrogenase complex, subunit A, flavoprotein (SDHA), and beta-2-microglobulin (B2M)] every 7 d (1, 7, 14, 21, 28, and 35 d) by RT-qPCR. The data were analysed by the four main algorithms, ΔCt method, geNorm, NormFinder, and BestKeeper and ranked by the RefFinder method. We subdivided the samples into eight subgroups.

RESULTS

All of the data sets from clonogenic and osteogenic samples were analysed using the RefFinder algorithm. The final ranking showed RPLP0/TBP as the two most stable RGs and TUB/B2M as the two least stable RGs. Either the ΔCt method or NormFinder analysis showed TBP/RPLP0 as the two most stable genes. However, geNorm analysis showed RPLP0/EF1α in the first place. These algorithms’ two least stable RGs were B2M/GAPDH. For BestKeeper, ALAS1 was ranked as the most stable RG, and SDHA as the least stable RG. The pair RPLP0/TBP was detected in most subgroups as the most stable RGs, following the RefFinfer ranking.

CONCLUSION

For the first time, we show that RPLP0/TBP are the most stable RGs, whereas TUB/B2M are unstable RGs for long-term osteogenic differentiation of human DPSCs in traditional monolayers.

Identification of suitable reference genes for RT-qPCR studies in human parathyroid tissue glandular cells

Identifying a proper reference gene allows us to understand fundamental changes in many biological processes. Normalization during gene expression analyses is essential for every tissue/cell type, including parathyroid tissue glandular cells. Quantitative method of gene expression analyses via qRT-PCR method provides the accurate examination of every target gene. There are limited reports to present commonly used reference genes in human parathyroid tissues rather than for glandular cell types. This study aims to determine and compare the most stable to least stable genes for parathyroid tissue cells. 43 human parathyroid tissue obtained from primary and secondary hyperparathyroidism patients and glandular cells isolated enzymatically by the removal of extracellular matrix components. After extraction of the total RNA, cDNA synthesis was performed, then qRT-PCR evaluated 14 candidate reference genes. Stability was determined by RefFinder software (Delta ct, BestKeeper, Genorm, and NormFinder algorithms), and the outcome was evaluated for five groups. Even if assessed with different groups, the most stable genes were RPLP0 and GAPDH, while the CLTC and RNA 18S were the least stable. We have confirmed the comprehensive ranking of the most stable three genes alone with the NormFinder algorithm to understand intergroup variation and found out that RPLP0>GAPDH>PGK1. Lastly, comparisons of relative target gene (GCM2) expression revealed similar expression patterns for the most stable reference genes. The most stable reference gene is recommended for the stages where stability is evaluated using the results of four different approaches using RefFinder. We aspire for this study to assist future research to conduct thorough assessments of appropriate reference genes before engaging in gene expression analyses for parathyroid tissue.

Impact of culture medium on the interpretation of qRT-PCR data in HepG2 incubated with lactobacilli

Recently, an increasing number of studies have investigated the mechanism of action of lactobacilli in the treatment of non-alcoholic fatty liver disease. Using four computational tools (NormFinder, geNorm, Delta Ct, and BestKeeper), six potential reference genes (RGs) were analyzed in the human liver cell line HepG2 cultivated 24 h in the presence of two strains of heat-killed lactobacilli, Limosilactobacillus reuteri E and Lactiplantibacillus plantarum KG4, respectively, in different cultivation media [Dulbecco´s Modified Eagle´s Medium (DMEM) high glucose or Roswell Park Memorial Institute (RPMI)]. The analysis revealed that the suitability of RG was similar between the two lactobacilli but quite different between the two media. The commonly used RGs, 18S rRNA and glyceraldehyde-3-phosphate dehydrogenase were the most unstable in DMEM high glucose. Normalization of the mRNA expression of the target gene encoding sterol regulatory element-binding protein 1c (SREBP-1c) to different RGs resulted in different expression profiles. This demonstrates that validation of candidate RGs under specific experimental conditions is crucial for the correct interpretation of quantitative polymerase chain reaction data. In addition, the choice of media has a profound impact on the effect of lactobacilli on lipogenesis at the gene expression level, as shown by the transcription factor SREBP-1c.

Characterization and modulation of the pro-inflammatory effects of immune cells in the canine intervertebral disk

Background

Intervertebral disk (IVD) degeneration affects both humans and canines and is a major cause of low back pain (LBP). Mast cell (MC) and macrophage (MØ) infiltration has been identified in the pathogenesis of IVD degeneration (IVDD) in the human and rodent model but remains understudied in the canine. MC degranulation in the IVD leads to a pro-inflammatory cascade and activates protease activated receptor 2 (PAR2) on IVD cells. The objectives of the present study are to: (1) highlight the pathophysiological changes observed in the degenerate canine IVD, (2) further characterize the inflammatory effect of MCs co-cultured with canine nucleus pulposus (NP) cells, (3) evaluate the effect of construct stiffness on NP and MCs, and (4) identify potential therapeutics to mitigate pathologic changes in the IVD microenvironment.

Methods

Canine IVD tissue was isolated from healthy autopsy research dogs (beagle) and pet dogs undergoing laminectomy for IVD herniation. Morphology, protein content, and inflammatory markers were assessed. NP cells isolated from healthy autopsy (Mongrel hounds) tissue were co-cultured with canine MCs within agarose constructs and treated with cromolyn sodium (CS) and PAR2 antagonist (PAR2A). Gene expression, sulfated glycosaminoglycan content, and stiffness of constructs were assessed.

Results

CD 31+ blood vessels, mast cell tryptase, and macrophage CD 163+ were increased in the degenerate surgical canine tissue compared to healthy autopsy. Pro-inflammatory genes were upregulated when canine NP cells were co-cultured with MCs and the stiffer microenvironment enhanced these effects. Treatment with CS and PAR2 inhibitors mediated key pro-inflammatory markers in canine NP cells.

Conclusion

There is increased MC, MØs, and vascular ingrowth in the degenerate canine IVD tissue, similar to observations in the clinical population with IVDD and LBP. MCs co-cultured with canine NP cells drive inflammation, and CS and PAR2A are potential therapeutics that may mitigate the pathophysiology of IVDD in vitro.

Cloning, Expression Characteristics of Farnesyl Pyrophosphate Synthase Gene from Platycodon grandiflorus and Functional Identification in Triterpenoid Synthesis

Platycodon grandiflorus is a medicinal plant whose main component is platycodins, which have a variety of pharmacological effects and nutritional values. The farnesyl pyrophosphate synthase (FPS) is a key enzyme in the isoprenoid biosynthesis pathway, which catalyzes the synthesis of farnesyl diphosphate (FPP). In this study, we cloned the FPS gene from P. grandiflorus (PgFPS) with an ORF of 1260 bp, encoding 419 amino acids with a deduced molecular weight and theoretical pI of 46,200.98 Da and 6.52, respectively. The squalene content of overexpressed PgFPS in tobacco leaves and yeast cells extract was 1.88-fold and 1.21-fold higher than that of the control group, respectively, and the total saponin content was also increased by 1.15 times in yeast cells extract, which verified the biological function of PgFPS in terpenoid synthesis. After 48 h of MeJA treatment and 6 h of ethephon treatment, the expression of the PgFPS gene in roots and stems reached its peak, showing a 3.125-fold and 3.236-fold increase compared to the untreated group, respectively. Interestingly, the expression of the PgFPS gene in leaves showed a decreasing trend after exogenous elicitors treatment. The discovery of this enzyme will provide a novel perspective for enhancing the efficient synthesis of platycodins.

Identifying the stability of housekeeping genes to be used for the quantitative real-time PCR normalization in retinal tissue of streptozotocin-induced diabetic rats

AIM

To investigate the stability of the seven housekeeping genes: beta-actin (ActB), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), 18s ribosomal unit 5 (18s), cyclophilin A (CycA), hypoxanthine-guanine phosphoribosyl transferase (HPRT), ribosomal protein large P0 (36B4) and terminal uridylyl transferase 1 (U6) in the diabetic retinal tissue of rat model.

METHODS

The expression of these seven genes in rat retinal tissues was determined using real-time quantitative reverse transcription polymerase chain reaction (RT-qPCR) in two groups; normal control rats and streptozotocin-induced diabetic rats. The stability analysis of gene expression was investigated using geNorm, NormFinder, BestKeeper, and comparative delta-Ct (ΔCt) algorithms.

RESULTS

The 36B4 gene was stably expressed in the retinal tissues of normal control animals; however, it was less stable in diabetic retinas. The 18s gene was expressed consistently in both normal control and diabetic rats' retinal tissue. That this gene was the best reference for data normalisation in RT-qPCR studies that used the retinal tissue of streptozotocin-induced diabetic rats. Furthermore, there was no ideal gene stably expressed for use in all experimental settings.

CONCLUSION

Identifying relevant genes is a need for achieving RT-qPCR validity and reliability and must be appropriately achieved based on a specific experimental setting.

Time- and Region-Specific Selection of Reference Genes in the Rat Brain in the Lithium-Pilocarpine Model of Acquired Temporal Lobe Epilepsy

Reverse transcription followed by quantitative polymerase chain reaction (RT-qPCR) is a commonly used tool for gene expression analysis. The selection of stably expressed reference genes is required for accurate normalization. The aim of this study was to identify the optimal reference genes for RT-qPCR normalization in various brain regions of rats at different stages of the lithium-pilocarpine model of acquired epilepsy. We tested the expression stability of nine housekeeping genes commonly used as reference genes in brain research: Actb, Gapdh, B2m, Rpl13a, Sdha, Ppia, Hprt1, Pgk1, and Ywhaz. Based on four standard algorithms (geNorm, NormFinder, BestKeeper, and comparative delta-Ct), we found that after pilocarpine-induced status epilepticus, the stability of the tested reference genes varied significantly between brain regions and depended on time after epileptogenesis induction (3 and 7 days in the latent phase, and 2 months in the chronic phase of the model). Pgk1 and Ywhaz were the most stable, while Actb, Sdha, and B2m demonstrated the lowest stability in the analyzed brain areas. We revealed time- and region-specific changes in the mRNA expression of the housekeeping genes B2m, Actb, Sdha, Rpl13a, Gapdh, Hprt1, and Sdha. These changes were more pronounced in the hippocampal region during the latent phase of the model and are thought to be related to epileptogenesis. Thus, RT-qPCR analysis of mRNA expression in acquired epilepsy models requires careful selection of reference genes depending on the brain region and time of analysis. For the time course study of epileptogenesis in the rat lithium-pilocarpine model, we recommend the use of the Pgk1 and Ywhaz genes.

Reference genes to study the sex-biased expression of genes regulating Drosophila metabolism

Sex is an important variable in biology. Notable differences have been observed between male and female Drosophila in regulation of metabolism, in response to nutritional challenges, and in phenotypes relevant for obesity and metabolic disorders. The differences between males and females can be expected to result from differences in gene expression. We observed that expression levels of reference genes commonly used for normalization of qRT-PCR results such as GAPDH, β-actin, and 18SrRNA, show prominent sexual dimorphism. Since this will impact relative expression and conclusions related to that, we performed a systematic analysis of candidate reference genes with the objective of identifying reference genes with stable expression in male and female Drosophila. These reference genes (LamCa, βTub60D and βTub97EF) were then used to assess sex-specific differences in expression of metabolism associated genes. Additionally, we evaluated the utility of these reference genes following a nutritional challenge and showed that LamCa and βtub97EF are stably expressed between sexes and under different nutritional conditions and are thus suitable as reference genes. Our results highlight the importance of evaluating the stability of reference genes when sex-specific differences in gene expression are studied, and identify structural genes as a category worth exploring as reference genes in other species. Finally, we also uncovered hitherto unknown sexually dimorphic expression of a number of metabolism-associated genes, information of interest to others working in the field of metabolic disorders.

Validation of qPCR reference genes in the endangered annual killifish Austrolebias charrua considering different tissues, gender and environmental conditions

The annual killifish Austrolebias charrua is an endangered species, endemic to the southern region of South America, which inhabits temporary ponds that emerges in the rainy season. The main anthropogenic threat driving the extinction of A. charrua stems from extensive agriculture, primarily due to the widrespread use of glyphosate-based herbicides near their habitats. Annual killifishes have been used as models for ecotoxicological studies but, up to now, there are no studies about reference genes in any Austrolebias species. This represents an obstacle to the use of qPCR-based technologies, the standard method for gene expression quantification. The present study aimed to select and validate potential reference genes for qPCR normalization in the annual killifish Austrolebias charrua considering different tissues, gender and environmental conditions. The candidate reference genes 18 s, actb, gapdh, ef1a, shox, eif3g, and the control gene atp1a1 were evaluated in male and female individuals in three different tissues (brain, liver, and gills) under two experimental conditions (control and acute exposition to Roundup Transorb®). The collected tissues were submitted to RNA extraction, followed by cDNA synthesis, cloning, sequencing, and qPCR. Overall, 18 s was the most stable reference gene, and 18 s and ef1a were the most stable combination. Otherwise, considering all variables, gapdh and shox were the least stable candidate genes. Foremost, suitable reference genes were validated in A. charrua, facilitating accurate mRNA quantification in this species, which might be useful for developing molecular tools of ecotoxicological assessment based on gene expression analysis for environmental monitoring of annual killifish.

Antennal Transcriptome Evaluation and Analysis for Odorant-Binding Proteins, Chemosensory Proteins, and Suitable Reference Genes in the Leaf Beetle Pest Diorhabda rybakowi Weise (Coleoptera: Chrysomelidae)

Diorhabda rybakowi Weise is one of the dominant pests feeding on Nitraria spp., a pioneer plant used for windbreaking and sand fixation purposes, and poses a threat to local livestock and ecosystems. To clarify the key olfactory genes of D. rybakowi and provide a theoretical basis for attractant and repellent development, the optimal reference genes under two different conditions (tissue and sex) were identified, and the bioinformatics and characterization of the tissue expression profiles of two categories of soluble olfactory proteins (OBPs and CSPs) were investigated. The results showed that the best reference genes were RPL13a and RPS18 for comparison among tissues, and RPL19 and RPS18 for comparison between sexes. Strong expressions of DrybOBP3, DrybOBP6, DrybOBP7, DrybOBP10, DrybOBP11, DrybCSP2, and DrybCSP5 were found in antennae, the most important olfactory organ for D. rybakowi. These findings not only provide a basis for further in-depth research on the olfactory molecular mechanisms of host-specialized pests but also provide a theoretical basis for the future development of new chemical attractants or repellents using volatiles to control D. rybakowi.

EndoGeneAnalyzer: A tool for selection and validation of reference genes

The selection of proper reference genes is critical for accurate gene expression analysis in all fields of biological and medical research, mainly because there are many distinctions between different tissues and specimens. Given this variability, even in known classic reference genes, demands of a comprehensive analysis platform is needed to identify the most suitable genes for each study. For this purpose, we present an analysis tool for assisting in decision-making in the analysis of reverse transcription-quantitative polymerase chain reaction (RT-qPCR) data. EndoGeneAnalyzer, an open-source web tool for reference gene analysis in RT-qPCR studies, was used to compare the groups/conditions under investigation. This interactive application offers an easy-to-use interface that allows efficient exploration of datasets. Through statistical and stability analyses, EndoGeneAnalyzer assists in the select of the most appropriate reference gene or set of genes for each condition. It also allows researchers to identify and remove unwanted outliers. Moreover, EndoGeneAnalyzer provides a graphical interface to compare the evaluated groups, providing a visually informative differential analysis.

Improving the quality of quantitative polymerase chain reaction experiments: 15 years of MIQE

The quantitative polymerase chain reaction (qPCR) is fundamental to molecular biology. It is not just a laboratory technique, qPCR is a bridge between research and clinical practice. Its theoretical foundations guide the design of experiments, while its practical implications extend to diagnostics, treatment, and research advancements in the life sciences, human and veterinary medicine, agriculture, and forensics. However, the accuracy, reliability and reproducibility of qPCR data face challenges arising from various factors associated with experimental design, execution, data analysis and inadequate reporting details. Addressing these concerns, the Minimum Information for the Publication of Quantitative Real-Time PCR Experiments (MIQE) guidelines have emerged as a cohesive framework offering a standardised set of recommendations that describe the essential information required for assessing qPCR experiments. By emphasising the importance of methodological rigour, the MIQE guidelines have made a major contribution to improving the trustworthiness, consistency, and transparency of many published qPCR results. However, major challenges related to awareness, resources, and publication pressures continue to affect their consistent application.

Selection and Validation of qRT-PCR Internal Reference Genes to Study Flower Color Formation in Camellia impressinervis

Real-time quantitative PCR (qRT-PCR) is a pivotal technique for gene expression analysis. To ensure reliable and accurate results, the internal reference genes must exhibit stable expression across varied experimental conditions. Currently, no internal reference genes for Camellia impressinervis have been established. This study aimed to identify stable internal reference genes from eight candidates derived from different developmental stages of C. impressinervis flowers. We employed geNorm, NormFinder, and BestKeeper to evaluate the expression stability of these candidates, which was followed by a comprehensive stability analysis. The results indicated that CiTUB, a tubulin gene, exhibited the most stable expression among the eight reference gene candidates in the petals. Subsequently, CiTUB was utilized as an internal reference for the qRT-PCR analysis of six genes implicated in the petal pigment synthesis pathway of C. impressinervis. The qRT-PCR results were corroborated by transcriptome sequencing data, affirming the stability and suitability of CiTUB as a reference gene. This study marks the first identification of stable internal reference genes within the entire genome of C. impressinervis, establishing a foundation for future gene expression and functional studies. Identifying such stable reference genes is crucial for advancing molecular research on C. impressinervis.

Assessing the feasibility of free DNA for disaster victim identification and forensic applications

In tropical disaster victim identification (DVI) scenarios, challenging environmental conditions lead to accelerated DNA degradation in remains. To further enhance the utilization of leached DNA from tissue in the preservative solution (termed "free DNA") as an alternative source, we incorporated new results by assessing its integrity in postmortem and decomposing cadavers preserved in DNA/RNA Shield™ and modified TENT, with silica-based purification (QIAquick®) for faster processing. The psoas muscle tissues of one decomposed and ten cadavers were preserved in each solution at 25 °C and 35 °C for 3 months. Free DNA efficiency was compared with individual reference samples for reliable results in quantity, quality, and STR profiles. The findings revealed that DNA/RNA Shield™ effectively preserves free DNA integrity for extended storage, while modified TENT is more suitable for short-term storage due to higher degradation levels. Moreover, the use of free DNA samples with massive parallel sequencing displays potential for forensic DNA analysis. Successful amplification of the mtDNA control region enables variant calling and heteroplasmy analysis while also serving as quality control using ACTB and enabling differentiation within the 16S rRNA region for microbiome analysis. The simplicity of handling free DNA for PCR-based forensic analysis adds to its potential for various applications, including DVI and field-based analysis of biological evidence.

Assessment of Reference Genes Stability in Cortical Bone of Obese and Diabetic Mice

Introduction

Bone, a pivotal structural organ, is susceptible to disorders with profound health implications. The investigation of gene expression in bone tissue is imperative, particularly within the context of metabolic diseases such as obesity and diabetes that augment the susceptibility to bone fractures. The objective of this study is to identify a set of internal control genes for the analysis of gene expression.

Methods

This study employs reverse transcription-quantitative real-time polymerase chain reaction (RT-qPCR) to assess gene expression in bone tissue. We selected fourteen housekeeping genes and assessed their stability in the cortical bone of mouse models for obesity and diabetes using four well-established algorithms (GeNorm, BestKeeper, NormFinder, and the comparative Delta Ct method).

Results and Conclusion

We identified Rpl13a as the mostly stably expressed reference gene in cortical bone tissue from mouse models of obesity and diabetes (db/db), while Gapdh was found to be the most stable reference gene in another diabetes model, KKAy mice. Additionally, Ef1a, Ppia, Rplp0, and Rpl22 were identified as alternative genes suitable for normalizing gene expression in cortical bone from obesity and diabetes mouse models. These findings enhance RT-qPCR accuracy and reliability, offering a strategic guide to select reference gene for studying bone tissue gene expression in metabolic disorders.

Reference Gene Validation in the Embryonic and Postnatal Brain in the Rat Hyperhomocysteinemia Model

Maternal hyperhomocysteinemia (HCY) induced by genetic defects in methionine cycle enzymes or vitamin imbalance is known to be a pathologic factor that can impair embryonal brain development and cause long-term consequences in the postnatal brain development as well as changes in the expression of neuronal genes. Studies of the gene expression on this model requires the selection of optimal housekeeping genes. This work aimed to analyze the expression stability of housekeeping genes in offspring brain. Pregnant female Wistar rats were treated daily with a 0.15% L-methionine solution in the period starting on the 4th day of pregnancy until delivery, to cause the increase in the homocysteine level in fetus blood and brain. Housekeeping gene expression was assessed by RT-qPCR on whole embryonic brain and selected rat brain areas at P20 and P90. The amplification curves were analyzed, and raw means Cq data were imported to the RefFinder online tool to assess the reference genes stability. Most of the analyzed genes showed high stability of mRNA expression in the fetal brain at both periods of analysis (E14 and E20). However, the most stably expressed genes at different age points differed. Actb, Ppia, Rpl13a are the most stably expressed on E14, Ywhaz, Pgk1, Hprt1 - on E20 and P20, Hprt1, Actb, and Pgk1 - on P90. Gapdh gene used as a reference in various studies demonstrates high stability only in the hippocampus and cannot be recommended as the optimal reference gene on HCY model. Hprt1 and Pgk1 genes were found to be the most stably expressed in the brain of rat subjected to HCY. These two genes showed high stability in the brain on E20 and in various areas of the brain on the P20 and P90. On E14, the preferred genes for normalization are Actb, Ppia, Rpl13a.

Runx2 Regulates Galnt3 and Fgf23 Expressions and Galnt3 Decelerates Osteoid Mineralization by Stabilizing Fgf23

Runx2 (runt related transcription factor 2) is an essential transcription factor for osteoblast proliferation and differentiation. Uridine diphosphate (UDP)-N-acetylgalactosamine (GalNAc): polypeptide GalNAc-transferase 3 (Galnt3) prevents proteolytic processing of fibroblast growth factor 23 (Fgf23), which is a hormone that regulates the serum level of phosphorus. Runx2 and Galnt3 were expressed in osteoblasts and osteocytes, and Fgf23 expression was restricted to osteocytes in bone. Overexpression and knock-down of Runx2 upregulated and downregulated, respectively, the expressions of Galnt3 and Fgf23, and Runx2 directly regulated the transcriptional activity of Galnt3 in reporter assays. The expressions of Galnt3 and Fgf23 in osteoblast-specific Runx2 knockout (Runx2fl/flCre) mice were about half those in Runx2fl/fl mice. However, the serum levels of phosphorus and intact Fgf23 in Runx2fl/flCre mice were similar to those in Runx2fl/fl mice. The trabecular bone volume was increased during aging in both male and female Galnt3-/- mice, but the osteoid was reduced. The markers for bone formation and resorption in Galnt3-/- mice were similar to the control in both sexes. Galnt3-/- mice exhibited hyperphosphatemia and hypercalcemia, and the intact Fgf23 was about 40% that of wild-type mice. These findings indicated that Runx2 regulates the expressions of Galnt3 and Fgf23 and that Galnt3 decelerates the mineralization of osteoid by stabilizing Fgf23.

Role of BAL and Serum Krebs von den Lungen-6 (KL-6) in Patients with Pulmonary Fibrosis

Background: Interstitial lung diseases (ILDs) encompass a diverse group of disorders affecting the lung interstitium, leading to inflammation, fibrosis, and impaired respiratory function. Currently, the identification of new diagnostic and prognostic biomarkers for ILDs turns out to be necessary. Several studies show the role of KL-6 in various types of interstitial lung disease and suggest that serum KL-6 levels can be used as a prognostic marker of disease. The aim of this study was to analyze KL-6 expression either in serum or bronchoalveolar lavage samples in order to: (i) make a serum vs. BAL comparison; (ii) better understand the local behavior of fibrosis vs. the systemic one; and (iii) evaluate any differences in patients with progressive fibrosis (PPF) versus patients with non-progressive fibrosis (nPPF). Methods: We used qRT-PCR to detect KL-6 expression both in serum and BAL samples. Mann-Whitney's U test was used to compare the differential expression between groups. Results: In serum, KL-6 is more highly expressed in PPF than in non-progressive fibrosis (p = 0.0295). This difference is even more significant in BAL (p < 0.001). Therefore, it is clear that KL-6 values are related to disease progression. Significant differences were found by making a comparison between BAL and serum. KL-6 was markedly higher in serum than BAL (p = 0.0146). Conclusions: This study identifies KL-6 as a promising biomarker for the severity of the fibrosing process and disease progression in ILDs, with significantly higher levels observed in PPF compared to nPPF. Moreover, the marked difference in KL-6 levels between serum and BAL emphasizes its potential diagnostic and prognostic relevance, providing enlightening insights into both the local and systemic aspects of ILDs.

Validation of Endogenous Control Genes by Real-Time Quantitative Reverse Transcriptase Polymerase Chain Reaction for Acute Leukemia Gene Expression Studies

Reference genes are used as internal reaction controls for gene expression analysis, and for this reason, they are considered reliable and must meet several important criteria. In view of the absence of studies regarding the best reference gene for the analysis of acute leukemia patients, a panel of genes commonly used as endogenous controls was selected from the literature for stability analysis: Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), Abelson murine leukemia viral oncogene human homolog 1 (ABL), Hypoxanthine phosphoribosyl-transferase 1 (HPRT1), Ribosomal protein lateral stalk subunit P0 (RPLP0), β-actin (ACTB) and TATA box binding protein (TBP). The stability of candidate reference genes was analyzed according to three statistical methods of assessment, namely, NormFinder, GeNorm and R software (version 4.0.3). From this study's analysis, it was possible to identify that the endogenous set composed of ACTB, ABL, TBP and RPLP0 demonstrated good performances and stable expressions between the analyzed groups. In addition to that, the GAPDH and HPRT genes could not be classified as good reference genes, considering that they presented a high standard deviation and great variability between groups, indicating low stability. Given these findings, this study suggests the main endogenous gene set for use as a control/reference for the gene expression in peripheral blood and bone marrow samples from patients with acute leukemias is composed of the ACTB, ABL, TBP and RPLP0 genes. Researchers may choose two to three of these housekeeping genes to perform data normalization.

Identification and validation of the optimal reference genes for standardizing the gene expression profiling diagnostic panel of Ph-like B-lineage acute lymphoblastic leukemia

Gene expression profiling is the criterion standard for recognizing Ph-like ALL signatures among B-ALLs. The prerequisite of GEP is the accurate normalization of target genes with stable expression of housekeeping genes in a quantitative PCR. HKGs exhibit differential expression in the different experimental conditions and affect the target genes' expression, leading to imprecise qPCR results. The selection of stable HKGs is crucial in GEP experiments, especially in identifying high-risk Ph-like ALL cases. We have evaluated the expression stability of nine HKGs (GAPDH, ACTB, GUSB, RNA18S, EEF2, PGK1, B2M, TBP and ABL1) in identified Ph-like ALLs and Ph-negative (n = 23 each) using six algorithms, 4 traditional softwares; geNorm, BestKeeper, NormFinder, Delta Cq value method, and two algorithms, RefFinderTM and ComprFinder. Further, we have validated the expression of 8 overexpressed normalized genes in Ph-like ALL cases (JCHAIN, CA6, MUC4, SPATS2L, BMPR1B, CRLF2, ADGRF1 and NRXN3). GeNorm, BestKeeper, NormFinder, Delta Cq value method, RefFinderTM and ComprFinder algorithm analysis revealed that EEF2, GAPDH, and PGK1 form the best representative HKGs in Ph-like ALL cases, while RNA18s, ß-actin, and ABL1 in Ph-negative ALLs. Lastly, we performed a correlation analysis and found that the combination of EEF2, GAPDH, and PGK1 represents the best combination with a very high correlation in Ph-like ALL cases. This is the first report that shows EEF2, GAPDH, and PGK1 are the best HKG genes and can be used in the diagnostic panel of Ph-like ALL cases using qPCR at baseline diagnosis.

Stat1 is an inducible transcriptional repressor of neural stem cells self-renewal program during neuroinflammation

A central issue in regenerative medicine is understanding the mechanisms that regulate the self-renewal of endogenous stem cells in response to injury and disease. Interferons increase hematopoietic stem cells during infection by activating STAT1, but the mechanisms by which STAT1 regulates intrinsic programs in neural stem cells (NSCs) during neuroinflammation is less known. Here we explored the role of STAT1 on NSC self-renewal. We show that overexpressing Stat1 in NSCs derived from the subventricular zone (SVZ) decreases NSC self-renewal capacity while Stat1 deletion increases NSC self-renewal, neurogenesis, and oligodendrogenesis in isolated NSCs. Importantly, we find upregulation of STAT1 in NSCs in a mouse model of multiple sclerosis (MS) and an increase in pathological T cells expressing IFN-γ rather than interleukin 17 (IL-17) in the cerebrospinal fluid of affected mice. We find IFN-γ is superior to IL-17 in reducing proliferation and precipitating an abnormal NSC phenotype featuring increased STAT1 phosphorylation and Stat1 and p16ink4a gene expression. Notably, Stat1-/- NSCs were resistant to the effect of IFN-γ. Lastly, we identified a Stat1-dependent gene expression profile associated with an increase in the Sox9 transcription factor, a regulator of self-renewal. Stat1 binds and transcriptionally represses Sox9 in a transcriptional luciferase assay. We conclude that Stat1 serves as an inducible checkpoint for NSC self-renewal that is upregulated during chronic brain inflammation leading to decreased self-renewal. As such, Stat1 may be a potential target to modulate for next generation therapies to prevent progression and loss of repair function in NSCs/neural progenitors in MS.

Identification of upper thermal thresholds during development in the endangered Nechako white sturgeon with management implications for a regulated river

Climate change-induced warming effects are already evident in river ecosystems, and projected increases in temperature will continue to amplify stress on fish communities. In addition, many rivers globally are impacted by dams, which have many negative effects on fishes by altering flow, blocking fish passage, and changing sediment composition. However, in some systems, dams present an opportunity to manage river temperature through regulated releases of cooler water. For example, there is a government mandate for Kenney dam operators in the Nechako river, British Columbia, Canada, to maintain river temperature <20°C in July and August to protect migrating sockeye salmon (Oncorhynchus nerka). However, there is another endangered fish species inhabiting the same river, Nechako white sturgeon (Acipenser transmontanus), and it is unclear if these current temperature regulations, or timing of the regulations, are suitable for spawning and developing sturgeon. In this study, we aimed to identify upper thermal thresholds in white sturgeon embryos and larvae to investigate if exposure to current river temperatures are playing a role in recruitment failure. We incubated embryos and yolk-sac larvae in three environmentally relevant temperatures (14, 18 and 21°C) throughout development to identify thermal thresholds across different levels of biological organization. Our results demonstrate upper thermal thresholds at 21°C across physiological measurements in embryo and yolk-sac larvae white sturgeon. Before hatch, both embryo survival and metabolic rate were reduced at 21°C. After hatch, sublethal consequences continued at 21°C because larval sturgeon had decreased thermal plasticity and a dampened transcriptional response during development. In recent years, the Nechako river has reached 21°C by the end of June, and at this temperature, a decrease in sturgeon performance is evident in most of the traits measured. As such, the thermal thresholds identified here suggest current temperature regulations may not be suitable for developing white sturgeon and future recruitment.

Evaluation of reference genes for qRT-PCR studies in the colchicine producing Gloriosa superba L

The flame lily, Gloriosa superba L., is one of the two primary sources of the anti-inflammatory drug, colchicine. Previous studies have shown that a higher level of colchicine production occurs in the rhizomes than in leaves and roots. Earlier precursor feeding and transcriptome analysis of G. superba have provided a putative pathway and candidate genes involved in colchicine biosynthesis. Comparative analysis of expression levels of candidate pathway genes in different tissues of G. superba using quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR) can reveal highly expressed genes in the rhizome compared to other tissues which could suggest roles of the gene products in colchicine biosynthesis. Normalization is an important step in effectively analyzing differential gene expression by qRT-PCR with broader applications. The current study selected candidate reference genes from the transcriptome datasets and analyzed them to determine the most stable genes for normalization of colchicine biosynthesis-related genes. Using RefFinder, one stable reference gene, UBC22, was selected to normalize gene expression levels of candidate methyltransferase (MT) genes in the leaves, roots, and rhizomes of G. superba. With UBC22 as reference gene, the methyltransferases, GsOMT1, GsOMT3, and GsOMT4 showed significantly higher expression levels in the rhizome of G. superba, while MT31794 was more highly expressed in the roots. In conclusion, the current results showed a viable reference gene expression analysis system that could help elucidate colchicine biosynthesis and its exploitation for increased production of the drug in G. superba.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11816-023-00840-x.

Identification of quantitative polymerase chain reaction reference genes suitable for normalising gene expression in the brain of normal and dystrophic mice and dogs

Background: In addition to progressive, debilitating muscle degeneration, ~50% of patients with Duchenne muscular dystrophy (DMD) have associated cognitive and behavioural disorders secondary to deficiency of dystrophin protein in the brain. The brain expresses a variety of dystrophin isoforms (Dp427, Dp140 and Dp71) whose functions remain to be fully elucidated. Detailed comparative analysis of gene expression in healthy and dystrophin-deficient brain is fundamental to understanding the functions of each isoform, and the consequences of their deficiency, with animal models representing a key tool in this endeavour. Reverse transcription quantitative real-time PCR (RT-qPCR) is a widely used method to study gene expression. However, accurate quantitative assessment requires normalisation of expression data using validated reference genes. The aim of this study was to identify a panel of suitable reference genes that can be used to normalise gene expression in the brain of healthy and dystrophic dogs and mice. Methods: Using the DE50-MD dog and mdx mouse models of DMD we performed RT-qPCR from fresh frozen brain tissue and employed the geNorm, BestKeeper and Normfinder algorithms to determine the stability of expression of a panel of candidate reference genes across healthy and dystrophic animals, and across different brain regions. Results: We show that SDHA, UBC and YWHAZ are suitable reference genes for normalising gene expression in healthy and dystrophic canine brain, and GAPDH, RPL13A and CYC1 in healthy and dystrophic murine brain. Notably, there was no overlap in the highest performing reference genes between the two species. Conclusions: Our findings suggest that gene expression normalisation is possible across six regions of the canine brain, and three regions of the murine brain. Our results should facilitate future work to study gene expression in the brains of normal and dystrophic dogs and mice and thus decipher the transcriptional consequences of dystrophin deficiency in the brain.

Effect of HSPA8 gene on the proliferation, apoptosis and immune function of HD11 cells

HSPA8 (Heat shock 70 kDa protein 8) is a molecular chaperone involved in a variety of cellular processes. This gene may affect the proliferation, apoptosis and immune function of chicken macrophages, but the specific mechanism remains unclear. The purpose of this study was to explore the effect of the HSPA8 gene on the proliferation, apoptosis and immune function of chicken macrophages. In this study, a chicken HSPA8 overexpression plasmid, interference fragment and corresponding controls were transfected into HD11 cells, and then the expression of the HSPA8 gene, cell proliferation, cell cycle, apoptosis rate and immune function of each group were detected. The results showed that transfection of the HSPA8 overexpression plasmid significantly upregulated the level of HSPA8 expression in HD11 cells compared with the control; significantly promoted the proliferation of HD11 cells and the expression of PCNA, CCND1 and CCNB3; decreased the number of cells in the G1 phase and increased the number of cells in the S phase; decreased the rate of apoptosis and upregulated the expression of Bcl-2; and promoted the expression of the LPS-induced cytokines IL-1β, IL-6 and TNF-α. Transfection of the HSPA8 interference fragment significantly downregulated the level of HSPA8 expression in HD11 cells; significantly inhibited the proliferation of HD11 cells and the expression of PCNA, CCND1 and CDK1; increased the number of cells in the G1 phase and decreased the number of cells in the S phase; increased the rate of apoptosis, downregulated the expression of Bcl-2 and upregulated the expression levels of Fas and FasL; and inhibited the expression of the LPS-induced cytokines IL-1β and NF-κB. The results suggested that HSPA8 promotes the proliferation of and inhibits the apoptosis of HD11 cells and has a proinflammatory effect.

Identification of Optimal Reference Genes for qRT-PCR Normalization for Physical Activity Intervention and Omega-3 Fatty Acids Supplementation in Humans

The quantitative polymerase chain reaction (qRT-PCR) technique gives promising opportunities to detect and quantify RNA targets and is commonly used in many research fields. This study aimed to identify suitable reference genes for physical exercise and omega-3 fatty acids supplementation intervention. Forty healthy, physically active men were exposed to a 12-week eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) supplementation and standardized endurance training protocol. Blood samples were collected before and after the intervention and mRNA levels of six potential reference genes were tested in the leukocytes of 18 eligible participants using the qRT-PCR method: GAPDH (Glyceraldehyde-3-phosphate dehydrogenase), ACTB (Beta actin), TUBB (Tubulin Beta Class I), RPS18 (Ribosomal Protein S18), UBE2D2 (Ubiquitin-conjugating enzyme E2 D2), and HPRT1 (Hypoxanthine Phosphoribosyltransferase 1). The raw quantification cycle (Cq) values were then analyzed using RefFinder, an online tool that incorporates four different algorithms: NormFinder, geNorm, BestKeeper, and the comparative delta-Ct method. Delta-Ct, NormFinder, BestKeeper, and RefFinder comprehensive ranking have found GAPDH to be the most stably expressed gene. geNorm has identified TUBB and HPRT as the most stable genes. All algorithms have found ACTB to be the least stably expressed gene. A combination of the three most stably expressed genes, namely GAPDH, TUBB, and HPRT, is suggested for obtaining the most reliable results.

Validation of reference genes for the normalization of the RT-qPCR in peripheral blood mononuclear cells of septic patients

Objective

To screen and validate reference genes suitable for gene mRNA expression study in peripheral blood mononuclear cells (PBMCs) between septic patients and healthy controls (HC).

Methods

Total RNA in PBMCs was extracted and RT-qPCR was used to determine the mRNA expression profiles of 9 candidate genes, including ACTB, B2M, GAPDH, GUSB, HPRT1, PGK1, RPL13A, SDHA and YWHAZ. The genes expression stabilities were assessed by both geNorm and NormFinder software.

Results

YWHAZ was the most stable gene among the 9 candidate genes evaluated by both geNorm and NormFinder in mixed and sepsis groups. The most stable gene combination in mixed group analyzed by geNorm was the combination of GAPDH, PKG1 and YWHAZ, while that in sepsis group was the combination of ACTB, PKG1 and YWHAZ.

Conclusion

Our first systematic analysis of the reference genes in PBMC of septic patients suggested YWHAZ was the best candidate. The combination of ACTB, PKG1 and YWHAZ could improve RT-qPCR accuracy in septic patients. Our results identified the most stable reference genes to standardize RT-qPCR of sepsis patients, which can serve as a useful tool for gene function exploration in the future.

Study of the Genetic Expression of Antiretroviral Restriction Factors and Acute Phase Proteins in Cattle Infected with Bovine Leukemia Virus

The goal of this study was to analyze the genetic expression of antiretroviral restriction factors (ARF) and acute phase proteins (APP), as well as their correlation with proviral and viral loads in cattle with aleukemic (AL) and persistent lymphocytosis (PL). Complete blood samples were collected from a herd of dairy cows, and we extracted genetic material from peripheral blood leukocytes. Absolute quantification of the expression of ARF (APOBEC-Z1, Z2, and Z3; HEXIM-1, HEXIM-2, and BST2) and APP (haptoglobin (HP), and serum amyloid A (SAA)) was performed by qPCR. Statistical significance was observed in the expression of APOBEC-Z3 in BLV-infected animals. We only found positive correlations with a strong expression of the ARF genes in the AL group. The participation of APOBEC (Z1 and Z3), HEXIM-1, and HEXIM-2 was more frequently identified in BLV-infected animals. HEXIM-2 showed active gene expression in the AL group. Although the expression of ARF in early stages of infection (AL) maintains an important participation, in late stages (PL) it seems to have little relevance.

Quantitative evaluation of endogenous reference genes for ddPCR under salt stress using a moderate halophile

Droplet digital PCR (ddPCR) is being increasingly adopted for gene detection and quantification because of its higher sensitivity and specificity. According to previous observations and our laboratory data, it is essential to use endogenous reference genes (RGs) when investigating gene expression at the mRNA level under salt stress. This study aimed to select and validate suitable RGs for gene expression under salt stress using ddPCR. Six candidate RGs were selected based on the tandem mass tag (TMT)-labeled quantitative proteomics of Alkalicoccus halolimnae at four salinities. The expression stability of these candidate genes was evaluated using statistical algorithms (geNorm, NormFinder, BestKeeper and RefFinder). There was a small fluctuation in the cycle threshold (Ct) value and copy number of the pdp gene. Its expression stability was ranked in the vanguard of all algorithms and was the most suitable RG for quantification of expression by both qPCR and ddPCR of A. halolimnae under salt stress. Single RG pdp and RG combinations were used to normalize the expression of ectA, ectB, ectC and ectD under four salinities. The present study constitutes the first systematic analysis of endogenous RG selection for halophiles responding to salt stress. This work provides a valuable theory and an approach reference of internal control identification for ddPCR-based stress response models.

Genome-Wide Identification and Expression Profiling Analysis of the Trihelix Gene Family and response of PgGT1 under Abiotic Stresses in Platycodon grandiflorus

The trihelix gene family plays an important role in plant growth and abiotic stress responses. Through the analysis of genomic and transcriptome data, 35 trihelix family members were identified for the first time in Platycodon grandiflorus; they were classified into five subfamilies: GT-1, GT-2, SH4, GTγ, and SIP1. The gene structure, conserved motifs and evolutionary relationships were analyzed. Prediction of physicochemical properties of the 35 trihelix proteins founded, the number of amino acid molecules is between 93 and 960, theoretical isoelectric point is between 4.24 and 9.94, molecular weight is between 9829.77 and 107435.38, 4 proteins among them were stable, and all GRAVY is negative. The full-length cDNA sequence of the PgGT1 gene of the GT-1 subfamily was cloned by PCR. It is a 1165 bp ORF encoding a 387 amino acid protein, with a molecular weight of 43.54 kDa. The predicted subcellular localization of the protein in the nucleus was experimentally verified. After being treated with NaCl, PEG6000, MeJA, ABA, IAA, SA, and ethephon, the expression of PgGT1 gene showed an up-regulated trend except for the roots treated with NaCl and ABA. This study laid a bioinformatics foundation for the research of trihelix gene family and the cultivation of excellent germplasm of P. grandiflorus.

Reference gene expression stability within the rat brain under mild intermittent ketosis induced by supplementation with medium-chain triglycerides

Reverse transcription followed by quantitative (real-time) polymerase chain reaction (RT-qPCR) has become the gold standard in mRNA expression analysis. However, it requires an accurate choice of reference genes for adequate normalization. The aim of this study was to validate the reference genes for qPCR experiments in the brain of rats in the model of mild ketosis established through supplementation with medium-chain triglycerides (MCT) and intermittent fasting. This approach allows to reproduce certain neuroprotective effects of the classical ketogenic diet while avoiding its adverse effects. Ketogenic treatment targets multiple metabolic pathways, which may affect the reference gene expression. The standard chow of adult Wistar rats was supplemented with MCT (2 ml/kg orogastrically, during 6 h of fasting) or water (equivolume) for 1 month. The mRNA expression of 9 housekeeping genes (Actb, B2m, Gapdh, Hprt1, Pgk1, Ppia, Rpl13a, Sdha, Ywhaz) in the medial prefrontal cortex, dorsal and ventral hippocampus was measured by RT-qPCR. Using the RefFinder® online tool, we have found that the reference gene stability ranking strongly depended on the analyzed brain region. The most stably expressed reference genes were found to be Ppia, Actb, and Rpl13a in the medial prefrontal cortex; Rpl13a, Ywhaz, and Pgk1 in the dorsal hippocampus; Ywhaz, Sdha, and Ppia in the ventral hippocampus. The B2m was identified as an invalid reference gene in the ventral hippocampus, while Sdha, Actb, and Gapdh were unstable in the dorsal hippocampus. The stabilities of the examined reference genes were lower in the dorsal hippocampus compared to the ventral hippocampus and the medial prefrontal cortex. When normalized to the three most stably expressed reference genes, the Gapdh mRNA was upregulated, while the Sdha mRNA was downregulated in the medial prefrontal cortex of MCT-fed animals. Thus, the expression stability of reference genes strongly depends on the examined brain regions. The dorsal and ventral hippocampal areas differ in reference genes stability rankings, which should be taken into account in the RT-qPCR experimental design.

Invited review: Muscle protein breakdown and its assessment in periparturient dairy cows

Mobilization of body reserves including fat, protein, and glycogen is necessary to overcome phases of negative nutrient balance typical for high-yielding dairy cows during the periparturient period. Skeletal muscle, the largest internal organ in mammals, plays a crucial role in maintaining metabolic homeostasis. However, unlike in liver and adipose tissue, the metabolic and regulatory role of skeletal muscle in the adaptation of dairy cows to the physiological needs of pregnancy and lactation has not been studied extensively. The functional integrity and quality of skeletal muscle are maintained through a constant turnover of protein, resulting from both protein breakdown and protein synthesis. Thus, muscle protein breakdown (MPB) and synthesis are intimately connected and tightly controlled to ensure proper protein homeostasis. Understanding the regulation of MPB, the catabolic component of muscle turnover, and its assessment are therefore important considerations to provide information about the timing and extent of tissue mobilization in periparturient dairy cows. Based on animal models and human studies, it is now evident that MPB occurs via the integration of 3 main systems: autophagy-lysosomal, calpain Ca²⁺-dependent cysteine proteases, and the ubiquitin-proteasome system. These 3 main systems are interconnected and do not work separately, and the regulation is complex. The ubiquitin-proteasomal system is the most well-known cellular proteolytic system and plays a fundamental role in muscle physiology. Complete degradation of a protein often requires a combination of the systems, depending on the physiological situation. Determination of MPB in dairy cows is technically challenging, resulting in a relative dearth of information. The methods for assessing MPB can be divided into either direct or indirect measurements, both having their strengths and limitations. Available information on the direct measures of MPB primarily comes from stable isotopic tracer methods and those of indirect measurements from assessing expression and activity measures of the components of the 3 MPB systems in muscle biopsy samples. Other indirect approaches (i.e., potential indicators of MPB), including ultrasound imaging and measuring metabolites from muscle degradation (i.e., 3-methylhistidine and creatinine), seem to be applicable methods and can provide useful information about the extent and timing of MPB. This review presents our current understanding, including methodological considerations, of the process of MPB in periparturient dairy cows.

Validation of the Reference Genes for Expression Analysis in the Hippocampus after Transient Ischemia/Reperfusion Injury in Gerbil Brain

Transient brain ischemia in gerbils is a common model to study the mechanisms of neuronal changes in the hippocampus. In cornu ammonnis 2-3, dentate gyrus (CA2-3,DG) regions of the hippocampus, neurons are resistant to 5-min ischemia/reperfusion (I/R) insult, while cornu ammonnis 1 (CA1) is found to be I/R-vulnerable. The quantitative polymerase chain reaction (qRT-PCR) is widely used to study the expression of genes involved in these phenomena. It requires stable and reliable genes for normalization, which is crucial for comparable and reproducible analyses of expression changes of the genes of interest. The aim of this study was to determine the best housekeeping gene for the I/R gerbil model in two parts of the hippocampus in controls and at 3, 48, and 72 h after recanalization. We selected and tested six reference genes frequently used in central nervous system studies: Gapdh, Actb, 18S rRNA, Hprt1, Hmbs, Ywhaz, and additionally Bud23, using RefFinder, a comprehensive tool based on four commonly used algorithms: delta cycle threshold (Ct), BestKeeper, NormFinder, and geNorm, while Hprt1 and Hmbs were the most stable ones in CA2-3,DG. Hmbs was the most stable in the whole hippocampal formation. This indicates that the general use of Hmbs, especially in combination with Gapdh, a highly expressed reference gene, seems to be suitable for qRT-PCR normalization in all hippocampal regions in this model.

Selection and Evaluation of mRNA and miRNA Reference Genes for Expression Studies (qPCR) in Archived Formalin-Fixed and Paraffin-Embedded (FFPE) Colon Samples of DSS-Induced Colitis Mouse Model

The choice of appropriate reference genes is essential for correctly interpreting qPCR data and results. However, the majority of animal studies use a single reference gene without any prior evaluation. Therefore, many qPCR results from rodent studies can be misleading, affecting not only reproducibility but also translatability. In this study, the expression stability of reference genes for mRNA and miRNA in archived FFPE samples of 117 C57BL/6JOlaHsd mice (males and females) from 9 colitis experiments (dextran sulfate sodium; DSS) were evaluated and their expression analysis was performed. In addition, we investigated whether normalization reduced/neutralized the influence of inter/intra-experimental factors which we systematically included in the study. Two statistical algorithms (NormFinder and Bestkeeper) were used to determine the stability of reference genes. Multivariate analysis was made to evaluate the influence of normalization with different reference genes on target gene expression in regard to inter/intra-experimental factors. Results show that archived FFPE samples are a reliable source of RNA and imply that the FFPE procedure does not change the ranking of stability of reference genes obtained in fresh tissues. Multivariate analysis showed that the histological picture is an important factor affecting the expression levels of target genes.

A home-made pipette droplet microfluidics rapid prototyping and training kit for digital PCR, microorganism/cell encapsulation and controlled microgel synthesis

Droplet microfluidics offers a platform from which new digital molecular assay, disease screening, wound healing and material synthesis technologies have been proposed. However, the current commercial droplet generation, assembly and imaging technologies are too expensive and rigid to permit rapid and broad-range tuning of droplet features/cargoes. This rapid prototyping bottleneck has limited further expansion of its application. Herein, an inexpensive home-made pipette droplet microfluidics kit is introduced. This kit includes elliptical pipette tips that can be fabricated with a simple DIY (Do-It-Yourself) tool, a unique tape-based or 3D printed shallow-center imaging chip that allows rapid monolayer droplet assembly/immobilization and imaging with a smart-phone camera or miniature microscope. The droplets are generated by manual or automatic pipetting without expensive and lab-bound microfluidic pumps. The droplet size and fluid viscosity/surface tension can be varied significantly because of our particular droplet generation, assembly and imaging designs. The versatility of this rapid prototyping kit is demonstrated with three representative applications that can benefit from a droplet microfluidic platform: (1) Droplets as microreactors for PCR reaction with reverse transcription to detect and quantify target RNAs. (2) Droplets as microcompartments for spirulina culturing and the optical color/turbidity changes in droplets with spirulina confirm successful photosynthetic culturing. (3) Droplets as templates/molds for controlled synthesis of gold-capped polyacrylamide/gold composite Janus microgels. The easily fabricated and user-friendly portable kit is hence ideally suited for design, training and educational labs.

Identification of reference genes for the normalization of retinal mRNA expression by RT-qPCR in oxygen induced retinopathy, anemia, and erythropoietin administration

BACKGROUND

Anemia and retinopathy of prematurity (ROP) are common comorbidities experienced by preterm infants, yet the role of anemia on the pathogenesis of ROP remains unclear. Reverse-transcriptase quantitative polymerase chain reaction (RT-qPCR) is a sensitive technique for estimating the gene expression changes at the transcript level but requires identification of stably expressed reference genes for accurate data interpretation. This is particularly important for oxygen induced retinopathy studies given that some commonly used reference genes are sensitive to oxygen. This study aimed to identify stably expressed reference genes among eight commonly used reference genes in the neonatal rat pups' retina upon exposure to cyclic hyperoxia-hypoxia, anemia, and erythropoietin administration at two age groups (P14.5 and P20) using Bestkeeper, geNorm, and Normfinder, three publicly available, free algorithms, and comparing their results to the in-silico prediction program, RefFinder.

RESULTS

The most stable reference gene across both developmental stages was Rpp30, as predicted by Genorm, Bestkeeper, and Normfinder. RefFinder predicted Tbp to be the most stable across both developmental stages. At P14.5, stability varied by prediction program; at P20, RPP30 and MAPK1 were the most stable reference genes. Gapdh, 18S, Rplp0, and HPRT were predicted as the least stable reference genes by at least one of the prediction algorithms.

CONCLUSION

Expression of Rpp30 is the least affected by experimental conditions of oxygen induced retinopathy, phlebotomy induced anemia and erythropoietin administration at both timepoints of P14.5 and P20.

Precise measurement of gene expression changes in mouse brain areas denervated by injury

Quantitative PCR (qPCR) is a widely used method to study gene expression changes following brain injury. The accuracy of this method depends on the tissue harvested, the time course analyzed and, in particular on the choice of appropriate internal controls, i.e., reference genes (RGs). In the present study we have developed and validated an algorithm for the accurate normalization of qPCR data using laser microdissected tissue from the mouse dentate gyrus after entorhinal denervation at 0, 1, 3, 7, 14 and 28 days postlesion. The expression stabilities of ten candidate RGs were evaluated in the denervated granule cell layer (gcl) and outer molecular layer (oml) of the dentate gyrus. Advanced software algorithms demonstrated differences in stability for single RGs in the two layers at several time points postlesion. In comparison, a normalization index of several stable RGs covered the entire post-lesional time course and showed high stability. Using these RGs, we validated our findings and quantified glial fibrillary acidic protein (Gfap) mRNA and allograft inflammatory factor 1 (Aif1/Iba1) mRNA in the denervated oml. We compared the use of single RGs for normalization with the normalization index and found that single RGs yield variable results. In contrast, the normalization index gave stable results. In sum, our study shows that qPCR can yield precise, reliable, and reproducible datasets even under such complex conditions as brain injury or denervation, provided appropriate RGs for the model are used. The algorithm reported here can easily be adapted and transferred to any other brain injury model.

Selection of Ideal Reference Genes for Gene Expression Analysis in COVID-19 and Mucormycosis

Selection of reference genes during real-time quantitative PCR (qRT-PCR) is critical to determine accurate and reliable mRNA expression. Nonetheless, not a single study has investigated the expression stability of candidate reference genes to determine their suitability as internal controls in SARS-CoV-2 infection or COVID-19-associated mucormycosis (CAM). Using qRT-PCR, we determined expression stability of the nine most commonly used housekeeping genes, namely, TATA-box binding protein (TBP), cyclophilin (CypA), β-2-microglobulin (B2M), 18S rRNA (18S), peroxisome proliferator-activated receptor gamma (PPARG) coactivator 1 alpha (PGC-1α), glucuronidase beta (GUSB), hypoxanthine phosphoribosyltransferase 1 (HPRT-1), β-ACTIN, and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) in patients with COVID-19 of various severities (asymptomatic, mild, moderate, and severe) and those with CAM. We used statistical algorithms (delta-CT [threshold cycle], NormFinder, BestKeeper, GeNorm, and RefFinder) to select the most appropriate reference gene and observed that clinical severity profoundly influences expression stability of reference genes. CypA demonstrated the most consistent expression irrespective of disease severity and emerged as the most suitable reference gene in COVID-19 and CAM. Incidentally, GAPDH, the most commonly used reference gene, showed the maximum variations in expression and emerged as the least suitable. Next, we determined expression of nuclear factor erythroid 2-related factor 2 (NRF2), interleukin-6 (IL-6), and IL-15 using CypA and GAPDH as internal controls and show that CypA-normalized expression matches well with the RNA sequencing-based expression of these genes. Further, IL-6 expression correlated well with the plasma levels of IL-6 and C-reactive protein, a marker of inflammation. In conclusion, GAPDH emerged as the least suitable and CypA as the most suitable reference gene in COVID-19 and CAM. The results highlight the expression variability of housekeeping genes due to disease severity and provide a strong rationale for identification of appropriate reference genes in other chronic conditions as well. IMPORTANCE Gene expression studies are critical to develop new diagnostics, therapeutics, and prognostic modalities. However, accurate determination of expression requires data normalization with a reference gene, whose expression does not vary across different disease stages. Misidentification of a reference gene can produce inaccurate results. Unfortunately, despite the global impact of COVID-19 and an urgent unmet need for better treatment, not a single study has investigated the expression stability of housekeeping genes across the disease spectrum to determine their suitability as internal controls. Our study identifies CypA and then TBP as the two most suitable reference genes for COVID-19 and CAM. Further, GAPDH, the most commonly used reference gene in COVID-19 studies, turned out to be the least suitable. This work fills an important gap in the field and promises to facilitate determination of an accurate expression of genes to catalyze development of novel molecular diagnostics and therapeutics for improved patient care.

Identification of Suitable Reference Genes for Lowlanders Exposed to High Altitude and Ladakhi Highlanders

Sharma, Vandana, Rajeev Varshney, and Niroj Kumar Sethy. Identification of suitable reference genes for lowlanders exposed to high altitude and Ladakhi highlanders. High Alt Med Biol. 23:319-329, 2022. Background: Identifying a stable and reliable reference gene (RG) is a prerequisite for the unbiased and unambiguous analysis of gene expression data. It has become evident that conventionally used housekeeping genes such as beta-actin (ACTB), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and peptidylprolyl Isomerase A (PPIA) exhibit varied expression patterns under hypoxia. Hence, the identification of stable RGs for humans exposed to hypobaric hypoxia can enhance the accuracy of gene expression studies by limiting the negligent use of random housekeeping genes. Methods: Using TaqMan™ array-based quantitative real-time quantitative polymerase chain reaction, we evaluated the expression of 32 commonly used human RGs among lowlanders at Delhi (altitude 216 m, SL), lowlanders at Leh (altitude 3,524 m) after 1 day (HA-D1) and 7 days (HA-D7), as well as indigenous Ladakhi highlanders at the same altitude. The expression stability of the RGs was evaluated using geNorm, NormFinder, BestKeeper, Delta CT method, and RefFinder algorithms. Results: Our studies identify TATA-box binding protein (TBP), proteasome 26S subunit, ATPase 4 (PSMC4), and tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein zeta (YWHAZ) as the most stable human RGs for normalizing human gene expression under hypobaric hypoxia. In addition, we report the combination of TBP and cyclin-dependent kinase inhibitor 1B (CDKN1B) as the most stable RG for studying lowlander gene expression during high-altitude exposure. In contrast, RPL30 and 18S exhibited maximum variation across study groups and were identified as the least stable RGs.

Selection of species specific panel of reference genes in peripheral blood mononuclear cells of native livestock species adapted to trans-Himalayan region of Leh-Ladakh

The identification of appropriate references genes is an integral component of any gene expression-based study for getting accuracy and reliability in data interpretation. In this study, we evaluated the expression stability of 10 candidate reference genes (GAPDH, RPL4, EEF1A1, RPS9, HPRT1, UXT, RPS23, B2M, RPS15, ACTB) in peripheral blood mononuclear cells of livestock species that are adapted to high altitude hypoxia conditions of Leh-Ladakh. A total of 37 PBMCs samples from six native livestock species of Leh-Ladakh region such as Ladakhi cattle, Ladakhi yak, Ladakhi donkey, Chanthangi goat, Double hump cattle and Zanskar ponies were included in this study. The commonly used statistical algorithms such as geNorm, Normfinder, BestKeeper and RefFinder were employed to assess the stability of these RGs in all the livestock species. Our study has identified different panel of reference genes in each species; for example, EEF1A1, RPL4 in Ladakhi cattle; GAPDH, RPS9, ACTB in Ladakhi yak; HPRT1, B2M, ACTB in Ladakhi donkey; HPRT1, B2M, ACTB in Double hump camel, RPS9, HPRT1 in Changthangi goat, HPRT1 and ACTB in Zanskar ponies. To the best of our knowledge, this is the first systematic attempt to identify panel of RGs across different livestock species types adapted to high altitude hypoxia conditions. In future, the findings of the present study would be quite helpful in conducting any transcriptional studies to understand the molecular basis of high altitude adaptation of native livestock population of Leh-Ladakh.

Understanding quantitative polymerase chain reaction bioanalysis issues before validation planning: Japan Bioanalysis Forum discussion group

Investigating the biodistribution of cell and gene therapy products may play an important role in evaluating their safety and pharmacology. As quantitative polymerase chain reaction (qPCR) is often used for these analyses, it is essential to improve the reliability of bioanalysis performed using qPCR. In this report, the authors discuss the use of qPCR in nonclinical studies, as it can be used to detect target DNA/RNA and it is quantitative and applicable for long-term analysis. The authors also discuss points to consider during bioanalysis using qPCR and present appropriate validation items and their criteria. The authors anticipate the discussion provided herein to contribute to the development of validation and sample analysis for pharmaceuticals analyzed using qPCR.

Reference gene validation for the relative quantification of cannabinoid receptor expression in human odontoblasts via quantitative polymerase chain reaction

Objective

The aim of this study was to identify and validate the reference genes in cultured human odontoblasts to quantify their cannabinoid receptor transcripts.

Methods

The most stably transcribed genes in cultured human odontoblast cells were identified using the RefGenes tool and were selected for real-time polymerase chain reaction (PCR) amplification. Human odontoblast cells were differentiated from mesenchymal stem cells using a transforming growth factor-β-supplemented differentiation medium, and total RNA was purified. Reverse transcription-quantitative PCR and relative quantification analyses were performed using the Schefe's method. The relative expression dataset was analyzed to select the most stable genes.

Results

The analysis showed that the transcripts of cholinergic receptor nicotinic beta 2 subunit, LIM homeobox transcription factor 1 beta, and family with sequence similarity 223 member B presented the lowest standard deviation (SD) in expression (SD: 0.2, 0.17, and 0.16, respectively). These genes showed similar expression levels as the target genes (cannabinoid receptors). Significant differences were found in the relative expression levels of cannabinoid receptors using the selected genes compared to those calculated using beta actin transcripts as references (p < 0.05).

Conclusions

The strategy reported here for searching and verifying new reference genes will aid in the accurate and reliable expression of cannabinoid receptors in human odontoblast cells.

Selection and Validation of the Most Suitable Reference Genes for Quantitative Real-Time PCR Normalization in Salvia rosmarinus under In Vitro Conditions

Salvia rosmarinus L. (rosemary) is known to have a wide range of pharmacological effects including antidiabetic, anticarcinogenic, and antitumorigenic properties owing to its secondary metabolites. Studies aiming to elevate these metabolites have utilized various elicitors and stresses under in vitro conditions, although underlying molecular mechanisms remain unexplored. Gene expression studies using RT-qPCR might provide valuable information regarding how plant and plant cells interact and perceive various treatments and elicitors. However, despite being able to calculate accurate fold changes, the accuracy of the RT-qPCR data highly depends on the expression of reference genes. To the best of our knowledge, there is no information available on the stable reference genes in rosemary under in vitro conditions. Thus, in this paper, we assessed the stability of seven commonly used reference genes under different elicitor and stress conditions using RT-qPCR. Thereafter, the five most commonly used software and algorithms (comparative ΔCt, BestKeeper, NormFinder, geNorm, and RefFinder) were used to rank the candidates based on their expression stabilities. In conclusion, we recommend using a combination of F1-ATPase, ATP synthase and ACCase to normalize the gene expression experiments in rosemary under in vitro conditions. The selected reference genes were verified using 4-coumarate-CoA ligase, a pharmacologically important gene, whose expression might alter under nanoparticle treatment. Additionally, reference genes for several plant tissues, elicitors, and stresses are also proposed. The conclusions obtained from this current study will accelerate the future molecular work in S. rosmarinus and other related species.

The uncommon function and mechanism of the common enzyme glyceraldehyde-3-phosphate dehydrogenase in the metamorphosis of Helicoverpa armigera

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), a key enzyme in glycolysis, is commonly used as an internal reference gene in humans, mice, and insects. However, the function of GAPDH in insect development, especially in metamorphosis, has not been reported. In the present study, Helicoverpa armigera and Spodoptera frugiperda ovarian cell lines (Sf9 cells) were used as materials to study the function and molecular mechanism of GAPDH in larval metamorphosis. The results showed that HaGAPDH was more closely related to GAPDH of S. frugiperda and Spodoptera litura. The transcript peaks of HaGAPDH in sixth instar larvae were 6L-3 (epidermal and midgut) and 6L-1 (fat body) days, and 20E and methoprene significantly upregulated the transcripts of HaGAPDH of larvae in qRT-PCR. HaGAPDH–GFP–His was specifically localized in mitochondria in Sf9 cells. Knockdown of HaGAPDH by RNA interference (RNAi) in sixth instar larvae resulted in weight loss, increased mortality, and decreases in the pupation rate and emergence rates. HaGAPDH is directly bound to soluble trehalase (HaTreh1) physically and under 20E treatment in yeast two-hybrid, coimmunoprecipitation, and colocalization experiments. In addition, knockdown of HaGAPDH increased the Treh1 activity, which in turn decreased the trehalose content but increased the glucose content in larvae. Therefore, these data demonstrated that GAPDH controlled the glucose content within the normal range to ensure glucose metabolism and metamorphosis by directly binding with HaTreh1.

Identification of suitable reference genes for normalization of reverse transcription quantitative real-time PCR (RT-qPCR) in the fibrotic phase of the bleomycin mouse model of pulmonary fibrosis

Idiopathic pulmonary fibrosis (IPF) is a severe lung disease with a poor prognosis and few treatment options. In the most widely used experimental model for this disease, bleomycin is administered into the lungs of mice, causing a reaction of inflammation and consequent fibrosis that resembles the progression of human IPF. The inflammation and fibrosis together induce changes in gene expression that can be analyzed with reverse transcription quantitative real-time PCR (RT-qPCR), in which accurate normalization with a set of stably expressed reference genes is critical for obtaining reliable results. This work compares ten commonly used candidate reference genes in the late, fibrotic phase of bleomycin-induced pulmonary fibrosis and ranks them from the most to the least stable using NormFinder and geNorm. Sdha, Polr2a and Hprt were identified as the best performing and least variable reference genes when alternating between normal and fibrotic conditions. In order to validate the findings, we investigated the expression of Tnf and Col1a1, representing the hallmarks of inflammation and fibrotic changes, respectively. With the best three genes as references, both were found to be upregulated relative to untreated controls, unlike the situation when analyzed solely with Gapdh, a commonly used reference gene. We therefore recommend Sdha, Polr2a and Hprt as reference genes for RT-qPCR in the 4-week bleomycin challenge that represents the late fibrotic phase.

Copy Number Variations in the MICALL2 and MOGAT2 Genes Are Associated with Ashidan Yak Growth Traits

Copy number variations (CNVs) are a result of genomic rearrangement affecting DNA regions over 1 kb in length, and can include inversions, translocations, deletions, and duplications. The molecule interacting with CasL-like protein 2 (MICALL2) gene is primarily associated with mitochondrial protein targeting and exhibits predicted stress fiber colocalization. The monoacylglycerol O-acyltransferase 2 (MOGAT2) gene encodes an enzyme responsible for catalyzing diacylglycerol synthesis from 2-monoacylglycerol and fatty acyl-CoA. For this study, blood samples were obtained from 315 yaks, and the body weight, body length, withers height, and chest girth of these animals were measured at 6, 12, 18, and 30 months of age. Genomic DNA was harvested from the collected blood samples, and CNVs in these samples were detected by qPCR. The resultant data were compared using ANOVAs, revealing significant associations between MICALL2 gene CNVs and body weight at 6 months of age (p < 0.05), body length and chest girth at 30 months of age (p < 0.05), and withers height at 18 months of age (p < 0.01) in Ashidan yaks. Similarly, MOGAT2 CNVs were significantly associated with body weight at 6 and 30 months of age (p < 0.05), and with withers height at 18 months of age (p < 0.01) in these Ashidan yaks. MICALL2 and MOGAT2 gene expression was further analyzed in yak tissue samples, revealing that MICALL2 was most highly expressed in the adipose tissue, whereas MOGAT2 was most highly expressed in the lung. These results thus confirmed the relationship between CNVs in the MICALL2 and MOGAT2 genes and Ashidan yak growth traits, providing a valuable gene locus that can be leveraged for future marker-assisted yak breeding efforts.

Identification of qPCR reference genes suitable for normalising gene expression in the developing mouse embryo

Background: Progression through mammalian embryogenesis involves many interacting cell types and multiple differentiating cell lineages. Quantitative polymerase chain reaction (qPCR) analysis of gene expression in the developing embryo is a valuable tool for deciphering these processes, but normalisation to stably-expressed reference genes is essential for such analyses. Gene expression patterns change globally and dramatically as embryonic development proceeds, rendering identification of consistently appropriate reference genes challenging. Methods: We have investigated expression stability in mouse embryos from mid to late gestation (E11.5-E18.5), both at the whole-embryo level, and within the head and forelimb specifically, using 15 candidate reference genes ( ACTB, 18S, SDHA, GAPDH, HTATSF1, CDC40, RPL13A, CSNK2A2, AP3D1, HPRT1, CYC1, EIF4A, UBC, B2M and PAK1IP1), and four complementary algorithms (geNorm, Normfinder, Bestkeeper and deltaCt). Results: Unexpectedly, all methods suggest that many genes within our candidate panel are acceptable references, though AP3D1, RPL13A and PAK1IP1 are the strongest performing genes overall (scoring highly in whole embryos, heads or forelimbs alone, and in all samples collectively). HPRT1 and B2M are conversely poor choices, and show strong developmental regulation. We further show that normalisation using our three highest-scoring references can reveal subtle patterns of developmental expression even in genes ostensibly ranked as acceptably stable ( CDC40, HTATSF1). Conclusion: AP3D1, RPL13A and PAK1IP1 represent universally suitable reference genes for expression studies in the E11.5-E18.5 mouse embryo.