Gene expression mapping in neuropathic pain : molecular plasticity in nociceptors and superficial dorsal horn

RESUME : La douleur neuropathique est le résultat d'une lésion ou d'un dysfonctionnement du système nerveux. Les symptômes qui suivent la douleur neuropathique sont sévères et leur traitement inefficace. Une meilleure approche thérapeutique peut être proposée en se basant sur les mécanismes pathologiques de la douleur neuropathique. Lors d'une lésion périphérique une douleur neuropathique peut se développer et affecter le territoire des nerfs lésés mais aussi les territoires adjacents des nerfs non-lésés. Une hyperexcitabilité des neurones apparaît au niveau des ganglions spinaux (DRG) et de la corne dorsale (DH) de la moelle épinière. Le but de ce travail consiste à mettre en évidence les modifications moléculaires associées aux nocicepteurs lésés et non-lésés au niveau des DRG et des laminae I et II de la corne dorsale, là où l'information nociceptive est intégrée. Pour étudier les changements moléculaires liés à la douleur neuropathique nous utilisons le modèle animal d'épargne du nerf sural (spared nerve injury model, SNI) une semaine après la lésion. Pour la sélection du tissu d'intérêt nous avons employé la technique de la microdissection au laser, afin de sélectionner une sous-population spécifique de cellules (notamment les nocicepteurs lésés ou non-lésés) mais également de prélever le tissu correspondant dans les laminae superficielles. Ce travail est couplé à l'analyse à large spectre du transcriptome par puce ADN (microarray). Par ailleurs, nous avons étudié les courants électriques et les propriétés biophysiques des canaux sodiques (Na,,ls) dans les neurones lésés et non-lésés des DRG. Aussi bien dans le système nerveux périphérique, entre les neurones lésés et non-lésés, qu'au niveau central avec les aires recevant les projections des nocicepteurs lésés ou non-lésés, l'analyse du transcriptome montre des différences de profil d'expression. En effet, nous avons constaté des changements transcriptionnels importants dans les nocicepteurs lésés (1561 gènes, > 1.5x et pairwise comparaison > 77%) ainsi que dans les laminae correspondantes (618 gènes), alors que ces modifications transcriptionelles sont mineures au niveau des nocicepteurs non-lésés (60 gènes), mais important dans leurs laminae de projection (459 gènes). Au niveau des nocicepteurs, en utilisant la classification par groupes fonctionnels (Gene Ontology), nous avons observé que plusieurs processus biologiques sont modifiés. Ainsi des fonctions telles que la traduction des signaux cellulaires, l'organisation du cytosquelette ainsi que les mécanismes de réponse au stress sont affectés. Par contre dans les neurones non-lésés seuls les processus biologiques liés au métabolisme et au développement sont modifiés. Au niveau de la corne dorsale de la moelle, nous avons observé des modifications importantes des processus immuno-inflammatoires dans l'aire affectée par les nerfs lésés et des changements associés à l'organisation et la transmission synaptique au niveau de l'aire des nerfs non-lésés. L'analyse approfondie des canaux sodiques a démontré plusieurs changements d'expression, principalement dans les neurones lésés. Les analyses fonctionnelles n'indiquent aucune différence entre les densités de courant tétrodotoxine-sensible (TTX-S) dans les neurones lésés et non-lésés même si les niveaux d'expression des ARNm des sous-unités TTX-S sont modifiés dans les neurones lésés. L'inactivation basale dépendante du voltage des canaux tétrodotoxine-insensible (TTX-R) est déplacée vers des potentiels positifs dans les cellules lésées et non-lésées. En revanche la vitesse de récupération des courants TTX-S et TTX-R après inactivation est accélérée dans les neurones lésés. Ces changements pourraient être à l'origine de l'altération de l'activité électrique des neurones sensoriels dans le contexte des douleurs neuropathiques. En résumé, ces résultats suggèrent l'existence de mécanismes différenciés affectant les neurones lésés et les neurones adjacents non-lésés lors de la mise en place la douleur neuropathique. De plus, les changements centraux au niveau de la moelle épinière qui surviennent après lésion sont probablement intégrés différemment selon la perception de signaux des neurones périphériques lésés ou non-lésés. En conclusion, ces modulations complexes et distinctes sont probablement des acteurs essentiels impliqués dans la genèse et la persistance des douleurs neuropathiques. ABSTRACT : Neuropathic pain (NP) results from damage or dysfunction of the peripheral or central nervous system. Symptoms associated with NP are severe and difficult to treat. Targeting NP mechanisms and their translation into symptoms may offer a better therapeutic approach.Hyperexcitability of the peripheral and central nervous system occurs in the dorsal root ganglia (DRG) and the dorsal horn (DH) of the spinal cord. We aimed to identify transcriptional variations in injured and in adjacent non-injured nociceptors as well as in corresponding laminae I and II of DH receiving their inputs.We investigated changes one week after the injury induced by the spared nerve injury model of NP. We employed the laser capture microdissection (LCM) for the procurement of specific cell-types (enrichment in nociceptors of injured/non-injured neurons) and laminae in combination with transcriptional analysis by microarray. In addition, we studied functionál properties and currents of sodium channels (Nav1s) in injured and neighboring non-injured DRG neurons.Microarray analysis at the periphery between injured and non-injured DRG neurons and centrally between the area of central projections from injured and non-injured neurons show significant and differential expression patterns. We reported changes in injured nociceptors (1561 genes, > 1.5 fold, >77% pairwise comparison) and in corresponding DH laminae (618 genes), while less modifications occurred in non-injured nociceptors (60 genes) and in corresponding DH laminae (459 genes). At the periphery, we observed by Gene Ontology the involvement of multiple biological processes in injured neurons such as signal transduction, cytoskeleton organization or stress responses. On contrast, functional overrepresentations in non-injured neurons were noted only in metabolic or developmentally related mechanisms. At the level of superficial laminae of the dorsal horn, we reported changes of immune and inflammatory processes in injured-related DH and changes associated with synaptic organization and transmission in DH corresponding to non-injured neurons. Further transcriptional analysis of Nav1s indicated several changes in injured neurons. Functional analyses of Nav1s have established no difference in tetrodotoxin-sensitive (TTX-S) current densities in both injured and non-injured neurons, despite changes in TTX-S Nav1s subunit mRNA levels. The tetrodotoxin-resistant (TTX-R) voltage dependence of steady state inactivation was shifted to more positive potentials in both injured and non-injured neurons, and the rate of recovery from inactivation of TTX-S and TTX-R currents was accelerated in injured neurons. These changes may lead to alterations in neuronal electrogenesis. Taken together, these findings suggest different mechanisms occurring in the injured neurons and the adjacent non-injured ones. Moreover, central changes after injury are probably driven in a different manner if they receive inputs from injured or non-injured neurons. Together, these distinct and complex modulations may contribute to NP.

