The number and timing of time -dependent covariate measurements for the Cox regression model

The problem of analyzing data with updated measurements in the time-dependent proportional hazards model arises frequently in practice. One available option is to reduce the number of intervals (or updated measurements) to be included in the Cox regression model. We empirically investigated the bias of the estimator of the time-dependent covariate while varying the effect of failure rate, sample size, true values of the parameters and the number of intervals. We also evaluated how often a time-dependent covariate needs to be collected and assessed the effect of sample size and failure rate on the power of testing a time-dependent effect.^ A time-dependent proportional hazards model with two binary covariates was considered. The time axis was partitioned into k intervals. The baseline hazard was assumed to be 1 so that the failure times were exponentially distributed in the ith interval. A type II censoring model was adopted to characterize the failure rate. The factors of interest were sample size (500, 1000), type II censoring with failure rates of 0.05, 0.10, and 0.20, and three values for each of the non-time-dependent and time-dependent covariates (1/4,1/2,3/4).^ The mean of the bias of the estimator of the coefficient of the time-dependent covariate decreased as sample size and number of intervals increased whereas the mean of the bias increased as failure rate and true values of the covariates increased. The mean of the bias of the estimator of the coefficient was smallest when all of the updated measurements were used in the model compared with two models that used selected measurements of the time-dependent covariate. For the model that included all the measurements, the coverage rates of the estimator of the coefficient of the time-dependent covariate was in most cases 90% or more except when the failure rate was high (0.20). The power associated with testing a time-dependent effect was highest when all of the measurements of the time-dependent covariate were used. An example from the Systolic Hypertension in the Elderly Program Cooperative Research Group is presented. ^

THE ROLE OF A2B ADENOSINE RECEPTOR SIGNALING IN ADENOSINE DEPENDENT LUNG DISEASE

Chronic lung diseases and acute lung injuries are two distinctive pulmonary disorders that result in significant morbidity and mortality. Adenosine is a signaling nucleoside generated in response to injury and can serve both protective and destructive functions in tissues and cells through interaction with four G-protein coupled adenosine receptors: A1R, A2AR, A2BR, and A3R. However, the relationship between these factors is poorly understood. Recent findings suggest the A2BR has been implicated in the regulation of both chronic lung disease and acute lung injury. The work presented in this dissertation utilized the adenosine deaminase-deficient mouse model and the bleomycin-induced pulmonary injury model to determine the distinctive roles of the A2BR at different stages of the disease. Results demonstrate that the A2BR plays a protective role in attenuating vascular leakage in acute lung injuries and a detrimental role at chronic stages of the disease. In addition, tissues from patients with chronic obstructive pulmonary disease and idiopathic pulmonary fibrosis were utilized to examine adenosine metabolism and signaling in chronic lung diseases. Results demonstrate that components of adenosine metabolism and signaling are altered in a manner that promotes adenosine production and signaling in the lungs of these patients. Furthermore, this study provides the first evidence that A2BR signaling can promote the production of inflammatory and fibrotic mediators in patients with these disorders. Taken together, these findings suggest that the A2BR may have a bi-phasic effect at different stages of lung disease. It is protective in acute injury, whereas pro-inflammatory and pro-fibrotic at the chronic stage. Patients with acute lung injury or chronic lung disease may both benefit from adenosine and A2BR-based therapeutics.

A model of the roles of essential kinases in the induction and expression of late long-term potentiation.

The induction of late long-term potentiation (L-LTP) involves complex interactions among second-messenger cascades. To gain insights into these interactions, a mathematical model was developed for L-LTP induction in the CA1 region of the hippocampus. The differential equation-based model represents actions of protein kinase A (PKA), MAP kinase (MAPK), and CaM kinase II (CAMKII) in the vicinity of the synapse, and activation of transcription by CaM kinase IV (CAMKIV) and MAPK. L-LTP is represented by increases in a synaptic weight. Simulations suggest that steep, supralinear stimulus-response relationships between stimuli (e.g., elevations in [Ca(2+)]) and kinase activation are essential for translating brief stimuli into long-lasting gene activation and synaptic weight increases. Convergence of multiple kinase activities to induce L-LTP helps to generate a threshold whereby the amount of L-LTP varies steeply with the number of brief (tetanic) electrical stimuli. The model simulates tetanic, -burst, pairing-induced, and chemical L-LTP, as well as L-LTP due to synaptic tagging. The model also simulates inhibition of L-LTP by inhibition of MAPK, CAMKII, PKA, or CAMKIV. The model predicts results of experiments to delineate mechanisms underlying L-LTP induction and expression. For example, the cAMP antagonist RpcAMPs, which inhibits L-LTP induction, is predicted to inhibit ERK activation. The model also appears useful to clarify similarities and differences between hippocampal L-LTP and long-term synaptic strengthening in other systems.

