Housekeeping gene expression is affected by antidepressant treatment in a mouse fibroblast cell line

Selection and validation of reference genes for the normalization of quantitative real-time PCR in different muscle tissues of rabbits

Background

In molecular biology studies, the selection of optimal reference genes is of vital importance for accurately quantifying gene expression. The purpose of the present study was to screen the most stable reference genes in different muscle tissues of New Zealand white rabbits and Yufeng yellow rabbits.

Methods and results

Results indicated that the most stable reference genes in the muscle tissues of New Zealand white rabbits were HPRT1, ACTB and PPIC, while HPRT1, PPIC, and RPL13A were the most stable reference genes in muscle tissues of Yufeng yellow rabbits. However, in the longissimus dorsi muscle and the abdominal wall muscle of both varieties, the most stable reference genes were HPRT1, RPL13A, and SDHA. In the quadriceps femoris muscle, the most stable reference genes were ACTB, HPRT1, and SDHA. Furthermore, the relative abundance of MYOG, MYH3 and MSTN was used to confirm the suitability and reliability of the selected most stable reference genes and the most unstable reference gene. Results revealed the same expression patterns of these myogenic genes when normalized according to the most stable genes, while normalization against the unstable reference gene altered the observed expression patterns.

Conclusions

Taken together, our results demonstrated that the most stable reference genes varied among different muscle tissues and different breeds of rabbits. However, HPRT1, PPIC and SDHA presented high stability among all examined reference genes; thus, the combined analysis of HPRT1/ PPIC/ SDHA gene provides the best reference for RT-qPCR in muscle tissues of New Zealand white rabbits and Yufeng yellow rabbits, while HPRT1 is a better choice than other reference genes when using a single reference gene to assess target gene expression. Our results provide basic data for better measuring target gene expression profiles in muscle tissues of rabbits.

Validation of reference genes for use in untreated bovine fibroblasts

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

Identification of reference genes for quantitative PCR during C3H10T1/2 chondrogenic differentiation

C3H10T1/2, a mouse mesenchymal stem cell line, is a well-known in vitro model of chondrogenesis that can be easily employed to recapitulate some of the mechanisms intervening in this process. Moreover, these cells can be used to validate the effect of candidate molecules identified by high throughput screening approaches applied to the development of targeted therapy for human disorders in which chondrogenic differentiation may be involved, as in conditions characterized by heterotopic endochondral bone formation. Chondrogenic differentiation of C3H10T1/2 cells can be monitored by applying quantitative polymerase chain reaction (qPCR), one of the most sensitive methods that allows detection of small dynamic changes in gene expression between samples obtained under different experimental conditions. In this work, we have used qPCR to monitor the expression of specific markers during chondrogenic differentiation of C3H10T1/2 cells in micromass cultures. Then we have applied the geNorm approach to identify the most stable reference genes suitable to get a robust normalization of the obtained expression data. Among 12 candidate reference genes (Ap3d1, Csnk2a2, Cdc40, Fbxw2, Fbxo38, Htatsf1, Mon2, Pak1ip1, Zfp91, 18S, ActB, GAPDH) we identified Mon2 and Ap3d1 as the most stable ones during chondrogenesis. ActB, GAPDH and 18S, the most commonly used in the literature, resulted to have an expression level too high compared to the differentiation markers (Sox9, Collagen type 2a1, Collagen type 10a1 and Collagen type 1a1), therefore are actually less recommended for these experimental conditions. In conclusion, we identified nine reference genes that can be equally used to obtain a robust normalization of the gene expression variation during the C3H10T1/2 chondrogenic differentiation.

Differential Peripheral Proteomic Biosignature of Fluoxetine Response in a Mouse Model of Anxiety/Depression