A modular approach for integrative analysis of large-scale gene-expression and drug-response data

High-throughput technologies are now used to generate more than one type of data from the same biological samples. To properly integrate such data, we propose using co-modules, which describe coherent patterns across paired data sets, and conceive several modular methods for their identification. We first test these methods using in silico data, demonstrating that the integrative scheme of our Ping-Pong Algorithm uncovers drug-gene associations more accurately when considering noisy or complex data. Second, we provide an extensive comparative study using the gene-expression and drug-response data from the NCI-60 cell lines. Using information from the DrugBank and the Connectivity Map databases we show that the Ping-Pong Algorithm predicts drug-gene associations significantly better than other methods. Co-modules provide insights into possible mechanisms of action for a wide range of drugs and suggest new targets for therapy

Comparative study and identification of potent eukaryotic transcriptional repressors in gene switch systems.

In mammalian cells, proper gene regulation is achieved by the complex interplay of transcription factors that activate or repress gene expression by binding to the regulatory regions of target promoters. While transcriptional activators have been extensively characterised and classified into functional groups, relatively little is known about the comparative strength and cell type-specificity of transcriptional repressors. Here, we have compared the ability of a series of eukaryotic repression domains to silence basal and activated transcription. A series of the most potent repression domains was further tested in the context of a gene therapy gene-switch system in various cell types. The results indicate that the analysed repression domains exert varying silencing activities in different promoter contexts. Furthermore, their potential for gene silencing varies also depending on the cellular context. When multimerised within one chimeric repressor protein, particular combinations of repressor domains were found to display synergistic repressing effects and efficient repression in a panel of cell lines. This approach thus allowed the identification of transcriptional repressors that are both potent and versatile in terms of cellular specificity as a basis for gene switch systems.

Post-genomic update on a classical candidate gene for coronary artery disease: ESR1.

