Reverse transcription quantitative real-time polymerase chain reaction reference genes in the spared nerve injury model of neuropathic pain : validation and literature search

La douleur neuropathique est définie comme une douleur causée par une lésion du système nerveux somato-sensoriel. Elle se caractérise par des douleurs exagérées, spontanées, ou déclenchées par des stimuli normalement non douloureux (allodynie) ou douloureux (hyperalgésie). Bien qu'elle concerne 7% de la population, ses mécanismes biologiques ne sont pas encore élucidés. L'étude des variations d'expressions géniques dans les tissus-clés des voies sensorielles (notamment le ganglion spinal et la corne dorsale de la moelle épinière) à différents moments après une lésion nerveuse périphérique permettrait de mettre en évidence de nouvelles cibles thérapeutiques. Elles se détectent de manière sensible par reverse transcription quantitative real-time polymerase chain reaction (RT- qPCR). Pour garantir des résultats fiables, des guidelines ont récemment recommandé la validation des gènes de référence utilisés pour la normalisation des données ("Minimum information for publication of quantitative real-time PCR experiments", Bustin et al 2009). Après recherche dans la littérature des gènes de référence fréquemment utilisés dans notre modèle de douleur neuropathique périphérique SNI (spared nerve injury) et dans le tissu nerveux en général, nous avons établi une liste de potentiels bons candidats: Actin beta (Actb), Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal proteins 18S (18S), L13a (RPL13a) et L29 (RPL29), hypoxanthine phosphoribosyltransferase 1 (HPRT1) et hydroxymethyl-bilane synthase (HMBS). Nous avons évalué la stabilité d'expression de ces gènes dans le ganglion spinal et dans la corne dorsale à différents moments après la lésion nerveuse (SNI) en calculant des coefficients de variation et utilisant l'algorithme geNorm qui compare les niveaux d'expression entre les différents candidats et détermine la paire de gènes restante la plus stable. Il a aussi été possible de classer les gènes selon leur stabilité et d'identifier le nombre de gènes nécessaires pour une normalisation la plus précise. Les gènes les plus cités comme référence dans le modèle SNI ont été GAPDH, HMBS, Actb, HPRT1 et 18S. Seuls HPRT1 and 18S ont été précédemment validés dans des arrays de RT-qPCR. Dans notre étude, tous les gènes testés dans le ganglion spinal et dans la corne dorsale satisfont au critère de stabilité exprimé par une M-value inférieure à 1. Par contre avec un coefficient de variation (CV) supérieur à 50% dans le ganglion spinal, 18S ne peut être retenu. La paire de gènes la plus stable dans le ganglion spinal est HPRT1 et Actb et dans la corne dorsale il s'agit de RPL29 et RPL13a. L'utilisation de 2 gènes de référence stables suffit pour une normalisation fiable. Nous avons donc classé et validé Actb, RPL29, RPL13a, HMBS, GAPDH, HPRT1 et 18S comme gènes de référence utilisables dans la corne dorsale pour le modèle SNI chez le rat. Dans le ganglion spinal 18S n'a pas rempli nos critères. Nous avons aussi déterminé que la combinaison de deux gènes de référence stables suffit pour une normalisation précise. Les variations d'expression génique de potentiels gènes d'intérêts dans des conditions expérimentales identiques (SNI, tissu et timepoints post SNI) vont pouvoir se mesurer sur la base d'une normalisation fiable. Non seulement il sera possible d'identifier des régulations potentiellement importantes dans la genèse de la douleur neuropathique mais aussi d'observer les différents phénotypes évoluant au cours du temps après lésion nerveuse.

Reverse transcription quantitative real-time polymerase chain reaction reference genes in the spared nerve injury model of neuropathic pain: validation and literature search.