Learning reward timing in cortex through reward dependent expression of synaptic plasticity.

The ability to represent time is an essential component of cognition but its neural basis is unknown. Although extensively studied both behaviorally and electrophysiologically, a general theoretical framework describing the elementary neural mechanisms used by the brain to learn temporal representations is lacking. It is commonly believed that the underlying cellular mechanisms reside in high order cortical regions but recent studies show sustained neural activity in primary sensory cortices that can represent the timing of expected reward. Here, we show that local cortical networks can learn temporal representations through a simple framework predicated on reward dependent expression of synaptic plasticity. We assert that temporal representations are stored in the lateral synaptic connections between neurons and demonstrate that reward-modulated plasticity is sufficient to learn these representations. We implement our model numerically to explain reward-time learning in the primary visual cortex (V1), demonstrate experimental support, and suggest additional experimentally verifiable predictions.

Spike timing dependent plasticity: a consequence of more fundamental learning rules.

Spike timing dependent plasticity (STDP) is a phenomenon in which the precise timing of spikes affects the sign and magnitude of changes in synaptic strength. STDP is often interpreted as the comprehensive learning rule for a synapse - the "first law" of synaptic plasticity. This interpretation is made explicit in theoretical models in which the total plasticity produced by complex spike patterns results from a superposition of the effects of all spike pairs. Although such models are appealing for their simplicity, they can fail dramatically. For example, the measured single-spike learning rule between hippocampal CA3 and CA1 pyramidal neurons does not predict the existence of long-term potentiation one of the best-known forms of synaptic plasticity. Layers of complexity have been added to the basic STDP model to repair predictive failures, but they have been outstripped by experimental data. We propose an alternate first law: neural activity triggers changes in key biochemical intermediates, which act as a more direct trigger of plasticity mechanisms. One particularly successful model uses intracellular calcium as the intermediate and can account for many observed properties of bidirectional plasticity. In this formulation, STDP is not itself the basis for explaining other forms of plasticity, but is instead a consequence of changes in the biochemical intermediate, calcium. Eventually a mechanism-based framework for learning rules should include other messengers, discrete change at individual synapses, spread of plasticity among neighboring synapses, and priming of hidden processes that change a synapse's susceptibility to future change. Mechanism-based models provide a rich framework for the computational representation of synaptic plasticity.

A biophysical model of synaptic plasticity and metaplasticity can account for the dynamics of the backward shift of hippocampal place fields.

Hippocampal place cells in the rat undergo experience-dependent changes when the rat runs stereotyped routes. One such change, the backward shift of the place field center of mass, has been linked by previous modeling efforts to spike-timing-dependent plasticity (STDP). However, these models did not account for the termination of the place field shift and they were based on an abstract implementation of STDP that ignores many of the features found in cortical plasticity. Here, instead of the abstract STDP model, we use a calcium-dependent plasticity (CaDP) learning rule that can account for many of the observed properties of cortical plasticity. We use the CaDP learning rule in combination with a model of metaplasticity to simulate place field dynamics. Without any major changes to the parameters of the original model, the present simulations account both for the initial rapid place field shift and for the subsequent slowing down of this shift. These results suggest that the CaDP model captures the essence of a general cortical mechanism of synaptic plasticity, which may underlie numerous forms of synaptic plasticity observed both in vivo and in vitro.

Detection and sequence analysis of a nickel-induced mutation in a retroviral model system