BACKGROUND: After age, sex is the most important risk factor for coronary artery disease (CAD). The mechanism through which women are protected from CAD is still largely unknown, but the observed sex difference suggests the involvement of the reproductive steroid hormone signaling system. Genetic association studies of the gene-encoding Estrogen Receptor α (ESR1) have shown conflicting results, although only a limited range of variation in the gene has been investigated. METHODS AND RESULTS: We exploited information made available by advanced new methods and resources in complex disease genetics to revisit the question of ESR1's role in risk of CAD. We performed a meta-analysis of 14 genome-wide association studies (CARDIoGRAM discovery analysis, N=≈87,000) to search for population-wide and sex-specific associations between CAD risk and common genetic variants throughout the coding, noncoding, and flanking regions of ESR1. In addition to samples from the MIGen (N=≈6000), WTCCC (N=≈7400), and Framingham (N=≈3700) studies, we extended this search to a larger number of common and uncommon variants by imputation into a panel of haplotypes constructed using data from the 1000 Genomes Project. Despite the widespread expression of ERα in vascular tissues, we found no evidence for involvement of common or low-frequency genetic variation throughout the ESR1 gene in modifying risk of CAD, either in the general population or as a function of sex. CONCLUSIONS: We suggest that future research on the genetic basis of sex-related differences in CAD risk should initially prioritize other genes in the reproductive steroid hormone biosynthesis system.

A mucin like gene different from the previously reported members of the mucin like gene families is transcribed in Trypanosoma cruzi but not in Trypanosoma rangeli

Trypanosoma cruzi expresses mucin like glycoproteins encoded by a complex multigene family. In this work, we report the transcription in T. cruzi but not in T. rangeli of a mucin type gene automatically annotated by the T. cruzi genome project. The gene showed no nucleotide similarities with the previously reported T. cruzi mucin like genes, although the computational analysis of the deduced protein showed that it has the characteristic features of mucins: a signal peptide sequence, O-glycosylation sites, and glycosylphosphatidylinositol (GPI) anchor sequence. The presence in this gene of N- terminal and C- terminal coding sequences common to other annotated mucin like genes suggests the existence of a new mucin like gene family.

Unraveling the complex network of cuticular structure and function.

A hydrophobic cuticle is deposited at the outermost extracellular matrix of the epidermis in primary tissues of terrestrial plants. Besides forming a protective shield against the environment, the cuticle is potentially involved in several developmental processes during plant growth. A high degree of variation in cuticle composition and structure exists between different plant species and tissues. Lots of progress has been made recently in understanding the different steps of biosynthesis, transport, and deposition of cuticular components. However, the molecular mechanisms that underlie cuticular function remain largely elusive.

Characterization of the gene expression related to the process of DNA damage tolerance in Schistosoma mansoni

In the course of its complex life cycle, the parasite Schistosoma mansoni need to adapt to distinct environments, and consequently is exposed to various DNA damaging agents. The Schistosoma genome sequencing initiative has uncovered sequences from genes and transcripts related to the process of DNA damage tolerance as the enzymes UBC13, MMS2, and RAD6. In the present work, we evaluate the importance of this process in different stages of the life cycle of this parasite. The importance is evidenced by expression and phylogenetic profiles, which show the conservation of this pathway from protozoa to mammalians on evolution.

Functional and transcriptional induction of aquaporin-1 gene by hypoxia; analysis of promoter and role of Hif-1α.

Aquaporin-1 (AQP1) is a water channel that is highly expressed in tissues with rapid O(2) transport. It has been reported that this protein contributes to gas permeation (CO(2), NO and O(2)) through the plasma membrane. We show that hypoxia increases Aqp1 mRNA and protein levels in tissues, namely mouse brain and lung, and in cultured cells, the 9L glioma cell line. Stopped-flow light-scattering experiments confirmed an increase in the water permeability of 9L cells exposed to hypoxia, supporting the view that hypoxic Aqp1 up-regulation has a functional role. To investigate the molecular mechanisms underlying this regulatory process, transcriptional regulation was studied by transient transfections of mouse endothelial cells with a 1297 bp 5' proximal Aqp1 promoter-luciferase construct. Incubation in hypoxia produced a dose- and time-dependent induction of luciferase activity that was also obtained after treatments with hypoxia mimetics (DMOG and CoCl(2)) and by overexpressing stabilized mutated forms of HIF-1α. Single mutations or full deletions of the three putative HIF binding domains present in the Aqp1 promoter partially reduced its responsiveness to hypoxia, and transfection with Hif-1α siRNA decreased the in vitro hypoxia induction of Aqp1 mRNA and protein levels. Our results indicate that HIF-1α participates in the hypoxic induction of AQP1. However, we also demonstrate that the activation of Aqp1 promoter by hypoxia is complex and multifactorial and suggest that besides HIF-1α other transcription factors might contribute to this regulatory process. These data provide a conceptual framework to support future research on the involvement of AQP1 in a range of pathophysiological conditions, including edema, tumor growth, and respiratory diseases.