The incorporation of peripheral biomarkers in the treatment of major depressive disorders (MDD) could improve the efficiency of treatments and increase remission rate. Peripheral blood mononuclear cells (PBMCs) represent an attractive biological substrate allowing the identification of a drug response signature. Using a proteomic approach with high-resolution mass spectrometry, the present study aimed to identify a biosignature of antidepressant response (fluoxetine, a Selective Serotonin Reuptake Inhibitor) in PBMCs in a mouse model of anxiety/depression. Following determination of an emotionality score, using complementary behavioral analysis of anxiety/depression across three different tests (Elevated Plus Maze, Novelty Suppressed Feeding, Splash Test), we showed that a 4-week corticosterone treatment (35 μg/ml, CORT model) in C57BL/6NTac male mice induced an anxiety/depressive-like behavior. Then, chronic fluoxetine treatment (18 mg/kg/day for 28 days in the drinking water) reduced corticosterone-induced increase in emotional behavior. However, among 46 fluoxetine-treated mice, only 30 of them presented a 50% decrease in emotionality score, defining fluoxetine responders (CORT/Flx-R). To determine a peripheral biological signature of fluoxetine response, proteomic analysis was performed from PBMCs isolated from the “most” affected corticosterone/vehicle (CORT/V), corticosterone/fluoxetine responders and non-responders (CORT/Flx-NR) animals. In comparison to CORT/V, a total of 263 proteins were differently expressed after fluoxetine exposure. Expression profile of these proteins showed a strong similarity between CORT/Flx-R and CORT/Flx-NR (R = 0.827, p < 1e^-7). Direct comparison of CORT/Flx-R and CORT/Flx-NR groups revealed 100 differently expressed proteins, representing a combination of markers associated either with the maintenance of animals in a refractory state, or associated with behavioral improvement. Finally, 19 proteins showed a differential direction of expression between CORT/Flx-R and CORT/Flx-NR that drove them away from the CORT-treated profile. Among them, eight upregulated proteins (RPN2, HSPA9, NPTN, AP2B1, UQCRC2, RACK-1, TOLLIP) and one downregulated protein, TLN2, were previously associated with MDD or antidepressant drug response in the literature. Future preclinical studies will be required to validate whether proteomic changes observed in PBMCs from CORT/Flx-R mice mirror biological changes in brain tissues.

Human dermal fibroblasts in psychiatry research

In order to decipher the disease etiology, progression and treatment of multifactorial human brain diseases we utilize a host of different experimental models. Recently, patient-derived human dermal fibroblast (HDF) cultures have re-emerged as promising in vitro functional system for examining various cellular, molecular, metabolic and (patho)physiological states and traits of psychiatric disorders. HDF studies serve as a powerful complement to postmortem and animal studies, and often appear to be informative about the altered homeostasis in neural tissue. Studies of HDFs from patients with schizophrenia (SZ), depression, bipolar disorder (BD), autism, attention deficit and hyperactivity disorder and other psychiatric disorders have significantly advanced our understanding of these devastating diseases. These reports unequivocally prove that signal transduction, redox homeostasis, circadian rhythms and gene*environment (G*E) interactions are all amenable for assessment by the HDF model. Furthermore, the reported findings suggest that this underutilized patient biomaterial, combined with modern molecular biology techniques, may have both diagnostic and prognostic value, including prediction of response to therapeutic agents.

Identification of stable reference genes for quantitative PCR in cells derived from chicken lymphoid organs

Quantitative polymerase chain reaction (qPCR) is a powerful technique for quantification of gene expression, especially genes involved in immune responses. Although qPCR is a very efficient and sensitive tool, variations in the enzymatic efficiency, quality of RNA and the presence of inhibitors can lead to errors. Therefore, qPCR needs to be normalised to obtain reliable results and allow comparison. The most common approach is to use reference genes as internal controls in qPCR analyses. In this study, expression of seven genes, including β-actin (ACTB), β-2-microglobulin (B2M), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), β-glucuronidase (GUSB), TATA box binding protein (TBP), α-tubulin (TUBAT) and 28S ribosomal RNA (r28S), was determined in cells isolated from chicken lymphoid tissues and stimulated with three different mitogens. The stability of the genes was measured using geNorm, NormFinder and BestKeeper software. The results from both geNorm and NormFinder were that the three most stably expressed genes in this panel were TBP, GAPDH and r28S. BestKeeper did not generate clear answers because of the highly heterogeneous sample set. Based on these data we will include TBP in future qPCR normalisation. The study shows the importance of appropriate reference gene normalisation in other tissues before qPCR analysis.