We have developed a novel way to assess the mutagenicity of environmentally important metal carcinogens, such as nickel, by creating a positive selection system based upon the conditional expression of a retroviral transforming gene. The target gene is the v-mos gene in MuSVts110, a murine retrovirus possessing a growth temperature dependent defect in expression of the transforming gene due to viral RNA splicing. In normal rat kidney cells infected with MuSVts110 (6m2 cells), splicing of the MuSVts110 RNA to form the mRNA from which the transforming protein, p85$\sp{\rm gag-mos}$, is translated is growth-temperature dependent, occurring at 33 C and below but not at 39 C and above. This splicing "defect" is mediated by cis-acting viral sequences. Nickel chloride treatment of 6m2 cells followed by growth at 39 C, allowed the selection of "revertant" cells which constitutively express p85$\sp{\rm gag-mos}$ due to stable changes in the viral RNA splicing phenotype, suggesting that nickel, a carcinogen whose mutagenicity has not been well established, could induce mutations in mammalian genes. We also show by direct sequencing of PCR-amplified integrated MuSVts110 DNA from a 6m2 nickel-revertant cell line that the nickel-induced mutation affecting the splicing phenotype is a cis-acting 70-base duplication of a region of the viral DNA surrounding the 3$\sp\prime$ splice site. These findings provide the first example of the molecular basis for a nickel-induced DNA lesion and establish the mutagenicity of this potent carcinogen. ^

Early stages of experimental colon carcinogenesis inhibit mast cell-dependent mucosal immune function of the distal colon

Inhibition of local host immune reactions is one mechanism contributing to tumor progression. To determine if alterations in local immune functioning occur during colon carcinogenesis, a model mucosal immune response, type I hypersensitivity against the intestinal parasite Trichinella spiralis, was first characterized in normal mice and then examined during experimental colon carcinogenesis. Segments of sensitized colon mounted in Ussing chambers and challenged with T. spiralis-derived antigen resulted in a rise in short-circuit current ($\rm\Delta I\sb{sc}$) that was antigen-specific and inhibited by furosemide, implicating epithelial Cl$\sp-$ secretion as the ionic mechanism. The immune-regulated Cl$\sp-$ secretion by colonic epithelial cells required the presence of mast cells with surface IgE. Inhibition of potential anaphylactic mediators with various pharmacological agents in vitro implicated prostaglandins and leukotrienes as the principal mediators of the antigen-induced $\rm\Delta I\sb{sc}$, with 5-hydroxytryptamine also playing a role. Distal colon from immune mice fed an aspirin-containing diet (800 mg/kg powdered diet) ad libitum for 6 wk had a decreased response to antigen, confirming the major role of prostaglandins in generating the colonic I$\sb{\rm sc}$. To determine the effects of early stages of colon carcinogenesis on this mucosal immune response, mice were immunized with T. spiralis 1 day after or 8 wk prior to the first of 6 weekly injections of the procarcinogen 1,2-dimethylhydrazine (DMH). Responsiveness to antigenic challenge was suppressed in the distal colon 4-6 wk after the final injection of DMH. One injection of DMH was not sufficient to inhibit antigen responsiveness. The colonic epithelium remained sensitive to direct stimulation by exogenous Cl$\sp-$ secretagogues. Decreased antigen-induced $\rm\Delta I\sb{sc}$ in the distal colon was not due to systemic immune suppression by DMH, as the proximal colon and jejunum maintained responsiveness to antigen. Also, rejection of a secondary T. spiralis infection from the small intestine was not altered. Tumors eventually developed 25-30 wk after the final injection of DMH only in the distal portions of the colon. These results suggest that early stages of DMH-induced colon carcinogenesis manipulate the microenvironment such that mucosal immune function, as measured by immune-regulated Cl$\sp-$ secretion, is suppressed in the distal colon, but not in other regions of the gut. Future elucidation of the mechanisms by which this localized inhibition of immune-mediated ion transport occurs may provide possible clues to the microenvironmental changes necessary for tumor progression in the distal colon. ^

The differential expression of the $\beta\sb2$-adrenergic receptor modifies agonist-dependent adenylylcyclase activation

