The Renilla luciferase gene as a reference gene for normalization of gene expression in transiently transfected cells

Background: The importance of appropriate normalization controls in quantitative real-time polymerase chain reaction (qPCR) experiments has become more apparent as the number of biological studies using this methodology has increased. In developing a system to study gene expression from transiently transfected plasmids, it became clear that normalization using chromosomally encoded genes is not ideal, at it does not take into account the transfection efficiency and the significantly lower expression levels of the plasmids. We have developed and validated a normalization method for qPCR using a co-transfected plasmid.Results: The best chromosomal gene for normalization in the presence of the transcriptional activators used in this study, cadmium, dexamethasone, forskolin and phorbol-12-myristate 13-acetate was first identified. qPCR data was analyzed using geNorm, Normfinder and BestKeeper. Each software application was found to rank the normalization controls differently with no clear correlation. Including a co-transfected plasmid encoding the Renilla luciferase gene (Rluc) in this analysis showed that its calculated stability was not as good as the optimised chromosomal genes, most likely as a result of the lower expression levels and transfection variability. Finally, we validated these analyses by testing two chromosomal genes (B2M and ActB) and a co-transfected gene (Rluc) under biological conditions. When analyzing co-transfected plasmids, Rluc normalization gave the smallest errors compared to the chromosomal reference genes.Conclusions: Our data demonstrates that transfected Rluc is the most appropriate normalization reference gene for transient transfection qPCR analysis; it significantly reduces the standard deviation within biological experiments as it takes into account the transfection efficiencies and has easily controllable expression levels. This improves reproducibility, data validity and most importantly, enables accurate interpretation of qPCR data. © 2010 Jiwaji et al; licensee BioMed Central Ltd.

Mechanical stimulation alters pleiotrophin and aggrecan expression by human intervertebral disc cells and influences their capacity to stimulate endothelial migration

Study Design. The influence of mechanical load on pleiotrophin (PTM) and aggrecan expression by intervertebral disc (IVD) cells, and the effects of disc cell conditioned medium on endothelial cell migration was investigated. Objective. To examine possible interactions of mechanical loads and known pro- and antiangiogenic factors, which may regulate disc angiogenesis during degeneration. Summary of Background Data. Pleiotrophin expression can be influenced by mechanical stimulation and has been associated with disc vascularization. Disc aggrecan inhibits endothelial cell migration, suggesting an antiangiogenic role. A possible interplay between these factors is unknown. Methods. The influence of the respective predominant load (cyclic strain for anulus fibrosus and hydrostatic pressure for nucleus pulposus cells) on PTN and aggrecan expression by IVD cells was determined by real-time RT-PCR and Western blotting (PTN only). The effects of IVD cell conditioned medium on endothelial cell migration were analyzed in a bioassay using human microvascular endothelial (HMEC-1) cells. Results. Application of both mechanical loads resulted in significant alterations of gene expression of PTN (+67%, P = 0.004 in anulus cells; +29%, P = 0.03 in nucleus cells) and aggrecan (+42%, P = 0.03 in anulus cells, -25%, P = 0.03 in nucleus cells). These effects depended on the cell type, the applied load, and timescale. Conditioned media of nucleus pulposus cells enhanced HMEC-1 migration, but this effect was diminished after 2.5 MPa hydrostatic pressure, when aggrecan expression was diminished, but not 0.25 MPa, when expression levels were unchanged. Conclusion. Mechanical loading influences PTN expression by human IVD cells. Conditioned media from nucleus pulposus cell cultures stimulated HMEC-1 endothelial cell migration. This study demonstrates that the influence of mechanical loads on vascularization of the human IVD is likely to be complex and does not correlate simply with altered expression of known pro- and antiangiogenic factors.

Regulation of beta-cell viability and gene expression by distinct agonist fragments of adiponectin

