A network control theory pipeline for studying the dynamics of the structural connectome.

Network control theory (NCT) is a simple and powerful tool for studying how network topology informs and constrains the dynamics of a system. Compared to other structure-function coupling approaches, the strength of NCT lies in its capacity to predict the patterns of external control signals that may alter the dynamics of a system in a desired way. An interesting development for NCT in the neuroscience field is its application to study behavior and mental health symptoms. To date, NCT has been validated to study different aspects of the human structural connectome. NCT outputs can be monitored throughout developmental stages to study the effects of connectome topology on neural dynamics and, separately, to test the coherence of empirical datasets with brain function and stimulation. Here, we provide a comprehensive pipeline for applying NCT to structural connectomes by following two procedures. The main procedure focuses on computing the control energy associated with the transitions between specific neural activity states. The second procedure focuses on computing average controllability, which indexes nodes' general capacity to control the dynamics of the system. We provide recommendations for comparing NCT outputs against null network models, and we further support this approach with a Python-based software package called 'network control theory for python'. The procedures in this protocol are appropriate for users with a background in network neuroscience and experience in dynamical systems theory.

LRRK2 kinase inhibition reverses G2019S mutation-dependent effects on tau pathology progression.

BACKGROUND: Mutations in leucine-rich repeat kinase 2 (LRRK2) are the most common cause of familial Parkinson's disease (PD). These mutations elevate the LRRK2 kinase activity, making LRRK2 kinase inhibitors an attractive therapeutic. LRRK2 kinase activity has been consistently linked to specific cell signaling pathways, mostly related to organelle trafficking and homeostasis, but its relationship to PD pathogenesis has been more difficult to define. LRRK2-PD patients consistently present with loss of dopaminergic neurons in the substantia nigra but show variable development of Lewy body or tau tangle pathology. Animal models carrying LRRK2 mutations do not develop robust PD-related phenotypes spontaneously, hampering the assessment of the efficacy of LRRK2 inhibitors against disease processes. We hypothesized that mutations in LRRK2 may not be directly related to a single disease pathway, but instead may elevate the susceptibility to multiple disease processes, depending on the disease trigger. To test this hypothesis, we have previously evaluated progression of α-synuclein and tau pathologies following injection of proteopathic seeds. We demonstrated that transgenic mice overexpressing mutant LRRK2 show alterations in the brain-wide progression of pathology, especially at older ages. METHODS: Here, we assess tau pathology progression in relation to long-term LRRK2 kinase inhibition. Wild-type or LRRK2G2019S knock-in mice were injected with tau fibrils and treated with control diet or diet containing LRRK2 kinase inhibitor MLi-2 targeting the IC50 or IC90 of LRRK2 for 3-6 months. Mice were evaluated for tau pathology by brain-wide quantitative pathology in 844 brain regions and subsequent linear diffusion modeling of progression. RESULTS: Consistent with our previous work, we found systemic alterations in the progression of tau pathology in LRRK2G2019S mice, which were most pronounced at 6 months. Importantly, LRRK2 kinase inhibition reversed these effects in LRRK2G2019S mice, but had minimal effect in wild-type mice, suggesting that LRRK2 kinase inhibition is likely to reverse specific disease processes in G2019S mutation carriers. Additional work may be necessary to determine the potential effect in non-carriers. CONCLUSIONS: This work supports a protective role of LRRK2 kinase inhibition in G2019S carriers and provides a rational workflow for systematic evaluation of brain-wide phenotypes in therapeutic development.

Neural Network Connectivity Following Opioid Dependence is Altered by a Common Genetic Variant in the µ-Opioid Receptor, OPRM1 A118G.

