Validation of endogenous reference genes for qRT-PCR analysis of human visceral adipose samples

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 genetic loci associated with different responses to high-fat diet-induced obesity in C57BL/6N and C57BL/6J substrains

We have recently demonstrated that C57BL/6NTac and C57BL/6JRj substrains are significantly different in their response to high-fat diet-induced obesity (DIO). The C57BL/6JRj substrain seems to be protected from DIO and genetic differences between C57BL/6J and C57BL/6N substrains at 11 single nucleotide polymorphism (SNP) loci have been identified. To define genetic variants as well as differences in parameters of glucose homeostasis and insulin sensitivity between C57BL/6NTac and C57BL/6JRj substrains that may explain the different response to DIO, we analyzed 208 first backcross (BC1) hybrids of C57BL/6NTac and C57BL/6JRj [(C57BL/6NTac × C57BL/6JRj)F1 × C57BL/6NTac] mice. Body weight, epigonadal and subcutaneous fat mass, circulating leptin, as well as parameters of glucose metabolism were measured after 10 wk of high-fat diet (HFD). Genetic profiling of BC1 hybrids were performed using TaqMan SNP genotyping assays. Furthermore, to assess whether SNP polymorphisms could affect mRNA level, we carried out gene expression analysis in murine liver samples. Human subcutaneous adipose tissue was used to verify murine data of SNAP29. We identified four sex-specific variants that are associated with the extent of HFD-induced weight gain and fat depot mass. BC1 hybrids carrying the combination of risk or beneficial alleles exhibit the phenotypical extremes of the parental strains. Murine and human SC expression analysis revealed Snap29 as strongest candidate. Our data indicate an important role of these loci in responsiveness to HFD-induced obesity and suggest genes of the synaptic vesicle release system such as Snap29 being involved in the regulation of high-fat DIO.

Validation of reference genes for normalization gene expression in reverse transcription quantitative PCR in human normal thyroid and goiter tissue

Reverse transcription quantitative polymerase chain reaction (RT-qPCR) has been recognized as the most accurate method for quantifying mRNA transcripts, but normalization of samples is a prerequisite for correct data interpretation. So, this study aimed to evaluate the most stable reference gene for RT-qPCR in human normal thyroid and goiter tissues. Beta-actin (ACTB); glyceraldehyde-3-phosphate dehydrogenase (GAPDH); succinate dehydrogenase, subunit A, flavoprotein (Fp) (SDHA); hypoxanthine phosphoribosyltransferase I (HPRTI); tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide (YWHAZ); and beta-2-microglobulin (B2M) were evaluated in 14 thyroid tissue samples (7 normal and 7 goiter tissues) by RT-qPCR. The mean Cq and the maximum fold change (MFC) and NormFinder software were used to assess the stability of the genes. As a result, ACTB gene was more stable than GAPDH, SDHA, HPRTI, YWHAZ, and B2M. In conclusion, ACTB could be used to normalize RT-qPCR data in normal thyroid and goiter tissues.

The expression of LEP, LEPR, IGF1 and IL10 in obesity and the relationship with microRNAs

Obesity is a multifactorial disease, with epigenetic alterations. Have been described modifications in the expression of some microRNAs, and some proteins related to obesity. The objective was to determine and correlate, in obese patients, the gene expression of LEP, LEPR, IGF1, IL10 and of miR-27a, miR-27b, miR-143 and miR-145. RNA was extracted from biopsies of subcutaneous fat, liver and visceral fat of 15 obese subjects submitted to bariatric surgery and of 15 non-obese subjects submitted to cholecystectomy for cDNA synthesis and for RT-PCR. The microRNAs were chosen using the TargetScan software. An increased expression of LEP and IGF1 was detected in the subcutaneous fat of the obese group compared to control, while the expression of IGF1 was higher in the control group than in the obese one. MiRNA-27a had a higher expression in the omentum of the obese patients and there was also a correlation in the expression of miRNA-145 and LEPR in the omentum of this group.

Identification of housekeeping genes suitable for gene expression analysis in the pearl mussel, Hyriopsis cumingii, during biomineralization

Quantitative real-time polymerase chain reaction is a sensitive technique for quantifying gene expression levels. One or more appropriate reference genes must be selected to accurately compare mRNA transcript levels across different samples and tissues. The freshwater pearl, Hyriopsis cumingii (Lea), is an important economic species cultured in China. To date, no reference genes for gene expression analysis in this species have been validated. This study aimed to compare the relative expression of seven housekeeping genes across different tissue types and in the mantle or pearl sac during three biomineralization processes: seasonal shell growth, shell healing and pearl-sac formation in H. cumingii. Three programs evaluated the expression stabilities of the seven genes: BestKeeper, geNorm and NormFinder. The beta actin gene (ACTB), commonly used as a housekeeping gene in many studies, was the least stable. The expressions of Ubiquitin (Ubi) and Ribosomal protein L18 (Rpl18) and Elongation factor 1-alpha (EF1α) were more stable than the remaining four genes. Therefore, we suggest that Ubi, Rpl18 and EF1α are suitable reference genes. The three selected reference genes are expected to facilitate analysis of gene expressions during shell or pearl formation in H. cumingii.

A time course study demonstrating mRNA, microRNA, 18S rRNA, and U6 snRNA changes to estimate PMI in deceased rat's spleen