Obesity is an established risk factor for type 2 diabetes. Activation of the adiponectin receptors has a clear role in improving insulin resistance although conflicting evidence exists for its effects on pancreatic beta-cells. Previous reports have identified both adiponectin receptors (ADR-1 and ADR-2) in the beta-cell. Recent evidence has suggested that two distinct regions of the adiponectin molecule, the globular domain and a small N-terminal region, have agonist properties. This study investigates the effects of two agonist regions of adiponectin on insulin secretion, gene expression, cell viability and cell signalling in the rat beta-cell line BRIN-BD11, as well as investigating the expression levels of adiponectin receptors (ADRs) in these cells. Cells were treated with globular adiponectin and adiponectin (15-36) +/-leptin to investigate cell viability, expression of key beta-cell genes and ERK1/2 activation. Both globular adiponectin and adiponectin (15-36) caused significant ERK1/2 dependent increases in cell viability. Leptin co-incubation attenuated adiponectin (15-36) but not globular adiponectin induced cell viability. Globular adiponectin, but not adiponectin (15-36), caused a significant 450% increase in PDX-1 expression and a 45% decrease in LPL expression. ADR-1 was expressed at a higher level than ADR-2, and ADR mRNA levels were differentially regulated by non-esterified fatty acids and peroxisome-proliferator-activated receptor agonists. These data provide evidence of roles for two distinct adiponectin agonist domains in the beta-cell and confirm the potentially important role of adiponectin receptor agonism in maintaining beta-cell mass.

Decreased NADPH oxidase expression and antioxidant activity in cachectic skeletal muscle

Background - Cancer cachexia is the progressive loss of skeletal muscle protein that contributes significantly to cancer morbidity and mortality. Evidence of antioxidant attenuation and the presence of oxidised proteins in patients with cancer cachexia indicate a role for oxidative stress. The level of oxidative stress in tissues is determined by an imbalance between reactive oxygen species production and antioxidant activity. This study aimed to investigate the superoxide generating NADPH oxidase (NOX) enzyme and antioxidant enzyme systems in murine adenocarcinoma tumour-bearing cachectic mice. Methods - Superoxide levels, mRNA levels of NOX enzyme subunits and the antioxidant enzymes superoxide dismutase (SOD), glutathione peroxidise (GPx) and catalase was measured in the skeletal muscle of mice with cancer and cancer cachexia. Protein expression levels of NOX enzyme subunits and antioxidant enzyme activity was also measured in the same muscle samples. Results - Superoxide levels increased 1.4-fold in the muscle of mice with cancer cachexia, and this was associated with a decrease in mRNA of NOX enzyme subunits, NOX2, p40phox and p67phox along with the antioxidant enzymes SOD1, SOD2 and GPx. Cancer cachexia was also associated with a 1.3-fold decrease in SOD1 and 2.0-fold decrease in GPx enzyme activity. Conclusion - Despite increased superoxide levels in cachectic skeletal muscle, NOX enzyme subunits, NOX2, p40phox and p67phox, were downregulated along with the expression and activity of the antioxidant enzymes. Therefore, the increased superoxide levels in cachectic skeletal muscle may be attributed to the reduction in the activity of endogenous antioxidant enzymes.

Microarray analysis of expression of cell death-associated genes in spinal cord cells with cyclic tensile strain

Previous studies have described alterations in gene expression following spinal cord injury, but this response to mechanical stimuli is difficult to investigate in vivo. Therefore, we have investigated the effect of cyclic tensile strain on cultured spinal cord cells from E15 Sprague-Dawley rats. Microarray analysis of gene expression and categorization of identified genes were performed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) systems. The application of cyclic tensile strain reduced the viability of cultured spinal cord cells significantly in a dose- and time-dependent manner. GO analysis identified candidate genes related to apoptosis (44) and to response to stimulus (17). KEGG analysis identified changes in the expression levels of 12 genes of the mitogen-activated protein kinase (MAPK) signaling pathway, which were confirmed to be upregulated and validated by RT-PCR analysis. Spinal cord cells undergo cell death in response to cyclic tensile strain, which were dose- and time-dependent, with upregulation of various genes, in particular of the MAPK pathway.

