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

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.

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.

Genome-Wide Screening and Stability Verification of the Robust Internal Control Genes for RT-qPCR in Filamentous Fungi

In real-time quantitative PCR (RT-qPCR), internal control genes (ICGs) are crucial for normalization. This study screened 6 novel ICGs: Pre-mRNA-splicing factor cwc15 (Cwf15); ER associated DnaJ chaperone (DnaJ); E3 ubiquitin-protein ligase NEDD4 (HUL4); ATP-binding cassette, subfamily B (MDR/TAP), member 1 (VAMP); Exosome complex exonuclease DIS3/RRP44 (RNB); V-type H⁺-transporting ATPase sub-unit A (V-ATP) from the 22-transcriptome data of 8 filamentous fungi. The six novel ICGs are all involved in the basic biological process of cells and share the different transcription levels from high to low. In order to further verify the stability of ICGs candidates, the six novel ICGs as well as three traditional housekeeping genes: β-actin (ACTB); β-tubulin (β-TUB); glyceraldehyde-3-phosphate dehydrogenase gene (GAPDH) and the previously screened reference genes: SPRY-domain-containing protein (SPRYp); Ras-2 protein (Ras); Vacuolar protein sorting protein 26 (Vps26) were evaluated by geNorm and NormFinder statistical algorithms. RT-qPCR of 12 ICGs were performed at different developmental stages in Flammulina filiformis and under different treatment conditions in Neurospora crassa. The consistent results of the two algorithms suggested that the novel genes, RNB, V-ATP, and VAMP, showed the highest stability in F. filiformis and N. crassa. RNB, V-ATP, and VAMP have high expression stability and universal applicability and therefore have great potential as ICGs for standardized calculation in filamentous fungi. The results also provide a novel guidance for the screening stable reference genes in RT-qPCR and a wide application in gene expression analysis of filamentous fungi.

Reference gene identification for normalisation of RT-qPCR analysis in plasma samples of the rat middle cerebral artery occlusion model

OBJECTIVE

In quantitative reverse transcription-polymerase chain reaction (RT-qPCR) studies, the selection and validation of reference genes are crucial for the accurate analysis of MicroRNAs (miRNAs) expression. In this work, the optimal reference genes for RT-qPCR normalisation in plasma samples of rat middle cerebral artery occlusion (MCAO) models were identified.

METHODS

Six rat MCAO models were established. Blood samples were collected before modelling and approximately 16-24 h after modelling. Two commonly used reference genes (U6 and 5S) and three miRNAs (miR-24, miR-122 and miR-9a) were selected as candidate reference genes, and the expression of these genes was detected with RT-qPCR. The acquired data were analysed using geNorm, Normfinder, BestKeeper, RefFinder and comparative delta threshold cycle statistical models.

RESULTS

The analysed results consistently showed that miR-24 was the most stably expressed reference gene. The 'optimal combination' calculated by geNorm was miR-24, U6 and5S. The expression level of the target gene miR124 was similar when the most stable reference gene miR-24 or the 'optimal combination' was used as a reference gene. However, compared with miR24 or the 'optimal combination', the less stable reference genes influenced the fold change and the data accuracy with a large standard deviation.

CONCLUSION

These results confirmed the importance of selecting suitable reference genes for normalisation to obtain reliable results in RT-qPCR studies and demonstrated that the identified reference gene miR-24 or the 'optimal combination' could be used as an internal control for gene expression analysis in the rat MCAO model.

Validation of reference genes as an internal control for studying Avena sativa-Puccinia coronata interaction by RT-qPCR

In this study we evaluated eleven candidate reference genes in Avena sativa during compatible and incompatible interactions with two different pathotypes of Puccinia coronata f. sp. avenae in six time points post-inoculation. The identification of genes with high expression stability was performed by four algorithms (geNorm, NormFinder, BestKeeper and ΔCt method). The results obtained confirmed that the combination of two genes would be sufficient for reliable normalization of the expression data. In general, the most stable in the tested plant-pathogen system were HNR (heterogeneous nuclear ribonucleoprotein 27C) and EF1A (elongation factor 1-alpha). ARF (ADP-ribosylation factor) and EIF4A (eukaryotic initiation factor 4A-3) could also be considered as exhibiting high expression stability. CYP (cyclophilin) was shown by all assessment methods to be the worst candidate for normalization in this dataset. To date, this is the first report of reference genes selection in A. sativa–P. coronata interaction system. Identified reference genes enable reliable and comprehensive RT-qPCR analysis of oat gene expression in response to crown rust infection. Understanding the molecular mechanisms involved in the host–pathogen interactions may expand knowledge of durable resistance strategies beneficial to modern oat breeding.

Reference Genes across Nine Brain Areas of Wild Type and Prader-Willi Syndrome Mice: Assessing Differences in Igfbp7, Pcsk1, Nhlh2 and Nlgn3 Expression

Prader−Willi syndrome (PWS) is a complex neurodevelopmental disorder caused by the deletion or inactivation of paternally expressed imprinted genes at the chromosomal region 15q11−q13. The PWS-critical region (PWScr) harbors tandemly repeated non-protein coding IPW-A exons hosting the intronic SNORD116 snoRNA gene array that is predominantly expressed in brain. Paternal deletion of PWScr is associated with key PWS symptoms in humans and growth retardation in mice (PWScr model). Dysregulation of the hypothalamic−pituitary axis (HPA) is thought to be causally involved in the PWS phenotype. Here we performed a comprehensive reverse transcription quantitative PCR (RT-qPCR) analysis across nine different brain regions of wild-type (WT) and PWScr mice to identify stably expressed reference genes. Four methods (Delta Ct, BestKeeper, Normfinder and Genorm) were applied to rank 11 selected reference gene candidates according to their expression stability. The resulting panel consists of the top three most stably expressed genes suitable for gene-expression profiling and comparative transcriptome analysis of WT and/or PWScr mouse brain regions. Using these reference genes, we revealed significant differences in the expression patterns of Igfbp7, Nlgn3 and three HPA associated genes: Pcsk1, Pcsk2 and Nhlh2 across investigated brain regions of wild-type and PWScr mice. Our results raise a reasonable doubt on the involvement of the Snord116 in posttranscriptional regulation of Nlgn3 and Nhlh2 genes. We provide a valuable tool for expression analysis of specific genes across different areas of the mouse brain and for comparative investigation of PWScr mouse models to discover and verify different regulatory pathways affecting this complex disorder.

microRNA signatures associated with fetal growth restriction: a systematic review

Placental-origin microRNA (miRNA) profiles can be useful toward early diagnosis and management of fetal growth restriction (FGR) and associated complications. We conducted a systematic review to identify case-control studies that have examined miRNA signatures associated with human FGR. We systematically searched PubMed and ScienceDirect databases for relevant articles and manually searched reference lists of the relevant articles till May 18th, 2021. Of the 2133 studies identified, 21 were included. FGR-associated upregulation of miR-210 and miR-424 and downregulation of a placenta-specific miRNA cluster miRNA located on C19MC (miR-518b, miR-519d) and miR-221-3p was reported by >1 included studies. Analysis of the target genes of these miRNA as well as pathway analysis pointed to the involvement of angiogenesis and growth signaling pathways, such as the phosphatidylinositol 3-kinase- protein kinase B (PI3K-Akt) pathway. Only 3 out of the 21 included studies reported FGR-associated miRNAs in matched placental and maternal blood samples. We conclude that FGR-associated placental miRNAs could be utilized to inform clinical practice towards early diagnosis of FGR, provided enough evidence from studies on matched placental and maternal blood samples become available.Prospective Register of Systematic Reviews (PROSPERO) registration number: CRD42019136762.

Metabolic pathways for biosynthesis and degradation of starch in Tetraselmis chui during nitrogen deprivation and recovery

Tetraselmis chui is known to accumulate starch when subjected to stress. This phenomenon is widely studied for the purpose of industrial production and process development. Yet, knowledge about the metabolic pathways involved is still immature. Hence, in this study, transcription of 27 starch-related genes was monitored under nitrogen deprivation and resupply in 25 L tubular photobioreactors. T. chui proved to be an efficient starch producer under nitrogen deprivation, accumulating starch up to 56% of relative biomass content. The prolonged absence of nitrogen led to an overall down-regulation of the tested genes, in most instances maintained even after nitrogen replenishment when starch was actively degraded. These gene expression patterns suggest post-transcriptional regulatory mechanisms play a key role in T. chui under nutrient stress. Finally, the high productivity combined with an efficient recovery after nitrogen restitution makes this species a suitable candidate for industrial production of high-starch biomass.

