A decade of improvements in quantification of gene expression and internal standard selection

Reference genes for gene expression profiling in mouse models of Listeria monocytogenes infection

RT-qPCR dissects transcription-based processes but requires reference genes (RGs) for data normalization. This study prospected RGs for mouse macrophages (pMØ) and spleen infected with Listeria monocytogenes. The pMØ were infected in vitro with L. monocytogenes or vehicle for 4 h. Mice were injected with L. monocytogenes (or vehicle) and euthanized 24 h post-injection. The RGs came from a multispecies primer set, from the literature or designed here. The RG ranking relied on GeNorm, NormFinder, BestKeeper, Delta-CT and RefFinder. B2m-H3f3a-Ppia were the most stable RGs for pMØ, albeit RG indexes fine-tuned estimations of cytokine relative expression. Actβ-Ubc-Ppia were the best RGs for spleen but modestly impacted the cytokine relative expression. Hence, mouse models of L. monocytogenes require context-specific RGs for RT-qPCR, thus reinforcing its paramount contribution to accurate gene expression profiling.

Transcriptome-Based Identification of the Optimal Reference Genes for Quantitative Real-Time Polymerase Chain Reaction Analyses of Lingonberry Fruits throughout the Growth Cycle

(1) Background: Vaccinium vitis-idaea is a nutritionally and economically valuable natural wild plant species that produces berries useful for treating various diseases. There is growing interest in lingonberry, but there is limited information regarding lingonberry reference genes suitable for gene expression analyses of different tissues under various abiotic stress conditions. The objective of this study was to identify stable reference genes suitable for different lingonberry tissues in response to abiotic stress. (2) Methods: The delta Ct method and the GeNorm v3.5 and NormFinder v20 programs were used to comprehensively analyze gene expression stability. (3) Results: Actin Unigene23839 was the best reference gene for analyzing different cultivars, whereas Actin CL5740.Contig2 was the most suitable reference gene for analyzing different tissues and alkali stress. In contrast, 18S rRNA CL5051.Contig1 was the most stable reference gene under drought conditions. (4) Conclusions: These suitable reference genes may be used in future qRT-PCR analyses of different lingonberry tissues and the effects of abiotic stresses. Furthermore, the study data may be useful for functional genomics studies and the molecular breeding of lingonberry. In summary, internal reference genes or internal reference gene combinations should be carefully selected according to the experimental conditions to ensure that the generated gene expression data are accurate.

Selection of reference genes for RT-qPCR analysis in developing chicken embryonic ovary

BACKGROUND

Normalization of the expression profiling of target genes, in a tissue-specific manner and under different experimental conditions, requires stably expressed gene(s) to be used as internal reference(s). However, to study the molecular regulation of oocyte meiosis initiation during ovary development in chicken embryos, stable reference gene(s) still need to be compared and confirmed.

METHODS AND RESULTS

Six candidate genes previously used as internal references for the chicken embryo (Actb, Cvh, Dazl, Eef1a, Gapdh and Rpl15) were chosen, and their expression profiles in left ovaries dissected at five chicken embryonic days (E12.5, E15.5, E17.5, E18.5 and E20.5) were evaluated, respectively. Separately, GeNorm, NormFinder, BestKeeper and Comparative ΔCt methods were used to assess the stability of candidate reference genes, and all results were combined to give the final rank by RefFinder. All methods identified that Eef1a and Rpl15 were the two most stable internal reference genes, whereas Cvh is the most unstable one. Moreover, expression levels of three marker genes for chicken oocyte meiosis entry (Stra8, Scp3 and Dmc1) were normalized, based on Eef1a, Rpl15, or their combinations, respectively.

CONCLUSION

Our findings provide the most suitable internal reference genes (Eef1a and Rpl15), to investigate further molecular regulation of ovary development and oocyte meiosis initiation in chicken embryos.

Resistance Exercise-Induced Increases in Muscle Myostatin mRNA and Protein Expression Are Subsequently Decreased in Circulation in the Presence of Increased Levels of the Extracellular Matrix Stabilizing Protein Decorin

Resistance exercise (RE) activates cell signaling pathways associated with myostatin. Decorin is located in the extracellular matrix (ECM) and can block the inhibitory effect of myostatin. This study sought to determine the impact of low-load (LL) and high-load (HL) RE on myostatin mRNA and protein expression along with changes in muscle decorin and circulating follistatin. Ten resistance-trained men performed a LL (50% 1RM) and HL (80% 1RM) RE session using the angled leg press and leg extension with load and volume equated. Venous blood samples and muscle biopsies were obtained prior to and at 3h and 24h following each RE session. Muscle myostatin mRNA expression was increased at 24h post-exercise (p = 0.032) in LL and at 3h (p = 0.044) and 24h (p = 0.003) post-exercise in HL. Muscle decorin was increased at 24h post-exercise (p < 0.001) in LL and HL; however, muscle myostatin was increased at 24h post-exercise (p < 0.001) only in HL. For muscle Smad 2/3, no significant differences were observed (p > 0.05). Serum follistatin was increased and myostatin decreased at 24h post-exercise (p < 0.001) in LL and HL. Muscle myostatin gene and protein expression increased in response to HL RE. However, serum myostatin was decreased in the presence of increases in decorin in muscle and follistatin in circulation. Therefore, our data suggest a possible mechanism may exist where decorin within the ECM is able to bind to, and decrease, myostatin that might otherwise enter the circulation for activin IIB (ACTIIB) receptor binding and subsequent canonical signaling through Smad 2/3.

Reference Gene Validation for RT-qPCR in PBMCs from Asthmatic Patients with or without Obesity

Obesity is known to impair the efficacy of glucocorticoid medications for asthma control. Glucocorticoid-induced gene expression studies may be useful to discriminate those obese asthmatic patients who present a poor response to glucocorticoids. The expression of genes of interest is normalized with respect to reference genes (RGs). Ideally, RGs have a stable expression in different samples and are not affected by experimental conditions. The objective of this work was to analyze suitable RGs to study the role of glucocorticoid-induced genes in obese asthmatic patients in further research. The gene expression of eight potential RGs (GUSB, B2M, POLR2A, PPIA, ACTB, GAPDH, HPRT1, and TBP) was assessed with reverse transcription-quantitative polymerase chain reaction in peripheral blood mononuclear cells (PBMCs) from asthmatic, obese asthmatic, and healthy individuals. Their stability was analyzed using four different algorithms-BestKeeper, ΔCt, geNorm, and NormFinder. geNorm analysis recommended the use of a minimum of three genes for normalization. Moreover, intergroup variation due to the treatment was calculated by NormFinder, which found that B2M was the gene that was least affected by different treatments. Comprehensive rankings indicated GUSB and HPRT1 as the best RGs for qPCR in PBMCs from healthy and asthmatic subjects, while B2M and PPIA were the best for obese asthmatic subjects. Finally, our results demonstrated that B2M and HPRT1 were the most stable RGs among all groups, whereas ACTB, TBP, and GAPDH were the worst shared ones.

Identification of Suitable Reference Genes for qRT-PCR Normalization in Kiwifruit

Reference genes are used for the correction of qRT-PCR data, and it is necessary to investigate the optimum reference gene under certain conditions. The expression levels of seven traditional reference genes ACT1, ACT2, GAPDH, 18S rRNA, UBQ, TUB and CYP were analyzed using qRT-PCR in different varieties, tissues, developmental stages and hormone (or pollen polysaccharide) treatments in kiwifruit. Gene expression stability was assessed with the help of three common software (geNorm, NormFinder, BestKeeper), and the minimum number of reference genes necessary for normalization was also determined. GAPDH, ACT1 and ACT2 were selected as reference genes for different genotypes of kiwifruit. GAPDH and UBQ were the best combinations of reference genes for root, stem, leaf, flower and fruit. GAPDH and ACT1 could be the preferred reference genes for normalization of qRT-PCR data during fruit development. The pairing of ACT1 and UBQ constituted the optimal combination of reference genes in kiwifruit treated with different hormones (or pollen polysaccharide). This study provides a new and reliable option for the use of reference genes in the analysis of gene expression patterns of interest in kiwifruit.