Development of Gene Expression Markers of Acute Heat-Light Stress in Reef-Building Corals of the Genus Porites

Coral reefs are declining worldwide due to increased incidence of climate-induced coral bleaching, which will have widespread biodiversity and economic impacts. A simple method to measure the sub-bleaching level of heat-light stress experienced by corals would greatly inform reef management practices by making it possible to assess the distribution of bleaching risks among individual reef sites. Gene expression analysis based on quantitative PCR (qPCR) can be used as a diagnostic tool to determine coral condition in situ. We evaluated the expression of 13 candidate genes during heat-light stress in a common Caribbean coral Porites astreoides, and observed strong and consistent changes in gene expression in two independent experiments. Furthermore, we found that the apparent return to baseline expression levels during a recovery phase was rapid, despite visible signs of colony bleaching. We show that the response to acute heat-light stress in P. astreoides can be monitored by measuring the difference in expression of only two genes: Hsp16 and actin. We demonstrate that this assay discriminates between corals sampled from two field sites experiencing different temperatures. We also show that the assay is applicable to an Indo-Pacific congener, P. lobata, and therefore could potentially be used to diagnose acute heat-light stress on coral reefs worldwide.

XBP1s levels are implicated in the biology and outcome of myeloma mediating different clinical outcomes to thalidomide-based treatments.

Immunoglobulin production by myeloma plasma cells depends on the unfolded protein response for protein production and folding. Recent studies have highlighted the importance of IRE1alpha and X box binding protein 1 (XBP1), key members of this pathway, in normal B-plasma cell development. We have determined the gene expression levels of IRE1alpha, XBP1, XBP1UNSPLICED (XBP1u), and XBP1SPLICED (XBP1s) in a series of patients with myeloma and correlated findings with clinical outcome. We show that IRE1alpha and XBP1 are highly expressed and that patients with low XBP1s/u ratios have a significantly better overall survival. XBP1s is an independent prognostic marker and can be used with beta2 microglobulin and t(4;14) to identify a group of patients with a poor outcome. Furthermore, we show the beneficial therapeutic effects of thalidomide in patients with low XBP1s/u ratios. This study highlights the importance of XBP1 in myeloma and its significance as an independent prognostic marker and as a predictor of thalidomide response.

Genetic variation of temperature-regulated curd induction in cauliflower: elucidation of floral transition by genome-wide association mapping and gene expression analysis

Cauliflower (Brassica oleracea var. botrytis) is a vernalization-responsive crop. High ambient temperatures delay harvest time. The elucidation of the genetic regulation of floral transition is highly interesting for a precise harvest scheduling and to ensure stable market supply. This study aims at genetic dissection of temperature-dependent curd induction in cauliflower by genome-wide association studies and gene expression analysis. To assess temperature dependent curd induction, two greenhouse trials under distinct temperature regimes were conducted on a diversity panel consisting of 111 cauliflower commercial parent lines, genotyped with 14,385 SNPs. Broad phenotypic variation and high heritability (0.93) were observed for temperature-related curd induction within the cauliflower population. GWA mapping identified a total of 18 QTL localized on chromosomes O1, O2, O3, O4, O6, O8, and O9 for curding time under two distinct temperature regimes. Among those, several QTL are localized within regions of promising candidate flowering genes. Inferring population structure and genetic relatedness among the diversity set assigned three main genetic clusters. Linkage disequilibrium (LD) patterns estimated global LD extent of r(2) = 0.06 and a maximum physical distance of 400 kb for genetic linkage. Transcriptional profiling of flowering genes FLOWERING LOCUS C (BoFLC) and VERNALIZATION 2 (BoVRN2) was performed, showing increased expression levels of BoVRN2 in genotypes with faster curding. However, functional relevance of BoVRN2 and BoFLC2 could not consistently be supported, which probably suggests to act facultative and/or might evidence for BoVRN2/BoFLC-independent mechanisms in temperature regulated floral transition in cauliflower. Genetic insights in temperature-regulated curd induction can underpin genetically informed phenology models and benefit molecular breeding strategies toward the development of thermo-tolerant cultivars.