Determining the postmortem interval (PMI) is important in criminal, civil, and forensic cases. We examined the feasibility of using the transcript abundances of mRNAs, 18S rRNA, U6 snRNA, and microRNAs as a means to estimate the PMI. We removed spleen tissues from rats at different PMIs under 4°C or 25°C and examined gene transcript abundances in these samples by RT-qPCR. Using the algorithm geNorm, we found that microRNAs to be appropriate control markers because they were less affected by PMI and temperature. We also characterized relationships between observed PMI and the transcript levels of the above-mentioned RNAs. GAPDH1 and ACTB1 fluctuated slightly like cubic curves, while GAPDH2 and ACTB2 decreased rapidly. 18S rRNA transcript level exhibited a parabolic-like trend at 25°C and exponential growth at 4°C, while U6 transcript level exhibited exponential decay at 25°C and a parabolic-like trend at 4°C. Following validation, we conclude that GAPDH2, ACTB2, and 18S rRNA are suitable makers in the accurate determination of PMI.

Rosiglitazone and fenofibrate exacerbate liver steatosis in a mouse model of obesity and hyperlipidemia. A transcriptomic and metabolomic study

Peroxisome proliferator-activated receptors (PPAR) play an important role in the regulation of lipid and glucose metabolism, inflammatory, and vascular responses. We show the effect of treatment with two PPAR agonists, fenofibrate (FF) and rosiglitazone (RSG), on ob/ob and LDLR-double deficient mice, by combined gene-expression and metabolomic analyses. Male mice were daily treated for 12 weeks with RSG (10 mg·kg(1-)·day(-1) per os (p.o.), n = 8) and FF (50 mg·kg(1-)·day(-1) p.o., n = 8). Twelve untreated ob/ob and LDLR-double deficient mice were used as controls. To integrate the transcriptomic and metabolomic results, we designed a hierarchical algorithm, based on the average linkage method in clustering. Data were also interpreted with the Ingenuity Pathway Analysis program. FF and RSG treatments significantly increased the hepatic triglyceride content in the liver when compared with the control group, and the treatments induced an increase in the number and size of hepatic lipid droplets. Both drugs simultaneously activate pro-steatotic and antisteatotic metabolic pathways with a well-ordered result of aggravation of the hepatic lipid accumulation. The present study is a cautionary note not only to researchers on the basic mechanism of the action of PPAR activators but also to the use of these compounds in clinical practice.

Expression of inflammation-related genes is associated with adipose tissue location in horses

Background

In humans, adipose tissue (AT) originating from different depots shows varying gene expression profiles. In horses, the risk of certain metabolic disorders may also be influenced by the impact of specific AT depots. Macrophage infiltration in human and rat AT is considered to be a source of inflammatory changes. In horses, this relationship has not been extensively studied yet. Reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR), a useful method to evaluate differences in mRNA expression across different tissues, can be used to evaluate differences between equine AT depots. For a correct interpretation of the RT-qPCR results, expression data have to be normalized by the use of validated reference genes. The main objectives of this study were to compare mRNA expression of inflammation-related genes, as well as adipocyte morphology and number between different equine AT depots; and in addition, to investigate the presence of antigen presenting cells in equine AT and any potential relationship with adipokine mRNA expression.

Results

In this study, the mRNA expression of inflammation-related genes (leptin, chemokine ligand 5, interleukin 1β, interleukin 6, interleukin 10, adiponectin, matrix metalloproteinase 2, and superoxide dismutase 2) and candidate reference gene stability was investigated in 8 different AT depots collected from the nuchal, abdominal (mesenteric, retroperitoneal, and peri-renal) and subcutaneous (tail head and loin) AT region. By using GeNorm analysis, HPRT1, RPL32, and GAPDH were found to be the most stable genes in equine AT. The mRNA expression of leptin, chemokine ligand 5, interleukin 10, interleukin 1β, adiponectin, and matrix metalloproteinase 2 significantly differed across AT depots (P < 0.05). No significant AT depot effect was found for interleukin 6 and superoxide dismutase 2 (P > 0.05). Adipocyte area and number of antigen presenting cells per adipocyte significantly differed between AT depots (P < 0.05).

Conclusions

Adipose tissue location was associated with differences in mRNA expression of inflammation-related genes. This depot-specific difference in mRNA expression suggests that the overall inflammatory status of horses could be partially determined by the relative proportion of the different AT depots.

Mitochondrial Ca(2+) uniporter (MCU)-dependent and MCU-independent Ca(2+) channels coexist in the inner mitochondrial membrane

A protein referred to as CCDC109A and then renamed to mitochondrial calcium uniporter (MCU) has recently been shown to accomplish mitochondrial Ca²⁺ uptake in different cell types. In this study, we investigated whole-mitoplast inward cation currents and single Ca²⁺ channel activities in mitoplasts prepared from stable MCU knockdown HeLa cells using the patch-clamp technique. In whole-mitoplast configuration, diminution of MCU considerably reduced inward Ca²⁺ and Na⁺ currents. This was accompanied by a decrease in occurrence of single channel activity of the intermediate conductance mitochondrial Ca²⁺ current (i-MCC). However, ablation of MCU yielded a compensatory 2.3-fold elevation in the occurrence of the extra large conductance mitochondrial Ca²⁺ current (xl-MCC), while the occurrence of bursting currents (b-MCC) remained unaltered. These data reveal i-MCC as MCU-dependent current while xl-MCC and b-MCC seem to be rather MCU-independent, thus, pointing to the engagement of at least two molecularly distinct mitochondrial Ca²⁺ channels.