Transcriptome-Based Selection and Validation of Reference Genes for Gene Expression Analysis of Alicyclobacillus acidoterrestris Under Acid Stress

Alicyclobacillus acidoterrestris is a major concern in fruit juice industry due to its spoilage potential of acidic fruit juice. Quantifying the expression levels of functional genes by real-time quantitative polymerase chain reaction (RT-qPCR) is necessary to elucidate the response mechanisms of A. acidoterrestris to acid stress. However, appropriate reference genes (RGs) for data normalization are required to obtain reliable RT-qPCR results. In this study, eight novel candidate RGs were screened based on transcriptome datasets of A. acidoterrestris under acid stress. The expression stability of eight new RGs and commonly used RG 16s rRNA was assessed using geNorm, NormFinder, and BestKeeper algorithms. Moreover, the comprehensive analysis using the RefFinder program and the validation using target gene ctsR showed that dnaG and dnaN were the optimal multiple RGs for normalization at pH 4.0; ytvI, dnaG, and 16s rRNA at pH 3.5; icd and dnaG at pH 3.0; and ytvI, dnaG, and spoVE at pH 2.5. This study revealed for the first time that A. acidoterrestris had different suitable RGs under different acid conditions, with implications for further deciphering the acid response mechanisms of this spoilage-causing bacterium.

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.

Efficiency Correction Is Required for Accurate Quantitative PCR Analysis and Reporting

BACKGROUND

Quantitative PCR (qPCR) aims to measure the DNA or RNA concentration in diagnostic and biological samples based on the quantification cycle (Cq) value observed in the amplification curves. Results of qPCR experiments are regularly calculated as if all assays are 100% efficient or reported as just Cq, ΔCq, or ΔΔCq values.

CONTENTS

When the reaction shows specific amplification, it should be deemed to be positive, regardless of the observed Cq. Because the Cq is highly dependent on amplification efficiency that can vary among targets and samples, accurate calculation of the target quantity and relative gene expression requires that the actual amplification efficiency be taken into account in the analysis and reports. PCR efficiency is frequently derived from standard curves, but this approach is affected by dilution errors and hampered by properties of the standard and the diluent. These factors affect accurate quantification of clinical and biological samples used in diagnostic applications and collected in challenging conditions. PCR efficiencies determined from individual amplification curves avoid these confounders. To obtain unbiased efficiency-corrected results, we recommend absolute quantification with a single undiluted calibrator with a known target concentration and efficiency values derived from the amplification curves of the calibrator and the unknown samples.

SUMMARY

For meaningful diagnostics or biological interpretation, the reported results of qPCR experiments should be efficiency corrected. To avoid ambiguity, the Minimal Information for Publications on Quantitative Real-Time PCR Experiments (MIQE) guidelines checklist should be extended to require the methods that were used (1) to determine the PCR efficiency and (2) to calculate the reported target quantity and relative gene expression value.

Identification and validation of new reference genes for accurate quantitative reverse transcriptase-PCR normalization in the Antarctic plant Colobanthus quitensis under abiotic stress conditions

The Antarctic ecotype of Colobanthus quitensis is a vascular plant highly adapted to the harsh environmental conditions of Maritime Antarctica which is now facing with the rapid local warming experienced in the Antarctic Peninsula during the last decades. Thus, the identification of the molecular mechanisms leading to the adaptation to this warming trend is a new target for modern cell physiology. The selection of suitable reference genes for quantification of key stress-responsive genes through quantitative Reverse Transcriptase-Polymerase Chain Reaction (qRT-PCR) is important to ensure accurate and reliable results. In this study, we evaluated the expression stability of eleven candidate genes in C. quitensis under different abiotic stress conditions using geNorm and RefFinder tools. The statistical analysis showed that the appropriate reference genes varied depending on the experimental conditions, even if EF1α and PP2Acs ranked as the most stable reference genes when all stress conditions were considered. To further validate the stability of the selected reference genes, the expression patterns of C. quitensis catalase gene (CqCAT) was analyzed. The reference genes validated in this study will be useful for improving the accuracy of qRT-PCR analysis for gene expression studies of the Antarctic ecotype of C. quitensis and could be extended to other ecotypes adapted to low temperatures.

Guidelines for Setting Up a mRNA Sequencing Experiment and Best Practices for Bioinformatic Data Analysis

RNA-sequencing, commonly referred to as RNA-seq, is the most recently developed method for the analysis of transcriptomes. It uses high-throughput next-generation sequencing technologies and has revolutionized our understanding of the complexity and dynamics of whole transcriptomes.In this chapter, we recall the key developments in transcriptome analysis and dissect the different steps of the general workflow that can be run by users to design and perform a mRNA-seq experiment as well as to process mRNA-seq data obtained by the Illumina technology. The chapter proposes guidelines for completing a mRNA-seq study properly and makes available recommendations for best practices based on recent literature and on the latest developments in technology and algorithms. We also remark the large number of choices available (especially for bioinformatic data analysis) in front of which the scientist may be in trouble.In the last part of the chapter we discuss the new frontiers of single-cell RNA-seq and isoform sequencing by long read technology.

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.

There is no magic bullet: the importance of testing reference gene stability in RT-qPCR experiments across multiple closely related species

Reverse Transcriptase quantitative Polymerase Chain Reaction (RT-qPCR) is the current gold standard tool for the study of gene expression. This technique is highly dependent on the validation of reference genes, which exhibit stable expression levels among experimental conditions. Often, reference genes are assumed to be stable a priori without a rigorous test of gene stability. However, such an oversight can easily lead to misinterpreting expression levels of target genes if the references genes are in fact not stable across experimental conditions. Even though most gene expression studies focus on just one species, comparative studies of gene expression among closely related species can be very informative from an evolutionary perspective. In our study, we have attempted to find stable reference genes for four closely related species of grasshoppers (Orthoptera: Acrididae) that together exhibit a spectrum of density-dependent phenotypic plasticity. Gene stability was assessed for eight reference genes in two tissues, two experimental conditions and all four species. We observed clear differences in the stability ranking of these reference genes, both between tissues and between species. Additionally, the choice of reference genes clearly influenced the results of a gene expression experiment. We offer suggestions for the use of reference genes in further studies using these four species, which should be taken as a cautionary tale for future studies involving RT-qPCR in a comparative framework.

Screening and Validation of Reference Genes for RT-qPCR Under Different Honey Bee Viral Infections and dsRNA Treatment

Honey bee viruses are one of the most important pathogens that have contributed to the decrease in honey bee colony health. To analyze the infection dynamics of honey bee viruses, quantification of viral gene expression by RT-qPCR is necessary. However, suitable reference genes have not been reported from viral and RNAi studies of honey bee. Here, we evaluated the expression of 11 common reference genes (ache2, rps18, β-actin, tbp, tif, rpl32, gadph, ubc, α-tubulin, rpl14, and rpsa) from Apis mellifera (Am) and Apis cerana (Ac) under Israeli acute paralysis virus (IAPV), chronic bee paralysis virus (CBPV), and Chinese sacbrood virus (CSBV) infection as well as dsRNA-PGRP-SA treatment, and we confirmed their validation by evaluating the levels of the defensin 1 and prophenoloxidase (ppo) genes during viral infection. Our results showed that the expression of selected genes varied under different viral infections. ache2, rps18, β-actin, tbp, and tif can be used to normalize expression levels in Apis mellifera under IAPV infection, while the combination of actin and tif is suitable for CBPV-infected experiments. The combination of rpl14, tif, rpsa, ubc, and ache2 as well as more reference genes is suitable for CSBV treatment in Apis cerana. Rpl14, tif, rps18, ubc, and α-tubulin were the most stable reference genes under dsRNA treatment in Apis mellifera. Furthermore, the geNorm and NormFinder algorithms showed that tif was the best suitable reference gene for these four treatments. This study screened and validated suitable reference genes for the quantification of viral levels in honey bee, as well as for RNAi experiments.

