Modeling gene expression regulatory networks with the sparse vector autoregressive model

Abstract Background To understand the molecular mechanisms underlying important biological processes, a detailed description of the gene products networks involved is required. In order to define and understand such molecular networks, some statistical methods are proposed in the literature to estimate gene regulatory networks from time-series microarray data. However, several problems still need to be overcome. Firstly, information flow need to be inferred, in addition to the correlation between genes. Secondly, we usually try to identify large networks from a large number of genes (parameters) originating from a smaller number of microarray experiments (samples). Due to this situation, which is rather frequent in Bioinformatics, it is difficult to perform statistical tests using methods that model large gene-gene networks. In addition, most of the models are based on dimension reduction using clustering techniques, therefore, the resulting network is not a gene-gene network but a module-module network. Here, we present the Sparse Vector Autoregressive model as a solution to these problems. Results We have applied the Sparse Vector Autoregressive model to estimate gene regulatory networks based on gene expression profiles obtained from time-series microarray experiments. Through extensive simulations, by applying the SVAR method to artificial regulatory networks, we show that SVAR can infer true positive edges even under conditions in which the number of samples is smaller than the number of genes. Moreover, it is possible to control for false positives, a significant advantage when compared to other methods described in the literature, which are based on ranks or score functions. By applying SVAR to actual HeLa cell cycle gene expression data, we were able to identify well known transcription factor targets. Conclusion The proposed SVAR method is able to model gene regulatory networks in frequent situations in which the number of samples is lower than the number of genes, making it possible to naturally infer partial Granger causalities without any a priori information. In addition, we present a statistical test to control the false discovery rate, which was not previously possible using other gene regulatory network models.

Characterization and optimization of ArtinM lectin expression in Escherichia coli

Background ArtinM is a D-mannose-specific lectin from Artocarpus integrifolia seeds that induces neutrophil migration and activation, degranulation of mast cells, acceleration of wound healing, induction of interleukin-12 production by macrophages and dendritic cells, and protective T helper 1 immune response against Leishmania major, Leishmania amazonensis and Paracoccidioides brasiliensis infections. Considering the important biological properties of ArtinM and its therapeutic applicability, this study was designed to produce high-level expression of active recombinant ArtinM (rArtinM) in Escherichia coli system. Results The ArtinM coding region was inserted in pET29a(+) vector and expressed in E. coli BL21(DE3)-Codon Plus-RP. The conditions for overexpression of soluble ArtinM were optimized testing different parameters: temperatures (20, 25, 30 or 37°C) and shaking speeds (130, 200 or 220 rpm) during induction, concentrations of the induction agent IPTG (0.01-4 mM) and periods of induction (1-19 h). BL21-CodonPlus(DE3)-RP cells induced under the optimized conditions (incubation at 20°C, at a shaking speed of 130 rpm, induction with 0.4 mM IPTG for 19 h) resulted in the accumulation of large amounts of soluble rArtinM. The culture provided 22.4 mg/L of rArtinM, which activity was determined by its one-step purification through affinity chromatography on immobilized D-mannose and glycoarray analysis. Gel filtration showed that rArtinM is monomeric, contrasting with the tetrameric form of the plant native protein (jArtinM). The analysis of intact rArtinM by mass spectrometry revealed a 16,099.5 Da molecular mass, and the peptide mass fingerprint and esi-cid-ms/ms of amino acid sequences of peptides from a tryptic digest covered 41% of the total ArtinM amino acid sequence. In addition, circular dichroism and fluorescence spectroscopy of rArtinM indicated that its global fold comprises β-sheet structure. Conclusions Overall, the optimized process to express rArtinM in E. coli provided high amounts of soluble, correctly folded and active recombinant protein, compatible with large scale production of the lectin.

Expression, purification and molecular analysis of the human ZNF706 protein

Background: The ZNF706 gene encodes a protein that belongs to the zinc finger family of proteins and was found to be highly expressed in laryngeal cancer, making the structure and function of ZNF706 worthy of investigation. In this study, we expressed and purified recombinant human ZNF706 that was suitable for structural analysis in Escherichia coli BL21(DH3). Findings: ZNF706 mRNA was extracted from a larynx tissue sample, and cDNA was ligated into a cloning vector using the TOPO method. ZNF706 protein was expressed according to the E. coli expression system procedures and was purified using a nickel-affinity column. The structural qualities of recombinant ZNF706 and quantification alpha, beta sheet, and other structures were obtained by spectroscopy of circular dichroism. ZNF706's structural modeling showed that it is composed of α-helices (28.3%), β-strands (19.4%), and turns (20.9%), in agreement with the spectral data from the dichroism analysis. Conclusions: We used circular dichroism and molecular modeling to examine the structure of ZNF706. The results suggest that human recombinant ZNF706 keeps its secondary structures and is appropriate for functional and structural studies. The method of expressing ZNF706 protein used in this study can be used to direct various functional and structural studies that will contribute to the understanding of its function as well as its relationship with other biological molecules and its putative role in carcinogenesis.

