The effect of transforming growth factor β1 gene polymorphisms in ankylosing spondylitis

Objectives. To determine whether genetic polymorphisms in or near the transforming growth factor β1 (TGFB1) locus were associated d with susceptibility to or severity of ankylosing spondylitis (AS). Methods. Five intragenic single-nucleotide polymorphisms (SNP) and three microsatellite markers flanking the TGFB1 locus were genotyped. Seven hundred and sixty-two individuals from 184 multiplex families were genotyped for the microsatellite markers and two of the promoter SNPs. One thousand and two individuals from 212 English and 170 Finnish families with AS were genotyped for all five intragenic SNPs. A structured questionnaire was used to assess the age of symptom onset, disease duration and disease severity scores, including the BASDAI (Bath Ankylosing Spondylitis Disease Activity Index) and BASFI (Bath Ankylosing Spondylitis Functional Index). Results. A weak association was noted between the rare TGFB1 + 1632 T allele and AS in the Finnish population (P = 0.04) and in the combined data set (P = 0.03). No association was noted between any other SNPs or SNP haplotype and AS, even among those families with positive non-parametric linkage scores. The TGFB1 +1632 polymorphism was also associated with a younger age of symptom onset (English population, allele 2 associated with age of onset greater by 4.2 yr, P = 0.05; combined data set, allele 2 associated with age of onset greater by 3.2 yr, P = 0.02). A haplotype of coding region SNPs (TGFB1 +869/ +915+1632 alleles 2/1/2) was associated with age of symptom onset in both the English parent-case trios and the combined data set (English data set, haplotype 2/1/2 associated with age of onset greater by 4.9 yr, P = 0.03; combined data set, haplotype 2/1/2 associated with greater age of onset by 4.2 yr, P = 0.006). Weak linkage with AS susceptibility was noted and the peak LOD score was 1.3 at distance 2 cM centromeric to the TGFB1 gene. No other linkage or association was found between quantitative traits and the markers. Conclusion. This study suggests that the polymorphisms within the TGFB1 gene play at most a small role in AS and that other genes encoded on chromosome 19 are involved in susceptibility to the disease.

Polymorphism in codon 17 of the CTLA-4 gene (+49 A/G) is not associated with susceptibility to rheumatoid arthritis in British Caucasians

The role of the CTLA-4 antigen in the development of autoimmune diseases is well documented, with several autoimmune disorders showing association or linkage with the CTLA-4 locus. Its role in the aetiology of rheumatoid arthritis (RA) however, remains unclear, as the functional studies of the B7-CTLA-4 pathway in mouse models of RA and genetic studies in humans have given contrasting results. We have studied the single nucleotide polymorphism at position +49 (A/G) of the CTLA-4 gene, in a cohort of 421 RA cases and 452 healthy controls from the UK. Despite the high statistical power to detect even a weak susceptibility effect, no significant association was found. We also analysed the distribution of the allele and genotype frequencies with respect to the presence of the shared epitope (a known RA susceptibility factor) and found no statistically significant differences. We conclude that, although the importance of the B7-CTLA-4 interaction in the development of RA can not be excluded, the CTLA-4 gene is unlikely to be a predisposing factor to this disease.

Polymorphisms of the CYP2D6 gene increase susceptibility to ankylosing spondylitis

Ankylosing spondylitis (AS) is a common and highly familial rheumatic disorder. The sibling recurrence risk ratio for the disease is 63 and heritability assessed in twins > 90%. Although MHC genes, including HLA-B27, contribute only 20-50% of the genetic risk for the disease, no non-MHC gene has yet been convincingly demonstrated to influence either susceptibility to the disease or its phenotypic expression. Previous linkage and association studies have suggested the presence of a susceptibility gene for AS close to, or within, the cytochrome P450 2D6 gene (CYP2D6, debrisoquine hydroxylase) located at chromosome 22q13.1. We performed a linkage study of chromosome 22 in 200 families with AS affected sibling-pairs. Association of alleles of the CYP2D6 gene was examined by both case-control and within-family means. For case-control studies, 617 unrelated individuals with AS (361 probands from sibling-pair and parent-case trio families and 256 unrelated non-familial sporadic cases) and 402 healthy ethnically matched controls were employed. For within-family association studies, 361 families including 161 parent-case trios and 200 affected sibling-pair families were employed. Homozygosity for poor metabolizer alleles was found to be associated with AS. Heterozygosity for the most frequent poor metabolizer allele (CYP2D6*4) was not associated with increased susceptibility to AS. Significant within-family association of CYP2D6*4 alleles and AS was demonstrated. Weak linkage was also demonstrated between CYP2D6 and AS. We postulate that altered metabolism of a natural toxin or antigen by the CYP2D6 gene may increase susceptibility to AS.

