Differential expression of microRNAs in mouse pain models

Background: MicroRNAs (miRNAs) are short non-coding RNAs that inhibit translation of target genes by binding to their mRNAs. The expression of numerous brain-specific miRNAs with a high degree of temporal and spatial specificity suggests that miRNAs play an important role in gene regulation in health and disease. Here we investigate the time course gene expression profile of miR-1, -16, and -206 in mouse dorsal root ganglion (DRG), and spinal cord dorsal horn under inflammatory and neuropathic pain conditions as well as following acute noxious stimulation. Results: Quantitative real-time polymerase chain reaction analyses showed that the mature form of miR-1, -16 and -206, is expressed in DRG and the dorsal horn of the spinal cord. Moreover, CFA-induced inflammation significantly reduced miRs-1 and -16 expression in DRG whereas miR-206 was downregulated in a time dependent manner. Conversely, in the spinal dorsal horn all three miRNAs monitored were upregulated. After sciatic nerve partial ligation, miR-1 and -206 were downregulated in DRG with no change in the spinal dorsal horn. On the other hand, axotomy increases the relative expression of miR-1, -16, and 206 in a time-dependent fashion while in the dorsal horn there was a significant downregulation of miR-1. Acute noxious stimulation with capsaicin also increased the expression of miR-1 and -16 in DRG cells but, on the other hand, in the spinal dorsal horn only a high dose of capsaicin was able to downregulate miR-206 expression. Conclusions: Our results indicate that miRNAs may participate in the regulatory mechanisms of genes associated with the pathophysiology of chronic pain as well as the nociceptive processing following acute noxious stimulation. We found substantial evidence that miRNAs are differentially regulated in DRG and the dorsal horn of the spinal cord under different pain states. Therefore, miRNA expression in the nociceptive system shows not only temporal and spatial specificity but is also stimulus-dependent.

The insulin signaling pathway in honey bee (Apis mellifera) caste development - differential expression of insulin-like peptides and insulin receptors in queen and worker larvae

The insulin/insulin-like signaling (IIS) pathway is an evolutionarily conserved module in the control of body size and correlated organ growth in metazoans. In the highly eusocial bees, the caste phenotypes differ not only in size and several structural features but also in individual fitness and life history. We investigated the developmental expression profiles of genes encoding the two insulin-like peptides (AmILP-1 and AmILP-2) and the two insulin receptors (AmInR-1 and AmInR-2) predicted in the honey bee genome. Quantitative PCR analysis for queen and worker larvae in critical stages of caste development showed that AmILP-2 is the predominantly transcribed ILP in both castes, with higher expression in workers than in queens. Expression of both InR genes sharply declined in fourth instar queen larvae, but showed little modulation in workers. On first sight, these findings are non-intuitive, considering the higher growth rates of queens, but they can be interpreted as possibly antagonistic crosstalk between the IIS module and juvenile hormone. Analyzing AmInR-1 and AmInR-2 expression in ovaries of queen and worker larvae revealed low transcript levels in queens and a sharp drop in AmInR-2 expression in fifth instar worker larvae, indicating relative independence in tissue-specific versus overall IIS pathway activity. (C) 2008 Elsevier Ltd. All rights reserved.

Differential expression of HDAC3, HDAC7 and HDAC9 is associated with prognosis and survival in childhood acute lymphoblastic leukaemia