Exogenous reference gene normalization for real-time reverse transcription-polymerase chain reaction analysis under dynamic endogenous transcription

Quantitative real-time reverse transcription-polymerase chain reaction (qPCR) is widely used to investigate transcriptional changes following experimental manipulations to the nervous system. Despite the widespread utilization of qPCR, the interpretation of results is marred by the lack of a suitable reference gene due to the dynamic nature of endogenous transcription. To address this inherent deficiency, we investigated the use of an exogenous spike-in mRNA, luciferase, as an internal reference gene for the 2(-∆∆Ct) normalization method. To induce dynamic transcription, we systemically administered capsaicin, a neurotoxin selective for C-type sensory neurons expressing the TRPV-1 receptor, to adult male Sprague-Dawley rats. We later isolated nodose ganglia for qPCR analysis with the reference being either exogenous luciferase mRNA or the commonly used endogenous reference β-III tubulin. The exogenous luciferase mRNA reference clearly demonstrated the dynamic expression of the endogenous reference. Furthermore, variability of the endogenous reference would lead to misinterpretation of other genes of interest. In conclusion, traditional reference genes are often unstable under physiologically normal situations, and certainly unstable following the damage to the nervous system. The use of exogenous spike-in reference provides a consistent and easily implemented alternative for the analysis of qPCR data.

Molecular signatures of mood stabilisers highlight the role of the transcription factor REST/NRSF

BACKGROUND

The purpose of this study was to address the affects of mood modifying drugs on the transcriptome, in a tissue culture model, using qPCR arrays as a cost effective approach to identifying regulatory networks and pathways that might coordinate the cell response to a specific drug.

METHODS

We addressed the gene expression profile of 90 plus genes associated with human mood disorders using the StellARray™ qPCR gene expression system in the human derived SH-SY5Y neuroblastoma cell line.

RESULTS

Global Pattern Recognition (GPR) analysis identified a total of 9 genes (DRD3(⁎), FOS(†), JUN(⁎), GAD1(⁎†), NRG1(⁎), PAFAH1B3(⁎), PER3(⁎), RELN(⁎) and RGS4(⁎)) to be significantly regulated in response to cellular challenge with the mood stabilisers sodium valproate ((⁎)) and lithium ((†)). Modulation of FOS and JUN highlights the importance of the activator protein 1 (AP-1) transcription factor pathway in the cell response. Enrichment analysis of transcriptional networks relating to this gene set also identified the transcription factor neuron restrictive silencing factor (NRSF) and the oestrogen receptor as an important regulatory mechanism.

LIMITATIONS

Cell line models offer a window of what might happen in vivo but have the benefit of being human derived and homogenous with regard to cell type.

CONCLUSIONS

This data highlights transcription factor pathways, acting synergistically or separately, in the modulation of specific neuronal gene networks in response to mood stabilising drugs. This model can be utilised in the comparison of the action of multiple drug regimes or for initial screening purposes to inform optimal drug design.

Mood-stabilizers differentially affect housekeeping gene expression in human cells

Recent studies have revealed that antidepressants affect the expression of constitutively expressed "housekeeping genes" commonly used as normalizing reference genes in quantitative polymerase chain reaction (qPCR) experiments. There has yet to be an investigation however on the effects of mood-stabilizers on housekeeping gene stability. The current study utilized lymphoblastoid cell lines (LCLs) derived from patients with mood disorders to investigate the effects of a range of doses of lithium (0, 1, 2 and 5 mM) and sodium valproate (0, 0.06, 0.03 and 0.6 mM) on the stability of 12 housekeeping genes. RNA was extracted from LCLs and qPCR was used to generate cycle threshold (Ct ) values which were input into RefFinder analyses. The study revealed drug-specific effects on housekeeping gene stability. The most stable housekeeping genes in LCLs treated: acutely with sodium valproate were ACTB and RPL13A; acutely with lithium were GAPDH and ATP5B; chronically with lithium were ATP5B and CYC1. The stability of GAPDH and B2M were particularly affected by duration of lithium treatment. The study adds to a growing literature that the selection of appropriate housekeeping genes is important for the accurate normalization of target gene expression in experiments investigating the molecular effects of mood disorder pharmacotherapies.