Expression of tumor-related Rac1b antagonizes B-Raf-induced senescence in colorectal cells

Mutations in the BRAF oncogene have been identified as a tumor-initiating genetic event in mainly melanoma, thyroid and colon cancer, resulting in an initial proliferative stimulus that is followed by a growth arrest period known as oncogene-induced senescence (OIS). It remains unknown what triggers subsequent escape from OIS to allow further tumor progression. A previous analysis revealed that overexpression of splice variant Rac1b occurs in around 80% of colorectal tumors carrying a mutation in BRAF. Using both BRaf-V600E-directed RNAi and overexpression we demonstrate that this mutation does not directly lead to Rac1b overexpression, indicating the latter as an independent event during tumor progression. Nonetheless, we observed that expression of oncogenic BRaf-V600E in non-transformed colonocytes (NCM460 cell line) increased both the transcript and protein levels of p14ARF, p15INK4b and p21CIP1 and led to increased expression of β-galactosidase, all indicators of OIS induction. Interestingly, whereas the protein levels of these markers were reduced upon Rac1b overexpression, the levels of their respective transcripts remained unchanged. Importantly, the co-expression of Rac1b with B-Raf-V600E reverted the OIS phenotype, reducing the expression levels of the cell-cycle inhibitors and β-galactosidase to those of control cells. These data identify increased Rac1b expression as one potential mechanism by which colorectal tumor cells can escape from B-Raf-induced OIS.

Transcriptional dynamics of the developing sweet cherry (Prunus avium L.) fruit: sequencing, annotation and expression profiling of exocarp-associated genes

The exocarp, or skin, of fleshy fruit is a specialized tissue that protects the fruit, attracts seed dispersing fruit eaters, and has large economical relevance for fruit quality. Development of the exocarp involves regulated activities of many genes. This research analyzed global gene expression in the exocarp of developing sweet cherry (Prunus avium L., 'Regina'), a fruit crop species with little public genomic resources. A catalog of transcript models (contigs) representing expressed genes was constructed from de novo assembled short complementary DNA (cDNA) sequences generated from developing fruit between flowering and maturity at 14 time points. Expression levels in each sample were estimated for 34 695 contigs from numbers of reads mapping to each contig. Contigs were annotated functionally based on BLAST, gene ontology and InterProScan analyses. Coregulated genes were detected using partitional clustering of expression patterns. The results are discussed with emphasis on genes putatively involved in cuticle deposition, cell wall metabolism and sugar transport. The high temporal resolution of the expression patterns presented here reveals finely tuned developmental specialization of individual members of gene families. Moreover, the de novo assembled sweet cherry fruit transcriptome with 7760 full-length protein coding sequences and over 20 000 other, annotated cDNA sequences together with their developmental expression patterns is expected to accelerate molecular research on this important tree fruit crop.

Synthetic regulation of eukaryotic gene expression by noncoding RNA

Synthetic biological systems promise to combine the spectacular diversity of biological functionality with engineering principles to design new life to address many pressing needs. As these engineered systems advance in sophistication, there is ever-greater need for customizable, situation-specific expression of desired genes. However, existing gene control platforms are generally not modular, or do not display performance requirements required for robust phenotypic responses to input signals. This work expands the capabilities of eukaryotic gene control in two important directions.

For development of greater modularity, we extend the use of synthetic self-cleaving ribozyme switches to detect changes in input protein levels and convey that information into programmed gene expression in eukaryotic cells. We demonstrate both up- and down-regulation of levels of an output transgene by more than 4-fold in response to rising input protein levels, with maximal output gene expression approaching the highest levels observed in yeast. In vitro experiments demonstrate protein-dependent ribozyme activity modulation. We further demonstrate the platform in mammalian cells. Our switch devices do not depend on special input protein activity, and can be tailored to respond to any input protein to which a suitable RNA aptamer can be developed. This platform can potentially be employed to regulate the expression of any transgene or any endogenous gene by 3’ UTR replacement, allowing for more complex cell state-specific reprogramming.