Opioid use disorder is a chronic, relapsing disease associated with persistent changes in brain plasticity. A common single nucleotide polymorphism (SNP) in the µ-opioid receptor gene, OPRM1 A118G, is associated with altered vulnerability to opioid addiction. Reconfiguration of neuronal connectivity may explain dependence risk in individuals with this SNP. Mice with the equivalent Oprm1 variant, A112G, demonstrate sex-specific alterations in the rewarding properties of morphine and heroin. To determine whether this SNP influences network-level changes in neuronal activity, we compared FOS expression in male and female mice that were opioid-naive or opioid-dependent. Network analyses identified significant differences between the AA and GG Oprm1 genotypes. Based on several graph theory metrics, including small-world analysis and degree centrality, we show that GG females in the opioid-dependent state exhibit distinct patterns of connectivity compared to other groups of the same genotype. Using a network control theory approach, we identified key cortical brain regions that drive the transition between opioid-naive and opioid-dependent brain states; however, these regions are less influential in GG females leading to sixfold higher average minimum energy needed to transition from the acute to the dependent state. In addition, we found that the opioid-dependent brain state is significantly less stable in GG females compared to other groups. Collectively, our findings demonstrate sex- and genotype-specific modifications in local, mesoscale, and global properties of functional brain networks following opioid exposure and provide a framework for identifying genotype differences in specific brain regions that play a role in opioid dependence.

Using network control theory to study the dynamics of the structural connectome.

Network control theory (NCT) is a simple and powerful tool for studying how network topology informs and constrains dynamics. Compared to other structure-function coupling approaches, the strength of NCT lies in its capacity to predict the patterns of external control signals that may alter dynamics in a desired way. We have extensively developed and validated the application of NCT to the human structural connectome. Through these efforts, we have studied (i) how different aspects of connectome topology affect neural dynamics, (ii) whether NCT outputs cohere with empirical data on brain function and stimulation, and (iii) how NCT outputs vary across development and correlate with behavior and mental health symptoms. In this protocol, we introduce a framework for applying NCT to structural connectomes following two main pathways. Our primary pathway focuses on computing the control energy associated with transitioning between specific neural activity states. Our second pathway focuses on computing average controllability, which indexes nodes' general capacity to control dynamics. We also provide recommendations for comparing NCT outputs against null network models. Finally, we support this protocol with a Python-based software package called network control theory for python (nctpy).

Network models to enhance the translational impact of cross-species studies.

Neuroscience studies are often carried out in animal models for the purpose of understanding specific aspects of the human condition. However, the translation of findings across species remains a substantial challenge. Network science approaches can enhance the translational impact of cross-species studies by providing a means of mapping small-scale cellular processes identified in animal model studies to larger-scale inter-regional circuits observed in humans. In this Review, we highlight the contributions of network science approaches to the development of cross-species translational research in neuroscience. We lay the foundation for our discussion by exploring the objectives of cross-species translational models. We then discuss how the development of new tools that enable the acquisition of whole-brain data in animal models with cellular resolution provides unprecedented opportunity for cross-species applications of network science approaches for understanding large-scale brain networks. We describe how these tools may support the translation of findings across species and imaging modalities and highlight future opportunities. Our overarching goal is to illustrate how the application of network science tools across human and animal model studies could deepen insight into the neurobiology that underlies phenomena observed with non-invasive neuroimaging methods and could simultaneously further our ability to translate findings across species.

Molecular and long-term behavioral consequences of neonatal opioid exposure and withdrawal in mice.