There have been multiple reports which indicate that variations in $\beta$AR expression affect the V$\sb{\rm max}$ observed for the agonist-dependent activation of adenylylcyclase. This observation has been ignored by most researchers when V$\sb{\rm max}$ values obtained for wild type and mutant receptors are compared. Such an imprecise analysis may lead to erroneous conclusions concerning the ability of a receptor to activate adenylylcyclase. Equations were derived from the Cassel-Selinger model of GTPase activity and Tolkovsky and Levitzki's Collision Coupling model which predict that the EC$\sb{50}$ and V$\sb{\rm max}$ for the activation of adenylylcyclase are a function of receptor number. Experimental results for L cell clones in which either hamster or human $\beta$AR were transfected at varying levels showed that EC$\sb{50}$ decreases and V$\sb{\rm max}$ increases as receptor number increases. Comparison of these results with simulations obtained from the equations describing EC$\sb{50}$ and V$\sb{\rm max}$ showed a close correlation. This documents that the kinetic parameters of adenylylcyclase activation change with the level of receptor expression and relates this phenomenon to a theoretical framework concerning the mechanisms involved in $\beta$AR signal transduction.^ One of the terms used in the equations which expressed the EC$\sb{50}$ and V$\sb{\rm max}$ as a function of receptor number is coupling efficiency, defined as $\rm k\sb1/k\sb{-1}$. Calculation of $\rm k\sb1/k\sb{-1}$ can be accomplished for wild type receptors with the easily measured experimental values of agonist K$\sb{\rm d}$, EC$\sb{50}$ and receptor number. This was demonstrated for hamster $\beta$AR which yielded a coupling efficiency of 0.15 $\pm$ 0.003 and human $\beta$AR which yielded a coupling efficiency of 0.90 $\pm$ 0.031. $\rm k\sb1/k\sb{-1}$ replaces the traditional qualitative evaluation of the ability to activate adenylylcyclase, which utilizes V$\sb{\rm max}$ without correction for variation in receptor number, with a quantitative definition that more accurately describes the ability of $\beta$AR to couple to G$\sb{\rm s}$.^ The equations which express the EC$\sb{50}$ and V$\sb{\rm max}$ for adenylylcyclase activation as a function of receptor number and coupling efficiency were tested to determine whether they could accurately simulate the changes seen in these parameters during desensitization. Data from original desensitization experiments and data from the literature (24,25,52,54,83) were compared to simulated changes in EC$\sb{50}$ and V$\sb{\rm max}$. In a variety of systems the predictions of the equations were consistent with the changes observed in EC$\sb{50}$ and V$\sb{\rm max}$. In addition reductions in the calculated value of $\rm k\sb1/k\sb{-1}$ was shown to correlate well with $\beta$AR phosphorylation and to be minimally affected by sequestration and down-regulation. ^

Analysis of recurrent failure times: A time-dependent Yule process approach

Analysis of recurrent events has been widely discussed in medical, health services, insurance, and engineering areas in recent years. This research proposes to use a nonhomogeneous Yule process with the proportional intensity assumption to model the hazard function on recurrent events data and the associated risk factors. This method assumes that repeated events occur for each individual, with given covariates, according to a nonhomogeneous Yule process with intensity function λx(t) = λ 0(t) · exp( x′β). One of the advantages of using a non-homogeneous Yule process for recurrent events is that it assumes that the recurrent rate is proportional to the number of events that occur up to time t. Maximum likelihood estimation is used to provide estimates of the parameters in the model, and a generalized scoring iterative procedure is applied in numerical computation. ^ Model comparisons between the proposed method and other existing recurrent models are addressed by simulation. One example concerning recurrent myocardial infarction events compared between two distinct populations, Mexican-American and Non-Hispanic Whites in the Corpus Christi Heart Project is examined. ^

Characterization of E2F1's oncogenic, apoptotic and tumor-suppressive activities utilizing the K5 mouse model