Analysis of differential gene expression based on Bayesian estimation of variance

Gene expression is arguably the most important indicator of biological function. Thus identifying differentially expressed genes is one of the main aims of high throughout studies that use microarray and RNAseq platforms to study deregulated cellular pathways. There are many tools for analysing differentia gene expression from transciptomic datasets. The major challenge of this topic is to estimate gene expression variance due to the high amount of ‘background noise’ that is generated from biological equipment and the lack of biological replicates. Bayesian inference has been widely used in the bioinformatics field. In this work, we reveal that the prior knowledge employed in the Bayesian framework also helps to improve the accuracy of differential gene expression analysis when using a small number of replicates. We have developed a differential analysis tool that uses Bayesian estimation of the variance of gene expression for use with small numbers of biological replicates. Our method is more consistent when compared to the widely used cyber-t tool that successfully introduced the Bayesian framework to differential analysis. We also provide a user-friendly web based Graphic User Interface for biologists to use with microarray and RNAseq data. Bayesian inference can compensate for the instability of variance caused when using a small number of biological replicates by using pseudo replicates as prior knowledge. We also show that our new strategy to select pseudo replicates will improve the performance of the analysis. - See more at: http://www.eurekaselect.com/node/138761/article#sthash.VeK9xl5k.dpuf

An Explicit Guidance Scheme For A Low-Thrust Terminal Rendezvous

An explicit near-optimal guidance scheme is developed for a terminal rendezvous of a spacecraft with a passive target in circular orbit around the earth. The thrust angle versus time profile for the continuous-thrust, constant-acceleration maneuver is derived, based on the assumption that the components of inertial acceleration due to relative position and velocity are negligible on account of the close proximity between the two spacecraft. The control law is obtained as a ''bilinear tangent law'' and an analytic solution to the state differential equations is obtained by expanding a portion of the integrand as an infinite series in time. A differential corrector method is proposed, to obtain real-time updates to the guidance parameters at regular time intervals. Simulation of the guidance scheme is carried out using the Clohessy-Wiltshire equations of relative motion as well as the inverse-square two-body equations of motion. Results for typical examples are presented.

Synthesis, X-ray structures, and spectroscopic and electrochemical properties of (mu-oxo)bis(mu-carboxylato)diruthenium complexes having six imidazole bases as terminal ligands