We also address an important concern with ribozyme switches, and riboswitch performance in general, their dynamic range. While riboswitches have generally allowed for versatile and modular regulation, so far their dynamic ranges of output gene modulation have been modest, generally at most 10-fold. We address this shortcoming by developing a modular genetic amplifier for near-digital control of eukaryotic gene expression. We combine ribozyme switch-mediated regulation of a synthetic TF with TF-mediated regulation of an output gene. The amplifier platform allows for as much as 20-fold regulation of output gene expression in response to input signal, with maximal expression approaching the highest levels observed in yeast, yet being tunable to intermediate and lower expression levels. EC50 values are more than 4 times lower than in previously best-performing non-amplifier ribozyme switches. The system design retains the modular-input architecture of the ribozyme switch platform, and the near-digital dynamic ranges of TF-based gene control.

Together, these developments suggest great potential for the wide applicability of these platforms for better-performing eukaryotic gene regulation, and more sophisticated, customizable reprogramming of cellular activity.

Arrangement of gene pairs, retrotransposon insertions, and regulation of gene expression in plants

Plant genomes are extremely complex. Myriad factors contribute to their evolution and organization, as well as to the expression and regulation of individual genes. Here we present investigations into several such factors and their influence on genome structure and gene expression: the arrangement of pairs of physically adjacent genes, retrotransposons closely associated with genes, and the effect of retrotransposons on gene pair evolution. All sequenced plant genomes contain a significant fraction of retrotransposons, including that of rice. We investigated the effects of retrotransposons within rice genes and within a 1 kb putative promoter region upstream of each gene. We found that approximately one-sixth of all rice genes are closely associated with retrotransposons. Insertions within a gene’s promoter region tend to block gene expression, while retrotransposons within genes promote the existence of alternative splicing forms. We also identified several other trends in retrotransposon insertion and its effects on gene expression. Several studies have previously noted a connection among genes between physical proximity and correlated expression profiles. To determine the degree to which this correlation depends on an exact physical arrangement, we studied the expression and interspecies conservation of convergent and divergent gene pairs in rice, Arabidopsis, and Populus trichocarpa. Correlated expression among gene pairs was quite common in all three species, yet conserved arrangement was rare. However, conservation of gene pair arrangement was significantly more common among pairs with strongly correlated expression levels. In order to uncover additional properties of gene pair conservation and rearrangement, we performed a comparative analysis of convergent, divergent, and tandem gene pairs in rice, sorghum, maize, and Brachypodium. We noted considerable differences between gene pair types and species. We also constructed a putative evolutionary history for each pair, which led to several interesting discoveries. To further elucidate the causes of gene pair conservation and rearrangement, we identified retrotransposon insertions in and near rice gene pairs. Retrotransposon-associated pairs are less likely to be conserved, although there are significant differences in the possible effect of different types and locations of retrotransposon insertions. The three types of gene pair also varied in their susceptibility to retrotransposon-associated evolutionary changes.

POPLAR GENE EXPRESSION DATA ANALYSIS PIPELINES