Factors influencing degradation kinetics of mRNAs and half-lives of microRNAs, circRNAs, lncRNAs in blood in vitro using quantitative PCR

RNAs are rapidly degraded in samples and during collection, processing and testing. In this study, we used the same method to explore the half-lives of different RNAs and the influencing factors, and compared the degradation kinetics and characteristics of different RNAs in whole blood and experimental samples. Fresh anticoagulant blood samples were incubated at room temperature for different durations, RNAs were extracted, and genes, including internal references, were amplified by real-time quantitative PCR. A linear half-life model was established according to cycle threshold (Ct) values. The effects of experimental operations on RNA degradation before and after RNA extraction were explored. Quantitative analysis of mRNA degradation in samples and during experimental processes were explored using an orthogonal experimental design. The storage duration of blood samples at room temperature had the greatest influence on RNA degradation. The half-lives of messenger RNAs (mRNAs) was 16.4 h. The half-lives of circular RNAs (circRNAs), long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) were 24.56 ± 5.2 h, 17.46 ± 3.0 h and 16.42 ± 4.2 h, respectively. RNA degradation occurred mainly in blood samples. The half-life of mRNAs was the shortest among the four kinds of RNAs. Quantitative experiments related to mRNAs should be completed within 2 h. The half-lives of circRNAs and lncRNAs were longer than those of the former two.

Identification and validation of suitable housekeeping genes for gene expression studies in BCR-ABL1 positive B-lineage acute lymphoblastic leukemia

BACKGROUND

The stability of the housekeeping gene (HKG) expression is an absolute prerequisite for accurate normalization of target gene expression in a quantitative real-time polymerase chain reaction (RQ-PCR). In RQ-PCR, the widely used normalization approach involves the standardization of target genes to the most stable HKG control genes. According to the recent literature, in different experimental conditions the HKGs exhibit either up or down-regulation and thus affecting the gene expression profiles of target genes which leads to erroneous results. This implies that it is very important to select the appropriate HKG and verify the expression stability of the HKG before quantification of the target gene.

METHODS AND RESULTS

The present study aims to analyze six different HKGs for their expression profiles and stability in BCR-ABL1 negative cases and validate them in BCR-ABL1 positive cases, detected by multiplex reverse transcribed polymerase chain reaction (RT-PCR). Six commonly used reference genes (GAPDH, ABL1, RNA18S, ACTB, GUSB, and EEF2) were selected in this study. RQ-PCR was performed on 24 BCR-ABL1 negative cases and the outcomes were validated on 24 BCR-ABL1 positive cases. RefFinder™, a web-based composite software was used to check the stability of HKG genes by different algorithms and comprehensive ranking of each HKG gene in BCR-ABL1 negative cases and finally validated in BCR-ABL1 positive cases.

CONCLUSIONS

It was found that RNA18S, ABL1 and GUSB are good stable HKG genes, which showed minimum variability in gene expression compared to GAPDH, EEF2, and ACTB, the most commonly used HKG.

RNA-Seq is not required to determine stable reference genes for qPCR normalization

Assessment of differential gene expression by qPCR is heavily influenced by the choice of reference genes. Although numerous statistical approaches have been proposed to determine the best reference genes, they can give rise to conflicting results depending on experimental conditions. Hence, recent studies propose the use of RNA-Seq to identify stable genes followed by the application of different statistical approaches to determine the best set of reference genes for qPCR data normalization. In this study, however, we demonstrate that the statistical approach to determine the best reference genes from commonly used conventional candidates is more important than the preselection of ‘stable’ candidates from RNA-Seq data. Using a qPCR data normalization workflow that we have previously established; we show that qPCR data normalization using conventional reference genes render the same results as stable reference genes selected from RNA-Seq data. We validated these observations in two distinct cross-sectional experimental conditions involving human iPSC derived microglial cells and mouse sciatic nerves. These results taken together show that given a robust statistical approach for reference gene selection, stable genes selected from RNA-Seq data do not offer any significant advantage over commonly used reference genes for normalizing qPCR assays.

RTqPCR is a powerful technique that is widely used to quantify gene expression in research and diagnostics of different diseases. The technique involves making multiple copies (amplification) of a specific target DNA. The amplified target DNA binds to a molecule that emits fluorescence upon binding. The extent of fluorescence correlates to the amount of DNA present. To precisely quantify this fluorescence (and thus the quantities of target DNA), internal control genes also called as reference genes need to be determined. Such genes, in principle, do not have varied expression across samples and would exhibit the same fluorescence in all samples. They can thus be used to normalize the expression of the Target DNA. Unfortunately, choosing the right reference gene is very tricky and poor choice of reference genes results in unreliable data both in research and in diagnostics. In this study, we validate a statistical approach to find stably expressed reference genes for any experimental setting using a given set of candidates. We compare our approach to RNA sequencing which quantifies the expression of thousands of genes at the same time. We highlight the advantages of our approach which is cost effective and saves a lot of time when compared to sequencing.

[Current views on methods of determining the age of injuries in the practice of the forensic expert]

The review objective is to summarize the current data on new methods development for detection and determination of the age of injuries (AI) and to analyze the prospects of their use in the practice of forensic medical experts. The injury healing processes in various human body tissues are described in detail, and data on biomarkers of healing and their role are provided. Three main diagnostic methods for AI were analyzed: immunohistochemical, molecular biological study, and biophysical objectivization. Their advantages and disadvantages, as well as ways of further improvement of these methods, are considered.

Reexamining assumptions about miRNA-guided gene silencing

MicroRNAs (miRNAs) are short endogenously expressed RNAs that have the potential to regulate the expression of any RNA. This potential has led to the publication of several thousand papers each year connecting miRNAs to many different genes and human diseases. By contrast, relatively few papers appear that investigate the molecular mechanism used by miRNAs. There is a disconnect between rigorous understanding of mechanism and the extraordinary diversity of reported roles for miRNAs. Consequences of this disconnect include confusion about the assumptions underlying the basic science of human miRNAs and slow development of therapeutics that target miRNAs. Here, we present an overview of investigations into miRNAs and their impact on gene expression. Progress in our understanding of miRNAs would be aided by a greater focus on the mechanism of miRNAs and a higher burden of evidence on researchers who seek to link expression of a particular miRNA to a biological phenotype.

Establishment of reference (housekeeping) genes via quantitative real-time PCR for investigation of the genomic basis of abiotic stress resistance in Psammochloa villosa (Poaceae)

Psammochloa villosa is a desert plant growing in Northwest China with considerable resistance to abiotic stress, including drought, cold, and salt. To facilitate future studies of stress resistance in Psammochloa villosa, we sought to establish a suite of reference (or housekeeping) genes for utilization within future gene expression studies. Specifically, we selected nine candidate genes based on prior studies and new transcriptomic data for P. villosa, and we evaluated their expression stability in three different tissues of P. villosa under different treatments simulating abiotic stress conditions using four different bioinformatics assessments. Our results showed that TIP41 (TIP41-like family protein) was the most stable reference gene in drought- and salt-stressed leaves and salt-stressed stems, ELF-1α (elongation factor 1-α) was the most stable in cold-stressed leaves and drought- and salt-stressed roots, ACT (actin) was the most stable in drought-stressed stems, TUA (α-tubulin) was the most stable in cold-stressed stems, and 18S rRNA (18S ribosomal RNA) was the most stable in cold-stressed roots. Additionally, we tested the utility of these candidate reference genes to detect the expression pattern of P5CS (Δ¹-pyrroline-5-carboxylate synthetase), which is a drought-related gene. This study is the first report on selecting and validating reference genes of P. villosa under various stress conditions and will benefit future investigations of the genomic mechanisms of stress resistance in this ecologically important species.

Effects of neonatal dopaminergic lesion on oral cocaine self-administration in rats: Higher female vulnerability to cocaine consumption