Reliable reference genes for the quantification of mRNA in human T-cells and PBMCs stimulated with live influenza virus

Background

Quantitative PCR (qPCR) is a powerful tool that is particularly well-suited to measure mRNA levels in clinical samples, especially those with relatively low cell counts. However, a caveat of this approach is that reliable, stably expressed reference (housekeeping) genes are vital in order to ensure reproducibility and appropriate biological inference. In this study, we evaluated the expression stability of six reference genes in peripheral blood mononuclear cells (PBMCs) and isolated CD3⁺ T-cells from young and old adults (n = 10), following ex vivo stimulation with mock (unstimulated) or live influenza virus. Our genes included: β-actin (ACTB), glyercaldehyde-3-phostphate dehydrogenase (GAPDH), ribosomal protein L13a (RPL13a), ribosomal protein S18 (RPS18), succinate dehydrogenase complex flavoprotein subunit A (SDHA), and ubiquitin-conjugating enzyme E2D2 (UBE2D2).

Results

Reference gene expression varied significantly depending on cell type and stimulation conditions, but not age. Using the comparative ΔCt method, and the previously published software BestKeeper, NormFinder, and geNorm, we show that in PBMCs and T-cells, UBE2D2 and RPS18 were the most stable reference genes, followed by ACTB; however, the expression of UBE2D2 and RPS18 was found to increase with viral stimulation in isolated T-cells, while ACTB expression did not change significantly. No age-related differences in stability were observed for any gene

Conclusions

This study suggests the use of a combination of UBE2D2, RPS18, and ACTB for the study of influenza responses in PBMCs and T-cells, although ACTB alone may be the most optimal choice if choosing to compare target gene expression before and after viral stimulation. Both GAPDH and RPL13a were found to be poor reference genes and should be avoided for studies of this nature.

Mineralocorticoid receptor excessive activation involved in glucocorticoid-related brain injury

The mechanisms involved in brain damage during chronic glucocorticoid exposure are poorly understood. Since mineralocorticoid receptor (MR) activation has been proven to be important in the pathophysiology of vascular damage and MRs are highly expressed in many brain regions, we hypothesized that the cerebral injury observed in subjects with Cushing syndrome is in part associated with the overactivation of MR. The aim of this study was to determine whether the cerebral injury observed in chronic hyperglucocorticoidemia animal models is related to excessive MR activation. Male SD rats were divided into five groups: vehicle, hydrocortisone (HC, 5 mg/kg/day, i.g.), HC + spironolactone (SL, 20 mg/kg/d in chow), dexamethasone (DXM, 0.25 mg/kg/day, i.g.), and DXM + SL (20 mg/kg/d in chow). Compared to the vehicle-treated group, HC-treated rats had higher blood pressure and higher levels of cerebral vascular fibrosis, cortical/hippocampal atrophy, reactive oxygen species (ROS) production and proinflammatory gene expression. However, in HC-treated animals, treatment with SL markedly alleviated ROS production, cerebral and cerebrovascular morphological changes and inflammation but failed to reduce blood pressure. In contrast, DXM induced no cerebral morphological changes except fibrosis in cerebral vessels, an effect that was not ameliorated by SL treatment. These findings demonstrate that the excessive MR activation observed following chronic hyperglucocorticoidemia exposure contributes to cerebrovascular fibrosis and remodeling and promotes neural apoptosis in the cerebral cortex/hippocampus.

Expression analysis of ethylene synthesis and signalling genes in kiwifruit stigmatic arms and their involvement in programmed cell death processes

Kiwifruit (Actinidia chinensis var. deliciosa (A. Chev) A. Chev.) is a widely cultivated crop due to the nutritional value of its fruits. Its commercialization is related to the fruit size, which is directly linked with the number of seeds and, consequently, with pollination. In this dioecious species pollination is dependent on a short effective pollination period which is related to a Programmed Cell Death (PCD) process. At the same time, this PCD process allows the growth of many pollen tubes. Several studies suggest that ethylene can play an important role in PCD in a number of systems. In this report, we determined the full sequence of the AcACS gene, encoding the enzyme that catalyses a rate-limiting step of the ethylene synthesis. Next, we monitored the expression pattern of this gene as well as of other genes involved in ethylene synthesis (ACO2-5) and signalling (AdERS1a, AdERS1b, AdETR1, AdETR2, AdETR3, AdCTR1, AdCTR2, AdEIL1) in pollinated and non-pollinated stigmatic arms of kiwifruit female flowers. The relative expression patterns observed for AcACS, ACOs and ethylene perception and signalling genes (AdERS1, AdETR1, AdCTR1 and AdEIL1) showed that they are expressed before anthesis. After anthesis, expression of the studied genes was detected earlier in pollinated than in non-pollinated stigmatic arms, as it was previously determined for PCD hallmarks. In addition, the expression pattern of the studied genes showed a clear relationship with the PCD hallmarks described in a previous report in the secretory tissue both in non-pollinated stigmatic arms (related to the short EPP in this species) and in pollinated ones (related to the growth of many pollen tubes during progamic phase). Overall, these results suggest an involvement of ethylene with PCD contributing to the high reproductive success of this species.

The role and robustness of the Gini coefficient as an unbiased tool for the selection of Gini genes for normalising expression profiling data

We recently introduced the Gini coefficient (GC) for assessing the expression variation of a particular gene in a dataset, as a means of selecting improved reference genes over the cohort ('housekeeping genes') typically used for normalisation in expression profiling studies. Those genes (transcripts) that we determined to be useable as reference genes differed greatly from previous suggestions based on hypothesis-driven approaches. A limitation of this initial study is that a single (albeit large) dataset was employed for both tissues and cell lines. We here extend this analysis to encompass seven other large datasets. Although their absolute values differ a little, the Gini values and median expression levels of the various genes are well correlated with each other between the various cell line datasets, implying that our original choice of the more ubiquitously expressed low-Gini-coefficient genes was indeed sound. In tissues, the Gini values and median expression levels of genes showed a greater variation, with the GC of genes changing with the number and types of tissues in the data sets. In all data sets, regardless of whether this was derived from tissues or cell lines, we also show that the GC is a robust measure of gene expression stability. Using the GC as a measure of expression stability we illustrate its utility to find tissue- and cell line-optimised housekeeping genes without any prior bias, that again include only a small number of previously reported housekeeping genes. We also independently confirmed this experimentally using RT-qPCR with 40 candidate GC genes in a panel of 10 cell lines. These were termed the Gini Genes. In many cases, the variation in the expression levels of classical reference genes is really quite huge (e.g. 44 fold for GAPDH in one data set), suggesting that the cure (of using them as normalising genes) may in some cases be worse than the disease (of not doing so). We recommend the present data-driven approach for the selection of reference genes by using the easy-to-calculate and robust GC.

Selection of reference genes for normalization of cranberry (Vaccinium macrocarpon Ait.) gene expression under different experimental conditions

Real-time fluorescent quantitative PCR (qRT-PCR) is often chosen as an effective experimental method for analyzing gene expression. However, an appropriate reference gene as a standard is needed to obtain accurate gene expression data. To date, no internal reference genes have been reported for research on cranberries. Expanding the selection of internal reference genes for cranberry will enable reliable gene expression analysis, and, at the same time, can also lay a solid foundation for revealing the biological characteristics of cranberry. Here, we selected ten candidate reference gene families and used three statistical software tools-geNorm, NormFinder and BestKeeper-to evaluate their expression stability under the influence of different experimental factors. The results showed that protein phosphatase 2A regulatory subunit (PP2A) or RNA helicase-like 8 (RH 8) was the best choice for an internal reference gene when analyzing different cranberry cultivars. In two sample sets comprising different cranberry organs and three abiotic stress treatments, sand family protein (SAND) was the best choice as a reference gene. In this study, we screened genes that are stably expressed under the influence of various experimental factors by qRT-PCR. Our results will guide future studies involving gene expression analysis of cranberry.

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

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