BACKGROUND: The reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) is a widely used, highly sensitive laboratory technique to rapidly and easily detect, identify and quantify gene expression. Reliable RT-qPCR data necessitates accurate normalization with validated control genes (reference genes) whose expression is constant in all studied conditions. This stability has to be demonstrated.We performed a literature search for studies using quantitative or semi-quantitative PCR in the rat spared nerve injury (SNI) model of neuropathic pain to verify whether any reference genes had previously been validated. We then analyzed the stability over time of 7 commonly used reference genes in the nervous system - specifically in the spinal cord dorsal horn and the dorsal root ganglion (DRG). These were: Actin beta (Actb), Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal proteins 18S (18S), L13a (RPL13a) and L29 (RPL29), hypoxanthine phosphoribosyltransferase 1 (HPRT1) and hydroxymethylbilane synthase (HMBS). We compared the candidate genes and established a stability ranking using the geNorm algorithm. Finally, we assessed the number of reference genes necessary for accurate normalization in this neuropathic pain model. RESULTS: We found GAPDH, HMBS, Actb, HPRT1 and 18S cited as reference genes in literature on studies using the SNI model. Only HPRT1 and 18S had been once previously demonstrated as stable in RT-qPCR arrays. All the genes tested in this study, using the geNorm algorithm, presented gene stability values (M-value) acceptable enough for them to qualify as potential reference genes in both DRG and spinal cord. Using the coefficient of variation, 18S failed the 50% cut-off with a value of 61% in the DRG. The two most stable genes in the dorsal horn were RPL29 and RPL13a; in the DRG they were HPRT1 and Actb. Using a 0.15 cut-off for pairwise variations we found that any pair of stable reference gene was sufficient for the normalization process. CONCLUSIONS: In the rat SNI model, we validated and ranked Actb, RPL29, RPL13a, HMBS, GAPDH, HPRT1 and 18S as good reference genes in the spinal cord. In the DRG, 18S did not fulfill stability criteria. The combination of any two stable reference genes was sufficient to provide an accurate normalization.

Telomerase activity and human telomerase reverse transcriptase mRNA expression in soft tissue tumors: correlation with grade, histology, and proliferative activity.

Telomerase activity (TA) is detected in most human cancers but, with few exceptions, not in normal somatic cells. Little is known about TA in soft tissue tumors. We have examined a series of benign and malignant soft tissue tumors for TA using the telomerase repeat amplification protocol assay. Analysis of the expression of the human telomerase reverse transcriptase was also carried out using RT-PCR. TA was undetectable in benign lesions (15 of 15) and low-grade sarcomas (6 of 6) and was detectable in 50% (19 of 38) of intermediate-/high-grade sarcomas. Although the presence of TA in soft tissue tumors is synonymous with malignancy, it is neither a reliable method in making the distinction between reactive/benign and malignant (especially low-grade) lesions nor a reliable marker of tumor aggressiveness. Leiomyosarcomas and storiform/pleomorphic malignant fibrous histiocytomas rarely showed TA, irrespective of their grade. A strong correlation between human telomerase reverse transcriptase mRNA expression and TA was observed, supporting the close relationship between both parameters. No significant relationship was observed between proliferative activity (as assessed by MIB-1 immunolabeling) and TA. We verified that the absence of telomerase expression was not due to the presence of telomerase inhibitors and therefore alternative mechanism(s) for cell immortalization, yet to be determined, seem to be involved in the development and/or maintenance of some soft tissue sarcomas.

Selecting control genes for RT-QPCR using public microarray data.

Background: Gene expression analysis has emerged as a major biological research area, with real-time quantitative reverse transcription PCR (RT-QPCR) being one of the most accurate and widely used techniques for expression profiling of selected genes. In order to obtain results that are comparable across assays, a stable normalization strategy is required. In general, the normalization of PCR measurements between different samples uses one to several control genes (e. g. housekeeping genes), from which a baseline reference level is constructed. Thus, the choice of the control genes is of utmost importance, yet there is not a generally accepted standard technique for screening a large number of candidates and identifying the best ones. Results: We propose a novel approach for scoring and ranking candidate genes for their suitability as control genes. Our approach relies on publicly available microarray data and allows the combination of multiple data sets originating from different platforms and/or representing different pathologies. The use of microarray data allows the screening of tens of thousands of genes, producing very comprehensive lists of candidates. We also provide two lists of candidate control genes: one which is breast cancer-specific and one with more general applicability. Two genes from the breast cancer list which had not been previously used as control genes are identified and validated by RT-QPCR. Open source R functions are available at http://www.isrec.isb-sib.ch/similar to vpopovic/research/ Conclusion: We proposed a new method for identifying candidate control genes for RT-QPCR which was able to rank thousands of genes according to some predefined suitability criteria and we applied it to the case of breast cancer. We also empirically showed that translating the results from microarray to PCR platform was achievable.