The dopaminergic system is associated with cocaine-seeking behaviors, being influenced by other neurotransmitters such as GABA and deregulated by chronic cocaine self-administration. Administration of 6-hydroxydopamine (6-OHDA) to neonatal rats produces a depletion of brain dopamine, mainly, that results in behavioral alterations in adulthood. This model can be applied to better understanding of the role of the dopaminergic system in cocaine use and how its behavioral effects can modulate drug intake. Though there are well-established sex differences in the pattern of drug use, there are no published studies investigating sex-dependent effects of neonatal lesions with 6-OHDA on cocaine self-administration nor regarding GABA_A receptor (GABA_AR) subunits expression. Herein, neurotoxic lesion was induced in male and female neonatal rats by intracisternal injection of 6-OHDA at PND 4, and locomotion was evaluated before and after cocaine self-administration. Cocaine was diluted in a sweet solution (sucrose 1.5%) and offered for 27 consecutive 3-h daily sessions via a dispenser for oral intake, in an operant chamber under a fixed-ratio 1 (FR1) schedule. The 6-OHDA lesion reduced oral cocaine self-administration in male and female rats. Female rats, independent of dopaminergic condition, consumed more cocaine-containing solution than sucrose-only solution. Furthermore, as expected, 6-OHDA-lesioned animals presented a higher basal locomotor activity when compared to sham rats. We evaluated GABA_AR subunit expression and found no statistically significant differences between rats that self-administered a sucrose-only solution and those that self-administered a cocaine-containing solution. Even when the reward system is depleted, some behavioral differences remain in females, providing more data that highlight the female vulnerability to cocaine consumption.

Evaluation of Reference Genes in Glenea cantor (Fabricius) by Using qRT-PCR

Kapok is the main host of Glenea cantor (Fabricius), which causes serious damage and is difficult to control. In severe cases, it often causes the kapok trees to die continuously, which seriously affects the results of urban landscaping. To provide reference for the functional research on related genes in G. cantor, we screened the stable expression of candidate reference genes at different developmental stages (i.e., eggs, larvae, pupae, and adults), in various adult tissues (i.e., head, thorax, abdomen, feet, antennae, and wings), and sexes (i.e., male pupae, female pupae, male adults, and female adults). In this study, 12 candidate reference genes (i.e., ACTINLIKE, ACTININ, TUB, RPL36, RPL32, RPS20, TBP, GAPDH, 18S rRNA, EF1A1, EF1A2, and UBQ) were evaluated using different adult tissues, developmental stages, and sexes. RefFinder, geNorm, NormFinder, and BestKeeper were used to evaluate and comprehensively analyze the stability of the expression of the candidate reference genes. The results show that RPL32 and EF1A1 were the most suitable reference genes in the different adult tissues, and RPL36 and EF1A1 were best at the different developmental stages. RPL36 and EF1A2 were the best fit for the qRT-PCR reference genes in the different sexes, while RPL36 and EF1A1 were the most appropriate qRT-PCR reference genes in all samples. Results from geNorm showed that the optimal number of reference genes was two. We also surveyed the expression of cellulase at the different developmental stages and in the different adult tissues. Results further verified the reliability of the reference genes, and confirmed the best reference genes under the different experimental conditions. This study provides a useful tool for molecular biological studies on G. cantor.

ACTB and GAPDH appear at multiple SDS-PAGE positions, thus not suitable as reference genes for determining protein loading in techniques like Western blotting

We performed polyacrylamide gel electrophoresis of human proteins with sodium dodecyl sulfate, isolated proteins at multiple positions, and then used liquid chromatography and tandem mass spectrometry (LC-MS/MS) to determine the protein identities. Although beta-actin (ACTB) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) are 41.7 and 36 kDa proteins, respectively, LC-MS/MS identified their peptides at all the positions studied. The National Center for Biotechnology Information (USA) database lists only one ACTB mRNA but five GAPDH mRNAs and one noncoding RNA. The five GAPDH mRNAs encode three protein isoforms, while our bioinformatics analysis identified a 17.6 kDa isoform encoded by the noncoding RNA. All LC-MS/MS-identified GAPDH peptides at all positions studied are unique, but some of the identified ACTB peptides are shared by ACTC1, ACTBL2, POTEF, POTEE, POTEI, and POTEJ. ACTC1 and ACTBL2 belong to the ACT family with significant similarities to ACTB in protein sequence, whereas the four POTEs are ACTB-containing chimeric genes with the C-terminus of their proteins highly similar to the ACTB. These data lead us to conclude that GAPDH and ACTB are poor reference genes for determining the protein loading in such techniques as Western blotting, a leading role these two genes have been playing for decades in biomedical research.

Modulating OPG and TGF-β1 mRNA expression via bioelectrical stimulation

BACKGROUND

Bone remodeling is a lifelong process that ranges from orthodontic tooth movement/alignment to bone damage/healing, to overall bone health. Osteoprotegerin (OPG) and transforming growth factor β1 (TGF-β1) are secreted by osteoblasts and participate in bone remodeling. OPG promotes bone remineralization and stabilization prominent in post-mechanical repositioning of the teeth in the dental alveolus. TGF-β1 participates in regulatory processes to promote osteoblast and osteoclast equilibrium. In the context of orthodontic tooth movement, post-treatment fixation requires additional, exogenous, stabilization support. Recent research showcases supplementary solutions, in conjunction to standard tooth fixation techniques, such as OPG injections into gum and periodontal tissues to accelerate tooth anchorage; however, injections are prone to post-procedure complications and discomfort. This study utilizes noninvasive bioelectric stimulation (BES) to modulate OPG and TGF-β1 as a novel solution to regulate bone remineralization specifically in the context of post-orthodontic tooth movement.

PURPOSE

The aim of this study was to investigate a spectrum of BES parameters that would modulate OPG and TGF-β1 expression in osteoblasts.

METHODS

Osteoblasts were cultured and stimulated using frequencies from 25 Hz to 3 MHz. RT-qPCR was used to quantify changes in OPG and TGFb-1 mRNA expression.

RESULTS

OPG mRNA expression was significantly increased at frequencies above 10,000 Hz with a maximum expression increase of 332 ± 8% at 100 kHz. Conversely, OPG mRNA expression was downregulated at frequencies lower than 1000 Hz. TGF-β1 mRNA expression increased throughout all stimulation frequencies with a peak of 332 ± 72% at 250 kHz. Alizarin Red tests for calcium, indicated that mineralization of stimulated osteoblasts in vitro increased 28% after 6 weeks in culture.

DISCUSSION

Results support the working hypothesis that OPG and TGF-β1 mRNA expression can be modulated through BES. Noninvasive BES approaches have the potential to accelerate bone remineralization by providing a novel tool to supplement the anchorage process, reduce complications, and promote patient compliance and reduce post-treatment relapse. Noninvasive BES may be applicable to other clinical applications as a novel therapeutic tool to modulate bone remodeling.

Reference gene screening for analyzing gene expression in the heart, liver, spleen, lung and kidney of forest musk deer

The screening of reference genes for real-time quantitative PCR (qPCR) in forest musk deer (FMD) tissue is of great significance to the basic research on FMD. However, there are few reports on the stability analysis of FMD reference genes so far. In this study, We used qPCR to detect the expression levels of 11 reference gene candidates (18S rRNA, beta-actin [ACTB], glyceraldehyde-3-phosphate dehydrogenase [GAPDH], TATA box-binding protein [TBP], hypoxanthine phosphoribosyltransferase 1 [HPRT1], tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein zeta polypeptide [YWHAZ], hydroxymethylbilane synthase [HMBS], eukaryotic translation elongation factor 1 alpha 1 [EEF1A1], succinate dehydrogenase complex flavoprotein subunit A [SDHA], peptidylprolyl isomerase B [PPIB], and ubiquitin C [UBC]) in heart, liver, spleen, lung and kidney of FMD. After removing 18S rRNA on account of its high expression level, geNorm, NormFinder, BestKeeper and ΔCt algorithms were used to evaluate the expression stability of the remaining genes in the five organs, and further comprehensive ranking was calculated by RefFinder. According to the results, the selected reference genes with the most stable expression in the heart of FMD are SDHA and YWHAZ, while in the liver are ACTB and SDHA; in the spleen and lung are YWHAZ and HPRT1; in the kidney are YWHAZ and PPIB. The use of common reference genes in all five organs is not recommended. The analyses showed that tissue is an important variability factor in genes expression stability. Meanwhile, the result can be used as a reference for the selection of reference genes for qPCR in further study.