Selection and validation of reference genes for gene expression analysis in apomictic and sexual Cenchrus ciliaris

Background

Apomixis is a naturally occurring asexual mode of seed reproduction resulting in offspring genetically identical to the maternal plant. Identifying differential gene expression patterns between apomictic and sexual plants is valuable to help deconstruct the trait. Quantitative RT-PCR (qRT-PCR) is a popular method for analyzing gene expression. Normalizing gene expression data using proper reference genes which show stable expression under investigated conditions is critical in qRT-PCR analysis. We used qRT-PCR to validate expression and stability of six potential reference genes (EF1alpha, EIF4A, UBCE, GAPDH, ACT2 and TUBA) in vegetative and reproductive tissues of B-2S and B-12-9 accessions of C. ciliaris.

Findings

Among tissue types evaluated, EF1alpha showed the highest level of expression while TUBA showed the lowest. When all tissue types were evaluated and compared between genotypes, EIF4A was the most stable reference gene. Gene expression stability for specific ovary stages of B-2S and B-12-9 was also determined. Except for TUBA, all other tested reference genes could be used for any stage-specific ovary tissue normalization, irrespective of the mode of reproduction.

Conclusion

Our gene expression stability assay using six reference genes, in sexual and apomictic accessions of C. ciliaris, suggests that EIF4A is the most stable gene across all tissue types analyzed. All other tested reference genes, with the exception of TUBA, could be used for gene expression comparison studies between sexual and apomictic ovaries over multiple developmental stages. This reference gene validation data in C. ciliaris will serve as an important base for future apomixis-related transcriptome data validation.

Reliable reference genes for normalization of gene expression in cucumber grown under different nitrogen nutrition

In plants, nitrogen is the most important nutritional factor limiting the yield of cultivated crops. Since nitrogen is essential for synthesis of nucleotides, amino acids and proteins, studies on gene expression in plants cultivated under different nitrogen availability require particularly careful selection of suitable reference genes which are not affected by nitrogen limitation. Therefore, the objective of this study was to select the most reliable reference genes for qPCR analysis of target cucumber genes under varying nitrogen source and availability. Among twelve candidate cucumber genes used in this study, five are highly homologous to the commonly used internal controls, whereas seven novel candidates were previously identified through the query of the cucumber genome. The expression of putative reference genes and the target CsNRT1.1 gene was analyzed in roots, stems and leaves of cucumbers grown under nitrogen deprivation, varying nitrate availability or different sources of nitrogen (glutamate, glutamine or NH3). The stability of candidate genes expression significantly varied depending on the tissue type and nitrogen supply. However, in most of the outputs genes encoding CACS, TIP41, F-box protein and EFα proved to be the most suitable for normalization of CsNRT1.1 expression. In addition, our results suggest the inclusion of 3 or 4 references to obtain highly reliable results of target genes expression in all cucumber organs under nitrogen-related stress.

Selection and validation of reference genes for quantitative gene expression studies by real-time PCR in eggplant (Solanum melongena L)

BACKGROUND

Analysis of gene expression patterns leads to functional understanding of biological processes. Quantitative real-time PCR has become the most commonly used technique for in-depth studies of gene expression. To quantify variation in specific gene expression, accurate and reliable normalization across different samples and tissues is necessary. This can be achieved by selecting one or more suitable reference genes to compare the target mRNA transcript levels. In the present work, we illustrate the first evaluation of potential internal control or reference genes across different developmental stages of eggplant for reliable quantification of transcripts by real-time PCR.

RESULTS

We have evaluated the stability in expression of six candidate reference genes (18s rRNA, APRT, GAPDH, Cyclophilin, Actin, and RuBP) in a set of tissues representing six developmental stages of eggplant. The candidate genes were cloned from cDNA and analysed by real-time PCR. The expression data analyzed by three statistical methods (geNorm, NormFinder and BestKeeper) identified 18s rRNA, Cyclophilin and APRT as the most stable and suitable reference genes in eggplant. This was further confirmed in four different varieties, two representative lines of transgenic eggplant as well as in nematode infected eggplant.

CONCLUSION

18s rRNA, Cyclophilin and APRT have been found to be appropriate for the normalization of real-time PCR data for gene expression studies in eggplant.

Identification and validation of reference genes for transcript normalization in strawberry (Fragaria × ananassa) defense responses

Strawberry (Fragaria spp) is an emerging model for the development of basic genomics and recombinant DNA studies among rosaceous crops. Functional genomic and molecular studies involve relative quantification of gene expression under experimental conditions of interest. Accuracy and reliability are dependent upon the choice of an optimal reference control transcript. There is no information available on validated endogenous reference genes for use in studies testing strawberry-pathogen interactions. Thirteen potential pre-selected strawberry reference genes were tested against different tissues, strawberry cultivars, biotic stresses, ripening and senescent conditions, and SA/JA treatments. Evaluation of reference candidate's suitability was analyzed by five different methodologies, and information was merged to identify best reference transcripts. A combination of all five methods was used for selective classification of reference genes. The resulting superior reference genes, FaRIB413, FaACTIN, FaEF1α and FaGAPDH2 are strongly recommended as control genes for relative quantification of gene expression in strawberry. This report constitutes the first systematic study to identify and validate optimal reference genes for accurate normalization of gene expression in strawberry plant defense response studies.