Analyzing large-scale gene expression data is a labor-intensive and time-consuming process. To make data analysis easier, we developed a set of pipelines for rapid processing and analysis poplar gene expression data for knowledge discovery. Of all pipelines developed, differentially expressed genes (DEGs) pipeline is the one designed to identify biologically important genes that are differentially expressed in one of multiple time points for conditions. Pathway analysis pipeline was designed to identify the differentially expression metabolic pathways. Protein domain enrichment pipeline can identify the enriched protein domains present in the DEGs. Finally, Gene Ontology (GO) enrichment analysis pipeline was developed to identify the enriched GO terms in the DEGs. Our pipeline tools can analyze both microarray gene data and high-throughput gene data. These two types of data are obtained by two different technologies. A microarray technology is to measure gene expression levels via microarray chips, a collection of microscopic DNA spots attached to a solid (glass) surface, whereas high throughput sequencing, also called as the next-generation sequencing, is a new technology to measure gene expression levels by directly sequencing mRNAs, and obtaining each mRNA’s copy numbers in cells or tissues. We also developed a web portal (http://sys.bio.mtu.edu/) to make all pipelines available to public to facilitate users to analyze their gene expression data. In addition to the analyses mentioned above, it can also perform GO hierarchy analysis, i.e. construct GO trees using a list of GO terms as an input.

Unconjugated bile acids influence expression of circadian genes: a potential mechanism for microbe-host crosstalk

Disruptions to circadian rhythm in mice and humans have been associated with an increased risk of obesity and metabolic syndrome. The gut microbiota is known to be essential for the maintenance of circadian rhythm in the host suggesting a role for microbe-host interactions in the regulation of the peripheral circadian clock. Previous work suggested a role for gut bacterial bile salt hydrolase (BSH) activity in the regulation of host circadian gene expression. Here we demonstrate that unconjugated bile acids, known to be generated through the BSH activity of the gut microbiota, are potentially chronobiological regulators of host circadian gene expression. We utilised a synchronised Caco-2 epithelial colorectal cell model and demonstrated that unconjugated bile acids, but not the equivalent tauro-conjugated bile salts, enhance the expression levels of genes involved in circadian rhythm. In addition oral administration of mice with unconjugated bile acids significantly altered expression levels of circadian clock genes in the ileum and colon as well as the liver with significant changes to expression of hepatic regulators of circadian rhythm (including Dbp) and associated genes (Per2, Per3 and Cry2). The data demonstrate a potential mechanism for microbe-host crosstalk that significantly impacts upon host circadian gene expression. Disruptions to circadian rhythm in mice and humans have been associated with an increased risk of obesity and metabolic syndrome. The gut microbiota is known to be essential for the maintenance of circadian rhythm in the host suggesting a role for microbe-host interactions in the regulation of the peripheral circadian clock. Previous work suggested a role for gut bacterial bile salt hydrolase (BSH) activity in the regulation of host circadian gene expression. Here we demonstrate that unconjugated bile acids, known to be generated through the BSH activity of the gut microbiota, are potentially chronobiological regulators of host circadian gene expression. We utilised a synchronised Caco-2 epithelial colorectal cell model and demonstrated that unconjugated bile acids, but not the equivalent tauro-conjugated bile salts, enhance the expression levels of genes involved in circadian rhythm. In addition oral administration of mice with unconjugated bile acids significantly altered expression levels of circadian clock genes in the ileum and colon as well as the liver with significant changes to expression of hepatic regulators of circadian rhythm (including Dbp) and associated genes (Per2, Per3 and Cry2). The data demonstrate a potential mechanism for microbe-host crosstalk that significantly impacts upon host circadian gene expression.

Selection And Validation Of Reference Genes For Functional Studies In The Calliphoridae Family.

The genera Cochliomyia and Chrysomya contain both obligate and saprophagous flies, which allows the comparison of different feeding habits between closely related species. Among the different strategies for comparing these habits is the use of qPCR to investigate the expression levels of candidate genes involved in feeding behavior. To ensure an accurate measure of the levels of gene expression, it is necessary to normalize the amount of the target gene with the amount of a reference gene having a stable expression across the compared species. Since there is no universal gene that can be used as a reference in functional studies, candidate genes for qPCR data normalization were selected and validated in three Calliphoridae (Diptera) species, Cochliomyia hominivorax Coquerel, Cochliomyia macellaria Fabricius, and Chrysomya albiceps Wiedemann . The expression stability of six genes ( Actin, Gapdh, Rp49, Rps17, α -tubulin, and GstD1) was evaluated among species within the same life stage and between life stages within each species. The expression levels of Actin, Gapdh, and Rp49 were the most stable among the selected genes. These genes can be used as reliable reference genes for functional studies in Calliphoridae using similar experimental settings.