Tumor necrosis factor and its targets in the inflammatory cytokine pathway are identified as putative transcriptomic biomarkers for escitalopram response

Converging evidence suggests that the activation of the inflammatory cytokine pathway is important in the pathophysiology of unipolar depression. Antidepressants have anti-inflammatory properties and evidence suggests that inter-individual variability in response to antidepressants may reflect genetic differences in the inflammatory cytokine pathway. In particular, protein levels of Tumor Necrosis Factor (TNF) and the SNPs rs1126757 in interleukin-11 (IL11), and rs7801617 in interleukin-6 (IL6), have previously been implicated in the clinical response to the selective serotonin reuptake inhibitor (SSRI) antidepressant escitalopram. This study investigated the transcription of TNF, IL11 and IL6 as well as genes in the wider inflammatory cytokine pathway both at baseline and after escitalopram treatment in depressed patients who were either clinical "responders" (n=25) or "non-responders" (n=21). Samples were obtained as a subset of the Genome-Based Therapeutic Drugs for Depression (GENDEP) project and response status is based on changes in the Montgomery-Asberg Depression Rating Scores over a 12 wk treatment period. Binary logistic regressions revealed significant expression differences at baseline between responders and non-responders in TNF, and after escitalopram treatment in TNF and IL11. Differences in IL11 after treatment were found to be driven by drug-induced allele-specific expression differences relating to rs1126757. Top hits in the wider inflammatory cytokine pathway at both baseline and after escitalopram treatment were found to be targets of TNF. The current study adds substantial support for the role of the inflammatory cytokine pathway in mediating response to the SSRI escitalopram, and is the first to identify TNF and its targets as putative transcriptomic predictors of clinical response.

Allele-specific expression of the serotonin transporter and its transcription factors following lamotrigine treatment in vitro

Lamotrigine, a mood stabilizer used clinically in the treatment of bipolar disorder, is thought to exert actions on the serotonin system. However lamotrigine's exact mechanism of action remains unclear. The current study investigated whether lamotrigine might exert its effects through altering the expression of the serotonin transporter (5-HTT) gene and its regulatory transcription factors Y box binding protein 1 (YB-1) and CCCTC-binding factor (CTCF). We further considered whether functional variable number tandem repeat (VNTR) polymorphisms in the promoter region of 5-HTT, (5-HTTLPR) and within intron 2 (Stin2) of the gene, moderated any putative gene expression changes. The study employed an in vitro design carried out in human lymphoblastoid cell lines (LCLs) to investigate the effects of lamotrigine treatment at 0.04, 0.2, and 0.4 mM doses for 24 hr on the mRNA expression of 5-HTT, YB-1, and CTCF. LCLs were selected based on combinations of haplotypes of the two VNTRs in the serotonin transporter gene; creating low-expressing and high-expressing LCL groups. Ubiquitin C (UBC) and topoisomerase I (TOP1) genes were found to be the most stably expressed housekeeping genes in drug-treated LCLs. Subsequently, quantitative PCR revealed that higher doses of lamotrigine significantly lowered 5-HTT expression and increased CTCF expression. Haplotype-specific differences in CTCF expression were found in response to lamotrigine, with strongest expression changes observed in the high-expressing LCLs. These data provide an allele-specific in vitro model for examining the molecular targets of lamotrigine, and support the important role of the serotonin transporter gene in its clinical mechanism of action.

PepGMV rep-protein expression in mammalian cells