Sequencing and validation of housekeeping genes for quantitative real-time PCR during the gonadotrophic cycle of Diploptera punctata

Background

Quantitative RT-PCR (q-RT-PCR) is a powerful tool that allows for the large scale analysis of small changes in gene expression. Accurate and reliable results depend on the use of stable reference genes for normalization. However, the expression of some widely used housekeeping genes can vary under different experimental setups. To our knowledge, no validation studies have been reported for reference genes in cockroaches. The aim of the current study is the identification and validation of a set of eight housekeeping genes during the first gonadotrophic cycle of the cockroach, Diploptera punctata. This study made use of two different algorithms (geNorm and Normfinder) to evaluate the stability of gene expression.

Results

Candidate housekeeping genes were sequenced: β-actin (Actin), elongation factor 1 alpha (EF1a), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), armadillo (Arm), ribosomal protein L32 (RpL32), succinate dehydrogenase (SDHa), annexin IX (AnnIX) and α-tubulin (Tub). The expression of these eight genes was analyzed in corpora allata (CA) and ovaries of adult female D. punctata. Both geNorm, as well as Normfinder characterized SDHa, EF1a and Arm as being the most stably expressed in the corpora allata. In the ovary, the geNorm calculation showed Tub, EF1a and RpL32 to be most stable, whereas Normfinder identified Tub, EF1a and Arm as the best. In ovary, the least stable gene was Actin, challenging its usefulness in normalization. As a proof of principle, the expression of follicle cell protein 3c and CYP15A1 was monitored during the first gonadotrophic cycle.

Conclusion

Arm and EF1a form the most stably expressed combination of two reference genes out of the eight candidates that were tested in the corpora allata. Our results show that the combined use of Tub, EF1a and RpL32 ensures an accurate normalization of gene expression levels in ovary of D. punctata. Our study has indicated that neither Actin nor AnnIX should be used for normalization of transcript levels when studying the first gonadotrophic cycle in CA or ovary of D. punctata. The results stress the necessity for validation of reference genes in q-RT-PCR studies in cockroaches.

Molecular signature of adipose tissue in patients with both non-alcoholic fatty liver disease (NAFLD) and polycystic ovarian syndrome (PCOS)

Background

Polycystic ovarian syndrome (PCOS) is one of the most common reproductive disorders with strong association with both insulin resistance and non-alcoholic fatty liver disease (NAFLD). To untangle the complex relationship between PCOS and NAFLD, we analyzed serum biomarkers of apoptosis, some adipokines and mRNA profiles in the visceral adipose tissue of obese patients with NAFLD who were also diagnosed with PCOS and compared to a group with NAFLD only.

Methods

We included patients with biopsy-proven NAFLD and PCOS (N = 12) and BMI-matched biopsy-proven NAFLD patients without PCOS (N = 12). Expression levels of individual mRNAs and soluble serum biomarkers were compared by non-parametric Mann–Whitney test. The analysis also included Spearman rank correlation tests and multiple regression analysis. For co-correlated genes, the factor analysis was performed.

Results

The total serum levels of apoptotic biomarker M30 were significantly elevated in PCOS patients with liver steatosis as compared to non-PCOS NAFLD controls (P < 0.02), pointing that androgen-dependent proapoptotic PCOS environment that may directly contribute to NAFLD progression in these patients. Similarly, hyperandrogenism may explain the observed PCOS-specific decrease (P < 0.04) in adipose LDLR mRNA expression that may be connected to the proneness of PCOS patients to NAFLD. The levels of mRNA encoding angiogenesis-associated GSK-3B interacting protein ninein were also significantly increased in the adipose tissue of NAFLD patients with PCOS (P < 0.007). Furthermore, the levels of resistin positively correlated with expression levels of LDLR and prothrombin time (PT).

Conclusion

An androgen-dependent proapoptotic PCOS environment may directly contribute to NAFLD progression in these patients. Hyperandrogenism may explain an observed decrease in adipose LDLR mRNA expression. An inflammation-associated increase in the release of resistin into circulation might contribute to the prothrombotic state observed under conditions associated with insulin resistance, including PCOS. The studies of larger cohorts of NAFLD with and without PCOS patients are needed to further assess these potential interactions.

Cross-platform microarray meta-analysis for the mouse jejunum selects novel reference genes with highly uniform levels of expression

Reference genes (RGs) with uniform expression are used for normalization of reverse transcription quantitative PCR (RT-qPCR) data. Their optimization for a specific biological context, e.g. a specific tissue, has been increasingly considered. In this article, we compare RGs identified by expression data meta-analysis restricted to the context tissue, the jejunum of Mus musculus domesticus, i) to traditional RGs, ii) to expressed interspersed repeated DNA elements, and iii) to RGs identified by meta-analysis of expression data from diverse tissues and conditions. To select the set of candidate RGs, we developed a novel protocol for the cross-platform meta-analysis of microarray data. The expression stability of twenty-four putative RGs was analysed by RT-qPCR in at least 14 jejunum samples of the mouse strains C57Bl/6N, CD1, and OF1. Across strains, the levels of expression of the novel RGs Plekha7, Zfx, and Ube2v1 as well as of Oaz1 varied less than two-fold irrespective of genotype, sex or their combination. The gene set consisting of Plekha7 and Oaz1 showed superior expression stability analysed with the tool RefFinder. The novel RGs are functionally diverse. This facilitates expression studies over a wide range of conditions. The highly uniform expression of the optimized RGs in the jejunum points towards their involvement in tightly regulated pathways in this tissue. We also applied our novel protocol of cross-microarray platform meta-analysis to the identification of RGs in the duodenum, the ileum and the entire small intestine. The selection of RGs with improved expression stability in a specific biological context can reduce the number of RGs for the normalization step of RT-qPCR expression analysis, thus reducing the number of samples and experimental costs.