Exploring regulatory networks in plants: transcription factors of starch metabolism

Biological networks are complex (non-linear), redundant (cyclic) and compartmentalized at the subcellular level. Rational manipulation of plant metabolism may have failed due to inherent difficulties of a comprehensive understanding of regulatory loops. We first need to identify key factors controlling the regulatory loops of primary metabolism. The paradigms of plant networks are revised in order to highlight the differences between metabolic and transcriptional networks. Comparison between animal and plant transcription factors (TFs) reveal some important differences. Plant transcriptional networks function at a lower hierarchy compared to animal regulatory networks. Plant genomes contain more TFs than animal genomes, but plant proteins are smaller and have less domains as animal proteins which are often multifunctional. We briefly summarize mutant analysis and co-expression results pinpointing some TFs regulating starch enzymes in plants. Detailed information is provided about biochemical reactions, TFs and cis regulatory motifs involved in sucrose-starch metabolism, in both source and sink tissues. Examples about coordinated responses to hormones and environmental cues in different tissues and species are listed. Further advancements require combined data from single-cell transcriptomic and metabolomic approaches. Cell fractionation and subcellular inspection may provide valuable insights. We propose that shuffling of promoter elements might be a promising strategy to improve in the near future starch content, crop yield or food quality.

Stable Reference Gene Selection for RT-qPCR Analysis in Synechococcus elongatus PCC 7942 under Abiotic Stresses

Synechococcus elongatus PCC 7942 (S. elongatus PCC 7942) is a model cyanobacteria species for circadian clock mechanism studies. It has also been widely used as a bioreactor to produce biofuels and other metabolic products. Quantitative real-time PCR (qPCR) technology is the most commonly used method for studying the expression of specific genes, in which the relative expression level of target genes is calibrated by stably expressed internal reference genes. In this work, we examined the expression of nine candidate reference genes in time-course samples of S. elongatus PCC 7942 under no treatment (control), NaCl-stress conditions, H₂O₂-stress conditions, and high light-stress conditions. Based on the qPCR amplification parameters, the stability ranking of these candidate reference genes was established by three statistical software programs, geNorm, NormFinder, and BestKeeper. Considering all the stress conditions or high light stress alone, the results showed that the combination of prs and secA was the best choice for the double reference gene calibration method by qPCR. The combination of secA and ppc, rimM and rnpA, rnpA, and ilvD was most stable under no treatment, NaCl-stress conditions, and H₂O₂-stress conditions, respectively. rimM was stable under only special conditions and should be carefully chosen. 16S and rnpB were not suitable as internal reference genes for S. elongatus PCC 7942 qPCR experiments under all experimental conditions. To validate the above results, a cyanobacterial core clock gene, kaiC, was used to evaluate the actual performance of the optimized reference genes by qPCR, as well as the worst reference genes under different stress conditions. The results indicated that the best reference gene yielded more accurate calibration results for qPCR experiments carried out in S. elongatus PCC 7942 time-course samples.

Underlying mechanisms and ecological context of variation in exploratory behavior of the Argentine ant, Linepithema humile

Uncovering how and why animals explore their environment is fundamental for understanding population dynamics, the spread of invasive species, species interactions, etc. In social animals, individuals within a group can vary in their exploratory behavior, and the behavioral composition of the group can determine its collective success. Workers of the invasive Argentine ant (Linepithema humile) exhibit individual variation in exploratory behavior, which affects the colony's collective nest selection behavior. Here, we examine the mechanisms underlying this behavioral variation in exploratory behavior and determine its implications for the ecology of this species. We first establish that individual variation in exploratory behavior is repeatable and consistent across situations. We then show a relationship between exploratory behavior and the expression of genes that have been previously linked with other behaviors in social insects. Specifically, we found a negative relationship between exploratory behavior and the expression of the foraging (Lhfor) gene. Finally, we determine how colonies allocate exploratory individuals in natural conditions. We found that ants from inside the nest are the least exploratory individuals, whereas workers on newly formed foraging trails are the most exploratory individuals. Furthermore, we found temporal differences throughout the year: in early-mid spring, when new resources emerge, workers are more exploratory than at the end of winter, potentially allowing the colony to find and exploit new resources. These findings reveal the importance of individual variation in behavior for the ecology of social animals.

Selection of reference genes for quantitative real-time PCR analysis in chicken ovary following silver nanoparticle treatment

Reliable results of quantitative real time PCR (qPCR) analysis require normalization of target gene expression level using reference genes (RGs). However housekeeping genes expression may vary under experimental conditions, so selection of the proper RGs is a crucial step in a qPCR analysis. Several algorithms have been developed to address this problem: geNorm, NormFinder and BestKeeper. In this study, we have used these three tools to evaluate the stability of RGs in the ovarian tissues of hens treated with silver nanoparticles. Eight genes were selected for the validation: HPRT, HMBS, VIM, SDHA, TBP, RPL13, GAPDH and 18S rRNA. According to geNorm the best combination of reference genes is SDHA and TPP. NormFinder also selected SDHA as the most suitable gene, but in combination with RPL13. Analysis in BestKeeper showed that SDHA, RPL13 might be the best choice in gene expression studies using the chicken ovary. In conclusion, the results obtained depend on the algorithm used and it arises from the diverse calculation strategies used in these programs. The outcome from the NormFinder is considered to be the most trustworthy and used in further qPCR analysis.

Selection of reference genes for qRT-PCR and expression analysis of high-altitude-related genes in grassland caterpillars (Lepidoptera: Erebidae: Gynaephora) along an altitude gradient

Changes in gene expression patterns can reflect the adaptation of organisms to divergent environments. Quantitative real‐time PCR (qRT‐PCR) is an important tool for ecological adaptation studies at the gene expression level. The quality of the results of qRT‐PCR analysis largely depends on the availability of reliable reference genes (RGs). To date, reliable RGs have not been determined for adaptive evolution studies in insects using a standard approach. Here, we evaluated the reliability of 17 candidate RGs for five Gynaephora populations inhabiting various altitudes of the Tibetan Plateau (TP) using four independent (geNorm, NormFinder, BestKeeper, and the deltaCt method) and one comprehensive (RefFinder) algorithms. Our results showed that EF1‐α, RPS15, and RPS13 were the top three most suitable RGs, and a combination of these three RGs was the most optimal for normalization. Conversely, RPS2,ACT, and RPL27 were the most unstable RGs. The expression profiles of two target genes (HSP70 and HSP90) were used to confirm the reliability of the chosen RGs. Additionally, the expression patterns of four other genes (GPI,HIF1A,HSP20, and USP) associated with adaptation to extreme environments were assessed to explore the adaptive mechanisms of TP Gynaephora species to divergent environments. Each of these six target genes showed discrepant expression patterns among the five populations, suggesting that the observed expression differences may be associated with the local adaptation of Gynaephora to divergent altitudinal environments. This study is a useful resource for studying the adaptive evolution of TP Gynaephora to divergent environments using qRT‐PCR, and it also acts as a guide for selecting suitable RGs for ecological and evolutionary studies in insects.

Evaluation of Candidate Reference Genes for Quantitative Gene Expression Analysis in Spodoptera exigu a after Long-time Exposure to Cadmium

Studies on the transcriptional control of gene expression play an important role in many areas of biology. Reference genes, which are often referred to as housekeeping genes, such as GAPDH, G3PDH, EF2, RpL7A, RpL10, TUBα and Actin, have traditionally been assumed to be stably expressed in all conditions, and they are frequently used to normalize mRNA levels between different samples in qPCR analysis. However, it is known that the expression of these genes is influenced by numerous factors, such as experimental conditions. The difference in gene expression underlies a range of biological processes, including development, reproduction and behavior. The aim of this study was to show the problems associated with using reference genes in the qPCR technique, in a study on inbred strains of Spodoptera exigua selected toward cadmium resistance. We present and discuss our results and observations, and give some recommendations concerning the use and limitations of housekeeping genes as internal standards, especially in research on insects. Our results suggest that holometabolism and poikilothermia, as well as time since metamorphosis and the level of exposure to the selective factor (cadmium in this case), have a significant effect on the expression of reference genes.