Complexes [Ru2O(O2CR)(2)(1-MeIm)(6)](ClO4)(2) (la-c), [Ru2O(O2CR)(2)(ImH)(6)](ClO4)(2) (2a,b), and [Ru2O(O2CR)(2)(4-MeImH)(6)](ClO4)(2) (3a,b) with a (mu-oxo)bis(mu-carboxylato)diruthenium(III) core have been prepared by reacting Ru2Cl(O2CR)(4) with the corresponding imidazole base, viz. 1-methylimidazole (1-MeIm), imidazole (ImH), and 4-methylimidazole (4-MeImH) in methanol, followed by treatment with NaClO4 in water (R: Me, a; C6H4-p-OMe, b; C6H4-p-Me, c). Diruthenium(III,IV) complexes [Ru2O(O2CR)(2)(1-MeIm)(6)](ClO4)(3) (R: Me, 4a; C6H4-p-OMe, 4b; C6H4-p-Me, 4c) have been prepared by one-electron oxidation of 1 in MeCN with K2S2O8 in water. Complexes la, 2a . 3H(2)O, and 4a . 1.5H(2)O have been structurally characterized. Crystal data for the complexes are as follows: la, orthorhombic, P2(1)2(1)2(1), a = 7.659(3) Angstrom, b = 22.366(3) Angstrom, c = 23.688(2) Angstrom, V = 4058(2) Angstrom(3), Z = 4, R = 0.0475, and R-w = 0.0467 for 2669 reflections with F-o > 2 sigma(F-o); 2a . 3H(2)O, triclinic, , a = 13.735(3) Angstrom, b = 14.428(4) Angstrom, c = 20.515(8) Angstrom, alpha = 87.13(3)degrees, beta = 87.61(3)degrees, gamma = 63.92(2)degrees, V = 3646(2) Angstrom(3), Z = 4, R = 0.0485 and R-w = 0.0583 for 10 594 reflections with F-o > 6 sigma(F-o); 4a . 1.5H(2)O triclinic, , a = 11.969(3) Angstrom, b = 12.090(6) Angstrom, c = 17.421(3) Angstrom, alpha = 108.93(2)degrees, beta = 84.42(2)degrees, gamma = 105.97(2)degrees, V = 2292(1) Angstrom(3), Z = 2, R = 0.0567, and R-w = 0.0705 for 6775 reflections with F-o > 6 sigma(F-o). The complexes have a diruthenium unit held by an oxo and two carboxylate ligands, and the imidazole ligands occupy the terminal sites of the core. The Ru-Ru distance and the Ru-O-oxo-Ru angle in la and 2a . 3H(2)O are 3.266(1), 3.272(1) Angstrom and 122.4(4), 120.5(2)degrees, while in 4a . 1.5H(2)O these values are 3.327(1) Angstrom and 133.6(2)degrees. The diruthenium(III) complexes 1-3 are blue in color and they exhibit an intense visible band in the range 560-575 nm. The absorption is charge transfer in nature involving the Ru(III)-d pi and O-oxo-p pi orbitals. The diruthenium(III,IV) complexes are red in color and show an intense band near 500 nm. The diruthenium(III) core readily gets oxidized with K2S2O8 forming quantitatively the diruthenium(III,IV) complex. The visible spectral record of the conversion shows an isosbestic point at 545 nm for 1 and at 535 nm for 2 and 3. Protonation of the oxide bridge by HClO4 in methanol yields the [Ru-2(mu-OH)(mu-O2CR)(2)](3+) core. The hydroxo species shows a visible band al 550 nm. The pK(a) value for la is 2.45. The protonated species are unstable. The 1-MeIm species converts to the diruthenium(III,IV) core, while the imidazole complex converts to [Ru(ImH)(6)](3+) and some uncharacterized products. Complex [Ru(ImH)(6)](ClO4)(3) has been structurally characterized. The diruthenium(III) complexes are essentially diamagnetic and show characteristic H-1 NMR spectra indicating the presence of the dimeric structure in solution. The diruthenium(III,IV) complexes are paramagnetic and display rhombic EPR spectral features. Complexes 1-3 are redox active. Complex 1 shows the one-electron reversible Ru-2(III)/(RuRuIV)-Ru-III, one-electron quasireversible (RuRuIV)-Ru-III/Ru-2(IV), and two-electron quasireversible Ru-2(III)/Ru-2(II) couples near 0.4, 1.5, and -1.0 V vs SCE In MeCN-0.1 M TBAP, respectively, in the cyclic and differential pulse voltammetric studies. Complexes 2 and 3 exhibit only reversible Ru-2(III)/(RuRuIV)-Ru-III and the quasireversible (RuRuIV)-Ru-III/Ru-2(IV) couples near 0.4 and 1.6 V vs SCE, respectively, The observation of a quasireversible one-step two-electron transfer reduction process in 1 is significant considering its relevance to the rapid and reversible Fe-2(III)/Fe-2(II) redox process known for the tribridged diiron core in the oxy and deoxy forms of hemerythrin.

Role of TATATAA element in the regulation of tRNA(1)(Gly) gene expression in Bombyx mori is position dependent

Transcription of tRNA genes by RNA polymerase III is controlled by the internal conserved sequences within the coding region and the immediate upstream flanking sequences. A highly transcribed copy of glycyl tRNA gene tRNA1(Gly)-1 from Bombyx mori is down regulated by sequences located much farther upstream in the region -150 to -300 nucleotides (nt), with respect to the +1 nt of tRNA. The negative regulatory effect has been narrowed down to a sequence motif 'TATATAA', a perfect consensus recognised by the TATA binding protein, TBP. This sequence element, when brought closer to the transcription start point, on the other hand, exerts a positive effect by promoting transcription of the gene devoid of other cis regulatory elements. The identity of the nuclear protein interacting with this 'TATATAA' element to TBP has been established by antibody and mutagenesis studies. The 'TATATAA' element thus influences the transcription of tRNA genes positively or negatively in a position-dependent manner either by recruitment or sequestration of TBP from the transcription machinery.

Genome-Wide Gene Expression Analysis Reveals a Dynamic Interplay between Luteotropic and Luteolytic Factors in the Regulation of Corpus Luteum Function in the Bonnet Monkey (Macaca radiata)