Altered expression of histone deacetylases (HDACs) is a common feature in several human malignancies and may represent an interesting target for cancer treatment, including haematological malignancies. We evaluated the mRNA gene expression profile of 12 HDAC genes by quantitative real-time polymerase chain reaction in 94 consecutive childhood acute lymphoblastic leukaemia (ALL) samples and its association with clinical/biological features and survival. ALL samples showed higher expression levels of HDAC2, HDAC3, HDAC8, HDAC6 and HDAC7 when compared to normal bone marrow samples. HDAC1 and HDAC4 showed high expression in T-ALL and HDAC5 was highly expressed in B-lineage ALL. Higher than median expression levels of HDAC3 were associated with a significantly lower 5-year event-free survival (EFS) in the overall group of patients (P = 0.03) and in T-ALL patients (P = 0.01). HDAC7 and HADC9 expression levels higher than median were associated with a lower 5-year EFS in the overall group (P = 0.04 and P = 0.003, respectively) and in B-lineage CD10-positive patients (P = 0.009 and P = 0.005, respectively). Our data suggest that higher expression of HDAC7 and HDAC9 is associated with poor prognosis in childhood ALL and could be promising therapeutic targets for the treatment of refractory childhood ALL.

Differential expression of CDC25 phosphatases splice variants in human breast cancer cells

Background: CDC25 phosphatases control cell cycle progression by activating cyclin dependent kinases. The three CDC25 isoforms encoding genes are submitted to alternative splicing events which generate at least two variants for CDC25A and five for both CDC25B and CDC25C. An over-expression of CDC25 was reported in several types of cancer, including breast cancer, and is often associated with a poor prognosis. Nevertheless, most of the previous studies did not address the expression of CDC25 splice variants. Here, we evaluated CDC25 spliced transcripts expression in anti-cancerous drug-sensitive and resistant breast cancer cell lines in order to identify potential breast cancer biomarkers. Methods: CDC25 splice variants mRNA levels were evaluated by semi-quantitative RT-PCR and by an original real-time RT-PCR assay. Results: CDC25 spliced transcripts are differentially expres-sed in the breast cancer cell lines studied. An up-regulation of CDC25A2 variant and an increase of the CDC25C5/C1 ratio are associated to the multidrug-resistance in VCREMS and DOXOR breast cancer cells, compared to their sensitive counterpart cell line MCF-7. Additionally, CDC25B2 tran-script is exclusively over-expressed in VCREMS resistant cells and could therefore be involved in the development of certain type of drug resistance. Conclusions: CDC25 splice variants could represent interesting potential breast cancer prognostic biomarkers.

Differential expression of peroxisome proliferator-activated receptors (PPARs): tissue distribution of PPAR-alpha, -beta, and -gamma in the adult rat.

Peroxisome proliferator-activated receptors (PPARs) are members of the nuclear hormone receptor superfamily that can be activated by various xenobiotics and natural fatty acids. These transcription factors primarily regulate genes involved in lipid metabolism and also play a role in adipocyte differentiation. We present the expression patterns of the PPAR subtypes in the adult rat, determined by in situ hybridization using specific probes for PPAR-alpha, -beta and -gamma, and by immunohistochemistry using a polyclonal antibody that recognizes the three rat PPAR subtypes. In numerous cell types from either ectodermal, mesodermal, or endodermal origin, PPARs are coexpressed, with relative levels varying between them from one cell type to the other. PPAR-alpha is highly expressed in hepatocytes, cardiomyocytes, enterocytes, and the proximal tubule cells of kidney. PPAR-beta is expressed ubiquitously and often at higher levels than PPAR-alpha and -gamma. PPAR-gamma is expressed predominantly in adipose tissue and the immune system. Our results suggest new potential directions to investigate the functions of the different PPAR subtypes.

A comparison of methods for differential expression analysis of RNA-seq data.