microRNAs as reference genes for quantitative PCR in cotton

Cotton (Gossypium hirsutum) is the most important non-food plant in the world. Studies concerning the fiber quality and plant fitness of cotton at molecular level depend on high sensitive and reproducible gene-expression assays. However, only a few reports have described genes suitable for normalizing gene expression data. In this study, we report for the first time that microRNAs (miRNAs) are reliable reference genes (RGs) for cotton gene expression normalization in quantitative real-time reverse transcription (RT)-PCR. The stability of cotton miRNAs was assayed in root, stem, leaf and flower samples from three different cultivars [FiberMax (FM966), Delta Opal (DO) and Cedro] and under conditions of biotic stress caused by infection with Cotton leafroll dwarf virus (CLRDV). The stability of mRNAs already described as reference genes in cotton was also assessed. The geNorm, NormFinder, BestKeeper and ΔCt algorithms were used to select the best reference genes. In 8 of the 12 sets tested, miRNAs (miR172, 168 and 390) were found to be the best RGs. To validate the best selected RGs, miR159, miR164, miR2118, miR2910, miR3476, GhDCL2 and GhDCL4 expression levels were evaluated under biotic stress conditions, and miR164 and a putative myo-inositol oxigenase gene (GhMIOX) were assessed in leaves and flowers. The RGs selected in this work proved to be excellent reference genes in the two cases studied. Our results support the use of miRNAs as reference genes for miRNA and protein-coding genes.

Contemporary nucleic acid-based molecular techniques for detection, identification, and characterization of Bifidobacterium

Bifidobacteria are one of the most important bacterial groups found in the gastrointestinal tract of humans. Medical and food industry researchers have focused on bifidobacteria because of their health-promoting properties. Researchers have historically relied on classic phenotypic approaches (culture and biochemical tests) for detection and identification of bifidobacteria. Those approaches still have values for the identification and detection of some bifidobacterial species, but they are often labor-intensive and time-consuming and can be problematic in differentiating closely related species. Rapid, accurate, and reliable methods for detection, identification, and characterization of bifidobacteria in a mixed bacterial population have become a major challenge. The advent of nucleic acid-based molecular techniques has significantly advanced isolation and detection of bifidobacteria. Diverse nucleic acid-based molecular techniques have been employed, including hybridization, target amplification, and fingerprinting. Certain techniques enable the detection, characterization, and identification at genus-, species-, and strains-levels, whereas others allow typing of species or strains of bifidobacteria. In this review, an overview of methodological principle, technique complexity, and application of various nucleic acid-based molecular techniques for detection, identification, and characterization of bifidobacteria is presented. Advantages and limitations of each technique are discussed, and significant findings based on particular techniques are also highlighted.

Excision Repair Cross-complementation Group 1 is a Prognostic Biomarker in Patients with Colorectal Cancer Receiving Chemotherapy

Background:

Conflicting results about the association between expression level of excision repair cross-complementation group 1 (ERCC1) and clinical outcome in patients with colorectal cancer (CRC) receiving chemotherapy have been reported. Thus, we searched the available articles and performed the meta-analysis to elucidate the prognostic role of ERCC1 expression in patients with CRC.

Methods:

A thorough literature search using PubMed (Medline), Embase, Cochrane Library, Web of Science databases, and Chinese Science Citation Database was conducted to obtain the relevant studies. Pooled hazard ratios (HRs) or odds ratios (ORs) with 95% confidence intervals (CIs) were calculated to estimate the results.

Results:

A total of 11 studies were finally enrolled in this meta-analysis. Compared with patients with lower ERCC1 expression, patients with higher ERCC1 expression tended to have unfavorable overall survival (OS) (HR = 2.325, 95% CI: 1.720–3.143, P < 0.001), progression-free survival (PFS) (HR = 1.917, 95% CI: 1.366–2.691, P < 0.001) and poor response to chemotherapy (OR = 0.491, 95% CI: 0.243–0.990, P = 0.047). Subgroup analyses by treatment setting, ethnicity, HR extraction, detection methods, survival analysis, and study design demonstrated that our results were robust.

Conclusions:

ERCC1 expression may be taken as an effective prognostic factor predicting the response to chemotherapy, OS, and PFS. Further studies with better study design and longer follow-up are warranted in order to gain a deeper understanding of ERCC1's prognostic value.

Selection and Validation of Appropriate Reference Genes for Quantitative Real-Time PCR Normalization in Staminate and Perfect Flowers of Andromonoecious Taihangia rupestris

Quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR) is the most commonly used and powerful method for gene expression analysis due to its high sensitivity, specificity, and high throughput, and the accuracy of this approach depends on the stability of reference genes used for normalization. Taihangia rupestris Yu and Li (Rosaceae), an andromonoecious plant, produces both bisexual flowers and unisexual male flowers within the same individual. Using qRT-PCR technique, investigation of the gene expression profiling in staminate and perfect flowers would improve our understanding of the molecular mechanism in regulation of flower formation and sex differentiation in andromonoecious T. rupestris. To accurate normalize the gene expression level in Taihangia flower, 16 candidate reference genes, including 10 traditional housekeeping genes, and 6 newly stable genes, were selected based on transcriptome sequence data and previous studies. The expressions of these genes were assessed by qRT-PCR analysis in 51 samples, including 30 staminate and perfect flower samples across developmental stages and 21 different floral tissue samples from mature flowers. By using geNorm, NormFinder, BestKeeper, and comprehensive RefFinder algorithms, ADF3 combined with UFD1 were identified as the optimal reference genes for staminate flowers, while the combination of HIS3/ADF3 was the most accurate reference genes for perfect floral samples. For floral tissues, HIS3, UFD1, and TMP50 were the most suitable reference genes. Furthermore, two target genes, TruPI, and TruFBP24, involved in floral organ identity were selected to validate the most and least stable reference genes in staminate flowers, perfect flowers, and different floral tissues, indicating that the use of inappropriate reference genes for normalization will lead to the adverse results. The reference genes identified in this study will improve the accuracy of qRT-PCR quantification of target gene expression in andromonoecious T. rupestris flowers, and will facilitate the functional genomics studies on flower development and sex differentiation in the future.

Selection of Suitable Reference Genes for Quantitative Real-time PCR in Sapium sebiferum

Chinese tallow (Sapium sebiferum L.) is a promising landscape and bioenergy plant. Measuring gene expression by quantitative real-time polymerase chain reaction (qRT-PCR) can provide valuable information on gene function. Stably expressed reference genes for normalization are a prerequisite for ensuring the accuracy of the target gene expression level among different samples. However, the reference genes in Chinese tallow have not been systematically validated. In this study, 12 candidate reference genes (18S, GAPDH, UBQ, RPS15, SAND, TIP41, 60S, ACT7, PDF2, APT, TBP, and TUB) were investigated with qRT-PCR in 18 samples, including those from different tissues, from plants treated with sucrose and cold stresses. The data were calculated with four common algorithms, geNorm, BestKeeper, NormFinder, and the delta cycle threshold (ΔCt). TIP41 and GAPDH were the most stable for the tissue-specific experiment, GAPDH and 60S for cold treatment, and GAPDH and UBQ for sucrose stresses, while the least stable genes were 60S, TIP41, and 18S respectively. The comprehensive results showed APT, GAPDH, and UBQ to be the top-ranked stable genes across all the samples. The stability of 60S was the lowest during all experiments. These selected reference genes were further validated by comparing the expression profiles of the chalcone synthase gene in Chinese tallow in different samples. The results will help to improve the accuracy of gene expression studies in Chinese tallow.

Selection of reliable reference genes for normalization of quantitative RT-PCR from different developmental stages and tissues in amphioxus