The Geminiviruses genome is a small, single strand DNA that replicates in the plant cell nucleus. Analogous to animal DNA viruses, Geminiviruses depend on the host replication machinery to amplify their genomes and only supply the factors required to initiate their replication. Consequently, Geminiviruses remove the cell-cycle arrest and induce the host replication machinery using an endocycle process. They encode proteins, such as the conserved replication-associated proteins (Rep) that interact with retinoblastoma-like proteins in plants and alter the cell division cycle in yeasts. Therefore, the aim of this work is to analyze the impact of Pepper Golden Mosaic Virus (PepGMV) Rep protein in mammalian cells. Results indicate that the pTracer-SV40:Rep construction obtained in this work can be used to analyze the Rep protein effect in mammalian cells in order to compare the cell cycle regulation mechanisms in plants and animals.

Gene-expression analysis reveals that embryonic stem cells cultured under osteogenic conditions produce mineral non-specifically compared to marrow stromal cells or osteoblasts

Pluripotent cells, such as embryonic stem cells (ESCs), divide indefinitely and can differentiate to form mineralised nodules in response to osteogenic supplements. This suggests that they may be used as a cell source for bone replacement strategies. Here, we related the expression of osteogenic and chondrogenic genes in cultures of murine ESCs, marrow stromal cells (MSCs) and calvarial osteoblasts (OBs) cultured under osteogenic conditions to the biochemical composition and quantity of mineral formed. Mineralisation, measured by calcium sequestration, was >2-fold greater in ESC cultures than in either MSCs or OBs. Micro-Raman spectroscopy and spectral mapping revealed a lower mineral-to-matrix ratio and confirmed a more diffuse pattern of mineralisation in ESCs compared to MSCs and OBs. Baseline expression of chondrogenic and osteogenic genes was between 1 and 4 orders of magnitude greater in MSCs and OBs than in ESCs. Osteogenic culture of MSCs and OBs was accompanied by increases in osteogenic gene expression by factors of ~100 compared to only ~10 in ESCs. Consequentially, peak expression of osteogenic and chondrogenic genes was greater in MSCs and OBs than ESCs by factors of 100-1000, despite the fact that mineralisation was more extensive in ESCs than either MSCs or OBs. We also observed significant cell death in ESC nodules. We conclude that the mineralised material observed in cultures of murine ESCs during osteogenic differentiation may accumulate non-specifically, perhaps in necrotic cell layers, and that thorough characterisation of the tissue formed by ESCs must be achieved before these cells can be considered as a cell source for clinical applications.

Evaluation of reference genes for quantitative PCR analysis of mouse lymphocytes

Normalization to a reference gene is the method of choice for quantitative PCR analysis. The stability of reference genes is critical for accurate gene expression analysis, as significant variations in reference gene expression can alter experimental results and conclusions. In this study, we evaluated the expression stability of five commonly used reference genes found in mouse lymphocytes. Using NormFinder and BestKeeper algorithms, we consistently show that ubiquitin C (Ubc) is the optimal reference gene for normalizing qPCR data obtained from mouse lymphocytes, whereas beta-actin (Actb) is not a suitable reference gene due to its extensive variability in expression. Our findings emphasize the importance of validating reference genes for qPCR analyses. We provide a shortlist of reference genes to use for normalization and recommend freely available software programs as a rapid approach to validate potential reference genes.

Murine 3T3-L1 adipocyte cell differentiation model: validated reference genes for qPCR gene expression analysis

BACKGROUND

Analysis of gene expression at the mRNA level, using real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR), mandatorily requires reference genes (RGs) as internal controls. However, increasing evidences have shown that RG expression may vary considerably under experimental conditions. We sought for an appropriate panel of RGs to be used in the 3T3-L1 cell line model during their terminal differentiation into adipocytes. To this end, the expression levels of a panel of seven widely used RG mRNAs were measured by qRT-PCR. The 7 RGs evaluated were ß-actin (ACTB), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), hypoxanthine phosphoribosyl-transferase I (HPRT), ATP synthase H+ transporting mitochondrial F1 complex beta subunit (ATP-5b), tyrosine 3-monooxygenase/tryptophan 5- monooxygenase activation protein, zeta polypeptide (Ywhaz), Non-POU-domain containing octamer binding protein (NoNo), and large ribosomal protein L13a (RPL).

METHODOLOGY/PRINCIPAL FINDINGS