BACKGROUND: Finding genes that are differentially expressed between conditions is an integral part of understanding the molecular basis of phenotypic variation. In the past decades, DNA microarrays have been used extensively to quantify the abundance of mRNA corresponding to different genes, and more recently high-throughput sequencing of cDNA (RNA-seq) has emerged as a powerful competitor. As the cost of sequencing decreases, it is conceivable that the use of RNA-seq for differential expression analysis will increase rapidly. To exploit the possibilities and address the challenges posed by this relatively new type of data, a number of software packages have been developed especially for differential expression analysis of RNA-seq data. RESULTS: We conducted an extensive comparison of eleven methods for differential expression analysis of RNA-seq data. All methods are freely available within the R framework and take as input a matrix of counts, i.e. the number of reads mapping to each genomic feature of interest in each of a number of samples. We evaluate the methods based on both simulated data and real RNA-seq data. CONCLUSIONS: Very small sample sizes, which are still common in RNA-seq experiments, impose problems for all evaluated methods and any results obtained under such conditions should be interpreted with caution. For larger sample sizes, the methods combining a variance-stabilizing transformation with the 'limma' method for differential expression analysis perform well under many different conditions, as does the nonparametric SAMseq method.

Differential expression of peroxisome proliferator-activated receptor-alpha, -beta, and -gamma during rat embryonic development.

The expression patterns of the three different peroxisome proliferator-activated receptor (PPAR) isotypes have been determined during rat embryonic development by in situ hybridization. The expression of PPARalpha starts late in development, with increasing levels in organs such as liver, kidney, intestine, and pancreas, in which it will also be present later in adulthood to regulate its specific target genes. PPARalpha is also transiently expressed in the embryonic epidermis and central nervous system. PPARgamma presents a very restricted pattern of expression, being strongly expressed in brown adipose tissue, in which differentiation it has been shown to participate. Like PPARalpha, it is also expressed transiently in the central nervous system. Interestingly, PPARalpha, -beta and -gamma are coexpressed at high levels in brown adipose tissue. Finally, the high and ubiquitous expression of PPARbeta suggests some fundamental role(s) that this receptor might play throughout development.

Differential expression and hypoxic regulation of adenosine receptors and TRPC channels in the developing heart.

Objectives: To characterize the modifications of gene expression of adenosine receptors (AR), TRPC channels, HIF-1α and iNOS during the early cardiogenesis in response to chronic hypoxia exposure. Methods: 4-day-old chick embryos were subjected in ovo to 6H, 12H and 24H of hypoxia (10% O2). The mRNA expression was quantified by RT-qPCR. Results: The targeted genes were found to be expressed at mRNA level with a differential expression pattern within the heart. Hypoxia has no significant effect on mRNA expression of ARs, TRPCs channels and iNOS within the heart. By contrast, HIF-1α mRNA expression shows a tendency to be down-regulated by hypoxia. Conclusion: These results suggest that an intrauterine oxygen lack does not significantly affect expression of genes involved in adenosine signaling and in calcium handling by store operated channels (TRPC).

Differential expression of NLRP3 among hematopoietic cells.

Although the importance of the NOD-like receptor family, pyrin domain containing 3 (NLRP3) inflammasome in health and disease is well appreciated, a precise characterization of NLRP3 expression is yet undetermined. To this purpose, we generated a knock-in mouse in which the Nlrp3 coding sequence was substituted for the GFP (enhanced GFP [egfp]) gene. In this way, the expression of eGFP is driven by the endogenous regulatory elements of the Nlrp3 gene. In this study, we show that eGFP expression indeed mirrors that of NLRP3. Interestingly, splenic neutrophils, macrophages, and, in particular, monocytes and conventional dendritic cells showed robust eGFP fluorescence, whereas lymphoid subsets, eosinophils, and plasmacytoid dendritic cells showed negligible eGFP levels. NLRP3 expression was highly inducible in macrophages, both by MyD88- and Trif-dependent pathways. In vivo, when mice were challenged with diverse inflammatory stimuli, differences in both the number of eGFP-expressing cells and fluorescence intensity were observed in the draining lymph node. Thus, NLRP3 levels at the site of adaptive response initiation are controlled by recruitment of NLRP3-expressing cells and by NLRP3 induction.

Differential Expression of ADAM23, CDKN2A (P16), MMP14 and VIM Associated with Giant Cell Tumor of Bone