Amphioxus is a closest living proxy to the ancestor of cephalochordates with vertebrates, and key animal for novel understanding in the evolutionary origin of vertebrate body plan, genome, tissues and immune system. Reliable analyses using quantitative real-time PCR (qRT-PCR) for answering these scientific questions is heavily dependent on reliable reference genes (RGs). In this study, we evaluated stability of thirteen candidate RGs in qRT-PCR for different developmental stages and tissues of amphioxus by four independent (geNorm, NormFinder, BestKeeper and deltaCt) and one comparative algorithms (RefFinder). The results showed that the top two stable RGs were the following: (1) S20 and 18 S in thirteen developmental stages, (2) EF1A and ACT in seven normal tissues, (3) S20 and L13 in both intestine and hepatic caecum challenged with lipopolysaccharide (LPS), and (4) S20 and EF1A in gill challenged with LPS. The expression profiles of two target genes (EYA and HHEX) in thirteen developmental stages were used to confirm the reliability of chosen RGs. This study identified optimal RGs that can be used to accurately measure gene expression under these conditions, which will benefit evolutionary and functional genomics studies in amphioxus.

The Challenge of Stability in High-Throughput Gene Expression Analysis: Comprehensive Selection and Evaluation of Reference Genes for BALB/c Mice Spleen Samples in the Leishmania infantum Infection Model

The interaction of Leishmania with BALB/c mice induces dramatic changes in transcriptome patterns in the parasite, but also in the target organs (spleen, liver…) due to its response against infection. Real-time quantitative PCR (qPCR) is an interesting approach to analyze these changes and understand the immunological pathways that lead to protection or progression of disease. However, qPCR results need to be normalized against one or more reference genes (RG) to correct for non-specific experimental variation. The development of technical platforms for high-throughput qPCR analysis, and powerful software for analysis of qPCR data, have acknowledged the problem that some reference genes widely used due to their known or suspected “housekeeping” roles, should be avoided due to high expression variability across different tissues or experimental conditions. In this paper we evaluated the stability of 112 genes using three different algorithms: geNorm, NormFinder and RefFinder in spleen samples from BALB/c mice under different experimental conditions (control and Leishmania infantum-infected mice). Despite minor discrepancies in the stability ranking shown by the three methods, most genes show very similar performance as RG (either good or poor) across this massive data set. Our results show that some of the genes traditionally used as RG in this model (i.e. B2m, Polr2a and Tbp) are clearly outperformed by others. In particular, the combination of Il2rg + Itgb2 was identified among the best scoring candidate RG for every group of mice and every algorithm used in this experimental model. Finally, we have demonstrated that using “traditional” vs rationally-selected RG for normalization of gene expression data may lead to loss of statistical significance of gene expression changes when using large-scale platforms, and therefore misinterpretation of results. Taken together, our results highlight the need for a comprehensive, high-throughput search for the most stable reference genes in each particular experimental model.

The absence of pleiotrophin modulates gene expression in the hippocampus in vivo and in cerebellar granule cells in vitro

Pleiotrophin (PTN) is a secreted growth factor recently proposed to act as a neuromodulatory peptide in the Central Nervous System. PTN appears to be involved in neurodegenerative diseases and neural disorders, and it has also been implicated in learning and memory. Specifically, PTN-deficient mice exhibit a lower threshold for LTP induction in the hippocampus, which is attenuated in mice overexpressing PTN. However, there is little information about the signaling systems recruited by PTN to modulate neural activity. To address this issue, the gene expression profile in hippocampus of mice lacking PTN was analyzed using microarrays of 22,000 genes. In addition, we corroborated the effect of the absence of PTN on the expression of these genes by silencing this growth factor in primary neuronal cultures in vitro. The microarray analysis identified 102 genes that are differentially expressed (z-score>3.0) in PTN null mice, and the expression of eight of those modified in the hippocampus of KO mice was also modified in vitro after silencing PTN in cultured neurons with siRNAs. The data obtained indicate that the absence of PTN affects AKT pathway response and modulates the expression of genes related with neuroprotection (Mgst3 and Estrogen receptor 1, Ers 1) and cell differentiation (Caspase 6, Nestin, and Odz4), both in vivo and in vitro.

Identification and validation of reference genes for accurate normalization of real-time quantitative PCR data in kiwifruit

Identification and validation of reference genes are required for the normalization of qPCR data. We studied the expression stability produced by eight primer pairs amplifying four common genes used as references for normalization. Samples representing different tissues, organs and developmental stages in kiwifruit (Actinidia chinensis var. deliciosa (A. Chev.) A. Chev.) were used. A total of 117 kiwifruit samples were divided into five sample sets (mature leaves, axillary buds, stigmatic arms, fruit flesh and seeds). All samples were also analysed as a single set. The expression stability of the candidate primer pairs was tested using three algorithms (geNorm, NormFinder and BestKeeper). The minimum number of reference genes necessary for normalization was also determined. A unique primer pair was selected for amplifying the 18S rRNA gene. The primer pair selected for amplifying the ACTIN gene was different depending on the sample set. 18S 2 and ACT 2 were the candidate primer pairs selected for normalization in the three sample sets (mature leaves, fruit flesh and stigmatic arms). 18S 2 and ACT 3 were the primer pairs selected for normalization in axillary buds. No primer pair could be selected for use as the reference for the seed sample set. The analysis of all samples in a single set did not produce the selection of any stably expressing primer pair. Considering data previously reported in the literature, we validated the selected primer pairs amplifying the FLOWERING LOCUS T gene for use in the normalization of gene expression in kiwifruit.

Selection of Reference Genes for Real-Time Quantitative PCR in Pinus massoniana Post Nematode Inoculation

Pinus massoniaia Lamb has gained more and more attention as the most important tree species for timber and forestation in South China. Gene expression studies are of great importance to identify new and elite cultivars. Real-time quantitative PCR, a highly sensitive and specific method, is commonly used in the analysis of gene expression. The appropriate reference genes must be employed to normalize the calculation program for ascertaining repeatable and significant results. Herein, eleven housekeeping genes were evaluated during different stages of P. massoniana post nematode inoculation in this study. Three statistical approaches such as geNorm, NormFinder and BestKeeper were selected to analyze the stability of candidate genes. The results indicated that U2af and β-TUB were the most stable reference genes. These two genes could be used for the normalization in most of the experiments of P. massoniana, while Histone and AK were the least stable ones. In addition, EF expressed at the lowest average Ct value was the most abundant candidate gene. As an important gene associated with defense mechanisms, ABC transporter was analyzed by qRT-PCR, and the results were used to confirm the reliability of two genes. The selected reference genes in the present study will be conducive to future gene expression normalized by qRT-PCR in P. massoniana.

Identification of Suitable Reference Genes for Gene Expression Normalization in the Quantitative Real-Time PCR Analysis of Sweet Osmanthus (Osmanthus fragrans Lour.)

Quantitative real-time PCR (RT-qPCR), a sensitive technique for quantifying gene expression, depends on the stability of the reference gene(s) used for data normalization. Several studies examining the selection of reference genes have been performed in ornamental plants but none in sweet osmanthus (Osmanthus fragrans Lour.). Based on transcriptomic sequencing data from O. fragrans buds at four developmental stages, six reference genes (OfACT, OfEF1α, OfIDH, OfRAN1, OfTUB, and OfUBC2) with stable expression (0.5 to 2 fold change in expression levels between any two developmental stages), as well as the commonly used reference gene Of18S, were selected as candidates for gene expression normalization in the RT-qPCR analysis of O. fragrans. For the normalization of RT-qPCR with two dyes, SYBR Green and EvaGreen, the expressional stability of seven candidate reference genes in 43 O. fragrans samples was analyzed using geNorm, NormFinder and BestKeeper. For RT-qPCR using SYBR Green, OfRAN1 and OfUBC2 were the optimal reference genes for all samples and different cultivars, OfACT and OfEF1α were suitable for different floral developmental stages, and OfACT was the optimal reference gene for different temperature treatments. The geometric mean values of the optimal reference gene pairs for the normalization of RT-qPCR are recommended to be used for all samples, different cultivars and different floral developmental stages in O. fragrans. For RT-qPCR using EvaGreen, OfUBC2 was the optimal reference gene for all samples and different cultivars, and OfACT was the optimal reference gene for different floral developmental stages and different temperature treatments. As the worst reference gene, Of18S should not be used as a reference gene in O. fragrans in the future. Our results provide a reference gene application guideline for O. fragrans gene expression characterization using RT-qPCR.