Expression levels of obesity-related genes are associated with weight change in kidney transplant recipients

Background

The aim of this study was to investigate the association of gene expression profiles in subcutaneous adipose tissue with weight change in kidney transplant recipients and to gain insights into the underlying mechanisms of weight gain.

Methodology/Principal Findings

A secondary data analysis was done on a subgroup (n = 26) of existing clinical and gene expression data from a larger prospective longitudinal study examining factors contributing to weight gain in transplant recipients. Measurements taken included adipose tissue gene expression profiles at time of transplant, baseline and six-month weight, and demographic data. Using multivariate linear regression analysis controlled for race and gender, expression levels of 1553 genes were significantly (p<0.05) associated with weight change. Functional analysis using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes classifications identified metabolic pathways that were enriched in this dataset. Furthermore, GeneIndexer literature mining analysis identified a subset of genes that are highly associated with obesity in the literature and Ingenuity pathway analysis revealed several significant gene networks associated with metabolism and endocrine function. Polymorphisms in several of these genes have previously been linked to obesity.

Conclusions/Significance

We have successfully identified a set of molecular pathways that taken together may provide insights into the mechanisms of weight gain in kidney transplant recipients. Future work will be done to determine how these pathways may contribute to weight gain.

Nicotinamide nucleotide transhydrogenase mRNA expression is related to human obesity

OBJECTIVE

A spontaneous deletion in the nicotinamide nucleotide transhydrogenase (Nnt) gene eliminating exons 7-11 in C57BL/6J (B6J) mice is associated with reduced glucose-stimulated insulin secretion in vitro, impaired glucose tolerance, higher epigonadal fat mass, and altered susceptibility to diet induced obesity of male B6J mice was proposed. A potential implication for NNT in human adipose tissue distribution has not been investigated so far.

DESIGN AND METHODS

Therefore, NNT mRNA expression in paired human samples of visceral (vis) and subcutaneous (sc) adipose tissue from 221 subjects with a wide range of body mass index (BMI), insulin sensitivity, and glucose tolerance was analyzed.

RESULTS

NNT mRNA expression is significantly higher in visceral fat of obese patients and correlates with body weight, BMI, % body fat, visceral and sc fat area, waist and hip circumference, and fasting plasma insulin (FPI). Multivariate linear regression analysis revealed visceral NNT expression as age and gender independent predictor of BMI, waist circumference, visceral fat area, and % body fat, but not FPI and 2 h OGTT glucose.

CONCLUSION

In conclusion, a functional relevance of NNT in the development of human obesity and visceral fat distribution was suggested here.

Validation of internal reference genes for real-time quantitative polymerase chain reaction studies in the tick, Ixodes scapularis (Acari: Ixodidae)

Obtaining reliable gene expression data using real-time quantitative polymerase chain reaction (qPCR) is highly dependent on the choice of normalization method. We tested the expression stability of multiple candidate genes in the salivary glands (SG) and synganglia (SYN) of female Ixodes scapularis (Say) ticks in multiple blood-feeding phases. We found that the amount of total RNA in both the SG and SYN increases dramatically during tick feeding, with 34x and 5.8x increases from 62 and 7.1 ng of unfed tick, respectively. We tested candidate genes that were predicted from I. scapularis genome data to encode glyceraldehyde 3-phosphate dehydrogenase (gapdh), ribosomal protein L13A (l13a), TATA box-binding protein (tbp), ribosomal protein S4 (rps4), glucose 6-phosphate dehydrogenase (gpdh), and beta-glucuronidase (gusb). The geNorm and NormFinder algorithms were used to analyze data from different feeding phases (i.e., daily samples from unfed to fully engorged females over a 7-d period in three replicate experiments). We found that the rps4 and l13a genes showed highly stable expression patterns over the feeding duration in both the SG and SYN. Furthermore, the highly expressed rps4 gene makes it useful as a normalization factor when we perform studies using minute amounts of dissected tissue for qPCR. We conclude that rps4 and l13a, whether individually or as a pair, serve as suitable internal reference genes for qRT-PCR studies in the SG and SYN of I. scapularis.

Histological evidence: housekeeping genes beta-actin and GAPDH are of limited value for normalization of gene expression

Housekeeping genes are widely used as internal controls for gene expression normalization for western blotting, northern blotting, RT-PCR, etc. They are generally thought to be expressed in all cells of the organism at similar levels because it is assumed that these genes are required for the maintenance of basic cellular function as constitutive genes. However, real-time RT-PCR experiments revealed that their expression may vary depending on the developmental stage, type of tissue examined, experimental condition, and so on. To date, no histological data on their expression are available for embryonic development. In the present study, we compared the histological expression profile of two commonly used housekeeping genes, GAPDH and beta-actin, in the developing chicken embryo by using section and whole mount in situ hybridization supplemented by RT-PCR. Our results show that neither GAPDH mRNA nor beta-actin mRNA is expressed in all cell types or tissues at high levels. Strikingly, expression levels are very low in some organs. Moreover, the two genes show partially complementary expression patterns in the liver, the vascular system and the digestive tract. For example, GAPDH is more strongly expressed in the liver than beta-actin, but at lower levels in the arteries. Vice versa, beta-actin is more strongly expressed in the gizzard than GAPDH, but it is almost absent from cardiac muscle cells. Researchers should consider these histological results when using GAPGD and beta-actin for gene expression normalization in their experiments.