Though benign, giant cell tumor of bone (GCTB) can become aggressive and can exhibit a high mitotic rate, necrosis and rarely vascular invasion and metastasis. GCTB has unique histologic characteristics, a high rate of multinucleated cells, a variable and unpredictable growth potential and uncertain biological behavior. In this study, we sought to identify genes differentially expressed in GCTB, thus building a molecular profile of this tumor. We performed quantitative real-time polymerase chain reaction (qPCR), immunohistochemistry and analyses of methylation to identify genes that are putatively associated with GCTB. The expression of the ADAM23 and CDKN2A genes was decreased in GCTB samples compared to normal bone tissue, measured by qPCR. Additionally, a high hypermethylation frequency of the promoter regions of ADAM23 and CDKN2A in GCTB was observed. The expression of the MAP2K3, MMP14, TIMP2 and VIM genes was significantly higher in GCTB than in normal bone tissue, a fact that was confirmed by qPCR and immunohistochemistry. The set of genes identified here furthers our understanding of the molecular basis of GCTB.

Differential expression of toll-like receptor mRNAs in recurrent aphthous ulceration

BACKGROUND: Toll-like receptors (TLR) are membrane proteins that recognize conserved molecules derived from bacterial, viral, fungal or host tissues. They are responsible for promoting the production of cytokines and chemokines, increasing the expression of costimulatory molecules and influencing the T Helper response (Th) toward either a Th1 or Th2 profile, thereby modulating the regulatory T cell response and controlling the integrity of the epithelial barrier. The key factors responsible for increased susceptibility to recurrent aphthous ulceration (RAU) are unclear, and because TLRs are involved in both immune regulation and control of the epithelial barrier, a deficiency in TLR activity is likely to cause increased susceptibility. METHODS: We investigated the gene expression of TLRs one through 10 in tissue samples and peripheral blood mononuclear cells (PBMC) of RAU patients in comparison to healthy controls using real-time quantitative reverse transcription PCR. RESULTS: The analysis of mRNA expression levels in oral lesion showed significant (P < 0.01) overexpression of the TLR2(similar to 6-fold) gene and decreased expression of the TLR3 (similar to 5-fold) and TLR5 (similar to 6-fold) genes in comparison with healthy oral mucosa. The analysis of mRNA expression in PBMC indicated a down-regulation of TLR5 gene expression in the cells from RAU patients (P < 0.05; similar to 2-fold). CONCLUSION: Our results support the hypothesis that a subset of RAU patients has fewer TLR expression that have been tentatively implicated in antiinflammatory effects. This derangement of TLR gene expression may cause an overlay exuberant inflammation reaction in situations where normal individuals are resistant. J Oral Pathol Med (2012) 41: 8085

Differential expression analysis for sequence count data via mixtures of negative binomials

The recent advent of Next-generation sequencing technologies has revolutionized the way of analyzing the genome. This innovation allows to get deeper information at a lower cost and in less time, and provides data that are discrete measurements. One of the most important applications with these data is the differential analysis, that is investigating if one gene exhibit a different expression level in correspondence of two (or more) biological conditions (such as disease states, treatments received and so on). As for the statistical analysis, the final aim will be statistical testing and for modeling these data the Negative Binomial distribution is considered the most adequate one especially because it allows for "over dispersion". However, the estimation of the dispersion parameter is a very delicate issue because few information are usually available for estimating it. Many strategies have been proposed, but they often result in procedures based on plug-in estimates, and in this thesis we show that this discrepancy between the estimation and the testing framework can lead to uncontrolled first-type errors. We propose a mixture model that allows each gene to share information with other genes that exhibit similar variability. Afterwards, three consistent statistical tests are developed for differential expression analysis. We show that the proposed method improves the sensitivity of detecting differentially expressed genes with respect to the common procedures, since it is the best one in reaching the nominal value for the first-type error, while keeping elevate power. The method is finally illustrated on prostate cancer RNA-seq data.