Though E2F1 is deregulated in most human cancers by mutations of the p16-cyclin D-Rb pathway, it also exhibits tumor suppressive activity. A transgenic mouse model overexpressing E2F1 under the control of the bovine keratin 5 (K5) promoter exhibits epidermal hyperplasia and spontaneously develops tumors in the skin and other epithelial tissues after one year of age. In a p53-deficient background, aberrant apoptosis in K5 E2F1 transgenic epidermis is reduced and tumorigenesis is accelerated. In sharp contrast, K5 E2F1 transgenic mice are resistant to papilloma formation in the DMBA/TPA two-stage carcinogenesis protocol. K5 E2F4 and K5 DP1 transgenic mice were also characterized and both display epidermal hyperplasia but do not develop spontaneous tumors even in cooperation with p53 deficiency. These transgenic mice do not have increased levels of apoptosis in their skin and are more susceptible to papilloma formation in the two-stage carcinogenesis model. These studies show that deregulated proliferation does not necessarily lead to tumor formation and that the ability to suppress skin carcinogenesis is unique to E2F1. E2F1 can also suppress skin carcinogenesis when okadaic acid is used as the tumor promoter and when a pre-initiated mouse model is used, demonstrating that E2F1's tumor suppressive activity is not specific for TPA and occurs at the promotion stage. E2F1 was thought to induce p53-dependent apoptosis through upregulation of p19ARF tumor suppressor, which inhibits mdm2-mediated p53 degradation. Consistent with in vitro studies, the overexpression of E2F1 in mouse skin results in the transcriptional activation of the p19ARF and the accumulation of p53. Inactivation of either p19ARF or p53 restores the sensitivity of K5 E2F1 transgenic mice to DMBA/TPA carcinogenesis, demonstrating that an intact p19ARF-p53 pathway is necessary for E2F1 to suppress carcinogenesis. Surprisingly, while p53 is required for E2F1 to induce apoptosis in mouse skin, p19ARF is not, and inactivation of p19ARF actually enhances E2F1-induced apoptosis and proliferation in transgenic epidermis. This indicates that ARF is important for E2F1-induced tumor suppression but not apoptosis. Senescence is another potential mechanism of tumor suppression that involves p53 and p19ARF. K5 E2F1 transgenic mice initiated with DMBA and treated with TPA show an increased number of senescence cells in their epidermis. These experiments demonstrate that E2F1's unique tumor suppressive activity in two-stage skin carcinogenesis can be genetically separated from E2F1-induced apoptosis and suggest that senescence utilizing the p19ARF-p53 pathway plays a role in tumor suppression by E2F1. ^

A physiological role of cytochromes P450 4Fs: Mouse model

Cytochromes P450 4Fs (CYP4F) are a subfamily of enzymes involved in arachidonic acid metabolism with highest catalytic activity towards leukotriene B 4 (LTB4), a potent chemoattractant involved in prompting inflammation. CYP4F-mediated metabolism of LTB4 leads to inactive ω-hydroxy products incapable of initiating chemotaxis and the inflammatory stimuli that result in the influx of inflammatory cells. Our hypothesis is based on the catalytic ability of CYP4Fs to inactivate pro-inflammatory LTB4 which assures these enzymes a pivotal role in the process of inflammation resolution. ^ To test this hypothesis and evaluate the changes in CYP4F expression under complex inflammatory conditions, we designed two mouse models, one challenged with lipopolysaccharide (LPS) as a sterile model of sepsis and the other challenged with a systemic live bacterial infection of Citrobacter rodentium, an equivalent of the human enterobacterium E. coli pathogen invasion. Based on the evidence that Peroxisome Proliferator Activated Receptors (PPARs) play an active role in inflammation regulation, we also examined PPARs as a regulation mechanism in CYP4F expression during inflammation using PPARα knockout mice under LPS challenge. Using the Citrobacter rodentium model of inflammation, we studied CYP4F levels to compare them to those in LPS challenged animals. LPS-triggered inflammation signal is mediated by Toll-like 4 (TLR4) receptors which specifically respond to LPS in association with several other proteins. Using TLR4 knockout mice challenged with Citrobacter rodentium we addressed possible mediation of CYP4F expression regulation via these receptors. ^ Our results show isoform- and tissue-specific CYP4F expression in all the tissues examined. The Citrobacter rodentium inflammation model revealed significant reduction in liver expression of CYP4F14 and CYP4F15 and an up-regulation of gene expression of CYP4F16 and CYP4F18. TLR4 knockout studies showed that the decrease in hepatic CYP4F15 expression is TLR4-dependent. CYP4F expression in kidney shows down-regulation of CYP4F14 and CYP4F15 and up-regulation of CYP4F18 expression. In the LPS inflammation model, we showed similar patterns of CYP4F changes as in Citrobacter rodentium -infected mice. The renal profile of CYP4Fs in PPARα knockout mice with LPS challenge showed CYP4F15 down-regulation to be PPARα dependent. Our study confirmed tissue- and isoform-specific regulation of CYP4F isoforms in the course of inflammation. ^

Lipid-dependent topogenesis and function of membrane proteins