18S rRNA is a reliable normalisation gene for real time PCR based on influenza virus infected cells

BACKGROUND

One requisite of quantitative reverse transcription PCR (qRT-PCR) is to normalise the data with an internal reference gene that is invariant regardless of treatment, such as virus infection. Several studies have found variability in the expression of commonly used housekeeping genes, such as beta-actin (ACTB) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH), under different experimental settings. However, ACTB and GAPDH remain widely used in the studies of host gene response to virus infections, including influenza viruses. To date no detailed study has been described that compares the suitability of commonly used housekeeping genes in influenza virus infections. The present study evaluated several commonly used housekeeping genes [ACTB, GAPDH, 18S ribosomal RNA (18S rRNA), ATP synthase, H+ transporting, mitochondrial F1 complex, beta polypeptide (ATP5B) and ATP synthase, H+ transporting, mitochondrial Fo complex, subunit C1 (subunit 9) (ATP5G1)] to identify the most stably expressed gene in human, pig, chicken and duck cells infected with a range of influenza A virus subtypes.

RESULTS

The relative expression stability of commonly used housekeeping genes were determined in primary human bronchial epithelial cells (HBECs), pig tracheal epithelial cells (PTECs), and chicken and duck primary lung-derived cells infected with five influenza A virus subtypes. Analysis of qRT-PCR data from virus and mock infected cells using NormFinder and BestKeeper software programmes found that 18S rRNA was the most stable gene in HBECs, PTECs and avian lung cells.

CONCLUSIONS

Based on the presented data from cell culture models (HBECs, PTECs, chicken and duck lung cells) infected with a range of influenza viruses, we found that 18S rRNA is the most stable reference gene for normalising qRT-PCR data. Expression levels of the other housekeeping genes evaluated in this study (including ACTB and GPADH) were highly affected by influenza virus infection and hence are not reliable as reference genes for RNA normalisation.

Impact of reference gene selection for type 2 cannabinoid receptor gene expression studies in human spermatozoa

Quantitative real-time polymerase chain reaction (qRT-PCR) has been employed to study the gene expression profiles in human spermatozoa, but accurate analysis is dependent upon normalisation of data against an endogenous control. β-Actin (ACTB) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) are the most commonly used reference genes for normalisation of gene expression in human spermatozoa, but the expression of these genes in many tissues has considerable variation under different physiological, pathological and experimental conditions which limits their effectiveness in normalisation. The expression stability of a panel of 12 reference genes was studied in normal and pathological human spermatozoa using geNorm and NormFinder software. Although there were some discrepancies in the ranking of reference gene stability, each software program ranked B2 M, ACTB, CYC1 and 18S RNA within the top 5 and recommended the combined use of at least two reference genes. Normalisation of qRT-PCR data for the cannabinoid receptor type 2 in spermatozoa using the different housekeeping genes demonstrated how, without validation, conflicting results are obtained. We recommend that the arbitrary use of reference genes should be avoided and the validation of reference gene stability should be undertaken prior to every study. For normalisation of CB2 expression, we would recommend using the geometric mean of B2 M and ACTB.

Validation of reference genes for expression analysis in the salivary gland and the intestine of Rhodnius prolixus (Hemiptera, Reduviidae) under different experimental conditions by quantitative real-time PCR

BACKGROUND

Rhodnius prolixus is a blood-feeding insect that can transmit Trypanosoma cruzi and Trypanosoma rangeli to vertebrate hosts. Recently, genomic resources for invertebrate vectors of human pathogens have increased significantly, and R. prolixus has been one of the main species studied among the triatomines. However, the paucity of information on many of the fundamental molecular aspects of this species limits the use of the available genomic information. The present study aimed to facilitate gene expression studies by identifying the most suitable reference genes for the normalization of mRNA expression data from qPCR.

RESULTS

The expression stability of five candidate reference genes (18S rRNA, GAPDH, β-actin, α-tubulin and ribosomal protein L26) was evaluated by qPCR in two tissues (salivary gland and intestine) and under different physiological conditions: before and after blood feeding and after infection with T. cruzi or T. rangeli. The results were analyzed with three software programs: geNorm, NormFinder and BestKeeper. All of the evaluated candidate genes proved to be acceptable as reference genes, but some were found to be more appropriate depending on the experimental conditions. 18S, GAPDH and α-tubulin showed acceptable stability for studies in all of the tissues and experimental conditions evaluated. β-actin, one of the most widely used reference genes, was confirmed to be one of the most suitable reference genes in studies with salivary glands, but it had the lowest expression stability in the intestine after insect blood feeding. L26 was identified as the poorest reference gene in the studies performed.

CONCLUSIONS

The expression stability of the genes varies in different tissue samples and under different experimental conditions. The results provided by three statistical packages emphasize the suitability of all five of the tested reference genes in both the crop and the salivary glands with a few exceptions. The results emphasise the importance of validating reference genes for qRT-PCR analysis in R. prolixus studies.

