Dataset Information

Custom selected reference genes outperform pre-defined reference genes in transcriptomic analysis.

ABSTRACT:

Background

RNA sequencing allows the measuring of gene expression at a resolution unmet by expression arrays or RT-qPCR. It is however necessary to normalize sequencing data by library size, transcript size and composition, among other factors, before comparing expression levels. The use of internal control genes or spike-ins is advocated in the literature for scaling read counts, but the methods for choosing reference genes are mostly targeted at RT-qPCR studies and require a set of pre-selected candidate controls or pre-selected target genes.

Results

Here, we report an R-based pipeline to select internal control genes based solely on read counts and gene sizes. This novel method first normalizes the read counts to Transcripts per Million (TPM) and then excludes weakly expressed genes using the DAFS script to calculate the cut-off. It then selects as references the genes with lowest TPM covariance. We used this method to pick custom reference genes for the differential expression analysis of three transcriptome sets from transgenic Arabidopsis plants expressing heterologous fungal effector proteins tagged with GFP (using GFP alone as the control). The custom reference genes showed lower covariance and fold change as well as a broader range of expression levels than commonly used reference genes. When analyzed with NormFinder, both typical and custom reference genes were considered suitable internal controls, but the custom selected genes were more stably expressed. geNorm produced a similar result in which most custom selected genes ranked higher (i.e. were more stably expressed) than commonly used reference genes.

Conclusions

The proposed method is innovative, rapid and simple. Since it does not depend on genome annotation, it can be used with any organism, and does not require pre-selected reference candidates or target genes that are not always available.

SUBMITTER: Dos Santos KCG

PROVIDER: S-EPMC6954607 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundAustism spectrum disorder (ASD) is a heterogeneous behavioral disorder or condition characterized by severe impairment of social engagement and the presence of repetitive activities. The molecular etiology of ASD is still largely unknown despite a strong genetic component. Part of the difficulty in turning genetics into disease mechanisms and potentially new therapeutics is the sheer number and diversity of the genes that have been associated with ASD and ASD symptoms. The goal of this work is to use shRNA-generated models of genetic defects proposed as causative for ASD to identify the common pathways that might explain how they produce a core clinical disability.MethodsTranscript levels of Mecp2, Mef2a, Mef2d, Fmr1, Nlgn1, Nlgn3, Pten, and Shank3 were knocked-down in mouse primary neuron cultures using shRNA constructs. Whole genome expression analysis was conducted for each of the knockdown cultures as well as a mock-transduced culture and a culture exposed to a lentivirus expressing an anti-luciferase shRNA. Gene set enrichment and a causal reasoning engine was employed to identify pathway level perturbations generated by the transcript knockdown.ResultsQuantification of the shRNA targets confirmed the successful knockdown at the transcript and protein levels of at least 75% for each of the genes. After subtracting out potential artifacts caused by viral infection, gene set enrichment and causal reasoning engine analysis showed that a significant number of gene expression changes mapped to pathways associated with neurogenesis, long-term potentiation, and synaptic activity.ConclusionsThis work demonstrates that despite the complex genetic nature of ASD, there are common molecular mechanisms that connect many of the best established autism candidate genes. By identifying the key regulatory checkpoints in the interlinking transcriptional networks underlying autism, we are better able to discover the ideal points of intervention that provide the broadest efficacy across the diverse population of autism patients.

Project description:BackgroundReal time RT-PCR (qPCR) is a useful and powerful tool for quantitative measurement of gene expression. The proper choice of internal standards such as reference genes is crucial for correct data evaluation. In female dogs, as in other species, the reproductive tract is continuously undergoing hormonal and cycle stage-dependent morphological changes, which are associated with altered gene expression. However, there have been few attempts published so far targeted to the dog aimed at determining optimal reference genes for the reproductive organs. Most of these approaches relied on genes previously described in other species. Large-scale transcriptome-based experiments are promising tools for defining potential candidate reference genes, but were never considered in this context in canine research.ResultsHere, using available microarray and RNA-seq datasets derived from reproductive organs (corpus luteum, placenta, healthy and diseased uteri) of dogs, we have performed multistudy analysis to identify the most stably expressed genes for expression studies, in each tissue separately and collectively for different tissues. The stability of newly identified reference genes (EIF4H, KDELR2, KDM4A and PTK2) has been determined and ranked relative to previously used reference genes, i.e., GAPDH, β-actin and cyclophillin A/PPIA, using RefFinder and NormFinder algorithms. Finally, expression of selected target genes (luteal IL-1b and MHCII, placental COX2 and VEGFA, and uterine IGF2 and LHR) was re-evaluated and normalized. All proposed candidate reference genes were more stable, ranked higher and introduced less variation than previously used genes.ConclusionsBased on our analyses, we recommend applying KDM4A and PTK2 for normalization of gene expression in the canine CL and placenta. The inclusion of a third reference gene, EIF4H, is suggested for healthy uteri. With this, the interpretation of qPCR data will be more reliable, allowing better understanding of canine reproductive physiology.