A comprehensive genome-wide analysis of long non-coding RNA and mRNA expression profiles of JAK2V617F-positive classical myeloproliferative neoplasms

Aberrant expression of long non-coding RNAs (lncRNAs) is involved in the progression of myeloid neoplasms, but the role of lncRNAs in the JAK2V617F-positive subtype of classical myeloproliferative neoplasms (cMPNs) remains unclear. This study was conducted to clarify the expression and regulation patterns of lncRNAs in JAK2V617F-positive cMPNs, and to explore new potential carcinogenic factors of cMPNs. Bioinformatics analysis of microarray detection and wet testing verification were performed to study the expression and regulation signature of differentially expressed lncRNAs (DELs) and related genes (DEGs) in cMPNs. The expression of lncRNAs and mRNAs were observed to significantly dysregulated in JAK2V617F-positive cMPN patients compared with the normal controls. Co-expression analysis indicated that there were significant differences of the co-expression pattern of lncRNAs and mRNAs in JAK2V617F-positive cMPN patients compared to normal controls. GO and KEGG pathway analysis of DEGs and DELs showed the involvement of several pathways previously reported to regulate the pathogenesis of leukemia and cMPNs. Cis- and trans-regulation analysis of lncRNAs showed that ZNF141, DHX29, NOC2L, MAS1L, AFAP1L1, and CPN2 were significantly cis-regulated by lncRNA ENST00000356347, ENST00000456816, hsa-mir-449c, NR_026874, TCONS_00012136, uc003lqp.2, and ENST00000456816, respectively, and DELs were mostly correlated with transcription factors including CTBP2, SUZ12, REST, STAT2, and GATA4 to jointly regulate multiple target genes. In summary, expression profiles of lncRNAs and mRNAs were significantly altered in JAK2V617F-positive cMPNs, the relative signaling pathway, co-expression, cis- and trans-regulation were regulated by dysregulation of lncRNAs and several important genes, such as ITGB3, which may act as a promising carcinogenic factor, warrant further investigation.

Genome-wide screening of novel RT-qPCR reference genes for study of GLRaV-3 infection in wine grapes and refinement of an RNA isolation protocol for grape berries

BACKGROUND

Grapevine, as an essential fruit crop with high economic values, has been the focus of molecular studies in diverse areas. Two challenges exist in the grapevine research field: (i) the lack of a rapid, user-friendly and effective RNA isolation protocol for mature dark-skinned berries and, (ii) the lack of validated reference genes that are stable for quantification of gene expression across desired experimental conditions. Successful isolation of RNA with sufficient yield and quality is essential for downstream analyses involving nucleic acids. However, ripe berries of dark-skinned grape cultivars are notoriously challenging in RNA isolation due to high contents of polyphenolics, polysaccharides, RNase and water.

RESULTS

We have optimized an RNA isolation protocol through modulating two factors at the lysis step that could impact results of RNA isolation - 2-ME concentration and berry mass. By finding the optimal combination among the two factors, our refined protocol was highly effective in isolating total RNA with high yield and quality from whole mature berries of an array of dark-skinned wine grape cultivars. Our protocol takes a much shorter time to complete, is highly effective, and eliminates the requirement for hazardous organic solvents. We have also shown that the resulting RNA preps were suitable for multiple downstream analyses, including the detection of viruses and amplification of grapevine genes using reverse transcription-polymerase chain reaction (RT-PCR), gene expression analysis via quantitative reverse transcription PCR (RT-qPCR), and RNA Sequencing (RNA-Seq). By using RNA-Seq data derived from Cabernet Franc, we have identified seven novel reference gene candidates (CYSP, NDUFS8, YLS8, EIF5A2, Gluc, GDT1, and EF-Hand) with stable expression across two tissue types, three developmental stages and status of infection with grapevine leafroll-associated virus 3 (GLRaV-3). We evaluated the stability of these candidate genes together with two conventional reference genes (actin and NAD5) using geNorm, NormFinder and BestKeeper. We found that the novel reference gene candidates outperformed both actin and NAD5. The three most stable reference genes were CYSP, NDUFS8 and YSL8, whereas actin and NAD5 were among the least stable. We further tested if there would be a difference in RT-qPCR quantification results when the most stable (CYSP) and the least stable (actin and NAD5) genes were used for normalization. We concluded that both actin and NAD5 led to erroneous RT-qPCR results in determining the statistical significance and fold-change values of gene expressional change.

CONCLUSIONS

We have formulated a rapid, safe and highly effective protocol for isolating RNA from recalcitrant berry tissue of wine grapes. The resulting RNA is of high quality and suitable for RT-qPCR and RNA-Seq. We have identified and validated a set of novel reference genes based on RNA-Seq dataset. We have shown that these new reference genes are superior over actin and NAD5, two of the conventional reference genes commonly used in early studies.

Case Study Using Recommended Reference Genes Actin and 18S for Reverse-Transcription Quantitative Real-Time PCR Analysis in Myzus persicae

Myzus persicae is a globally important pest with the ability to adjust to a wide range of environmental situations, and many molecular technologies have been developed and applied to understand the biology and/or control this pest insect directly. Reverse-transcription quantitative real-time PCR (RT-qPCR) is a primary molecular technology that is used to quantify gene expression. Choosing a stable reference gene is significantly important for precisely clarifying the expression level of the target gene. Actin and 18S have been recommended as stable compounds for real-time RT-qPCR in M. persicae under the tested biotic and abiotic conditions. In this study, we checked the stability of Actin and 18S by analyzing the relative expression levels of the cytochrome 450 monooxygenase family member genes CYP6CY3 and CYP6-1, carboxylesterase gene E4 and vacuolar protein sorting gene VPS11 via RT-qPCR under various conditions. The expression levels of these four target genes were normalized using both Actin and 18S individually and the combination of these two genes. Our results confirmed that Actin and 18S can be used as reference genes to normalize the expression of target genes under insecticide treatment and starvation in M. persicae. However, at the developmental stages of M. persicae, the expression of the four tested target genes was normalized stably by Actin but not 18S, with the latter presenting a problematic change with the developmental stages. Thus, the stability of reference genes in response to diverse biotic and abiotic factors should be evaluated before each RT-qPCR experiment.

Selection and Evaluation of Reference Genes for qRT-PCR in Spodoptera frugiperda (Lepidoptera: Noctuidae)

Simple Summary

The fall armyworm Spodoptera frugiperda (J. E. Smith) (Lepidoptera: Noctuidae) is an extremely important omnivorous agricultural pest, it poses a severe threat to food security and agricultural production. Quantitative real-time PCR (qRT-PCR) is an important molecular technology widely used for expression profile analyses of various target genes. It is essential to use reference genes as the benchmark to eliminate various errors and normalize the qRT-PCR analysis. In our study, 10 reference genes were evaluated under six experimental conditions, including developmental stages, tissues, mating status, hormones, diets, and temperatures. Finally, the expression profile of the target gene SfrOBP1 in various tissues of S. frugiperda was evaluated to verify the accuracy of the results. This study will provide a preliminary evaluation of reference genes of S. frugiperda, which can be beneficial to the further research of functional gene expression.

Abstract

As an accurate and convenient technique, the qRT-PCR is always used in the quantitative expression analysis of functional genes. Normalization of the data relies on stable reference genes. The fall armyworm Spodoptera frugiperda (J. E. Smith) is an important invasive and migratory pest that seriously threatens corn production around the world. In this paper, we selected 10 candidate reference genes (18S, AK, RPL10, RPS24, 28S, SOD, ATP, GAPDH, ACT, and a-TUB) and determined their expression levels under different conditions (different developmental stages, various tissues, mating status, hormones, diets, and temperatures). Subsequently, the stability of reference genes was evaluated by four algorithms (Delta Ct method, geNorm, NormFinder, BestKeeper). The optimal combination of reference genes for each treatment was obtained by geNorm. Finally, the comprehensive ranks were determined by the online tool RefFinder. Results showed that the most stable reference genes were SOD, RPL10, and RPS24 for developmental stages, α-TUB, RPL10, and ATP for different tissues, AK, RPL10, and 18S for mating status, 18S and AK under hormone treatment, 18S, RPL10, and SOD under diet treatment, RPL10, 18S, and RPS24 under temperature treatment. This study confirmed recent data on a few reference genes and provided an evaluation of a number of additional reference genes of S. frugiperda under various conditions.