Variation in stability of endogenous reference genes in fallopian tubes and endometrium from healthy and ectopic pregnant women

RT-qPCR is commonly employed in gene expression studies in ectopic pregnancy. Most use RN18S1, β-actin or GAPDH as internal controls without validation of their suitability as reference genes. A systematic study of the suitability of endogenous reference genes for gene expression studies in ectopic pregnancy is lacking. The aims of this study were therefore to evaluate the stability of 12 reference genes and suggest those that are stable for use as internal control genes in fallopian tubes and endometrium from ectopic pregnancy and healthy non-pregnant controls. Analysis of the results showed that the genes consistently ranked in the top six by geNorm and NormFinder algorithms, were UBC, GAPDH, CYC1 and EIF4A2 (fallopian tubes) and UBC and ATP5B (endometrium). mRNA expression of NAPE-PLD as a test gene of interest varied between the groups depending on which of the 12 reference genes was used as internal controls. This study demonstrates that arbitrary selection of reference genes for normalisation in RT-qPCR studies in ectopic pregnancy without validation, risk producing inaccurate data and should therefore be discouraged.

Evaluation of internal reference genes for quantitative expression analysis by real-time PCR in ovine whole blood

The use of reference genes is commonly accepted as the most reliable approach to normalize qRT-PCR and to reduce possible errors in the quantification of gene expression. The most suitable reference genes in sheep have been identified for a restricted range of tissues, but no specific data on whole blood are available. The aim of this study was to identify a set of reference genes for normalizing qRT-PCR from ovine whole blood. We designed 11 PCR assays for commonly employed reference genes belonging to various functional classes and then determined their expression stability in whole blood samples from control and disease-stressed sheep. SDHA and YWHAZ were considered the most suitable internal controls as they were stably expressed regardless of disease status according to both geNorm and NormFinder software; furthermore, geNorm indicated SDHA/HPRT, YWHAZ/GAPDH and SDHA/YWHAZ as the best reference gene combinations in control, disease-stressed and combined sheep groups, respectively. Our study provides a validated panel of optimal control genes which may be useful for the identification of genes differentially expressed by qRT-PCR in a readily accessible tissue, with potential for discovering new physiological and disease markers and as a tool to improve production traits (e.g., by identifying expression Quantitative Trait Loci). An additional outcome of the study is a set of intron-spanning primer sequences suitable for gene expression experiments employing SYBR Green chemistry on other ovine tissues and cells.

Selection of reference genes for qRT-PCR in high fat diet-induced hepatic steatosis mice model

With the epidemic proportions of obesity worldwide and the concurrent prevalence of hepatic steatosis, there is an urgent need for better understanding the intrinsic mechanism of hepatic steatosis, especially the changes of gene expression underlying the development of hepatic steatosis and its associated abnormal liver function. Quantitative real-time PCR (qRT-PCR) is a sensitive and highly reproducible technique of gene expression analysis. However, for accurate and reliable gene expression results, it is vital to have an internal control gene expressed at constant levels under all the experimental conditions being analyzed for. In this study, the authors validated candidate reference genes suitable for qRT-PCR profiling experiments using livers from control mice and high fat diet-induced obese mice. Cross-validation of expression stability of ten selected reference genes using three popular algorithms, GeNorm, NormFinder, and BestKeeper found HPRT1 and GAPDH as most stable reference genes. Thus, HPRT1 and GAPDH are recommended as stable reference genes most suitable for gene expression studies in the development of hepatic steatosis.

Housekeeping gene variability in the liver of alcoholic patients

BACKGROUND

Quantification of gene expression using reverse transcription-quantitative polymerase chain reaction (RT-qPCR) requires normalization to an endogenous reference gene termed housekeeping gene (HKG). Many of the commonly used HKGs are regulated and vary under experimental conditions and disease stages. Alcoholic liver disease (ALD) is associated with several different liver histological lesions that may modulate HKG expression. We investigated the variability of commonly used HGKs (18S, β-actin, glyceraldehyde-3-phosphate [GAPDH], and arginine/serine-rich splicing factor [SFRS4]) in the liver of patients with ALD.

METHODS

Fifty consecutive patients at different stages of ALD underwent liver biopsy. The stability of HKG was assessed according to liver histological lesions.

RESULTS

β-actin had the highest coefficient of dispersion (COD) (23.9). β-actin tended to decrease with steatosis and to increase with alcoholic hepatitis; β-actin also increased in patients with both alcoholic hepatitis and cirrhosis. GAPDH and SFRS4 COD were 2.8 and 2.1, respectively. GAPDH was decreased with steatosis and increased with alcoholic hepatitis and fibrosis. 18S had the lowest COD (1.4). Both 18S and SFRS4 levels were not significantly modified with respect to all alcohol-induced liver histological lesions.

CONCLUSIONS

In patients with ALD, the most constantly expressed HKGs are 18S and SFRS4. These genes are appropriate reference genes for normalization of RT-qPCR in the liver of patients with ALD. The use of other HKGs such as β-actin or GAPDH would lead to misinterpretation of the results.

Use of Maximum Likelihood-Mixed Models to select stable reference genes: a case of heat stress response in sheep

Background

Reference genes with stable expression are required to normalize expression differences of target genes in qPCR experiments. Several procedures and companion software have been proposed to find the most stable genes. Model based procedures are attractive because they provide a solid statistical framework. NormFinder, a widely used software, uses a model based method. The pairwise comparison procedure implemented in GeNorm is a simpler procedure but one of the most extensively used. In the present work a statistical approach based in Maximum Likelihood estimation under mixed models was tested and compared with NormFinder and geNorm softwares. Sixteen candidate genes were tested in whole blood samples from control and heat stressed sheep.