Bacterial reference genes for gene expression studies by RT-qPCR: survey and analysis

The appropriate choice of reference genes is essential for accurate normalization of gene expression data obtained by the method of reverse transcription quantitative real-time PCR (RT-qPCR). In 2009, a guideline called the Minimum Information for Publication of Quantitative Real-Time PCR Experiments (MIQE) highlighted the importance of the selection and validation of more than one suitable reference gene for obtaining reliable RT-qPCR results. Herein, we searched the recent literature in order to identify the bacterial reference genes that have been most commonly validated in gene expression studies by RT-qPCR (in the first 5 years following publication of the MIQE guidelines). Through a combination of different search parameters with the text mining tool MedlineRanker, we identified 145 unique bacterial genes that were recently tested as candidate reference genes. Of these, 45 genes were experimentally validated and, in most of the cases, their expression stabilities were verified using the software tools geNorm and NormFinder. It is noteworthy that only 10 of these reference genes had been validated in two or more of the studies evaluated. An enrichment analysis using Gene Ontology classifications demonstrated that genes belonging to the functional categories of DNA Replication (GO: 0006260) and Transcription (GO: 0006351) rendered a proportionally higher number of validated reference genes. Three genes in the former functional class were also among the top five most stable genes identified through an analysis of gene expression data obtained from the Pathosystems Resource Integration Center. These results may provide a guideline for the initial selection of candidate reference genes for RT-qPCR studies in several different bacterial species.

Tracking the best reference genes for RT-qPCR data normalization in filamentous fungi

BACKGROUND

A critical step in the RT-qPCR workflow for studying gene expression is data normalization, one of the strategies being the use of reference genes. This study aimed to identify and validate a selection of reference genes for relative quantification in Talaromyces versatilis, a relevant industrial filamentous fungus. Beyond T. versatilis, this study also aimed to propose reference genes that are applicable more widely for RT-qPCR data normalization in filamentous fungi.

RESULTS

A selection of stable, potential reference genes was carried out in silico from RNA-seq based transcriptomic data obtained from T. versatilis. A dozen functionally unrelated candidate genes were analysed by RT-qPCR assays over more than 30 relevant culture conditions. By using geNorm, we showed that most of these candidate genes had stable transcript levels in most of the conditions, from growth environments to conidial germination. The overall robustness of these genes was explored further by showing that any combination of 3 of them led to minimal normalization bias. To extend the relevance of the study beyond T. versatilis, we challenged their stability together with sixteen other classically used genes such as β-tubulin or actin, in a representative sample of about 100 RNA-seq datasets. These datasets were obtained from 18 phylogenetically distant filamentous fungi exposed to prevalent experimental conditions. Although this wide analysis demonstrated that each of the chosen genes exhibited sporadic up- or down-regulation, their hierarchical clustering allowed the identification of a promising group of 6 genes, which presented weak expression changes and no tendency to up- or down-regulation over the whole set of conditions. This group included ubcB, sac7, fis1 and sarA genes, as well as TFC1 and UBC6 that were previously validated for their use in S. cerevisiae.

CONCLUSIONS

We propose a set of 6 genes that can be used as reference genes in RT-qPCR data normalization in any field of fungal biology. However, we recommend that the uniform transcription of these genes is tested by systematic experimental validation and to use the geometric averaging of at least 3 of the best ones. This will minimize the bias in normalization and will support trustworthy biological conclusions.

Reference genes in real-time PCR

This paper aims to discuss various aspects of the use of reference genes in qPCR technique used in the thousands of present studies. Most frequently, these are housekeeping genes and they must meet several criteria so that they can lay claim to the name. Lots of papers report that in different conditions, for different organisms and even tissues the basic assumption—the constant level of the expression is not maintained for many genes that seem to be perfect candidates. Moreover, their transcription can not be affected by experimental factors. Sounds simple and clear but a great number of designed protocols and lack of consistency among them brings confusion on how to perform experiment properly. Since during selection of the most stable normalizing gene we can not use any reference gene, different ways and algorithms for their selection were developed. Such methods, including examples of best normalizing genes in some specific cases and possible mistakes are presented based on available sources. Numerous examples of reference genes applications, which are usually in too few numbers in relevant articles not allowing to make a solid fundament for a reader, will be shown along with instructive compilations to make an evidence for presented statements and an arrangement of future qPCR experiments. To include all the pitfalls and problems associated with the normalization methods there is no way not to begin from sample preparation and its storage going through candidate gene selection, primer design and statistical analysis. This is important because numerous short reviews available cover the topic only in lesser extent at the same time giving the reader false conviction of complete topic recognition.

Identification of two putative reference genes from grapevine suitable for gene expression analysis in berry and related tissues derived from RNA-Seq data

Background

Data normalization is a key step in gene expression analysis by qPCR. Endogenous control genes are used to estimate variations and experimental errors occurring during sample preparation and expression measurements. However, the transcription level of the most commonly used reference genes can vary considerably in samples obtained from different individuals, tissues, developmental stages and under variable physiological conditions, resulting in a misinterpretation of the performance of the target gene(s). This issue has been scarcely approached in woody species such as grapevine.

Results

A statistical criterion was applied to select a sub-set of 19 candidate reference genes from a total of 242 non-differentially expressed (NDE) genes derived from a RNA-Seq experiment comprising ca. 500 million reads obtained from 14 table-grape genotypes sampled at four phenological stages. From the 19 candidate reference genes, VvAIG1 (AvrRpt2-induced gene) and VvTCPB (T-complex 1 beta-like protein) were found to be the most stable ones after comparing the complete set of genotypes and phenological stages studied. This result was further validated by qPCR and geNorm analyses.

Conclusions

Based on the evidence presented in this work, we propose to use the grapevine genes VvAIG1 or VvTCPB or both as a reference tool to normalize RNA expression in qPCR assays or other quantitative method intended to measure gene expression in berries and other tissues of this fruit crop, sampled at different developmental stages and physiological conditions.

Label-free mass spectrometry exploits dozens of detected peptides to quantify lamins in wildtype and knockdown cells

Label-free quantitation and characterization of proteins by mass spectrometry (MS) is now feasible, especially for moderately expressed structural proteins such as lamins that typically yield dozens of tryptic peptides from tissue cells. Using standard cell culture samples, we describe general algorithms for quantitative analysis of peptides identified in liquid chromatography tandem mass spectrometry (LC-MS/MS). The algorithms were foundational to the discovery that the absolute stoichiometry of A-type to B-type lamins scales with tissue stiffness (Swift et al., Science 2013). Isoform dominance helps make sense of why mutations and changes with age of mechanosensitive lamin-A,C only affect "stiff" tissues such as heart, muscle, bone, or even fat, but not brain. A Peak Ratio Fingerprinting (PRF) algorithm is elaborated here through its application to lamin-A,C knockdown. After demonstrating the large dynamic range of PRF using calibrated mixtures of human and mouse lysates, we validate measurements of partial knockdown with standard cell biology analyses using quantitative immunofluorescence and immunoblotting. Optimal sets of MS-detected peptides as determined by PRF demonstrate that the strongest peptide signals are not necessarily the most reliable for quantitation. After lamin-A,C knockdown, PRF computes an invariant set of "housekeeping" proteins as part of a broader proteomic analysis that also shows the proteome of mesenchymal stem cells (MSCs) is more broadly perturbed than that of a human epithelial cancer line (A549s), with particular variation in nuclear and cytoskeletal proteins. These methods offer exciting prospects for basic and clinical studies of lamin-A,C as well as other MS-detectable proteins.