Although LH is essential for survival and function of the corpus luteum (CL) in higher primates, luteolysis occurs during nonfertile cycles without a discernible decrease in circulating LH levels. Using genome-wide expression analysis, several experiments were performed to examine the processes of luteolysis and rescue of luteal function in monkeys. Induced luteolysis with GnRH receptor antagonist (Cetrorelix) resulted in differential regulation of 3949 genes, whereas replacement with exogenous LH (Cetrorelix plus LH) led to regulation of 4434 genes (1563 down-regulation and 2871 up-regulation). A model system for prostaglandin (PG) F-2 alpha-induced luteolysis in the monkey was standardized and demonstrated that PGF(2 alpha) regulated expression of 2290 genes in the CL. Analysis of the LH-regulated luteal transcriptome revealed that 120 genes were regulated in an antagonistic fashion by PGF(2 alpha). Based on the microarray data, 25 genes were selected for validation by real-time RT-PCR analysis, and expression of these genes was also examined in the CL throughout the luteal phase and from monkeys treated with human chorionic gonadotropin (hCG) to mimic early pregnancy. The results indicated changes in expression of genes favorable to PGF(2 alpha) action during the late to very late luteal phase, and expressions of many of these genes were regulated in an opposite manner by exogenous hCG treatment. Collectively, the findings suggest that curtailment of expression of downstream LH-target genes possibly through PGF(2 alpha) action on the CL is among the mechanisms underlying cross talk between the luteotropic and luteolytic signaling pathways that result in the cessation of luteal function, but hCG is likely to abrogate the PGF(2 alpha)-responsive gene expression changes resulting in luteal rescue crucial for the maintenance of early pregnancy. (Endocrinology 150: 1473-1484, 2009)

Gene patent practice across plant and human genomes

The uses of genetic sequences to inform, enable or create products or services for human biomedicine are substantially different from their uses in crop-based agriculture. Here, we explore what similarities and differences may emerge in patent use and strategies, and map patent-disclosed sequences onto three important plant genomes: maize (corn), rice and soybean. We focus on those referenced in the granted patent claims to compare their uses to the approach used in human gene patenting.

The ownership question of plant gene and genome intellectual properties

The restructuring of the crop agriculture industry over the past two decades has enabled patent holders to exclude, prevent and deter others from using certain research tools and delay or block further follow-on inventions

Synthesis, structure and spectroscopic studies on hydroxo-bridged heterotrinuclear complexes of cobalt(m) with ethane-l,2-diamine and ammonia as terminal ligands

Hydroxo-bridged homo- and hetero-trinuclear cobalt(III) complexes of the type [MII(H2O)2{(OH)2CoIII(N4)}2]X2·nH2O [MII= a divalent metal ion such as CoII, NiII or ZnII; N4=(en)2(en = ethane-1,2-diamine) or (NH3)4; X = SO4 or (ClO4)2; n= 3 or 5] have been prepared and spectroscopically characterized. The structure of [Cu{(OH)2Co(en)2}2][SO4]2·2H2O 1 has been determined. The geometry around copper atom is a pseudo-square-based pyramid, with the basal sites occupied by four bridging hydroxide oxygens and the apical site is occupied by a weakly co-ordinated sulfate anion [Cu–O 2.516(4)Å]. The hydroxo groups bridge pairs of cobalt(III) ions which are in near-octahedral environments. The ethylenediamine chelate rings have the twist conformation. In the crystal structure of [Cu{(OH)2Co(en)2}2][ClO4]4·2H2O 2 the perchlorate ion is not co-ordinated and the en ligands have envelope conformations. The sulfate ion in [Cu{(OH)2Co(NH3)4}2][SO4]2·4H2O 3 is not co-ordinated to the central copper ion. Electronic, infrared and variable-temperature EPR spectral data are discussed.

Identification of a discrete intermediate in the assembly disassembly of physalis mottle tymovirus through mutational analysis