Results

A model including gene and treatment as fixed effects, sample (animal), gene by treatment, gene by sample and treatment by sample interactions as random effects with heteroskedastic residual variance in gene by treatment levels was selected using goodness of fit and predictive ability criteria among a variety of models. Mean Square Error obtained under the selected model was used as indicator of gene expression stability. Genes top and bottom ranked by the three approaches were similar; however, notable differences for the best pair of genes selected for each method and the remaining genes of the rankings were shown. Differences among the expression values of normalized targets for each statistical approach were also found.

Conclusions

Optimal statistical properties of Maximum Likelihood estimation joined to mixed model flexibility allow for more accurate estimation of expression stability of genes under many different situations. Accurate selection of reference genes has a direct impact over the normalized expression values of a given target gene. This may be critical when the aim of the study is to compare expression rate differences among samples under different environmental conditions, tissues, cell types or genotypes. To select reference genes not only statistical but also functional and biological criteria should be considered. Under the method here proposed SDHA/MDH1 have arisen as the best set of reference genes to be used in qPCR assays to study heat shock in ovine blood samples.

Effect of three common SNPs in 5'-flanking region of LEP and ADIPOQ genes on their expression in Polish obese children and adolescents

Genes encoding adipokines are considered as candidates for human obesity. In this study we analyzed the expression of leptin (LEP) and adiponectin (ADIPOQ) genes in relation to common 5′-flanking or 5′UTR variants: -2548G>A (LEP), 19A>G (LEP) and -11377C>G (ADIPOQ) in Polish obese children and adolescents. Relative transcription levels in the subcutaneous adipose tissue (real time RT–PCR) and serum protein concentrations (RIA) were measured in 48 obese subjects with known genotypes at three polymorphic sites and in five non-obese controls. None of the studied polymorphisms altered significantly the expression. Significantly elevated relative transcription levels of the LEP gene (P < 0.05) and serum leptin concentrations (P < 0.01) were recorded in obese patients, when compared with the non-obese controls, but such differences were not found for the ADIPOQ gene. Interestingly, the leptin to adiponectin protein concentration ratio (L/A) was approximately sevenfold higher in obese children and adolescents when compared with the non-obese controls (P < 0.001). Taking into consideration the observed relationship between the genotypes and the gene expression level we suggest that these SNPs are not conclusive markers for predisposition to obesity in Polish children and adolescents. On the other hand, we confirmed that the leptin to adiponectin gene expression ratio (L/A) is an informative index characterizing obesity.

Identification of internal control genes for quantitative expression analysis by real-time PCR in bovine peripheral lymphocytes

Gene expression studies in blood cells, particularly lymphocytes, are useful for monitoring potential exposure to toxicants or environmental pollutants in humans and livestock species. Quantitative PCR is the method of choice for obtaining accurate quantification of mRNA transcripts although variations in the amount of starting material, enzymatic efficiency, and the presence of inhibitors can lead to evaluation errors. As a result, normalization of data is of crucial importance. The most common approach is the use of endogenous reference genes as an internal control, whose expression should ideally not vary among individuals and under different experimental conditions. The accurate selection of reference genes is therefore an important step in interpreting quantitative PCR studies. Since no systematic investigation in bovine lymphocytes has been performed, the aim of the present study was to assess the expression stability of seven candidate reference genes in circulating lymphocytes collected from 15 dairy cows. Following the characterization by flow cytometric analysis of the cell populations obtained from blood through a density gradient procedure, three popular softwares were used to evaluate the gene expression data. The results showed that two genes are sufficient for normalization of quantitative PCR studies in cattle lymphocytes and that YWAHZ, S24 and PPIA are the most stable genes.

Impact of reference gene selection for target gene normalization on experimental outcome using real-time qRT-PCR in adipocytes

BACKGROUND

With the current rise in obesity-related morbidities, real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR) has become a widely used method for assessment of genes expressed and regulated by adipocytes. In order to measure accurate changes in relative gene expression and monitor intersample variability, normalization to endogenous control genes that do not change in relative expression is commonly used with qRT-PCR determinations. However, historical evidence has clearly demonstrated that the expression profiles of traditional control genes (e.g., β-actin, GAPDH, α-tubulin) are differentially regulated across multiple tissue types and experimental conditions.

METHODOLOGY/PRINCIPAL FINDINGS

Therefore, we validated six commonly used endogenous control genes under diverse experimental conditions of inflammatory stress, oxidative stress, synchronous cell cycle progression and cellular differentiation in 3T3-L1 adipocytes using TaqMan qRT-PCR. Under each study condition, we further evaluated the impact of reference gene selection on experimental outcome using examples of target genes relevant to adipocyte function and differentiation. We demonstrate that multiple reference genes are regulated in a condition-specific manner that is not suitable for use in target gene normalization.

CONCLUSION/SIGNIFICANCE

Data are presented demonstrating that inappropriate reference gene selection can have profound influence on study conclusions ranging from divergent statistical outcome to inaccurate data interpretation of significant magnitude. This study validated the use of endogenous controls in 3T3-L1 adipocytes and highlights the impact of inappropriate reference gene selection on data interpretation and study conclusions.