Constitutive Defense Strategy of Coffee Under Field Conditions: A Comparative Assessment of Resistant and Susceptible Cultivars to Rust

Coffea arabica is the most economically important coffee species worldwide. However, its production is severely limited by diseases such as rust. The mechanisms underlying constitutive defense responses in coffee are still poorly understood, compared with induced defense mechanisms. We aimed to characterize constitutive defense responses of thirteen cultivars of C. arabica. Cultivars were classified under field conditions according to the level of resistance to rust: resistant (R), moderately resistant (MR), and susceptible (S). Based on this classification, the stability of eight reference genes (RGs) was evaluated. The most stable RGs were EF1α, APT1, and 24S. We also evaluated the expression of CaWRKY1, CaPAL1, CaCAD1, and CaPOX1, and activities of PAL, CAD, and POX, which are involved in lignin biosynthesis, and leaf content of total phenolic compounds and lignin. Gene expression and enzymatic activity were not correlated with defense metabolites in the R cultivar group but showed a negative correlation with phenolic compounds in MR cultivars. Cultivar S showed positive correlations of gene expression and enzyme activity with phenolic compounds. These results may assist coffee breeding programs regarding selection of genotypes and in optimization of rust resistance.

Stability analysis of reference genes for RT-qPCR assays involving compatible and incompatible Ralstonia solanacearum-tomato 'Hawaii 7996' interactions

Reverse transcription-quantitative PCR (RT-qPCR) is an analytical tool for gene expression quantification. Reference genes are not yet available for gene expression analysis during interactions of Ralstonia solanacearum with ‘Hawaii 7996’ (the most stable source of resistance in tomato). Here, we carried out a multi-algorithm stability analysis of eight candidate reference genes during interactions of ‘Hawaii 7996’ with one incompatible/avirulent and two compatible/virulent (= resistance-breaking) bacterial isolates. Samples were taken at 24- and 96-h post-inoculation (HPI). Analyses were performed using the ∆∆Ct method and expression stability was estimated using BestKeeper, NormFinder, and geNorm algorithms. TIP41 and EF1α (with geNorm), TIP41 and ACT (with NormFinder), and UBI3 and TIP41 (with BestKeeper), were the best combinations for mRNA normalization in incompatible interactions at 24 HPI and 96 HPI. The most stable genes in global compatible and incompatible interactions at 24 HPI and 96 HPI were PDS and TIP41 (with geNorm), TIP41 and ACT (with NormFinder), and UBI3 and PDS/EXP (with BestKeeper). Global analyses on the basis of the three algorithms across 20 R. solanacearum-tomato experimental conditions identified UBI3, TIP41 and ACT as the best choices as reference tomato genes in this important pathosystem.

Housekeeping Genes for Parkinson's Disease in Humans and Mice

A critical aspect of real-time PCR is the presence of housekeeping genes (HKGs) as an internal control for the normalization of expression data for genes of interest. It is necessary to select correct HKGs in the investigation of various pathologies. Thereby, we analyzed the stability of expression of the HKGs in Parkinson’s disease (PD). The work was carried out in the peripheral blood of patients with PD and in the brain tissues and peripheral blood of mice with MPTP-induced PD. As a result, Aars was the most stably expressed HKG in the mouse brain as a whole. However, different genes were more stably expressed in different parts of the brain. Polr2f was the most stably expressed in the cortex, Psmd6 was the most stably expressed in the cerebellum, and Psmd7 was the most stably expressed in the striatum and substantia nigra. HKGs were different in similar tissues of the studied organisms. Polr2f was the most stably expressed HKG in the peripheral blood of mice, whereas PSMD6 was the most stably expressed gene in humans. Thus, there is no universal HKG both for different brain tissues of one organism and for similar tissues of different organisms. Furthermore, the identified most stably expressed HKGs can be considered as such only under conditions in PD.

Assessment of reference genes at six different developmental stages of Schistosoma mansoni for quantitative RT-PCR

Reverse-transcription quantitative real-time polymerase chain reaction (RT-qPCR) is the most used, fast, and reproducible method to confirm large-scale gene expression data. The use of stable reference genes for the normalization of RT-qPCR assays is recognized worldwide. No systematic study for selecting appropriate reference genes for usage in RT-qPCR experiments comparing gene expression levels at different Schistosoma mansoni life-cycle stages has been performed. Most studies rely on genes commonly used in other organisms, such as actin, tubulin, and GAPDH. Therefore, the present study focused on identifying reference genes suitable for RT-qPCR assays across six S. mansoni developmental stages. The expression levels of 25 novel candidates that we selected based on the analysis of public RNA-Seq datasets, along with eight commonly used reference genes, were systematically tested by RT-qPCR across six developmental stages of S. mansoni (eggs, miracidia, cercariae, schistosomula, adult males and adult females). The stability of genes was evaluated with geNorm, NormFinder and RefFinder algorithms. The least stable candidate reference genes tested were actin, tubulin and GAPDH. The two most stable reference genes suitable for RT-qPCR normalization were Smp_101310 (Histone H4 transcription factor) and Smp_196510 (Ubiquitin recognition factor in ER-associated degradation protein 1). Performance of these two genes as normalizers was successfully evaluated with females maintained unpaired or paired to males in culture for 8 days, or with worm pairs exposed for 16 days to double-stranded RNAs to silence a protein-coding gene. This study provides reliable reference genes for RT-qPCR analysis using samples from six different S. mansoni life-cycle stages.

Identification of Reference Genes for Circadian Studies on Brain Microvessels and Choroid Plexus Samples Isolated from Rats

Delivery of putative compounds of therapeutic value to the brain is limited by brain barriers: the blood–brain barrier located in the endothelium of the brain microvessels (BrMV) and the blood–cerebrospinal fluid barrier located in the epithelium of the choroid plexus (ChP). Understanding their function and modulation by the circadian clock may enhance the efficacy of brain-targeting therapies. The aim of the present study was to evaluate the stability of 10 reference genes in the BrMV and ChP, isolated from male and female rats at six time points (ZT1, 5, 9, 13, 17, and 21). Gene evaluations were performed by qPCR, analyzed by RefFinder tool, and verified by analyzing the expression of the brain and muscle ARNT-like 1 (Bmal1) using the qPCR and digital PCR methods. We identified as the most stable genes for circadian studies tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein zeta (Ywhaz) and apolipoprotein E (Apoe) for BrMV, and beta actin (Actb) and hypoxanthine-guanine phosphoribosyltransferase (Hprt1) for ChP. After verification, ribosomal protein (Rps18) was also included as a sufficient reference gene. Additionally, the observed gender difference in the Bmal1 oscillations in both BrMV and ChP suggests that separate studies for each gender are recommended.

Reference Gene Selection for Gene Expression Analyses in Mouse Models of Acute Lung Injury

qRT-PCR still remains the most widely used method for quantifying gene expression levels, although newer technologies such as next generation sequencing are becoming increasingly popular. A critical, yet often underappreciated, problem when analysing qRT-PCR data is the selection of suitable reference genes. This problem is compounded in situations where up to 25% of all genes may change (e.g., due to leukocyte invasion), as is typically the case in ARDS. Here, we examined 11 widely used reference genes for their suitability in commonly used models of acute lung injury (ALI): ventilator-induced lung injury (VILI), in vivo and ex vivo, lipopolysaccharide plus mechanical ventilation (MV), and hydrochloric acid plus MV. The stability of reference gene expression was determined using the NormFinder, BestKeeper, and geNorm algorithms. We then proceeded with the geNorm results because this is the only algorithm that provides the number of reference genes required to achieve normalisation. We chose interleukin-6 (Il-6) and C-X-C motif ligand 1 (Cxcl-1) as the genes of interest to analyse and demonstrate the impact of inappropriate normalisation. Reference gene stability differed between the ALI models and even within the subgroup of VILI models, no common reference gene index (RGI) could be determined. NormFinder, BestKeeper, and geNorm produced slightly different, but comparable results. Inappropriate normalisation of Il-6 and Cxcl1 gene expression resulted in significant misinterpretation in all four ALI settings. In conclusion, choosing an inappropriate normalisation strategy can introduce different kinds of bias such as gain or loss as well as under- or overestimation of effects, affecting the interpretation of gene expression data.

Reference genes for mesangial cell and podocyte qPCR gene expression studies under high-glucose and renin-angiotensin-system blocker conditions