Natural blood feeding and temperature shift modulate the global transcriptional profile of Rickettsia rickettsii infecting its tick vector

Rickettsia rickettsii is an obligate intracellular tick-borne bacterium that causes Rocky Mountain Spotted Fever (RMSF), the most lethal spotted fever rickettsiosis. When an infected starving tick begins blood feeding from a vertebrate host, R. rickettsii is exposed to a temperature elevation and to components in the blood meal. These two environmental stimuli have been previously associated with the reactivation of rickettsial virulence in ticks, but the factors responsible for this phenotype conversion have not been completely elucidated. Using customized oligonucleotide microarrays and high-throughput microfluidic qRT-PCR, we analyzed the effects of a 10 degrees C temperature elevation and of a blood meal on the transcriptional profile of R. rickettsii infecting the tick Amblyomma aureolatum. This is the first study of the transcriptome of a bacterium in the genus Rickettsia infecting a natural tick vector. Although both stimuli significantly increased bacterial load, blood feeding had a greater effect, modulating five-fold more genes than the temperature upshift. Certain components of the Type IV Secretion System (T4SS) were up-regulated by blood feeding. This suggests that this important bacterial transport system may be utilized to secrete effectors during the tick vector's blood meal. Blood feeding also up-regulated the expression of antioxidant enzymes, which might correspond to an attempt by R. rickettsii to protect itself against the deleterious effects of free radicals produced by fed ticks. The modulated genes identified in this study, including those encoding hypothetical proteins, require further functional analysis and may have potential as future targets for vaccine development.

New approaches to open problems in gene expression microarray data

In the past decade, the advent of efficient genome sequencing tools and high-throughput experimental biotechnology has lead to enormous progress in the life science. Among the most important innovations is the microarray tecnology. It allows to quantify the expression for thousands of genes simultaneously by measurin the hybridization from a tissue of interest to probes on a small glass or plastic slide. The characteristics of these data include a fair amount of random noise, a predictor dimension in the thousand, and a sample noise in the dozens. One of the most exciting areas to which microarray technology has been applied is the challenge of deciphering complex disease such as cancer. In these studies, samples are taken from two or more groups of individuals with heterogeneous phenotypes, pathologies, or clinical outcomes. these samples are hybridized to microarrays in an effort to find a small number of genes which are strongly correlated with the group of individuals. Eventhough today methods to analyse the data are welle developed and close to reach a standard organization (through the effort of preposed International project like Microarray Gene Expression Data -MGED- Society [1]) it is not unfrequant to stumble in a clinician's question that do not have a compelling statistical method that could permit to answer it.The contribution of this dissertation in deciphering disease regards the development of new approaches aiming at handle open problems posed by clinicians in handle specific experimental designs. In Chapter 1 starting from a biological necessary introduction, we revise the microarray tecnologies and all the important steps that involve an experiment from the production of the array, to the quality controls ending with preprocessing steps that will be used into the data analysis in the rest of the dissertation. While in Chapter 2 a critical review of standard analysis methods are provided stressing most of problems that In Chapter 3 is introduced a method to adress the issue of unbalanced design of miacroarray experiments. In microarray experiments, experimental design is a crucial starting-point for obtaining reasonable results. In a two-class problem, an equal or similar number of samples it should be collected between the two classes. However in some cases, e.g. rare pathologies, the approach to be taken is less evident. We propose to address this issue by applying a modified version of SAM [2]. MultiSAM consists in a reiterated application of a SAM analysis, comparing the less populated class (LPC) with 1,000 random samplings of the same size from the more populated class (MPC) A list of the differentially expressed genes is generated for each SAM application. After 1,000 reiterations, each single probe given a "score" ranging from 0 to 1,000 based on its recurrence in the 1,000 lists as differentially expressed. The performance of MultiSAM was compared to the performance of SAM and LIMMA [3] over two simulated data sets via beta and exponential distribution. The results of all three algorithms over low- noise data sets seems acceptable However, on a real unbalanced two-channel data set reagardin Chronic Lymphocitic Leukemia, LIMMA finds no significant probe, SAM finds 23 significantly changed probes but cannot separate the two classes, while MultiSAM finds 122 probes with score >300 and separates the data into two clusters by hierarchical clustering. We also report extra-assay validation in terms of differentially expressed genes Although standard algorithms perform well over low-noise simulated data sets, multi-SAM seems to be the only one able to reveal subtle differences in gene expression profiles on real unbalanced data. In Chapter 4 a method to adress similarities evaluation in a three-class prblem by means of Relevance Vector Machine [4] is described. In fact, looking at microarray data in a prognostic and diagnostic clinical framework, not only differences could have a crucial role. In some cases similarities can give useful and, sometimes even more, important information. The goal, given three classes, could be to establish, with a certain level of confidence, if the third one is similar to the first or the second one. In this work we show that Relevance Vector Machine (RVM) [2] could be a possible solutions to the limitation of standard supervised classification. In fact, RVM offers many advantages compared, for example, with his well-known precursor (Support Vector Machine - SVM [3]). Among these advantages, the estimate of posterior probability of class membership represents a key feature to address the similarity issue. This is a highly important, but often overlooked, option of any practical pattern recognition system. We focused on Tumor-Grade-three-class problem, so we have 67 samples of grade I (G1), 54 samples of grade 3 (G3) and 100 samples of grade 2 (G2). The goal is to find a model able to separate G1 from G3, then evaluate the third class G2 as test-set to obtain the probability for samples of G2 to be member of class G1 or class G3. The analysis showed that breast cancer samples of grade II have a molecular profile more similar to breast cancer samples of grade I. Looking at the literature this result have been guessed, but no measure of significance was gived before.