Introduction: Infants exposed to opioids in utero are at high risk of exhibiting Neonatal Opioid Withdrawal Syndrome (NOWS), a combination of somatic withdrawal symptoms including high pitched crying, sleeplessness, irritability, gastrointestinal distress, and in the worst cases, seizures. The heterogeneity of in utero opioid exposure, particularly exposure to polypharmacy, makes it difficult to investigate the underlying molecular mechanisms that could inform early diagnosis and treatment of NOWS, and challenging to investigate consequences later in life. Methods: To address these issues, we developed a mouse model of NOWS that includes gestational and post-natal morphine exposure that encompasses the developmental equivalent of all three human trimesters and assessed both behavior and transcriptome alterations. Results: Opioid exposure throughout all three human equivalent trimesters delayed developmental milestones and produced acute withdrawal phenotypes in mice reminiscent of those observed in infants. We also uncovered different patterns of gene expression depending on the duration and timing of opioid exposure (3-trimesters, in utero only, or the last trimester equivalent only). Opioid exposure and subsequent withdrawal affected social behavior and sleep in adulthood in a sex-dependent manner but did not affect adult behaviors related to anxiety, depression, or opioid response. Discussion: Despite marked withdrawal and delays in development, long-term deficits in behaviors typically associated with substance use disorders were modest. Remarkably, transcriptomic analysis revealed an enrichment for genes with altered expression in published datasets for Autism Spectrum Disorders, which correlate well with the deficits in social affiliation seen in our model. The number of differentially expressed genes between the NOWS and saline groups varied markedly based on exposure protocol and sex, but common pathways included synapse development, the GABAergic and myelin systems, and mitochondrial function.

Mindful attention promotes control of brain network dynamics for self-regulation and discontinues the past from the present.

Mindful attention is characterized by acknowledging the present experience as a transient mental event. Early stages of mindfulness practice may require greater neural effort for later efficiency. Early effort may self-regulate behavior and focalize the present, but this understanding lacks a computational explanation. Here we used network control theory as a model of how external control inputs-operationalizing effort-distribute changes in neural activity evoked during mindful attention across the white matter network. We hypothesized that individuals with greater network controllability, thereby efficiently distributing control inputs, effectively self-regulate behavior. We further hypothesized that brain regions that utilize greater control input exhibit shorter intrinsic timescales of neural activity. Shorter timescales characterize quickly discontinuing past processing to focalize the present. We tested these hypotheses in a randomized controlled study that primed participants to either mindfully respond or naturally react to alcohol cues during fMRI and administered text reminders and measurements of alcohol consumption during 4 wk postscan. We found that participants with greater network controllability moderated alcohol consumption. Mindful regulation of alcohol cues, compared to one's own natural reactions, reduced craving, but craving did not differ from the baseline group. Mindful regulation of alcohol cues, compared to the natural reactions of the baseline group, involved more-effortful control of neural dynamics across cognitive control and attention subnetworks. This effort persisted in the natural reactions of the mindful group compared to the baseline group. More-effortful neural states had shorter timescales than less effortful states, offering an explanation for how mindful attention promotes being present.

A common SNP in Chrna5 enhances morphine reward in female mice.

The single nucleotide polymorphism (SNP) D398N (rs16969968) in CHRNA5, the gene encoding the α5 subunit of the nicotinic acetylcholine receptors (nAChR), has been associated with both nicotine and opiate dependence in human populations. Expression of this SNP on presynaptic VTA dopaminergic (DA) neurons is known to cause a reduction in calcium signaling, leading to alterations in transmitter signaling and altered responses to drugs of abuse. To examine the impact of the Chrna5 SNP on opiate reward and underlying dopaminergic mechanisms, mice harboring two copies of the risk-associated allele (Chrna5 A/A) at a location equivalent to human rs16969968 were generated via CRISPR/cas9 genome editing. We sought to determine whether Chrna5 A/A mice show differences in sensitivity to rewarding properties of morphine using the conditioned place preference paradigm. When mice were tested two weeks after conditioning, female Chrna5 A/A mice showed significantly enhanced preference for the morphine-paired chamber relative to WT females, suggesting that this genotype may enhance opioid reward specifically in females. In contrast, Chrna5 genotype had no effect on locomotor sensitization in male or female mice. Relative to WT females, peak amplitude of ACh-gated currents recorded from VTA DA neurons in Chrna5 A/A females was potentiated 1 day after conditioning with morphine. Increased FOS expression was also observed in Chrna5 A/A mice relative to WT mice following exposure to the morphine CPP chamber. We propose that impaired α5 nAChR subunit function alters DA neuron response following repeated morphine exposures, and that this early cellular response could contribute to enhanced opiate reward two weeks after conditioning.