BACKGROUND

Real-time PCR remains currently the gold standard method for gene expression studies. Identification of the best reference gene is a key point in performing high-quality qPCR, providing strong support for results, and performing as a source of bias when inappropriately chosen. Mesangial cells and podocytes, as essential cell lines to study diabetic kidney disease (DKD) physiopathology, demand accurate analysis of the reference genes used thus far to enhance the validity of gene expression studies, especially regarding high glucose (HG) and DKD treatments, with angiotensin II receptor blockers (e.g., losartan) being the most commonly used. This study aimed to evaluate the suitability and define the most stable reference gene for mesangial cell and podocyte studies of an in vitro DKD model of disease and its treatment.

METHODS

Five software packages (RefFinder, NormFinder, GeNorm, Bestkeeper, and DataAssist) and the comparative ΔCt method were selected to analyze six different candidate genes: HPRT, ACTB, PGAM-1, GAPDH, PPIA, and B2M. RNA was extracted, and cDNA was synthesized from immortalized mouse mesangial cells and podocytes cultured in 4 groups: control (n = 5; 5 mM glucose), mannitol (n = 5; 30 mM, as osmotic control), HG (n = 5; 30 mM glucose), and HG + losartan (n = 5; 30 mM glucose and 10-4 mM losartan). Real-time PCR was performed according to MIQE guidelines.

RESULTS

We identified that the use of 2 genes was the best combination for qPCR normalization for both mesangial cells and podocytes. For mesangial cells, the combination of HPRT and ACTB presented higher stability values. For podocytes, HPRT and GAPDH showed the best results.

CONCLUSION

This analysis provides support for the use of HPRT and ACTB as reference genes in mouse mesangial cell studies of gene expression via real-time PCR, while for podocytes, HPRT and GAPDH should be chosen.

Establishment of suitable reference genes for studying relative gene expression during the transition from trophozoites to cyst-like stages and first evidences of stress-induced expression of meiotic genes in Trichomonas vaginalis

Trichomonas vaginalis is a parasite of the human urogenital tract and the causative agent of trichomoniasis, a sexually transmitted disease of worldwide importance. This parasite is usually found as a motile flagellated trophozoite. However, when subjected to stressful microenvironmental conditions, T. vaginalis trophozoites can differentiate into peculiar cyst-like stages, which exhibit notable physiological resistance to unfavourable conditions. Although well documented in morphological and proteomic terms, patterns of gene expression changes involved in the cellular differentiation into cyst-like stages are mostly unknown. The real-time reverse transcription polymerase chain reaction (RT-qPCR) is recognized as a sensitive and accurate method for quantification of gene expression, providing fluorescence-based data that are proportional to the amount of a target RNA. However, the reliability of relative expression studies depends on the validation of suitable reference genes, which RNAs exhibit a minimum of variation between tested conditions. Here, we attempt to determine suitable reference genes to be used as controls of invariant expression during cold-induced in vitro differentiation of T. vaginalis trophozoites into cyst-like forms. Furthermore, we reveal that the mRNA from the meiotic recombinase Dmc1 is upregulated during this process, indicating that cryptic sexual events may take place in cyst-like stages of T. vaginalis.

Identification of the most suitable reference gene for gene expression studies with development and abiotic stress response in Bromus sterilis

Bromus sterilis is an annual weedy grass, causing high yield losses in winter cereals. Frequent use of herbicides had led to the evolution of herbicide resistance in this species. Mechanisms underlying herbicide resistance in B. sterilis must be uncovered because this problem is becoming a global threat. qRT-PCR and the next-generation sequencing technologies can elucidate the resistance mechanisms. Although qRT-PCR can calculate precise fold changes, its preciseness depends on the expression of reference genes. Regardless of stable expression in any given condition, no gene can act as a universal reference gene. Hence, it is necessary to identify the suitable reference gene for each species. To our knowledge, there are no reports on the suitable reference gene in any brome species so far. Thus, in this paper, the stability of eight genes was evaluated using qRT-PCR experiments followed by expression stability ranking via five most commonly used software for reference gene selection. Our findings suggest using a combination of 18S rRNA and ACCase to normalise the qRT-PCR data in B. sterilis. Besides, reference genes are also recommended for different experimental conditions. The present study outcomes will facilitate future molecular work in B. sterilis and other related grass species.

miR-361-5p as a promising qRT-PCR internal control for tumor and normal breast tissues

BACKGROUND

One of the most widely used evaluation methods in miRNA experiments is qRT-PCR. However, selecting suitable internal controls (IC) is crucial for qRT-PCR experiments. Currently, there is no consensus on the ICs for miRNA qRT-PCR experiments in breast cancer. To this end, we tried to identify the most stable (the least expression alteration) and promising miRNAs in normal and tumor breast tissues by employing TCGA miRNA-Seq data and then experimentally validated them on fresh clinical samples.

METHODS

A multi-component scoring system was used which takes into account multiple expression stability criteria as well as correlation with clinical characteristics. Furthermore, we extended the scoring system for more than two biological sub-groups. TCGA BRCA samples were analyzed based on two grouping criteria: Tumor & Normal samples and Tumor subtypes. The top 10 most stable miRNAs were further investigated by differential expression and survival analysis. Then, we examined the expression level of the top scored miRNA (hsa-miR-361-5p) along with two commonly used ICs hsa-miR-16-5p and U48 on 34 pairs of Primary breast tumor and their adjacent normal tissues using qRT-PCR.

RESULTS

According to our multi-component scoring system, hsa-miR-361-5p had the highest stability score in both grouping criteria and hsa-miR-16-5p showed significantly lower scores. Based on our qRT-PCR assay, while U48 was the most abundant IC, hsa-miR-361-5p had lower standard deviation and also was the only IC capable of detecting a significant up-regulation of hsa-miR-21-5p in breast tumor tissue.

CONCLUSIONS

miRNA-Seq data is a great source to discover stable ICs. Our results demonstrated that hsa-miR-361-5p is a highly stable miRNA in tumor and non-tumor breast tissue and we recommend it as a suitable reference gene for miRNA expression studies in breast cancer. Additionally, although hsa-miR-16-5p is a commonly used IC, it's not a suitable one for breast cancer studies.

Identification of reliable reference genes for expression studies in the magnum of laying hens housed in cage and cage-free systems

Stress factors during poultry production can evoke changes in gene transcription and protein synthesis in the hen oviduct and could affect the internal and external egg quality. Studies of relative gene expression require the identification of the most stable reference genes for the quantitative polymerase chain reaction (qPCR) to investigate the reproductive tissues' response in laying hens kept in different production systems. The objective of this study was to determine the most stable reference genes of the magnum tissues of laying hens housed in two different production systems. Hy-Line Brown hens were reared under the same sanitary conditions until 15 weeks of age. Later on, they were transferred into two different production systems, conventional cage (CC) and cage free (CF), up to 82 weeks of age. At 50 and 60 weeks, a total of six hens from each production system were euthanized, and magnum samples were collected. The qPCR was used to determine the RNA transcription level of five reference genes, ACTB, 18S, GAPDH, MSX2 and HMBS. These genes were evaluated for transcript expression in magnum tissues by NormFinder, BestKeeper, geNorm and RefFinder software. The results indicated that the most stable gene in the CF housing system was HMBS in three of the algorithms and in the CC housing system was the 18S, and the best combination of reference genes was HMBS/GAPDH in CF and 18S/HMBS in CC. In conclusion, HMBS, 18S and GAPDH genes could be used together as reference genes for the normalization of the magnum tissues transcript expression of hens in CF and CC housing systems.

Validation of reference genes for use in untreated bovine fibroblasts

Proper normalization of RT-qPCR data is pivotal to the interpretation of results and accuracy of scientific conclusions. Though different approaches may be taken, normalization against multiple reference genes is now standard practice. Genes traditionally used and deemed constitutively expressed have demonstrated variability in expression under different experimental conditions, necessitating the proper validation of reference genes prior to utilization. Considering the wide use of fibroblasts in research and scientific applications, it is imperative that suitable reference genes for fibroblasts of different animal origins and conditions be elucidated. Previous studies on bovine fibroblasts have tested limited genes and/or samples. Herein, we present an extensive study investigating the expression stability of 16 candidate reference genes across 7 untreated bovine fibroblast cell lines subjected to controlled conditions. Data were analysed using various statistical tools and algorithms, including geNorm, NormFinder, BestKeeper, and RefFinder. A combined use of GUSB and RPL13A was determined to be the best approach for data normalization in untreated bovine fibroblasts.