Differential expression of TGFbeta-stimulated clone 22 in normal prostate and prostate cancer

The transforming growth factor-beta (TGFbeta) superfamily and its downstream effector genes are key regulators of epithelial homeostasis. Altered expression of these genes may be associated with malignant transformation of the prostate gland. The cDNA array analysis of differential expression of the TGFbeta superfamily and functionally related genes between patient-matched noncancerous prostate (NP) and prostate cancer (PC) bulk tissue specimens highlighted two genes, namely TGFbeta-stimulated clone-22 (TSC-22) and Id4. Verification of their mRNA expression by real-time PCR in patient-matched NP and PC bulk tissue, in laser-captured pure epithelial and cancer cells and in NP and PC cell lines confirmed TSC-22 underexpression, but not Id4 overexpression, in PC and in human PC cell lines. Immunohistochemical analysis showed that TSC-22 protein expression in NP is restricted to the basal cells and colocalizes with the basal cell marker cytokeratin 5. In contrast, all matched PC samples lack TSC-22 immunoreactivity. Likewise, PC cell lines do not show detectable TSC-22 protein expression as shown by immunoblotting. TSC-22 should be considered as a novel basal cell marker, potentially useful for studying lineage determination within the epithelial compartment of the prostate. Conversely, lack of TSC-22 seems to be a hallmark of malignant transformation of the prostate epithelium. Accordingly, TSC-22 immunohistochemistry may prove to be a diagnostic tool for discriminating benign lesions from malignant ones of the prostate. The suggested tumour suppressor function of TSC-22 warrants further investigation on its role in prostate carcinogenesis and on the TSC-22 pathway as a candidate therapeutic target in PC.

Differential Expression with the Bioconductor Project

A basic, yet challenging task in the analysis of microarray gene expression data is the identification of changes in gene expression that are associated with particular biological conditions. We discuss different approaches to this task and illustrate how they can be applied using software from the Bioconductor Project. A central problem is the high dimensionality of gene expression space, which prohibits a comprehensive statistical analysis without focusing on particular aspects of the joint distribution of the genes expression levels. Possible strategies are to do univariate gene-by-gene analysis, and to perform data-driven nonspecific filtering of genes before the actual statistical analysis. However, more focused strategies that make use of biologically relevant knowledge are more likely to increase our understanding of the data.

Differential expression of stx2 variants in Shiga toxin-producing Escherichia coli belonging to seropathotypes A and C

Only a subset of Shiga toxin (Stx)-producing Escherichia coli (STEC) are human pathogens, but the characteristics that account for differences in pathogenicity are not well understood. In this study, we investigated the distribution of the stx variants coding for Stx2 and its variants in highly virulent STEC of seropathotype A and low-pathogenic STEC of seropathotype C. We analysed and compared transcription of the corresponding genes, production of Shiga toxins, and stx-phage release in basal as well as in induced conditions. We found that the stx(2) variant was mainly associated with strains of seropathotype A, whereas most of the strains of seropathotype C possessed the stx(2-vhb) variant, which was frequently associated with stx(2), stx(2-vha) or stx(2c). Levels of stx(2) and stx(2)-related mRNA were higher in strains belonging to seropathotype A and in those strains of seropathotype C that express the stx(2) variant than in the remaining strains of seropathotype C. The stx(2-vhb) genes were the least expressed, in basal as well as in induced conditions, and in many cases did not seem to be carried by an inducible prophage. A clear correlation was observed between stx mRNA levels and stx-phage DNA in the culture supernatants, suggesting that most stx(2)-related genes are expressed only when they are carried by a phage. In conclusion, some relationship between stx(2)-related gene expression in vitro and the seropathotype of the STEC strains was observed. A higher expression of the stx(2) gene and a higher release of its product, in basal as well as in induced conditions, was observed in pathogenic strains of seropathotype A. A subset of strains of seropathotype C shows the same characteristics and could be a high risk to human health.