Linking the CHRNA5 SNP to drug abuse liability: From circuitry to cellular mechanisms.

Genetics are known to be a significant risk factor for drug abuse. In human populations, the single nucleotide polymorphism (SNP) D398N in the gene CHRNA5 has been associated with addiction to nicotine, opioids, cocaine, and alcohol. In this paper, we review findings from studies in humans, rodent models, and cell lines and provide evidence that collectively suggests that the Chrna5 SNP broadly influences the response to drugs of abuse in a manner that is not substance-specific. This finding has important implications for our understanding of the role of the cholinergic system in reward and addiction vulnerability. This article is part of the special issue on 'Vulnerabilities to Substance Abuse.'

Sex-specific transgenerational effects of morphine exposure on reward and affective behaviors.

Current estimates indicate that millions of people in the United States abuse opioid drugs, which may also affect their offspring. To determine whether parental exposure to morphine alters reward and affective behaviors in subsequent generations we exposed male and female C57BL/6NTac mice to morphine (75 mg) or placebo pellets for 4 weeks. Naïve mice were used as mating partners to create subsequent generations (F1 and F2). Adult male and female F1 and F2 mice were tested in the morphine conditioned place preference paradigm (CPP), marble burying (MB), acoustic startle response (ASR), and open field tests (OFT). Paternal morphine exposure resulted in significantly attenuated preference scores amongst F1 male offspring, but significantly higher preference scores amongst F1 female offspring at the lowest CPP dose tested (5 mg/kg). In contrast, maternal exposure to morphine did not affect morphine reward in the F1 generation; however, the F2 male offspring of morphine-exposed F0 females displayed significantly higher CPP preference scores. Preference scores in F2 females were not affected by F0 male or female morphine exposure. Sex-specific alterations in affective behaviors were observed only in the offspring of F0 males exposed to morphine with F1 males spending less time in the center of the open field and F1 females spending more time in the center of the open field. One generation later, affective behaviors were no longer altered in F2 males but F2 females from the F0 male morphine exposure buried more marbles in the MB test. In summary, early exposure to morphine in males and females causes lineage-specific inheritance of reward and affective behaviors.

Gene coexpression patterns predict opiate-induced brain-state transitions.

Opioid addiction is a chronic, relapsing disorder associated with persistent changes in brain plasticity. Reconfiguration of neuronal connectivity may explain heightened abuse liability in individuals with a history of chronic drug exposure. To characterize network-level changes in neuronal activity induced by chronic opiate exposure, we compared FOS expression in mice that are morphine-naïve, morphine-dependent, or have undergone 4 wk of withdrawal from chronic morphine exposure, relative to saline-exposed controls. Pairwise interregional correlations in FOS expression data were used to construct network models that reveal a persistent reduction in connectivity strength following opiate dependence. Further, we demonstrate that basal gene expression patterns are predictive of changes in FOS correlation networks in the morphine-dependent state. Finally, we determine that regions of the hippocampus, striatum, and midbrain are most influential in driving transitions between opiate-naïve and opiate-dependent brain states using a control theoretic approach. This study provides a framework for predicting the influence of specific therapeutic interventions on the state of the opiate-dependent brain.

Intrinsic differences in insular circuits moderate the negative association between nicotine dependence and cingulate-striatal connectivity strength.