Selection and validation of appropriate reference genes for RT-qPCR analysis of flowering stages and different genotypes of Iris germanica L

Iris germanica L. is a perennial herbaceous plant that has been widely cultivated worldwide and is popular for its elegant and vibrantly colorful flowers. Selection of appropriate reference genes is the prerequisite for accurate normalization of target gene expression by quantitative real-time PCR. However, to date, the most suitable reference genes for flowering stages have not been elucidated in I. germanica. In this study, eight candidate reference genes were examined for the normalization of RT-qPCR in three I. germanica cultivars, and their stability were evaluated by four different algorithms (GeNorm, NormFinder, BestKeeper, and Ref-finder). The results revealed that IgUBC and IgGAPDH were the most stable reference genes in '00246' and 'Elizabeth', and IgTUB and IgUBC showed stable expression in '2010200'. IgUBC and IgGAPDH were the most stable in all samples, while IgUBQ showed the least stability. Finally, to validate the reliability of the selected reference genes, the expression patterns of IgFT (Flowering Locus T gene) was analyzed and emphasized the importance of appropriate reference gene selection. This work presented the first systematic study of reference genes selection during flower bud development and provided guidance to research of the molecular mechanisms of flowering stages in I. germanica.

Defining stable reference genes in HIV latency reversal experiments

Quantification of cell associated HIV RNA (ca-RNA) is one of the most important and commonly used methods to evaluate the performance of latency-reversing agents (LRAs). Copies of HIV RNA measured by qPCR, are often normalized to the input RNA or cell number. However, these could be affected by biological variability and/or technical errors, which can be avoided by using an internal reference gene. To obtain reliable data, it is essential to select stable reference genes (RGs) of which the expression is not influenced by biological variability, the type of cells, or the LRAs used. However, to date, no study has carefully evaluated RG stability following LRA exposure. We analyzed the stability of six widely used RGs (GAPDH, TBP, YWHAZ, UBE2D2, HPRT1 and RPL27A) in human PBMC and CD4+ T cells. LRA exposure significantly influenced the stability of these RGs. Overall, TBP, UBE2D2, and RPL27A were the most stable RGs in all tested conditions. TBP was generally the most stable RG whereas GAPDH varied the most. Finally, we evaluated the impact of applying different RG normalizers to host genes and HIV ca-RNA data. Altered results were observed both in host and HIV gene expression when unstable RGs were used. Our data underline the importance of testing the stability of RGs utilized to evaluate LRA-induced HIV ca-RNA expression. To our knowledge, this is the first careful evaluation of the stability of RGs after LRA exposure and will significantly contribute to the quality of data analysis in regard to gene expression.IMPORTANCELatency-reversing agents (LRAs) are ubiquitously used in the "shock-and-kill" HIV cure strategy and their performance is often evaluated by ex-vivo quantification of cell associated HIV RNA. HIV RNA, measured by qPCR, is often normalized to internal reference genes, but the expression of these genes should not be influenced by the experimental settings. We found that treatment of human PBMC and CD4+ T cells with LRAs significantly altered the expression of several commonly used reference genes, such as GAPDH. Finally, we evaluate the impact of different reference genes on normalization of host genes and HIV cell associated RNA expression and demonstrated that using unstable reference genes dramatically altered experimental outcome. Our data highlight the importance of using reference genes that are unaffected by LRAs under study to correctly evaluate host gene and cell associated HIV RNA expression induced by latency-reversing agents.

Identification and Validation of Reference Genes in Clostridium beijerinckii NRRL B-598 for RT-qPCR Using RNA-Seq Data

Gene expression analysis through reverse transcription-quantitative real-time polymerase chain reaction (RT-qPCR) depends on correct data normalization by reference genes with stable expression. Although Clostridium beijerinckii NRRL B-598 is a promising Gram-positive bacterium for the industrial production of biobutanol, validated reference genes have not yet been reported. In this study, we selected 160 genes with stable expression based on an RNA sequencing (RNA-Seq) data analysis, and among them, seven genes (zmp, rpoB1, rsmB, greA, rpoB2, topB2, and rimO) were selected for experimental validation by RT-qPCR and gene ontology (GO) enrichment analysis. According to statistical analyses, zmp and greA were the most stable and suitable reference genes for RT-qPCR normalization. Furthermore, our methodology can be useful for selection of the reference genes in other strains of C. beijerinckii and it also suggests that the RNA-Seq data can be used for the initial selection of novel reference genes, however, their validation is required.

Selection and validation reference genes for qRT-PCR normalization in different cultivars during fruit ripening and softening of peach (Prunus persica)

Quantitative real-time PCR (qRT-PCR) has been emerged as an effective method to explore the gene function and regulatory mechanisms. However, selecting appropriate reference gene (s) is a prerequisite for obtaining accurate qRT-PCR results. Peach is one of important fruit in Rosaceae and is widely cultivated worldwide. In this study, to explore reliable reference gene (s) in peach with different types during fruit ripening and softening (S1-S4), nine candidate reference genes (EF-1α, GAPDH, TBP, UBC, eIF-4α, TUB-A, TUB-B, ACTIN, and HIS) were selected from the whole-genome data. Then, the expression levels of the nine selected genes were detected using qRT-PCR in three peach types, including 'Hakuho' (melting type), 'Xiacui' (stony hard type), 'Fantasia' and 'NJC108' (non-melting type) cultivars were detected using qRT-PCR. Four software (geNorm, NormFinder, BestKeeper and RefFinder) were applied to evaluate the expression stability of these candidate reference genes. Gene expression was characterized in different peach types during fruit ripening and softening stages. The overall performance of each candidate in all samples was evaluated. The Actin gene (ACTIN) was a suitable reference gene and displayed excellent stability in 'Total' set, 'Hakuho' samples, S3 and S4 fruit developmental stages. Ubiquitin C gene (UBC) showed the best stability in most independent samples, including 'Fantasia', 'NJC108', S2 sets. Elongation factor-1α gene (EF-1α) was the most unstable gene across the set of all samples, 'NJC108' and S2 sets, while showed the highest stability in 'Xiacui' samples. The stability of candidate reference genes was further verified by analyzing the relative expression level of ethylene synthase gene of Prunus persica (PpACS1) in fruit ripening and softening periods of 'Hakuho'. Taken together, the results from this study provide a basis for future research on the mining of important functional genes, expression patterns and regulatory mechanisms in peach.

Validation of Reference Genes via qRT-PCR in Multiple Conditions in Brandt's Voles, Lasiopodomys brandtii

Simple Summary

This study validated the stability of the expression profiles of nine common candidate reference genes (Gapdh, Hprt1, β-actin, PPIA, Rpl13a, Tbp, Sdha, Hmbs, and B2M) using qRT-PCR in different tissues, developmental stages, and photoperiods. None of these genes were suitable as optimal reference genes at 4 weeks postnatal in different tissues. Under different developmental stages in the hypothalamus, B2M for males and Rpl13a for females were suitable as reference genes. Under different photoperiods in the hypothalamus, none of the selected genes were suitable as reference genes at 6 weeks postnatal, β-actin and PPIA were the optimal reference genes at 12 weeks postnatal, while Hprt1, β-actin, PPIA, Hmbs, and B2M were excellent reference genes at 24 weeks postnatal.

Abstract

The choice of optimal reference gene is challenging owing to the varied expression of reference genes in different organs, development stages, and experimental treatments. Brandt’s vole (Lasiopodomys brandtii) is an ideal animal to explore the regulatory mechanism of seasonal breeding, and many studies on this vole involve gene expression analysis using quantitative real-time polymerase chain reaction (qRT-PCR). In this study, we used the method of the coefficient of variation and the NormFinder algorithm to evaluate the performance of nine commonly used reference genes Gapdh, Hprt1, β-actin, PPIA, Rpl13a, Tbp, Sdha, Hmbs, and B2M using qRT-PCR in eight different tissues, five developmental stages, and three different photoperiods. We found that all nine genes were not uniformly expressed among different tissues. B2M and Rpl13a were the optimal reference genes for different postnatal development stages in the hypothalamus for males and females, respectively. Under different photoperiods in the hypothalamus, none of the selected genes were suitable as reference genes at 6 weeks postnatal; β-actin and PPIA were the optimal reference genes at 12 weeks postnatal; Hprt1, β-actin, PPIA, Hmbs, and B2M were excellent reference genes at 24 weeks postnatal. The present study provides a useful basis for selecting the appropriate reference gene in Lasiopodomys brandtii.