Identification of a new cancer/testis gene family, CT47, among expressed multicopy genes on the human X chromosome.

Cancer/testis (CT) genes are normally expressed in germ cells only, yet are reactivated and expressed in some tumors. Of the approximately 40 CT genes or gene families identified to date, 20 are on the X chromosome and are present as multigene families, many with highly conserved members. This indicates that novel CT gene families may be identified by detecting duplicated expressed genes on chromosome X. By searching for transcript clusters that map to multiple locations on the chromosome, followed by in silico analysis of their gene expression profiles, we identified five novel gene families with testis-specific expression and >98% sequence identity among family members. The expression of these genes in normal tissues and various tumor cell lines and specimens was evaluated by qualitative and quantitative RT-PCR, and a novel CT gene family with at least 13 copies was identified on Xq24, designated as CT47. mRNA expression of CT47 was found mainly in the testes, with weak expression in the placenta. Brain tissue was the only positive somatic tissue tested, with an estimated CT47 transcript level 0.09% of that found in testis. Among the tumor specimens tested, CT47 expression was found in approximately 15% of lung cancer and esophageal cancer specimens, but not in colorectal cancer or breast cancer. The putative CT47 protein consists of 288 amino acid residues, with a C-terminus rich in alanine and glutamic acid. The only species other than human in which a gene homologous to CT47 has been detected is the chimpanzee, with the predicted protein showing approximately 80% identity in its carboxy terminal region.

Agenesis of human pancreas due to decreased half-life of insulin promoter factor 1.

Neonatal diabetes mellitus can be transient or permanent. The severe form of permanent neonatal diabetes mellitus can be associated with pancreas agenesis. Normal pancreas development is controlled by a cascade of transcription factors, where insulin promoter factor 1 (IPF1) plays a crucial role. Here, we describe two novel mutations in the IPF1 gene leading to pancreas agenesis. Direct sequence analysis of exons 1 and 2 of the IPF1 gene revealed two point mutations within the homeobox in exon 2. Genetic analysis of the parents showed that each mutation was inherited from one parent. Mutations localized in helices 1 and 2, respectively, of the homeodomain, decreased the protein half-life significantly, leading to intracellular IPF1 levels of 36% and 27% of wild-type levels. Both mutant forms of IPF1 were normally translocated to the nucleus, and their DNA binding activity on different known target promoters was similar to that of the wild-type protein. However, transcriptional activity of both mutant IPF1 proteins, alone or in combination with HNF3 beta/Foxa2, Pbx1, or the heterodimer E47-beta 2 was reduced, findings accounted for by decreased IPF1 steady state levels and not by impaired protein-protein interactions. We conclude that the IPF1 level is critical for human pancreas formation.

Ribonucleotide reductase genes of Bacillus prophages: a refuge to introns and intein coding sequences.

The ribonucleotide reductase gene tandem bnrdE/bnrdF in SPbeta-related prophages of different Bacillus spp. isolates presents different configurations of intervening sequences, comprising one to three of six non-homologous splicing elements. Insertion sites of group I introns and intein DNA are clustered in three relatively short segments encoding functionally important domains of the ribonucleotide reductase. Comparison of the bnrdE homologs reveals mutual exclusion of a group I intron and an intein coding sequence flanking the codon that specifies a conserved cysteine. In vivo splicing was demonstrated for all introns. However, for two of them a part of the mRNA precursor molecules remains unspliced. Intergenic bnrdE-bnrdF regions are unexpectedly long, comprising between 238 and 541 nt. The longest encodes a putative polypeptide related to HNH homing endonucleases.

Emergence of young human genes after a burst of retroposition in primates.