Recombinant expression of molluscan hemocyanin (KLH) substructures in a prokaryotic system: E. coli

The recombinant expression of 19 different substructures of KLH in the prokaryotic sys-tem E. coli has been successfully achieved: each one of the eight single FUs a to h of both isoforms, KLH1 and KLH2, two substructures consisting of two consecutive FUs (KLH1-bc and KLH1-gh) as well as a cDNA encompassing KLH1-abc. All recombinant proteins, fused to an N-terminal 6xHis tag, have successfully been detected by immuno precipitation using monoclonal α-His-antibodies and polyclonal α-KLH1- and α-KLH2-antibodies. One exception remained: SP-KLH2-a, which was not detected by the α-His-antibodies. This allows speculations as to whether the coexpressed signal peptide can lead, at one hand, to the secretion of the recombinant protein, and on the other to the simultaneous cut-off of the leader peptide, which results in the splitting off of even more N-terminal 6xHis tag, leading to failed recognition by the appropriate antibodies. The comparison of native KLH with recombinantly expressed prokaryotic (E. coli) and eukaryotic (Sf9 insect cells) KLH was done using FU-1h. The weak detection by the polyclonal α-KLH1-antibodies of both recombinantly expressed proteins showed that the native protein was the best recognized. For the prokaryotic one, both the denaturation applied for solubilisation of the bacterial inclusion bodies and the inability of bacterial cells to add N-linked glycosylation, are the reason for the poor hybridization. In contrast, KLH1-h expressed in eukaryotic insect cells is likely to be glycosylated. The incubation with the α-KLH1-antibodies resulting in the same weak detection, however, revealed that the linked carbohydrate side chains are not those expected. The establishment of SOE-PCR, together with further improvement, has enabled the generation of a clone encompassing the complete subunit KLH1-abcdefgh. The se-quence analysis compared to the original KLH1 sequence showed, however, that the resulting recombinant protein is defective in two histidines, required for the copper bind-ing sites in FU-1b and FU-1d and in three disulfide bridges (FU-1a, FU-1b and FU 1g). This is due to polymerase-related nucleotide exchanges, resulting in a changed amino acid sequence. Nevertheless, all eight potential N-glycosylation sites are present, leading to the speculation that the recombinant protein can in theory be fully glycosylated, which is the most important aspect for the clinical applicability of recombinant KLH as an im-munotherapeutic agent. The improvement of this method elaborated during the present work indicates bright prospects for the future generation of a correct cDNA sequence encoding for the complete KLH2 subunit.

Gene expression analysis in ‘Candidatus Phytoplasma mali’-resistant and -susceptible Malus genotypes