The development of brain-based biomarkers to assess nicotine dependence severity and treatment efficacy are essential to improve the current marginally effective treatment outcomes. Cross-sectional resting state functional connectivity (rsFC) studies in humans identified a circuit between the dorsal anterior cingulate cortex and the ventral striatum that negatively correlated with increased nicotine dependence severity but was unaffected by acute nicotine administration, suggesting a trait marker of addiction. However, whether this trait circuit dysregulation is predispositional to or resultant from nicotine dependence is unclear. Using a rat model of nicotine dependence with longitudinal fMRI measurements, we assessed the relationship between ACC-striatal rsFC and nicotine dependence severity. Data-driven modularity-based parcellation of the rat medial prefrontal cortex (mPFC) combined with seed-based connectivity analysis with the striatum recapitulated the cingulate-striatum relationship observed in humans. Furthermore, the relationship between cingulate-striatal brain circuits and nicotine dependence severity as indexed by the intensity of precipitated withdrawal, was fully statistically moderated by a predispositional insular-frontal cortical functional circuit. These data suggest that the identified trans-species ACC-striatal circuit relationship with nicotine dependence severity is dysregulated following chronic nicotine administration-induced dependence and may be biased by individual differences in predispositional insula-based striatal-frontal circuits, highlighting the circuit's potential as a biomarker of dependence severity.

Neurobehavioral effects of neonatal opioid exposure in mice: Influence of the OPRM1 SNP.

Opioid use among pregnant women is a growing public health concern in the United States. Infants exposed to opioids in utero are at risk of exhibiting neonatal opioid withdrawal syndrome (NOWS). The biological mechanisms underlying short and long-term consequences of in utero opioid exposure and NOWS are unknown. A potential genetic factor is a single-nucleotide polymorphism (SNP) in the mu-opioid receptor gene (OPRM1 A118G). Opioid exposed infants with the G-allele spend less time in hospitals after birth. To determine whether this SNP modulates the neurobehavioral effects of neonatal opioid exposure and withdrawal, we used mice possessing the equivalent Oprm1 SNP (A112G). Pups were treated chronically with saline or morphine from postnatal days (PNDs) 1 to 14, a developmental period equivalent to the third trimester of a human pregnancy and a sensitive period for opioid exposure in rodents. Morphine treatment produced significant developmental delays regardless of genotype and increased total ultrasonic vocalizations in males during spontaneous withdrawal. Animals were aged and tested for anxiety and drug response during adolescence and adulthood, respectively. AA morphine-treated animals showed reduced activity in the marble burying task compared with saline controls; however, this effect was absent in AG and GG animals. As adults, AA males exposed to morphine from PNDs 1 to 14 exhibited enhanced development of locomotor sensitization to morphine, whereas females showed reduced locomotor sensitization. These data suggest the involvement of the Oprm1 SNP for certain outcomes of neonatal opioid exposure and highlight the importance of considering sex and genetic variability for the prognosis of NOWS.

CREB deletion increases resilience to stress and downregulates inflammatory gene expression in the hippocampus.

The transcription factor CREB (cyclic AMP response element (CRE)-binding protein) is implicated in the pathophysiology and treatment of depression. Structural and functional studies in both animals and humans suggest that abnormalities of the hippocampus may play a role in depression. CREB regulates thousands of genes, yet to date, only a handful that mediate depression or antidepressant response have been identified as relevant CREB targets. In order to comprehensively identify genes regulated by CREB in the hippocampus, we employed translating ribosome affinity purification (TRAP) to detect actively translating mRNAs in wild type and CREB-deficient mice. Using CrebloxP/loxP; RosaLSL-GFP-L10a mice, we conducted whole genome sequencing to identify transcripts only in cells that lack CREB, as introduction of Cre-recombinase simultaneously deleted CREB and expressed GFP-tagged L10a ribosomes that enabled TRAP. We identified over 200 downregulated genes predominantly associated with inflammation and the immune system, including toll-like receptor 1 (TLR1). To determine if baseline disruption in gene expression in the hippocampus of CREB-deficient mice can modulate behavior, we used unpredictable chronic mild stress (UCMS) to produce a set of behavioral alterations with strong validity for depression. We found that CREB-deficient mice demonstrated resilience to the physiological effects of UCMS and also showed changes in affective behaviors specifically in the presence of stress. TLR1 expression was increased following UCMS in control but not in CREB-deficient mice. The results suggest that CREB-mediated regulation of immune system and inflammatory factors may provide additional targets for the treatment of depression.