Identification of Mycobacterium tuberculosis complex based on amplification and sequencing of the oxyR pseudogene from stored Ziehl-Neelsen-stained sputum smears in Brazil

A cross-sectional analysis of stored Ziehl-Neelsen (ZN)-stained sputum smear slides (SSS) obtained from two public tuberculosis referral laboratories located in Juiz de Fora, Minas Gerais, was carried out to distinguish Mycobacterium bovis from other members of the Mycobacterium tuberculosis complex (MTC). A two-step approach was used to distinguish M. bovis from other members of MTC: (i) oxyR pseudogene amplification to detect MTC and, subsequently, (ii) allele-specific sequencing based on the polymorphism at position 285 of this gene. The oxyR pseudogene was successfully amplified in 100 of 177 (56.5%) SSS available from 99 individuals. No molecular profile of M. bovis was found. Multivariate analysis indicated that acid-fast bacilli (AFB) results and the source laboratory were associated (p < 0.05) with oxyR pseudogene amplification. SSS that were AFB++ SSS showed more oxyR pseudogene amplification than those with AFB0, possibly due to the amount of DNA. One of the two source laboratories presented a greater chance of oxyR pseudogene amplification, suggesting that differences in sputum conservation between laboratories could have influenced the preservation of DNA. This study provides evidence that stored ZN-SSS can be used for the molecular detection of MTC.

Methods for the ex vivo characterization of human CD8+ T subsets based on gene expression and replicative history analysis.

The generation of an antigen-specific T-lymphocyte response is a complex multi-step process. Upon T-cell receptor-mediated recognition of antigen presented by activated dendritic cells, naive T-lymphocytes enter a program of proliferation and differentiation, during the course of which they acquire effector functions and may ultimately become memory T-cells. A major goal of modern immunology is to precisely identify and characterize effector and memory T-cell subpopulations that may be most efficient in disease protection. Sensitive methods are required to address these questions in exceedingly low numbers of antigen-specific lymphocytes recovered from clinical samples, and not manipulated in vitro. We have developed new techniques to dissect immune responses against viral or tumor antigens. These allow the isolation of various subsets of antigen-specific T-cells (with major histocompatibility complex [MHC]-peptide multimers and five-color FACS sorting) and the monitoring of gene expression in individual cells (by five-cell reverse transcription-polymerase chain reaction [RT-PCR]). We can also follow their proliferative life history by flow-fluorescence in situ hybridization (FISH) analysis of average telomere length. Recently, using these tools, we have identified subpopulations of CD8+ T-lymphocytes with distinct proliferative history and partial effector-like properties. Our data suggest that these subsets descend from recently activated T-cells and are committed to become differentiated effector T-lymphocytes.

How to keep the brain awake? The complex molecular pharmacogenetics of wake promotion.

Wake-promoting drugs are widely used to treat excessive daytime sleepiness. The neuronal pathways involved in wake promotion are multiple and often not well characterized. We tested d-amphetamine, modafinil, and YKP10A, a novel wake-promoting compound, in three inbred strains of mice. The wake duration induced by YKP10A and d-amphetamine depended similarly on genotype, whereas opposite strain differences were observed after modafinil. Electroencephalogram (EEG) analysis during drug-induced wakefulness revealed a transient approximately 2 Hz slowing of theta oscillations and an increase in beta-2 (20-35 Hz) activity only after YKP10A. Gamma activity (35-60 Hz) was induced by all drugs in a drug- and genotype-dependent manner. Brain transcriptome and clustering analyses indicated that the three drugs have both common and specific molecular signatures. The correlation between specific EEG and gene-expression signatures suggests that the neuronal pathways activated to stay awake vary among drugs and genetic background.

Molecular divergence in the timeless and cpr genes among three sympatric cryptic species of the Anopheles triannulatus complex