Using three Excel applications, GeNorm, NormFinder and BestKeeper, we observed that the number and the stability of potential RGs vary significantly during differentiation of 3T3-L1 cells into adipocytes. mRNA expression analyses using qRT-PCR revealed that during the entire differentiation program, only NoNo expression is relatively stable. Moreover, the RG sets that were acceptably stable were different depending on the phase of the overall differentiation process (i.e. mitotic clonal expansion versus the terminal differentiation phase). RPL, ACTB, and Ywhaz, are suitable for terminal differentiation, whereas ATP-5b and HPRT, are suitable during mitotic clonal expansion.

CONCLUSION

Our results demonstrate that special attention must be given to the choice of suitable RGs during the various well defined phases of adipogenesis to ensure accurate data analysis and that the use of several RGs is absolutely required. Consequently, our data show for the first time, that during mitotic clonal expansion, the most suitable RGs are ATP-5b, NoNo and HPRT, while during terminal differentiation the most suitable RGs are, NoNo, RPL, ACTB and Ywhaz.

The honey bee epigenomes: differential methylation of brain DNA in queens and workers

Using genome-wide methylation profiles in honey bee queen and worker brains to understand how contrasting organismal outputs are generated from the same genotype.

In honey bees (Apis mellifera) the behaviorally and reproductively distinct queen and worker female castes derive from the same genome as a result of differential intake of royal jelly and are implemented in concert with DNA methylation. To determine if these very different diet-controlled phenotypes correlate with unique brain methylomes, we conducted a study to determine the methyl cytosine (mC) distribution in the brains of queens and workers at single-base-pair resolution using shotgun bisulfite sequencing technology. The whole-genome sequencing was validated by deep 454 sequencing of selected amplicons representing eight methylated genes. We found that nearly all mCs are located in CpG dinucleotides in the exons of 5,854 genes showing greater sequence conservation than non-methylated genes. Over 550 genes show significant methylation differences between queens and workers, revealing the intricate dynamics of methylation patterns. The distinctiveness of the differentially methylated genes is underscored by their intermediate CpG densities relative to drastically CpG-depleted methylated genes and to CpG-richer non-methylated genes. We find a strong correlation between methylation patterns and splicing sites including those that have the potential to generate alternative exons. We validate our genome-wide analyses by a detailed examination of two transcript variants encoded by one of the differentially methylated genes. The link between methylation and splicing is further supported by the differential methylation of genes belonging to the histone gene family. We propose that modulation of alternative splicing is one mechanism by which DNA methylation could be linked to gene regulation in the honey bee. Our study describes a level of molecular diversity previously unknown in honey bees that might be important for generating phenotypic flexibility not only during development but also in the adult post-mitotic brain.

The queen honey bee and her worker sisters do not seem to have much in common. Workers are active and intelligent, skillfully navigating the outside world in search of food for the colony. They never reproduce; that task is left entirely to the much larger and longer-lived queen, who is permanently ensconced within the colony and uses a powerful chemical influence to exert control. Remarkably, these two female castes are generated from identical genomes. The key to each female's developmental destiny is her diet as a larva: future queens are raised on royal jelly. This specialized diet is thought to affect a particular chemical modification, methylation, of the bee's DNA, causing the same genome to be deployed differently. To document differences in this epigenomic setting and hypothesize about its effects on behavior, we performed high-resolution bisulphite sequencing of whole genomes from the brains of queen and worker honey bees. In contrast to the heavily methylated human genome, we found that only a small and specific fraction of the honey bee genome is methylated. Most methylation occurred within conserved genes that provide critical cellular functions. Over 550 genes showed significant methylation differences between the queen and the worker, which may contribute to the profound divergence in behavior. How DNA methylation works on these genes remains unclear, but it may change their accessibility to the cellular machinery that controls their expression. We found a tantalizing clue to a mechanism in the clustering of methylation within parts of genes where splicing occurs, suggesting that methylation could control which of several versions of a gene is expressed. Our study provides the first documentation of extensive molecular differences that may allow honey bees to generate different phenotypes from the same genome.