Intrinsic Insular-Frontal Networks Predict Future Nicotine Dependence Severity.

Although 60% of the US population have tried smoking cigarettes, only 16% smoke regularly. Identifying this susceptible subset of the population before the onset of nicotine dependence may encourage targeted early interventions to prevent regular smoking and/or minimize severity. While prospective neuroimaging in human populations can be challenging, preclinical neuroimaging models before chronic nicotine administration can help to develop translational biomarkers of disease risk. Chronic, intermittent nicotine (0, 1.2, or 4.8 mg/kg/d; N = 10-11/group) was administered to male Sprague Dawley rats for 14 d; dependence severity was quantified using precipitated withdrawal behaviors collected before, during, and following forced nicotine abstinence. Resting-state fMRI functional connectivity (FC) before drug administration was subjected to a graph theory analytical framework to form a predictive model of subsequent individual differences in nicotine dependence. Whole-brain modularity analysis identified five modules in the rat brain. A metric of intermodule connectivity, participation coefficient, of an identified insular-frontal cortical module predicted subsequent dependence severity, independent of nicotine dose. To better spatially isolate this effect, this module was subjected to a secondary exploratory modularity analysis, which segregated it into three submodules (frontal-motor, insular, and sensory). Higher FC among these three submodules and three of the five originally identified modules (striatal, frontal-executive, and sensory association) also predicted dependence severity. These data suggest that predispositional, intrinsic differences in circuit strength between insular-frontal-based brain networks before drug exposure may identify those at highest risk for the development of nicotine dependence.SIGNIFICANCE STATEMENT Developing biomarkers of individuals at high risk for addiction before the onset of this brain-based disease is essential for prevention, early intervention, and/or subsequent treatment decisions. Using a rodent model of nicotine dependence and a novel data-driven, network-based analysis of resting-state fMRI data collected before drug exposure, functional connections centered on an intrinsic insular-frontal module predicted the severity of nicotine dependence after drug exposure. The predictive capacity of baseline network measures was specific to inter-regional but not within-region connectivity. While insular and frontal regions have consistently been implicated in nicotine dependence, this is the first study to reveal that innate, individual differences in their circuit strength have the predictive capacity to identify those at greatest risk for and resilience to drug dependence.

Activation of AMPK by metformin improves withdrawal signs precipitated by nicotine withdrawal.

Cigarette smoking is the leading cause of preventable disease and death in the United States, with more persons dying from nicotine addiction than any other preventable cause of death. Even though smoking cessation incurs multiple health benefits, the abstinence rate remains low with current medications. Here we show that the AMP-activated protein kinase (AMPK) pathway in the hippocampus is activated following chronic nicotine use, an effect that is rapidly reversed by nicotine withdrawal. Increasing pAMPK levels and, consequently, downstream AMPK signaling pharmacologically attenuate anxiety-like behavior following nicotine withdrawal. We show that metformin, a known AMPK activator in the periphery, reduces withdrawal symptoms through a mechanism dependent on the presence of the AMPKα subunits within the hippocampus. This study provides evidence of a direct effect of AMPK modulation on nicotine withdrawal symptoms and suggests central AMPK activation as a therapeutic target for smoking cessation.

Physiological characterization of a robust survival rodent fMRI method.