Assembly intermediates of icosahedral viruses are usually transient and are difficult to identify. In the present investigation, site-specific and deletion mutants of the coat protein gene of physalis mottle tymovirus (PhMV) were used to delineate the role of specific amino acid residues in the assembly of the virus and to identify intermediates in this process. N-terminal 30, 34, 35 and 39 amino acid deletion and single C-terminal (N188) deletion mutant proteins of PhMV were expressed in Escherichia coli. Site-specific mutants H69A, C75A, W96A, D144N, D144N-T151A, K143E and N188A were also constructed and expressed. The mutant protein lacking 30 amino acid residues from the N terminus self-assembled to T = 3 particles in vivo while deletions of 34, 35 and 39 amino acid residues resulted in the mutant proteins that were insoluble. Interestingly, the coat protein (pR PhCP) expressed using pRSET B vector with an additional 41 amino acid residues at the N terminus also assembled into T = 3 particles that were more compact and had a smaller diameter. These results demonstrate that the amino-terminal segment is flexible and either the deletion or addition of amino acid residues at the N terminus does not affect T = 3 capsid assembly, in contrast, the deletion of even a single residue from the C terminus (PhN188 Delta 1) resulted in capsids that were unstable. These capsids disassembled to a discrete intermediate with a sedimentation coefficent of 19.4 S. However, the replacement of C-terminal asparagine 188 by alanine led to the formation of stable capsids. The C75A and D144N mutant proteins also assembled into capsids that were as stable as the pR PhCP, suggesting that C75A and D144 are not crucial for the T = 3 capsid assembly. pR PhW96A and pR PhD144N-T151A mutant proteins failed to form capsids and were present as heterogeneous aggregates. Interestingly, the pR PhK143E mutant protein behaved in a manner similar to the C-terminal deletion protein in forming unstable capsids. The intermediate with an s value of 19.4 S was the major assembly product of pR PhH69A mutant protein and could correspond to a 30mer. It is possible that the assembly or disassembly is arrested at a similar stage in pR PhN188 Delta 1, pR PhH69A and pR PhK143E mutant proteins.

Structures of three different neutral polysaccharides of Acinetobacter baumannii, NIPH190, NIPH201, and NIPH615, assigned to K30, K45, and K48 capsule types, respectively, based on capsule biosynthesis gene clusters

Neutral capsular polysaccharides (CPSs) were isolated from Acinetobacter baumannii NIPH190, NIPH201, and NIPH615. The CPSs were found to contain common monosaccharides only and to be branched with a side-chain 1→3-linked β-d-glucopyranose residue. Structures of the oligosaccharide repeat units (K units) of the CPSs were elucidated by 1D and 2D 1H and 13C NMR spectroscopy. Novel CPS biosynthesis gene clusters, designated KL30, KL45, and KL48, were found at the K locus in the genome sequences of NIPH190, NIPH201, and NIPH615, respectively. The genetic content of each gene cluster correlated with the structure of the CPS unit established, and therefore, the capsular types of the strains studied were designated as K30, K45, and K48, respectively. The initiating sugar of each K unit was predicted, and glycosyltransferases encoded by each gene cluster were assigned to the formation of the linkages between sugars in the corresponding K unit.

A linear programming approach for estimating the structure of a sparse linear genetic network from transcript profiling data

Background: A genetic network can be represented as a directed graph in which a node corresponds to a gene and a directed edge specifies the direction of influence of one gene on another. The reconstruction of such networks from transcript profiling data remains an important yet challenging endeavor. A transcript profile specifies the abundances of many genes in a biological sample of interest. Prevailing strategies for learning the structure of a genetic network from high-dimensional transcript profiling data assume sparsity and linearity. Many methods consider relatively small directed graphs, inferring graphs with up to a few hundred nodes. This work examines large undirected graphs representations of genetic networks, graphs with many thousands of nodes where an undirected edge between two nodes does not indicate the direction of influence, and the problem of estimating the structure of such a sparse linear genetic network (SLGN) from transcript profiling data. Results: The structure learning task is cast as a sparse linear regression problem which is then posed as a LASSO (l1-constrained fitting) problem and solved finally by formulating a Linear Program (LP). A bound on the Generalization Error of this approach is given in terms of the Leave-One-Out Error. The accuracy and utility of LP-SLGNs is assessed quantitatively and qualitatively using simulated and real data. The Dialogue for Reverse Engineering Assessments and Methods (DREAM) initiative provides gold standard data sets and evaluation metrics that enable and facilitate the comparison of algorithms for deducing the structure of networks. The structures of LP-SLGNs estimated from the INSILICO1, INSILICO2 and INSILICO3 simulated DREAM2 data sets are comparable to those proposed by the first and/or second ranked teams in the DREAM2 competition. The structures of LP-SLGNs estimated from two published Saccharomyces cerevisae cell cycle transcript profiling data sets capture known regulatory associations. In each S. cerevisiae LP-SLGN, the number of nodes with a particular degree follows an approximate power law suggesting that its degree distributions is similar to that observed in real-world networks. Inspection of these LP-SLGNs suggests biological hypotheses amenable to experimental verification. Conclusion: A statistically robust and computationally efficient LP-based method for estimating the topology of a large sparse undirected graph from high-dimensional data yields representations of genetic networks that are biologically plausible and useful abstractions of the structures of real genetic networks. Analysis of the statistical and topological properties of learned LP-SLGNs may have practical value; for example, genes with high random walk betweenness, a measure of the centrality of a node in a graph, are good candidates for intervention studies and hence integrated computational – experimental investigations designed to infer more realistic and sophisticated probabilistic directed graphical model representations of genetic networks. The LP-based solutions of the sparse linear regression problem described here may provide a method for learning the structure of transcription factor networks from transcript profiling and transcription factor binding motif data.