Apple proliferation (AP) disease is the most important graft-transmissible and vector-borne disease of apple in Europe. ‘Candidatus Phytoplasma mali’ (Ca. P. mali) is the causal agent of AP. Apple (Malus x domestica) and other Malus species are the only known woody hosts. In European apple orchards, the cultivars are mainly grafted on one rootstock, M. x domestica cv. M9. M9 like all other M. x domestica cultivars is susceptible to ‘Ca. P. mali’. Resistance to AP was found in the wild genotype Malus sieboldii (MS) and in MS-derived hybrids but they were characterised by poor agronomic value. The breeding of a new rootstock carrying the resistant and the agronomic traits was the major aim of a project of which this work is a part. The objective was to shed light into the unknown resistance mechanism. The plant-phytoplasma interaction was studied by analysing differences between the ‘Ca. P. mali’-resistant and -susceptible genotypes related to constitutively expressed genes or to induced genes during infection. The cDNA-Amplified Fragment Length Polymorphism (cDNA-AFLP) technique was employed in both approaches. Differences related to constitutively expressed genes were identified between two ‘Ca. P. mali’-resistant hybrid genotypes (4551 and H0909) and the ‘Ca. P. mali’-susceptible M9. 232 cDNA-AFLP bands present in the two resistant genotypes but absent in the susceptible one were isolated but several different products associated to each band were found. Therefore, two different macroarray hybridisation experiments were performed with the cDNA-AFLP fragments yielding 40 sequences encoding for genes of unknown function or a wide array of functions including plant defence. In the second approach, individuation and analysis of the induced genes was carried out exploiting an in vitro system in which healthy and ‘Ca. P. mali’-infected micropropagated plants were maintained under controlled conditions. Infection trials using in vitro grafting of ‘Ca. P. mali’ showed that the resistance phenotype could be reproduced in this system. In addition, ex vitro plants were generated as an independent control of the genes differentially expressed in the in vitro plants. The cDNA-AFLP analysis in in vitro plants yielded 63 bands characterised by over-expression in the infected state of both the H0909 and MS genotypes. The major part (37 %) of the associated sequences showed homology with products of unknown function. The other genes were involved in plant defence, energy transport/oxidative stress response, protein metabolism and cellular growth. Real-time qPCR analysis was employed to validate the differential expression of the genes individuated in the cDNA-AFLP analysis. Since no internal controls were available for the study of the gene expression in Malus, an analysis on housekeeping genes was performed. The most stably expressed genes were the elongation factor-1 α (EF1) and the eukaryotic translation initiation factor 4-A (eIF4A). Twelve out of 20 genes investigated through qPCR were significantly differentially expressed in at least one genotype either in in vitro plants or in ex vitro plants. Overall, about 20% of the genes confirmed their cDNA-AFLP expression pattern in M. sieboldii or H0909. On the contrary, 30 % of the genes showed down-regulation or were not differentially expressed. For the remaining 50 % of the genes a contrasting behaviour was observed. The qPCR data could be interpreted as follows: the phytoplasma infection unbalance photosynthetic activity and photorespiration down-regulating genes involved in photosynthesis and in the electron transfer chain. As result, and in contrast to M. x domestica genotypes, an up-regulation of genes of the general response against pathogens was found in MS. These genes involved the pathway of H2O2 and the production of secondary metabolites leading to the hypothesis that a response based on the accumulation of H2O2 in MS would be at the base of its resistance. This resembles a phenomenon known as “recovery” where the spontaneous remission of the symptoms is observed in old susceptible plants but occurring in a stochastic way while the resistance in MS is an inducible but stable feature. As additional product of this work three cDNA-AFLP-derived markers were developed which showed independent distribution among the seedlings of two breeding progenies and were associated to a genomic region characteristic of MS. These markers will contribute to the development of molecular markers for the resistance as well as to map the resistance on the Malus genome.

Untersuchungen zur Lokalisation und Expression des respiratorischen Proteins Neuroglobin bei Säugetieren