Anesthetics are commonly used in preclinical functional MRI studies. It is well-appreciated that proper choice of anesthetics is of critical importance for maintaining a physiologically normal range of autonomic functioning. A recent study, using a low dose of dexmedetomidine (active isomer of medetomidine) in combination with a low dose of isoflurane, suggested stable measurements across repeated fMRI experiments in individual animals with each session lasting up to several hours. The rat default mode network has been successfully identified using this preparation, indicating that this protocol minimally disturbs brain network functions. However, medetomidine is known to cause peripheral vasoconstriction, respiratory suppression, and bradycardia, each of which could independently confound the BOLD signal. The goal of this study was to systematically characterize physiological conditions for fMRI experiments under this anesthetic regimen. To this end, we acquired somatosensory stimulation "task-evoked" and resting-state fMRI to evaluate the integrity of neurovascular coupling and brain network function during three time windows (0-30min, 30-90min, and 90-150min) following dexmedetomidine initiation. Results demonstrate that both evoked BOLD response and resting-state fMRI signal remained stable during the 90-150min time window, while autonomic physiological parameters maintained near-normal conditions during this period. Our data suggest that using a spontaneously-inhaled, low dose of isoflurane in combination with a continuous low dose of dexmedetomidine is a viable option for longitudinal imaging studies in rats.

Constituents and functional implications of the rat default mode network.

The default mode network (DMN) has been suggested to support a variety of self-referential functions in humans and has been fractionated into subsystems based on distinct responses to cognitive tasks and functional connectivity architecture. Such subsystems are thought to reflect functional hierarchy and segregation within the network. Because preclinical models can inform translational studies of neuropsychiatric disorders, partitioning of the DMN in nonhuman species, which has previously not been reported, may inform both physiology and pathophysiology of the human DMN. In this study, we sought to identify constituents of the rat DMN using resting-state functional MRI (rs-fMRI) and diffusion tensor imaging. After identifying DMN using a group-level independent-component analysis on the rs-fMRI data, modularity analyses fractionated the DMN into an anterior and a posterior subsystem, which were further segregated into five modules. Diffusion tensor imaging tractography demonstrates a close relationship between fiber density and the functional connectivity between DMN regions, and provides anatomical evidence to support the detected DMN subsystems. Finally, distinct modulation was seen within and between these DMN subcomponents using a neurocognitive aging model. Taken together, these results suggest that, like the human DMN, the rat DMN can be partitioned into several subcomponents that may support distinct functions. These data encourage further investigation into the neurobiological mechanisms of DMN processing in preclinical models of both normal and disease states.

A novel method to induce nicotine dependence by intermittent drug delivery using osmotic minipumps.

Although osmotic minipumps are a reliable method for inducing nicotine dependence in rodents, continuous nicotine administration does not accurately model the intermittent pattern of nicotine intake in cigarette smokers. Our objectives, therefore, were to investigate whether intermittent nicotine delivery via osmotic minipumps could induce dependence in rats, and to compare the magnitude and duration of withdrawal following forced abstinence from intermittent nicotine to that induced by continuous nicotine administration. In order to administer nicotine intermittently, rats were surgically implanted with saline-filled osmotic minipumps attached to polyethylene tubing that contained hourly unit doses of nicotine alternating with mineral oil to mimic "injections". Three doses of nicotine (1.2, 2.4, and 4.8mg/kg/day) and saline were administered for 14days using this method. In order to compare our intermittent delivery method with the more traditional continuous nicotine delivery, a second group of rats was implanted with minipumps attached to tubing that delivered continuous nicotine for 14days. Rats were administered a 1.5mg/kg subcutaneous (SC) mecamylamine challenge and observed for somatic signs of withdrawal on days 7, 14, 21, and 28 following minipump implantation. Fifteen somatic withdrawal signs were summed within a 50-minute observation period to obtain a composite Dependence Score. A generalized linear mixed-effects model revealed a significant Day×Dose×Method interaction. Amongst continuously-treated rats, only 4.8mg/kg/d nicotine resulted in dependence scores significantly greater than those of controls at 14days of exposure. In contrast, all intermittent nicotine groups showed significantly higher scores beginning at 7days of exposure and persisting beyond 7days of abstinence. In general, intermittent delivery produced a more robust withdrawal syndrome than continuous delivery, and did so at a lower dose threshold and with greater persistence after forced abstinence.