Membranes are essential for the integrity and function of the cell. The collective property of the lipid bilayer is critical in providing an optimal functioning environment for membrane proteins. The simple yet well-characterized bacterium Escherichia coli serves an ideal model system to study the function of specific lipids since its lipid content can be easily manipulated. The most abundant lipid in E. coli membrane is phosphatidylethanolamine (PE, 70-80%). A PE-lacking E. coli mutant displays a complex mixture of deficient phenotypes, suggesting a profound role for PE in different aspects of cell function. A novel role of PE as a topological and functional determinant for membrane proteins has been established using lactose permease (LacY) as a model protein. PE is found to be required for energy-dependent uphill transport process of LacY. In PE-lacking membranes, LacY undergoes a dramatic conformational change, and the first half of the protein adopts an inverted topology with respect to the bilayer plane. ^ The work reported here was initiated to understand the molecular properties of lipids that enable their function as topological and functional determinants for membrane proteins. A glycolipid, monoglucosyldiacylglycerol (MGlcDAG) which shares physicochemical similarities with PE, was introduced to PE-lacking E. coli membranes. The introduction of MGlcDAG suppresses many of the PE-deficient phenotypes, and in particular supports the function and native topology of LacY. ^ The lipid-sensitive topogenic signals encoded in the amino acid sequence of LacY were also identified. Native LacY adopts an inverted topology when synthesized without PE, but mutation of specific acidic residues in the cytoplasmic extra-membrane domains can prevent this inversion and supports a native topological organization of LacY in PE-lacking membranes. These results suggest that it is the interplay between the collective charge properties of the lipid bilayer and extra-membrane loops of protein that determines the final orientation of transmembrane domains. By comparing the similarities as well as differences between these two lipids, we established how specific physical and chemical properties of lipids influence various cell functions and elucidated the molecular basis for the novel role of lipids in determining membrane protein topology. ^

Caenorhabditis elegans germinal center kinase (Gck-1) is a p-body component that inhibits precocious ectopic MAP kinase dependent meiotic maturation

The Caenorhabditis elegans germline is an excellent model system for studying meiosis, as the gonad contains germ cells in all stages of meiosis I prophase in a linear temporal and spatial pattern. To form healthy gametes, many events must be coordinated. Failure of any step in the process can reduce fertility. Here, we describe a C. elegans Germinal Center Kinase, GCK-1, that is essential for the accurate progression of germ cells through meiosis I prophase. In the absence of GCK-1, germ cells undergo precocious maturation due to the activation of a specific MAP kinase isoform. Furthermore, GCK-1 localizes to P-bodies, RNP particles that have been implicated in RNA degradation and translational control. Like two other components of C. elegans germline P-bodies, GCK-1 functions to limit physiological germ cell apoptosis. This is the first study to identify a role for a GCK-III kinase in metazoan germ cell development and to link P-body function with MAP kinase activation and germ cell maturation. ^

Blood-oxygen-level-dependent (BOLD) signal changes in total cortex and subcortex, pain, reward, and vigilance regions during mechanical pressure pain after morphine administration

Morphine is the most common clinical choice in the management of severe pain. Although the molecular mechanisms of morphine have already been characterized, the cerebral circuits by which it attenuates the sensation of pain have not yet been studied in humans. The objective of this two-arm (morphine versus placebo), between-subjects study was to examine whether morphine affects pain via pain-related cortical circuits, but also via reward regions that relate to the motivational state, as well as prefrontal regions that relate to vigilance as a result of morphine's sedative effects. Cortical activity was measured by the blood-oxygen-level-dependent (BOLD) signal changes using functional magnetic resonance imaging (fMRI). ^ The novelty of this study is at three levels: (i) to develop a methodology that will assess the average BOLD signal across subjects for the pain, reward, and vigilance cortical systems; (ii) to examine whether the reward and/or sedative effects of morphine are contributing factors to cortical regions associated with the motivational state and vigilance; and (iii) to propose a neuroanatomical model related to the opioid-sensitive effects of reward and sedation as a function of cortical activity related to pain in an effort to assess future analgesics. ^ Consistent with our hypotheses, our findings showed that the decrease in total pain-related volume activated between the post- and the pre-treatment morphine group was about 78%, while the post-treatment placebo group displayed only a 5% decrease when compared to pre-treatment levels of activation. The volume increase in reward regions was 451% in the post-treatment compared to the pre-treatment morphine condition. Finally, the volumetric decrease in vigilance regions was 63% in the posttreatment compared to the pre-treatment morphine condition. ^ These findings imply that changes in the blood flow of the reward and vigilance regions may be contributing factors in producing the analgesic effect under morphine administration. Future studies need to replicate this study in a higher resolution fMRI environment and to assess the proposed neuroanatomical model in patient populations. The necessity of pain research is apparent, since pain cuts across different diseases especially chronic ones, and thus, is recognized as a vital public health developing area. ^