In der vorliegenden Arbeit wurde Neuroglobin (Ngb), ein evolutiv altes und in Metazoen konserviertes respiratorisches Protein, funktionell untersucht. Mittels des induzierbaren Tet on / Tet off Systems wurde Ngb ektopisch in der murinen Leber und im Gehirn überexprimiert. Die Transkriptome von Leber und Gehirnregionen Ngb-transgener Mäuse wurden mittels Microarrays und RNA-Seq im Vergleich zum Wildtyp analysiert, um Auswirkungen der Ngb-Überexpression zu ermitteln. Die Transkriptom-Analyse in Leber und Gehirn zeigte eine nur geringe Anzahl differenziell regulierter Gene und Stoffwechselwege nach Ngb-Überexpression. Ngb transgene Mäuse wurden CCl4-induziertem ROS-Stress ausgesetzt und die Leberfunktion untersucht. Zudem wurden primäre Hepatozyten-Kulturen etabliert und in diesen in vitro die extrinsische Apoptose induziert. Die Stressversuche zeigten: (i) Die Ngb-Überexpression hat keine protektive Wirkung in der Leber in vivo. (ii) In Leberzellen in vitro hingegen verminderte eine Ngb-Überexpression effizient die Aktivierung der apoptotischen Kaskade. Eine protektive Wirkung von Ngb ist vermutlich von betrachtetem Gewebe und dem verwendeten Stressor abhängig und keine generelle, selektierte Funktion des Proteins.rnWeiterhin wurde eine Ngb-KnockOut-Mauslinie mit einem LacZ-KnockIn-Genotyp etabliert. Hierbei zeigten die KO-Mäuse keinen offensichtlichen Phänotyp in ihrer Entwicklung, Fortpflanzung und Retina-Funktion. Unter Verwendung des LacZ-Knockin-Konstrukts konnten kontrovers diskutierte Ngb-Expressionsorte im adulten Mausgehirn (Hippocampus, Cortex und Cerebellum) sowie in Testes experimentell bestätigt werden. Parallel wurden öffentlich verfügbare RNA-Seq Datensätze ausgewertet, um die regionale Ngb-Expression systematisch ohne Antikörper-assoziierte Spezifitätsprobleme zu charakterisieren. Eine basale Ngb-Expression (RPKM: ~1-5) wurde im Hippocampus, Cortex und Cerebellum, sowie in Retina und Testes ermittelt. Eine 20-40fach höhere, starke Expression (RPKM: ~160) wurde im Hypothalamus bzw. im Hirnstamm nachgewiesen. Die „digitale“ Expressionsuntersuchung wurde mittels qRT-PCR und Western Blot bestätigt. Dieses Expressionsprofil von Ngb in der Maus weist auf eine besondere funktionelle Bedeutung von Ngb im Hypothalamus hin. Eine Funktion von Ngb in der Sauerstoffversorgung der Retina und eine generelle Funktion von Ngb in der Protektion von Neuronen sind mit dem beobachteten Expressionsspektrum weniger gut vereinbar.

Targeting neuronal populations by AAV-mediated gene transfer for studying the endocannabinoid system

The cannabinoid type 1 (CB1) receptor is involved in a plethora of physiological functions and heterogeneously expressed on different neuronal populations. Several conditional loss-of-function studies revealed distinct effects of CB1 receptor signaling on glutamatergic and GABAergic neurons, respectively. To gain a comprehensive picture of CB1 receptor-mediated effects, the present study aimed at developing a gain-of-function approach, which complements conditional loss-of-function studies. Therefore, adeno-associated virus (AAV)-mediated gene delivery and Cre-mediated recombination were combined to recreate an innovative method, which ensures region- and cell type-specific transgene expression in the brain. This method was used to overexpress the CB1 receptor in glutamatergic pyramidal neurons of the mouse hippocampus. Enhanced CB1 receptor activity at glutamatergic terminals caused impairment in hippocampus-dependent memory performance. On the other hand, elevated CB1 receptor levels provoked an increased protection against kainic acid-induced seizures and against excitotoxic neuronal cell death. This finding indicates the protective role of CB1 receptor on hippocampal glutamatergic terminals as a molecular stout guard in controlling excessive neuronal network activity. Hence, CB1 receptor on glutamatergic hippocampal neurons may represent a target for novel agents to restrain excitotoxic events and to treat neurodegenerative diseases. Endocannabinoid synthesizing and degrading enzymes tightly regulate endocannabinoid signaling, and thus, represent a promising therapeutic target. To further elucidate the precise function of the 2-AG degrading enzyme monoacylglycerol lipase (MAGL), MAGL was overexpressed specifically in hippocampal pyramidal neurons. This genetic modification resulted in highly increased MAGL activity accompanied by a 50 % decrease in 2-AG levels without affecting the content of arachidonic acid and anandamide. Elevated MAGL protein levels at glutamatergic terminals eliminated depolarization-induced suppression of excitation (DSE), while depolarization-induced suppression of inhibition (DSI) was unchanged. This result indicates that the on-demand availability of the endocannabinoid 2-AG is crucial for short-term plasticity at glutamatergic synapses in the hippocampus. Mice overexpressing MAGL exhibited elevated corticosterone levels under basal conditions and an increase in anxiety-like behavior, but surprisingly, showed no changes in aversive memory formation and in seizure susceptibility. This finding suggests that 2 AG-mediated hippocampal DSE is essential for adapting to aversive situations, but is not required to form aversive memory and to protect against kainic acid-induced seizures. Thus, specific inhibition of MAGL expressed in hippocampal pyramidal neurons may represent a potential treatment strategy for anxiety and stress disorders. Finally, the method of AAV-mediated cell type-specific transgene expression was advanced to allow drug-inducible and reversible transgene expression. Therefore, elements of the tetracycline-controlled gene expression system were incorporated in our “conditional” AAV vector. This approach showed that transgene expression is switched on after drug application and that background activity in the uninduced state was only detectable in scattered cells of the hippocampus. Thus, this AAV vector will proof useful for future research applications and gene therapy approaches.