Anopheles triannulatus s.l. is a malaria vector with a wide geographic distribution, ranging from Argentina-Nicaragua and Trinidad. Here we analysed sequences of two genes, timeless and cpr, to assess the genetic variability and divergence among three sympatric cryptic species of this complex from Salobra, central-western Brazil. The timeless gene sequences did not conclusively differentiate Anopheles halophylus and An. triannulatus species "C". However, a partial separation has been observed between these species and An. triannulatus s.s. Importantly, the analysis of the cpr gene sequences revealed fixed differences, no shared polymorphisms and considerable genetic differentiation among the three species of the An. triannulatus complex. The results confirm that An. triannulatus s.s., An. halophylus and An. triannulatus species C are distinct taxa, with the latter two likely representing a more recent speciation event.

The classification of esterases: an important gene family involved in insecticide resistance - A review

The use of chemical insecticides continues to play a major role in the control of disease vector populations, which is leading to the global dissemination of insecticide resistance. A greater capacity to detoxify insecticides, due to an increase in the expression or activity of three major enzyme families, also known as metabolic resistance, is one major resistance mechanisms. The esterase family of enzymes hydrolyse ester bonds, which are present in a wide range of insecticides; therefore, these enzymes may be involved in resistance to the main chemicals employed in control programs. Historically, insecticide resistance has driven research on insect esterases and schemes for their classification. Currently, several different nomenclatures are used to describe the esterases of distinct species and a universal standard classification does not exist. The esterase gene family appears to be rapidly evolving and each insect species has a unique complement of detoxification genes with only a few orthologues across species. The examples listed in this review cover different aspects of their biochemical nature. However, they do not appear to contribute to reliably distinguish among the different resistance mechanisms. Presently, the phylogenetic criterion appears to be the best one for esterase classification. Joint genomic, biochemical and microarray studies will help unravel the classification of this complex gene family.

Amplicon DNA melting analysis for the simultaneous detection of Brucella spp and Mycobacterium tuberculosis complex. Potential use in rapid differential diagnosis between extrapulmonary tuberculosis and focal complications of brucellosis.

Some sites of extrapulmonary tuberculosis and focal complications of brucellosis are very difficult to differentiate clinically, radiologically, and even histopathologically. Conventional microbiological methods for the diagnosis of extrapulmonary tuberculosis and complicated brucellosis not only lack adequate sensitivity, they are also time consuming, which could lead to an unfavourable prognosis. The aim of this work was to develop a multiplex real-time PCR assay based on SYBR Green I to simultaneously detect Brucella spp and Mycobacterium tuberculosis complex and evaluate the efficacy of the technique with different candidate genes. The IS711, bcsp31 and omp2a genes were used for the identification of Brucella spp and the IS6110, senX3-regX3 and cfp31 genes were targeted for the detection of the M. tuberculosis complex. As a result of the different combinations of primers, nine different reactions were evaluated. A test was defined as positive only when the gene combinations were capable of co-amplifying both pathogens in a single reaction tube and showed distinguishable melting temperatures for each microorganism. According to the melting analysis, only three combinations of amplicons (senX3-regX3+bcsp31, senX3-regX3+IS711 and IS6110+IS711) were visible. Detection limits of senX3-regX3+bcsp31 and senX3-regX3+IS711 were of 2 and 3 genome equivalents for M. tuberculosis complex and Brucella while for IS6110+IS711 they were of 200 and 300 genome equivalents, respectively. The three assays correctly identified all the samples, showing negative results for the control patients. The presence of multicopy elements and GC content were the components most influencing the efficiency of the test; this should be taken into account when designing a multiplex-based SYBR Green I assay. In conclusion, multiplex real time PCR assays based on the targets senX3-regX3+bcsp31 and senX3-regX3+IS711 using SYBR Green I are highly sensitive and reproducible. This may therefore be a practical approach for the rapid differential diagnosis between extrapulmonary tuberculosis and complicated brucellosis.

Molecular analysis of an odorant-binding protein gene in two sympatric species of Lutzomyia longipalpis s.l.

Lutzomyia longipalpis s.l. is the main vector of American visceral leishmaniasis (AVL) and occurs as a species complex. DNA samples from two Brazilian sympatric species that differ in pheromone and courtship song production were used to analyse molecular polymorphisms in an odorant-binding protein ( obp29 ) gene. OBPs are proteins related to olfaction and are involved in activities fundamental to survival, such as foraging, mating and choice of oviposition site. In this study, the marker obp29 was found to be highly polymorphic in Lu. longipalpis s.l. , with no fixed differences observed between the two species. A pairwise fixation index test indicated a moderate level of genetic differentiation between the samples analysed.