Selection of reliable reference genes for analysis of gene expression in the rat placenta

The selection of suitable reference genes (RGs), especially the identification of the proper combination of RGs is the key to obtain reliable results of gene expression for quantitative real-time polymerase chain reaction (qRT-PCR). To date, there is no relevant study dealing with the stability of RGs in rat placenta. In this study, the geNorm, NormFinder, and BestKeeper software were used to analyze the expression stability of the candidate RGs in placenta under physiological and prenatal caffeine exposure (PCE) conditions. The expression of Tbp, Gapdh and Ywhaz in female and Polr2a, Gapdh and Ywhaz in male placenta were highly stable under physiological conditions, and there was no obvious gender difference. We further found that two RGs were sufficient for reliable normalization in female and male placenta and the combination of Ywhaz and Gapdh was the most suitable compound RGs under physiological conditions. Under PCE conditions, Ywhaz, Gapdh and Polr2a were the most stable genes in both female and male placenta. Among them, Ywhaz and Gapdh were chosen as the best paring. Finally, selected RGs were employed for normalization of the expression of a clear target gene and the results of standardization supported our choice. In conclusion, our study confirmed that Ywhaz/Gapdh combination was the most suitable RGs in rat placenta under physiological and PCE pathological conditions and provided a theoretical and experimental basis for physiological and pathological research of the rat placenta.

Selection of reference genes for quantitative PCR analysis in poultry red mite (Dermanyssus gallinae)

Poultry red mites (PRMs, Dermanyssus gallinae) are harmful ectoparasites that affect farmed chickens and cause serious economic losses in the poultry industry worldwide. Acaricides are used for PRM control; however, some PRMs have developed acaricide-resistant properties, which have indicated the need for different approaches for PRM control. Therefore, it is necessary to elucidate the biological status of PRMs to develop alternative PRM control strategies. Quantitative polymerase chain reaction (qPCR) allows analysis of the biological status at the transcript level. However, reference genes are preferable for accurate comparison of expression level changes given the large variation in the quality of the PRM samples collected in each farm. This study aimed to identify candidate reference genes with stable expression levels in the different blood feeding states and life stages of PRMs. First, we selected candidates based on the following criteria: sufficient expression intensity and no significant expression difference between fed and starved states. We selected and characterized seven candidate reference genes. Among them, we evaluated the gene expression stability between the starved and fed states using RefFinder; moreover, we compared their expression levels in each life-stage and identified two reference genes, Elongation factor 1-alpha (ELF1A)-like and apolipophorins-like. Finally, we evaluated the utility of the candidates as reference genes, and the use of ELF1A-like and apolipophorins-like successfully normalized ATP synthase subunit g -like gene expression. Thus, ELF1A-like and apolipophorins-like could be suitable reference genes in PRMs.

Differential gene expression and Ingenuity Pathway Analysis of bronchoalveolar lavage cells from horses with mild/moderate neutrophilic or mastocytic inflammation on BAL cytology

Mild to moderate equine asthma syndrome (mEAS) affects horses of all ages and breeds. To date, the etiology and pathophysiology of mEAS are active areas of research, and it remains incompletely understood whether mEAS horses with different immune cell 'signatures' on BAL cytology represent different phenotypes, distinct pathobiological mechanisms (endotypes), varied environmental conditions, disease severity, genetic predispositions, or all of the above. In this descriptive study, we compared gene expression data from BAL cells isolated from horses with normal BALF cytology (n = 5), to those isolated from horses with mild/moderate neutrophilic inflammation (n = 5), or mild/moderate mastocytic inflammation (n = 5). BAL cell protein lysates were analyzed for cytokine/chemokine levels using Multiplex Bead Immunoassay, and for select proteins using immunoblot. The transcriptome, determined by RNA-seq and analyzed with DEseq2, contained 20, 63, and 102 significantly differentially expressed genes in horses with normal vs. neutrophilic, normal vs. mastocytic, and neutrophilic vs. mastocytic BALF cytology, respectively. Pathway analyses revealed that BAL-isolated cells from horses with neutrophilic vs. normal cytology showed enrichment in inflammation pathways, and horses with mastocytic vs. normal cytology showed enrichment in pathways involved in fibrosis and allergic reaction. BAL cells from horses with mastocytic mEAS, compared to neutrophilic mEAS, showed enrichment in pathways involved in alteration of tissue structures. Cytokine analysis determined that IL-1β was significantly different in the lysates from horses with neutrophilic inflammation compared to those with normal or mastocytic BAL cytology. Immunoblot revealed significant difference in the relative level of MMP2 in horses with neutrophilic vs. mastocytic mEAS. Upregulation of mRNA transcripts involved in the IL-1 family cytokine signaling axis (IL1a, IL1b, and IL1R2) in neutrophilic mEAS, as well as KIT mRNA in mastocytic mEAS, are novel, potentially clinically relevant, findings of this study. These findings further inform our understanding of inflammatory cell subtypes in mEAS.

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.

Selection and validation of reference genes for quantitative real-time PCR in the green microalgae Tetraselmis chui

Quantitative real-time reverse transcription PCR (RT-qPCR) is a highly sensitive technique that can be applied to analyze how genes are modulated by culture conditions, but identification of appropriate reference genes for normalization is a critical factor to be considered. For this reason, the expression stability of 18 candidate reference genes was evaluated for the green microalgae Tetraselmis chui using the widely employed algorithms geNorm, NormFinder, BestKeeper, the comparative ΔCT method, and RefFinder. Microalgae samples were collected from large scale outdoor photobioreactors during the growing phase (OUT_GP), and during the semi-continuous phase at different times of the day (OUT_DC). Samples from standard indoor cultures under highly controlled conditions (IND) were also collected to complement the other data. Different rankings for the candidate reference genes were obtained depending on the culture conditions and the algorithm employed. After comparison of the achieved ranks with the different methods, the references genes selected for samples from specific culture conditions were ALD and EFL in OUT_GP, RPL32 and UBCE in OUT_DC, and cdkA and UBCE in IND. Moreover, the genes EFL and cdkA or EFL and UBCE appeared as appropriate combinations for pools generated from all samples (ALL). Examination in the OUT_DC cultures of genes encoding the large and small subunits of ADP-glucose pyrophosphorylase (AGPL and AGPS, respectively) confirmed the reliability of the identified reference genes, RPL32 and UBCE. The present study represents a useful contribution for studies of gene expression in T. chui, and also represents the first step to set-up an RT-qPCR platform for quality control of T. chui biomass production in industrial facilities.

Identification and validation of suitable reference genes for quantitative real-time PCR gene expression analysis in pregnant human myometrium

Accurate quantification of quantitative PCR (qPCR) data requires a set of stable reference genes (RGs) for normalisation. Despite its importance to mechanistic studies, no evaluation of RG stability has been conducted for pregnant human myometrium. A systematic search of the literature was performed to identify the most used RGs in human myometrial gene expression studies. The stability of these genes, and others, was then evaluated using geNorm and NormFinder algorithms, in samples of myometrium from singleton or twin pregnancies (n = 7 per group) delivering at term or preterm. The most frequently cited RGs were GAPDH, ACTB, B2M and 18s. There was strong agreement between algorithms on the most and least stable genes: Both indicated CYC1, YWHAZ and ATP5B were the most stably expressed. Despite being some of the most used RGs, B2M, 18s and ACTB expression was least stable and was too variable for use as accurate normalisation factors. Pairwise variation analysis determined that the optimal number of RGs for accurate normalisation is two. Validation of the choice of RGs by comparing relative expression of oxytocin receptors (OXTR) using the least stable 18s and B2M, with the most stable, CYC1 and YWHAZ, erroneously demonstrated significantly increased OXTR expression in myometrium in singleton pregnancies compared to twins. This study demonstrates the importance of appropriate RG selection for accurate quantification of relative expression in pregnant human myometrium qPCR studies. For normalisation, the geometric mean of CYC1 and YWHAZ or ATP5B is suggested. The use of ACTB, 18s and B2M, is not recommended.