Spastin in the human and mouse central nervous system with special reference to its expression in the hippocampus of mouse pilocarpine model of status epilepticus and temporal lobe epilepsy

In the present in situ hybridization and immunocytochemical studies in the mouse central nervous system (CNS), a strong expression of spastin mRNA and protein was found in Purkinje cells and dentate nucleus in the cerebellum, in hippocampal principal cells and hilar neurons, in amygdala, substantia nigra, striatum, in the motor nuclei of the cranial nerves and in different layers of the cerebral cortex except piriform and entorhinal cortices where only neurons in layer II were strongly stained. Spastin protein and mRNA were weakly expressed in most of the thalamic nuclei. In selected human brain regions such as the cerebral cortex, cerebellum, hippocampus, amygdala, substania nigra and striatum, similar results were obtained. Electron microscopy showed spastin immunopositive staining in the cytoplasma, dendrites, axon terminals and nucleus. In the mouse pilocarpine model of status epilepticus and subsequent temporal lobe epilepsy, spastin expression disappeared in hilar neurons as early as at 2h during pilocarpine induced status epilepticus, and never recovered. At 7 days and 2 months after pilocarpine induced status epilepticus, spastin expression was down-regulated in granule cells in the dentate gyrus, but induced expression was found in reactive astrocytes. The demonstration of widespread distribution of spastin in functionally different brain regions in the present study may provide neuroanatomical basis to explain why different neurological, psychological disorders and cognitive impairment occur in patients with spastin mutation. Down-regulation or loss of spastin expression in hilar neurons may be related to their degeneration and may therefore initiate epileptogenetic events, leading to temporal lobe epilepsy.

Expression and hypoxic regulation of the endothelin system in endocrine cells of human and rat pancreatic islets

CONTEXT: The success of pancreatic islet transplantation depends largely on the capacity of the islet graft to survive the initial phase immediately after transplantation until revascularization is completed. Endothelin-1 (ET-1) is a strong vasoconstrictor which has been involved in solid organ graft failure but is also known to be a potent mitogenic/anti-apoptotic factor which could also potentially enhance the survival of the transplanted islets. OBJECTIVE: Characterization of the endothelin system with regard to a potential endothelin agonist/antagonist treatment. DESIGN: Regulated expression of the endothelin system in human and rat pancreatic islets and beta-cell lines was assessed by means of immunohistochemistry, competition binding studies, western blot, RT-PCR, real-time PCR and transplant studies. RESULTS: ET-1, ETA- and ETB-receptor immunoreactivity was identified in the endocrine cells of human and rat pancreatic islets. The corresponding mRNA was detectable in rat beta-cell lines and isolated rat and human pancreatic islets. Competition binding studies on rat islets revealed binding sites for both receptor types. ET-1 stimulated the phosphorylation of mitogen-activated protein kinase, which was prevented by ETA- and ETB-receptor antagonists. After exposure to hypoxia equal to post-transplant environment oxygen tension, mRNA levels of ET-1 and ETB-receptor of human islets were robustly induced whereas ETA-receptor mRNA did not show significant changes. Immunostaining signals for ET-1 and ETA-receptor of transplanted rat islets were markedly decreased when compared to native pancreatic sections. CONCLUSIONS: In pancreatic islets, ET-1 and its receptors are differentially expressed by hypoxia and after transplantation. Our results provide the biological basis for the study of the potential use of endothelin agonists/antagonists to improve islet transplantation outcome.