The origin of new genes through gene duplication is fundamental to the evolution of lineage- or species-specific phenotypic traits. In this report, we estimate the number of functional retrogenes on the lineage leading to humans generated by the high rate of retroposition (retroduplication) in primates. Extensive comparative sequencing and expression studies coupled with evolutionary analyses and simulations suggest that a significant proportion of recent retrocopies represent bona fide human genes. We estimate that at least one new retrogene per million years emerged on the human lineage during the past approximately 63 million years of primate evolution. Detailed analysis of a subset of the data shows that the majority of retrogenes are specifically expressed in testis, whereas their parental genes show broad expression patterns. Consistently, most retrogenes evolved functional roles in spermatogenesis. Proteins encoded by X chromosome-derived retrogenes were strongly preserved by purifying selection following the duplication event, supporting the view that they may act as functional autosomal substitutes during X-inactivation of late spermatogenesis genes. Also, some retrogenes acquired a new or more adapted function driven by positive selection. We conclude that retroduplication significantly contributed to the formation of recent human genes and that most new retrogenes were progressively recruited during primate evolution by natural and/or sexual selection to enhance male germline function.

Gene expression profiling of isolated tumour cells from anaplastic large cell lymphomas: insights into its cellular origin, pathogenesis and relation to Hodgkin lymphoma.

Anaplastic large cell lymphoma (ALCL) is a main type of T-cell lymphomas and comprises three distinct entities: systemic anaplastic lymphoma kinase (ALK) positive, systemic ALK(-) and cutaneous ALK(-) ALCL (cALCL). Little is known about their pathogenesis and their cellular origin, and morphological and immunophenotypical overlap exists between ALK(-) ALCL and classical Hodgkin lymphoma (cHL). We conducted gene expression profiling of microdissected lymphoma cells of five ALK(+) and four ALK(-) systemic ALCL, seven cALCL and sixteen cHL, and of eight subsets of normal T and NK cells. The analysis supports a derivation of ALCL from activated T cells, but the lymphoma cells acquired a gene expression pattern hampering an assignment to a CD4(+), CD8(+) or CD30(+) T-cell origin. Indeed, ALCL display a down-modulation of many T-cell characteristic molecules. All ALCL types show significant expression of NFkappaB target genes and upregulation of genes involved in oncogenesis (e.g. EZH2). Surprisingly, few genes are differentially expressed between systemic and cALCL despite their different clinical behaviour, and between ALK(-) ALCL and cHL despite their different cellular origin. ALK(+) ALCL are characterized by expression of genes regulated by pathways constitutively activated by ALK. This study provides multiple novel insights into the molecular biology and pathogenesis of ALCL.

Phosphorylation regulates FOXC2-mediated transcription in lymphatic endothelial cells.

One of the key mechanisms linking cell signaling and control of gene expression is reversible phosphorylation of transcription factors. FOXC2 is a forkhead transcription factor that is mutated in the human vascular disease lymphedema-distichiasis and plays an essential role in lymphatic vascular development. However, the mechanisms regulating FOXC2 transcriptional activity are not well understood. We report here that FOXC2 is phosphorylated on eight evolutionarily conserved proline-directed serine/threonine residues. Loss of phosphorylation at these sites triggers substantial changes in the FOXC2 transcriptional program. Through genome-wide location analysis in lymphatic endothelial cells, we demonstrate that the changes are due to selective inhibition of FOXC2 recruitment to chromatin. The extent of the inhibition varied between individual binding sites, suggesting a novel rheostat-like mechanism by which expression of specific genes can be differentially regulated by FOXC2 phosphorylation. Furthermore, unlike the wild-type protein, the phosphorylation-deficient mutant of FOXC2 failed to induce vascular remodeling in vivo. Collectively, our results point to the pivotal role of phosphorylation in the regulation of FOXC2-mediated transcription in lymphatic endothelial cells and underscore the importance of FOXC2 phosphorylation in vascular development.

Identification of novel functional TBP-binding sites and general factor repertoires.

Our current knowledge of the general factor requirement in transcription by the three mammalian RNA polymerases is based on a small number of model promoters. Here, we present a comprehensive chromatin immunoprecipitation (ChIP)-on-chip analysis for 28 transcription factors on a large set of known and novel TATA-binding protein (TBP)-binding sites experimentally identified via ChIP cloning. A large fraction of identified TBP-binding sites is located in introns or lacks a gene/mRNA annotation and is found to direct transcription. Integrated analysis of the ChIP-on-chip data and functional studies revealed that TAF12 hitherto regarded as RNA polymerase II (RNAP II)-specific was found to be also involved in RNAP I transcription. Distinct profiles for general transcription factors and TAF-containing complexes were uncovered for RNAP II promoters located in CpG and non-CpG islands suggesting distinct transcription initiation pathways. Our study broadens the spectrum of general transcription factor function and uncovers a plethora of novel, functional TBP-binding sites in the human genome.