How to perform RT-qPCR accurately in plant species? A case study on flower colour gene expression in an azalea (Rhododendron simsii hybrids) mapping population

Background

Flower colour variation is one of the most crucial selection criteria in the breeding of a flowering pot plant, as is also the case for azalea (Rhododendron simsii hybrids). Flavonoid biosynthesis was studied intensively in several species. In azalea, flower colour can be described by means of a 3-gene model. However, this model does not clarify pink-coloration. The last decade gene expression studies have been implemented widely for studying flower colour. However, the methods used were often only semi-quantitative or quantification was not done according to the MIQE-guidelines. We aimed to develop an accurate protocol for RT-qPCR and to validate the protocol to study flower colour in an azalea mapping population.

Results

An accurate RT-qPCR protocol had to be established. RNA quality was evaluated in a combined approach by means of different techniques e.g. SPUD-assay and Experion-analysis. We demonstrated the importance of testing noRT-samples for all genes under study to detect contaminating DNA. In spite of the limited sequence information available, we prepared a set of 11 reference genes which was validated in flower petals; a combination of three reference genes was most optimal. Finally we also used plasmids for the construction of standard curves. This allowed us to calculate gene-specific PCR efficiencies for every gene to assure an accurate quantification. The validity of the protocol was demonstrated by means of the study of six genes of the flavonoid biosynthesis pathway. No correlations were found between flower colour and the individual expression profiles. However, the combination of early pathway genes (CHS, F3H, F3'H and FLS) is clearly related to co-pigmentation with flavonols. The late pathway genes DFR and ANS are to a minor extent involved in differentiating between coloured and white flowers. Concerning pink coloration, we could demonstrate that the lower intensity in this type of flowers is correlated to the expression of F3'H.

Conclusions

Currently in plant research, validated and qualitative RT-qPCR protocols are still rare. The protocol in this study can be implemented on all plant species to assure accurate quantification of gene expression. We have been able to correlate flower colour to the combined regulation of structural genes, both in the early and late branch of the pathway. This allowed us to differentiate between flower colours in a broader genetic background as was done so far in flower colour studies. These data will now be used for eQTL mapping to comprehend even more the regulation of this pathway.

Normalizing gene expression by quantitative PCR during somatic embryogenesis in two representative conifer species: Pinus pinaster and Picea abies

Suitable internal control genes to normalize qPCR data from different stages of embryo development and germination were identified in two representative conifer species. Clonal propagation by somatic embryogenesis has a great application potentiality in conifers. Quantitative PCR (qPCR) is widely used for gene expression analysis during somatic embryogenesis and embryo germination. No single reference gene is universal, so a systematic characterization of endogenous genes for concrete conditions is fundamental for accuracy. We identified suitable internal control genes to normalize qPCR data obtained at different steps of somatic embryogenesis (embryonal mass proliferation, embryo maturation and germination) in two representative conifer species, Pinus pinaster and Picea abies. Candidate genes included endogenous genes commonly used in conifers, genes previously tested in model plants, and genes with a lower variation of the expression along embryo development according to genome-wide transcript profiling studies. Three different algorithms were used to evaluate expression stability. The geometric average of the expression values of elongation factor-1α, α-tubulin and histone 3 in P. pinaster, and elongation factor-1α, α-tubulin, adenosine kinase and CAC in P. abies were adequate for expression studies throughout somatic embryogenesis. However, improved accuracy was achieved when using other gene combinations in experiments with samples at a single developmental stage. The importance of studies selecting reference genes to use in different tissues or developmental stages within one or close species, and the instability of commonly used reference genes, is highlighted.

Identification and testing of reference genes for Sesame gene expression analysis by quantitative real-time PCR

Sesame (Sesamum indicum L.) is an ancient and important oilseed crop. However, few sesame reference genes have been selected for quantitative real-time PCR until now. Screening and validating reference genes is a requisite for gene expression normalization in sesame functional genomics research. In this study, ten candidate reference genes, i.e., SiACT, SiUBQ6, SiTUB, Si18S rRNA, SiEF1α, SiCYP, SiHistone, SiDNAJ, SiAPT and SiGAPDH, were chosen and examined systematically in 32 sesame samples. Three qRT-PCR analysis methods, i.e., geNorm, NormFinder and BestKeeper, were evaluated systematically. Results indicated that all ten candidate reference genes could be used as reference genes in sesame. SiUBQ6 and SiAPT were the optimal reference genes for sesame plant development; SiTUB was suitable for sesame vegetative tissue development, SiDNAJ for pathogen treatment, SiHistone for abiotic stress, SiUBQ6 for bud development and SiACT for seed germination. As for hormone treatment and seed development, SiHistone, SiCYP, SiDNAJ or SiUBQ6, as well as SiACT, SiDNAJ, SiTUB or SiAPT, could be used as reference gene, respectively. To illustrate the suitability of these reference genes, we analyzed the expression variation of three functional sesame genes of SiSS, SiLEA and SiGH in different organs using the optimal qRT-PCR system for the first time. The stability levels of optimal and worst reference genes screened for seed development, anther sterility and plant development were validated in the qRT-PCR normalization. Our results provided a reference gene application guideline for sesame gene expression characterization using qRT-PCR system.

In silico estimation of translation efficiency in human cell lines: potential evidence for widespread translational control

Recently large scale transcriptome and proteome datasets for human cells have become available. A striking finding from these studies is that the level of an mRNA typically predicts no more than 40% of the abundance of protein. This correlation represents the overall figure for all genes. We present here a bioinformatic analysis of translation efficiency – the rate at which mRNA is translated into protein. We have analysed those human datasets that include genome wide mRNA and protein levels determined in the same study. The analysis comprises five distinct human cell lines that together provide comparable data for 8,170 genes. For each gene we have used levels of mRNA and protein combined with protein stability data from the HeLa cell line to estimate translation efficiency. This was possible for 3,990 genes in one or more cell lines and 1,807 genes in all five cell lines. Interestingly, our analysis and modelling shows that for many genes this estimated translation efficiency has considerable consistency between cell lines. Some deviations from this consistency likely result from the regulation of protein degradation. Others are likely due to known translational control mechanisms. These findings suggest it will be possible to build improved models for the interpretation of mRNA expression data. The results we present here provide a view of translation efficiency for many genes. We provide an online resource allowing the exploration of translation efficiency in genes of interest within different cell lines (http://bioanalysis.otago.ac.nz/TranslationEfficiency).

Quantification of yeast and bacterial gene transcripts in retail cheeses by reverse transcription-quantitative PCR

The cheese microbiota contributes to a large extent to the development of the typical color, flavor, and texture of the final product. Its composition is not well defined in most cases and varies from one cheese to another. The aim of the present study was to establish procedures for gene transcript quantification in cheeses by reverse transcription-quantitative PCR. Total RNA was extracted from five smear-ripened cheeses purchased on the retail market, using a method that does not involve prior separation of microbial cells. 16S rRNA and malate:quinone oxidoreductase gene transcripts of Corynebacterium casei, Brevibacterium aurantiacum, and Arthrobacter arilaitensis and 26S rRNA and beta tubulin gene transcripts of Geotrichum candidum and Debaryomyces hansenii could be detected and quantified in most of the samples. Three types of normalization were applied: against total RNA, against the amount of cheese, and against a reference gene. For the first two types of normalization, differences of reverse transcription efficiencies from one sample to another were taken into account by analysis of exogenous control mRNA. No good correlation was found between the abundances of target mRNA or rRNA transcripts and the viable cell concentration of the corresponding species. However, in most cases, no mRNA transcripts were detected for species that did not belong to the dominant species. The applications of gene expression measurement in cheeses containing an undefined microbiota, as well as issues concerning the strategy of normalization and the assessment of amplification specificity, are discussed.