Roles of forced nicotine exposure and Comt gene disruption in the development of addiction-related behavioural and neurochemical changes in mice

Activation of midbrain dopamine systems is thought to be critically involved in the addictive properties of abused substances. Drugs of abuse increase dopamine release in the nucleus accumbens and dorsal striatum, which are the target areas of mesolimbic and nigrostriatal dopamine pathways, respectively. Dopamine release in the nucleus accumbens is thought to mediate the attribution of incentive salience to rewards, and dorsal striatal dopamine release is involved in habit formation. In addition, changes in the function of prefrontal cortex (PFC), the target area of mesocortical dopamine pathway, may skew information processing and memory formation such that the addict pays an abnormal amount of attention to drug-related cues. In this study, we wanted to explore how long-term forced oral nicotine exposure or the lack of catechol-O-methyltransferase (COMT), one of the dopamine metabolizing enzymes, would affect the functioning of these pathways. We also wanted to find out how the forced nicotine exposure or the lack of COMT would affect the consumption of nicotine, alcohol, or cocaine. First, we studied the effect of forced chronic nicotine exposure on the sensitivity of dopamine D2-like autoreceptors in microdialysis and locomotor activity experiments. We found that the sensitivity of these receptors was unchanged after forced oral nicotine exposure, although an increase in the sensitivity was observed in mice treated with intermittent nicotine injections twice daily for 10 days. Thus, the effect of nicotine treatment on dopamine autoreceptor sensitivity depends on the route, frequency, and time course of drug administration. Second, we investigated whether the forced oral nicotine exposure would affect the reinforcing properties of nicotine injections. The chronic nicotine exposure did not significantly affect the development of conditioned place preference to nicotine. In the intravenous self-administration paradigm, however, the nicotine-exposed animals self-administered nicotine at a lower unit dose than the control animals, indicating that their sensitivity to the reinforcing effects of nicotine was enhanced. Next, we wanted to study whether the Comt gene knock-out animals would be a suitable model to study alcohol and cocaine consumption or addiction. Although previous work had shown male Comt knock-out mice to be less sensitive to the locomotor-activating effects of cocaine, the present study found that the lack of COMT did not affect the consumption of cocaine solutions or the development of cocaine-induced place preference. However, the present work did find that male Comt knock-out mice, but not female knock-out mice, consumed ethanol more avidly than their wild-type littermates. This finding suggests that COMT may be one of the factors, albeit not a primary one, contributing to the risk of alcoholism. Last, we explored the effect of COMT deficiency on dorsal striatal, accumbal, and prefrontal cortical dopamine metabolism under no-net-flux conditions and under levodopa load in freely-moving mice. The lack of COMT did not affect the extracellular dopamine concentrations under baseline conditions in any of the brain areas studied. In the prefrontal cortex, the dopamine levels remained high for a prolonged time after levodopa treatment in male, but not female, Comt knock-out mice. COMT deficiency induced accumulation of 3,4-dihydroxyphenylacetic acid, which increased further under levodopa load. Homovanillic acid was not detectable in Comt knock-out animals either under baseline conditions or after levodopa treatment. Taken together, the present results show that although forced chronic oral nicotine exposure affects the reinforcing properties of self-administered nicotine, it is not an addiction model itself. COMT seems to play a minor role in dopamine metabolism and in the development of addiction under baseline conditions, indicating that dopamine function in the brain is well-protected from perturbation. However, the role of COMT becomes more important when the dopaminergic system is challenged, such as by pharmacological manipulation.