Early onset collagen VI myopathies: Genetic and clinical correlations.

OBJECTIVE: Mutations in the genes encoding the extracellular matrix protein collagen VI (ColVI) cause a spectrum of disorders with variable inheritance including Ullrich congenital muscular dystrophy, Bethlem myopathy, and intermediate phenotypes. We extensively characterized, at the clinical, cellular, and molecular levels, 49 patients with onset in the first 2 years of life to investigate genotype-phenotype correlations. METHODS: Patients were classified into 3 groups: early-severe (18%), moderate-progressive (53%), and mild (29%). ColVI secretion was analyzed in patient-derived skin fibroblasts. Chain-specific transcript levels were quantified by quantitative reverse transcriptase polymerase chain reaction (qRT-PCR), and mutation identification was performed by sequencing of complementary DNA. RESULTS: ColVI secretion was altered in all fibroblast cultures studied. We identified 56 mutations, mostly novel and private. Dominant de novo mutations were detected in 61% of the cases. Importantly, mutations causing premature termination codons (PTCs) or in-frame insertions strikingly destabilized the corresponding transcripts. Homozygous PTC-causing mutations in the triple helix domains led to the most severe phenotypes (ambulation never achieved), whereas dominant de novo in-frame exon skipping and glycine missense mutations were identified in patients of the moderate-progressive group (loss of ambulation). INTERPRETATION: This work emphasizes that the diagnosis of early onset ColVI myopathies is arduous and time-consuming, and demonstrates that quantitative RT-PCR is a helpful tool for the identification of some mutation-bearing genes. Moreover, the clinical classification proposed allowed genotype-phenotype relationships to be explored, and may be useful in the design of future clinical trials.

Comparative gene finding in chicken indicates that we are closing in on the set of multi-exonic widely expressed human genes.

The recent availability of the chicken genome sequence poses the question of whether there are human protein-coding genes conserved in chicken that are currently not included in the human gene catalog. Here, we show, using comparative gene finding followed by experimental verification of exon pairs by RT-PCR, that the addition to the multi-exonic subset of this catalog could be as little as 0.2%, suggesting that we may be closing in on the human gene set. Our protocol, however, has two shortcomings: (i) the bioinformatic screening of the predicted genes, applied to filter out false positives, cannot handle intronless genes; and (ii) the experimental verification could fail to identify expression at a specific developmental time. This highlights the importance of developing methods that could provide a reliable estimate of the number of these two types of genes.

Involvement of microRNAs in the cytotoxic effects exerted by proinflammatory cytokines on pancreatic beta-cells.

OBJECTIVE: Pancreatic beta-cells exposed to proinflammatory cytokines display alterations in gene expression resulting in defective insulin secretion and apoptosis. MicroRNAs are small noncoding RNAs emerging as key regulators of gene expression. Here, we evaluated the contribution of microRNAs to cytokine-mediated beta-cell cytotoxicity. RESEARCH DESIGN AND METHODS: We used global microarray profiling and real-time PCR analysis to detect changes in microRNA expression in beta-cells exposed to cytokines and in islets of pre-diabetic NOD mice. We assessed the involvement of the microRNAs affected in cytokine-mediated beta-cell failure by modifying their expression in insulin-secreting MIN6 cells. RESULTS: We found that IL-1beta and TNF-alpha induce the expression of miR-21, miR-34a, and miR-146a both in MIN6 cells and human pancreatic islets. We further show an increase of these microRNAs in islets of NOD mice during development of pre-diabetic insulitis. Blocking miR-21, miR-34a, or miR-146a function using antisense molecules did not restore insulin-promoter activity but prevented the reduction in glucose-induced insulin secretion observed upon IL-1beta exposure. Moreover, anti-miR-34a and anti-miR-146a treatment protected MIN6 cells from cytokine-triggered cell death. CONCLUSIONS: Our data identify miR-21, miR-34a, and miR-146a as novel players in beta-cell failure elicited in vitro and in vivo by proinflammatory cytokines, notably during the development of peri-insulitis that precedes overt diabetes in NOD mice.