Multidimensional Spectroscopy of Nuclear Spin Clusters in Diamond.

Optically active spin defects in solids offer promising platforms to investigate nuclear spin clusters with high sensitivity and atomic-site resolution. To leverage near-surface defects for molecular structure analysis in chemical and biological contexts using nuclear magnetic resonance (NMR), further advances in spectroscopic characterization of nuclear environments are essential. Here, we report Fourier spectroscopy techniques to improve localization and mapping of the test bed ^{13}C nuclear spin environment of individual, shallow nitrogen-vacancy centers at room temperature. We use multidimensional spectroscopy, well-known from classical NMR, in combination with weak measurements of single-nuclear-spin precession. We demonstrate two examples of multidimensional NMR: (i) improved nuclear spin localization by separate encoding of the two hyperfine components along spectral dimensions and (ii) spectral editing of nuclear-spin pairs, including measurement of internuclear coupling constants. Our work adds important tools for the spectroscopic analysis of molecular structures by single-spin probes.

Estimation of species divergence times in presence of cross-species gene flow.

Cross-species introgression can have significant impacts on phylogenomic reconstruction of species divergence events. Here, we used simulations to show how the presence of even a small amount of introgression can bias divergence time estimates when gene flow is ignored in the analysis. Using advances in analytical methods under the multispecies coalescent (MSC) model, we demonstrate that by accounting for incomplete lineage sorting and introgression using large phylogenomic data sets this problem can be avoided. The multispecies-coalescent-with-introgression (MSci) model is capable of accurately estimating both divergence times and ancestral effective population sizes, even when only a single diploid individual per species is sampled. We characterize some general expectations for biases in divergence time estimation under three different scenarios: 1) introgression between sister species, 2) introgression between non-sister species, and 3) introgression from an unsampled (i.e., ghost) outgroup lineage. We also conducted simulations under the isolation-with-migration (IM) model and found that the MSci model assuming episodic gene flow was able to accurately estimate species divergence times despite high levels of continuous gene flow. We estimated divergence times under the MSC and MSci models from two published empirical datasets with previous evidence of introgression, one of 372 target-enrichment loci from baobabs (Adansonia), and another of 1000 transcriptome loci from 14 species of the tomato relative, Jaltomata. The empirical analyses not only confirm our findings from simulations, demonstrating that the MSci model can reliably estimate divergence times but also show that divergence time estimation under the MSC can be robust to the presence of small amounts of introgression in empirical datasets with extensive taxon sampling. [divergence time; gene flow; hybridization; introgression; MSci model; multispecies coalescent].

Multicone Diamond Waveguides for Nanoscale Quantum Sensing.

The long-lived electronic spin of the nitrogen-vacancy (NV) center in diamonds is a promising quantum sensor for detecting nanoscopic magnetic and electric fields in various environments. However, the poor signal-to-noise ratio (SNR) of prevalent optical spin-readout techniques presents a critical challenge in improving measurement sensitivity. Here, we address this limitation by coupling individual NVs to optimized diamond nanopillars, thereby enhancing the collection efficiency of fluorescence. Guided by near-field optical simulations, we predict improved performance for tall (≥5 µm) pillars with tapered sidewalls. This is subsequently verified by fabricating and characterizing a representative set of structures using a newly developed nanofabrication process. We observe increased SNR for optimized devices, owing to improved emission collimation and directionality. Promisingly, these devices are compatible with low-numerical-aperture collection optics and a reduced tip radius, reducing experimental overhead and facilitating improved spatial resolution for scanning applications.

[68Ga]Ga-DOTA-FAPI-04 PET/MR in patients with acute myocardial infarction: potential role of predicting left ventricular remodeling.

PURPOSE: To assess predictive value of 68Ga-labeled fibroblast activation protein inhibitor-04 ([68Ga]Ga-DOTA-FAPI-04) PET/MR for late left ventricular (LV) remodeling in patients with ST-segment elevated myocardial infarction (STEMI). METHODS: Twenty-six patients with STEMI were included in the study. [68Ga]Ga-DOTA-FAPI-04 PET/MR was performed at baseline and at average 12 months after STEMI. LV remodeling was defined as >10% increase in LV end-systolic volume (LVESV) from baseline to 12 months. RESULTS: The LV remodeling group demonstrated higher [68Ga]Ga-DOTA-FAPI-04 uptake volume (UV) at baseline than the non-LV remodeling group (p < 0.001). [68Ga]Ga-DOTA-FAPI-04 UV at baseline was a significant predictor (OR = 1.048, p = 0.011) for LV remodeling at 12 months after STEMI. Compared to clinical information, MR imaging and cardiac function parameters at baseline, [68Ga]Ga-DOTA-FAPI-04 UV demonstrated better predictive ability (AUC = 0.938, p < 0.001) for late LV remodeling, with sensitivity of 100.0% and specificity of 81.3%. CONCLUSIONS: [68Ga]Ga-DOTA-FAPI-04 PET/MR is an effective tool to non-invasively quantify myocardial fibroblasts activation, and baseline [68Ga]Ga-DOTA-FAPI-04 UV may have potential predictive value for late LV remodeling.

Single-Nitrogen-Vacancy NMR of Amine-Functionalized Diamond Surfaces.

Nuclear magnetic resonance (NMR) imaging with shallow nitrogen-vacancy (NV) centers in diamond offers an exciting route toward sensitive and localized chemical characterization at the nanoscale. Remarkable progress has been made to combat the degradation in coherence time and stability suffered by near-surface NV centers using suitable chemical surface termination. However, approaches that also enable robust control over adsorbed molecule density, orientation, and binding configuration are needed. We demonstrate a diamond surface preparation for mixed nitrogen- and oxygen-termination that simultaneously improves NV center coherence times for <10 nm-deep emitters and enables direct and recyclable chemical functionalization via amine-reactive cross-linking. Using this approach, we probe single NV centers embedded in nanopillar waveguides to perform 19F NMR sensing of covalently bound fluorinated molecules with detection on the order of 100 molecules. This work signifies an important step toward nuclear spin localization and structure interrogation at the single-molecule level.

Inflammatory Myofibroblastic Tumor of the Sigmoid Colon in an Infant: A Case Report and Literature Review.

Inflammatory myofibroblastic tumor (IMT) infrequently involves the sigmoid colon, and has not previously been described in an infant sigmoid colon.An inflammatory myofibroblastic tumor arose from the sigmoid colon of an 11-month-old boy, confirmed by anaplastic lymphoma kinase (ALK), smooth muscle actin (SMA) and desmin immunohistochemical staining. The patient recovered well after complete resection of the tumor.Sigmoid IMT can occur in infancy. This eighth case is the youngest so far. The child did well after surgical resection.

Complexity of the simplest species tree problem.

The multispecies coalescent model provides a natural framework for species tree estimation accounting for gene-tree conflicts. Although a number of species tree methods under the multispecies coalescent have been suggested and evaluated using simulation, their statistical properties remain poorly understood. Here, we use mathematical analysis aided by computer simulation to examine the identifiability, consistency, and efficiency of different species tree methods in the case of three species and three sequences under the molecular clock. We consider four major species-tree methods including concatenation, two-step, independent-sites maximum likelihood, and maximum likelihood. We develop approximations that predict that the probit transform of the species tree estimation error decreases linearly with the square root of the number of loci. Even in this simplest case, major differences exist among the methods. Full-likelihood methods are considerably more efficient than summary methods such as concatenation and two-step. They also provide estimates of important parameters such as species divergence times and ancestral population sizes,whereas these parameters are not identifiable by summary methods. Our results highlight the need to improve the statistical efficiency of summary methods and the computational efficiency of full likelihood methods of species tree estimation.

A simulation study to examine the impact of recombination on phylogenomic inferences under the multispecies coalescent model.

Phylogenomic analyses under the multispecies coalescent model assume no recombination within locus and free recombination among loci. Yet, in real data sets intralocus recombination causes different sites of the same locus to have different genealogical histories so that the model is misspecified. The impact of recombination on various coalescent-based phylogenomic analyses has not been systematically examined. Here, we conduct a computer simulation to examine the impact of recombination on several Bayesian analyses of multilocus sequence data, including species tree estimation, species delimitation (by Bayesian selection of delimitation models) and estimation of evolutionary parameters such as species divergence and introgression times, population sizes for modern and extinct species, and cross-species introgression probabilities. We found that recombination, at rates comparable to estimates from the human being, has little impact on coalescent-based species tree estimation, species delimitation and estimation of population parameters. At rates 10 times higher than the human rate, recombination may affect parameter estimation, causing positive biases in introgression times and ancestral population sizes, although species divergence times and cross-species introgression probabilities are estimated with little bias. Overall, the simulation suggests that phylogenomic inferences under the multispecies coalescent model are robust to realistic amounts of intralocus recombination.

The Asymptotic Behavior of Bootstrap Support Values in Molecular Phylogenetics.

The phylogenetic bootstrap is the most commonly used method for assessing statistical confidence in estimated phylogenies by non-Bayesian methods such as maximum parsimony and maximum likelihood (ML). It is observed that bootstrap support tends to be high in large genomic data sets whether or not the inferred trees and clades are correct. Here, we study the asymptotic behavior of bootstrap support for the ML tree in large data sets when the competing phylogenetic trees are equally right or equally wrong. We consider phylogenetic reconstruction as a problem of statistical model selection when the compared models are nonnested and misspecified. The bootstrap is found to have qualitatively different dynamics from Bayesian inference and does not exhibit the polarized behavior of posterior model probabilities, consistent with the empirical observation that the bootstrap is more conservative than Bayesian probabilities. Nevertheless, bootstrap support similarly shows fluctuations among large data sets, with no convergence to a point value, when the compared models are equally right or equally wrong. Thus, in large data sets strong support for wrong trees or models is likely to occur. Our analysis provides a partial explanation for the high bootstrap support values for incorrect clades observed in empirical data analysis. [Bootstrap; model selection; star-tree paradox; support value.].

Isolation and identification of plant growth-promoting rhizobacteria from tall fescue rhizosphere and their functions under salt stress.

Soil salinity has become one of the major factors that threaten tall fescue growth and turf quality. Plants recruit diverse microorganisms in the rhizosphere to cope with salinity stress. In this study, 15 plant growth-promoting rhizobacteria (PGPR) were isolated from the salt-treated rhizosphere of tall fescue and were annotated to 10 genera, including Agrobacterium, Fictibacillus, Rhizobium, Bhargavaea, Microbacterium, Paenarthrobacter, Pseudarthrobacter, Bacillus, Halomonas, and Paracoccus. All strains could produce indole-3-acetic acid (IAA). Additionally, eight strains exhibited the ability to solubilize phosphate and potassium. Most strains could grow on the medium containing 600 mM NaCl, such as Bacillus zanthoxyli and Bacillus altitudinis. Furthermore, Bacillus zanthoxyli and Bacillus altitudinis were inoculated with tall fescue seeds and seedlings to determine their growth-promoting effect. The results showed that Bacillus altitudinis and mixed culture significantly increased the germination rate of tall fescue seeds. Bacillus zanthoxyli can significantly increase the tillers number and leaf width of seedlings under salt conditions. Through the synergistic effect of FaSOS1, FaHKT1, and FaHAK1 genes, Bacillus zanthoxyli helps to expel the excess Na+ from aboveground parts and absorb more K+ in roots to maintain ion homeostasis in tall fescue. Unexpectedly, we found that Bacillus altitudinis displayed an inapparent growth-promoting effect on seedlings under salt stress. Interestingly, the mixed culture of the two strains was also able to alleviate, to some extent, the effects of salt stress on tall fescue. This study provides a preliminary understanding of tall fescue rhizobacteria and highlights the role of Bacillus zanthoxyli in tall fescue growth and salt tolerance.

A CMOS compatible thermoelectric device made of crystalline silicon membranes with nanopores.

Herein, we report the use of nanostructured crystalline silicon as a thermoelectric material and its integration into thermoelectric devices. The proof-of-concept relies on the partial suppression of lattice thermal conduction by introducing pores with dimensions scaling between the electron mean free path and the phonon mean free path. In other words, we artificially aimed at the well-known 'electron crystal and phonon glass' trade-off targeted in thermoelectricity. The devices were fabricated using CMOS-compatible processes and exhibited power generation up to 5.5 mW cm-2under a temperature difference of 280 K. These numbers demonstrate the capability to power autonomous devices with environmental heat sources using silicon chips of centimeter square dimensions. We also report the possibility of using the developed devices for integrated thermoelectric cooling.

Paclobutrazol Ameliorates Low-Light-Induced Damage by Improving Photosynthesis, Antioxidant Defense System, and Regulating Hormone Levels in Tall Fescue.

Paclobutrazol (PBZ) is a plant-growth regulator (PGR) in the triazole family that enhances plant tolerance to environmental stresses. Low-light (LL) intensity is a critical factor adversely affecting the growth of tall fescue (Festuca arundinacea Schreb.). Therefore, in this study, tall fescue seedlings were treated with PBZ under control and LL conditions to investigate the effects of PBZ on enhancing LL stress resistance by regulating the growth, photosynthesis, oxidative defense, and hormone levels. Our results reveal that LL stress reduced the total biomass, chlorophyll (Chl) content, photosynthetic capacity, and photochemical efficiency of photosystem II (PSII) but increased the membrane lipid peroxidation level and reactive oxygen species (ROS) generation. However, the application of PBZ increased the photosynthetic pigment contents, net photosynthetic rate (Pn), maximum quantum yield of PSII photochemistry (Fv/Fm), ribulose-1,5-bisphosphate carboxylase (RuBisCO) activity, and starch content. In addition, PBZ treatment activated the antioxidant enzyme activities, antioxidants contents, ascorbate acid-glutathione (AsA-GSH) cycle, and related gene expression, lessening the ROS burst (H2O2 and O2∙-). However, the gibberellic acid (GA) anabolism was remarkably decreased by PBZ treatment under LL stress, downregulating the transcript levels of kaurene oxidase (KO), kaurenoic acid oxidase (KAO), and GA 20-oxidases (GA20ox). At the same time, PBZ treatment up-regulated 9-cis-epoxycarotenoid dioxygenase (NCED) gene expression, significantly increasing the endogenous abscisic acid (ABA) concentration under LL stress. Thus, our study revealed that PBZ improves the antioxidation and photosynthetic capacity, meanwhile increasing the ABA concentration and decreasing GA concentration, which ultimately enhances the LL stress tolerance in tall fescue.

Melatonergic systems of AANAT, melatonin, and its receptor MT2 in the corpus luteum are essential for reproductive success in mammals†.

Corpus luteum (CL) plays a critical role in mammalian reproductive physiology. Its dysfunction will lead to infertility or habitual abortion. In the current study, by use of melatonin specific membrane receptor 2 (MT2) knocking out (KO) mice model combined with RNA-Seq, immunohistochemistry, and immunofluorescence analyses, the genes of melatonin synthetic enzyme arylalkylamine N-acetyltransferase (AANAT) and MT2 were identified to strongly express in the CL of sows and mice. KO MT2 significantly impaired the reproductive performance in mice indicated by the reduced litter sizes. Melatonin treatment elevated the progesterone production in sows suggesting the improved CL function. Mechanistic analysis showed that melatonin upregulated a set of progesterone synthesis-related genes including cytochrome P450 family 11 subfamily A member 1 (Cyp11a1), aldo-keto reductase family 1, member C18 (Akr1c18), isopentenyl-diphosphate delta isomerase 1 (Idi1), and luteinizing hormone/choriogonadotropin receptor (Lhcgr). The upregulation of these genes directly related to the increased progesterone production. The regulatory effects of melatonin on these gene expressions were mediated by MT2 and MT2KO diminished the effects of melatonin in this respect. Thus, the presence of melatonergic system of AANAT, melatonin, and its receptor MT2 in CL is essential for reproductive success in mammals.

Accounting for Uncertainty in the Evolutionary Timescale of Green Plants Through Clock-Partitioning and Fossil Calibration Strategies.

Establishing an accurate evolutionary timescale for green plants (Viridiplantae) is essential to understanding their interaction and coevolution with the Earth's climate and the many organisms that rely on green plants. Despite being the focus of numerous studies, the timing of the origin of green plants and the divergence of major clades within this group remain highly controversial. Here, we infer the evolutionary timescale of green plants by analyzing 81 protein-coding genes from 99 chloroplast genomes, using a core set of 21 fossil calibrations. We test the sensitivity of our divergence-time estimates to various components of Bayesian molecular dating, including the tree topology, clock models, clock-partitioning schemes, rate priors, and fossil calibrations. We find that the choice of clock model affects date estimation and that the independent-rates model provides a better fit to the data than the autocorrelated-rates model. Varying the rate prior and tree topology had little impact on age estimates, with far greater differences observed among calibration choices and clock-partitioning schemes. Our analyses yield date estimates ranging from the Paleoproterozoic to Mesoproterozoic for crown-group green plants, and from the Ediacaran to Middle Ordovician for crown-group land plants. We present divergence-time estimates of the major groups of green plants that take into account various sources of uncertainty. Our proposed timeline lays the foundation for further investigations into how green plants shaped the global climate and ecosystems, and how embryophytes became dominant in terrestrial environments.

Bayesian selection of misspecified models is overconfident and may cause spurious posterior probabilities for phylogenetic trees.

The Bayesian method is noted to produce spuriously high posterior probabilities for phylogenetic trees in analysis of large datasets, but the precise reasons for this overconfidence are unknown. In general, the performance of Bayesian selection of misspecified models is poorly understood, even though this is of great scientific interest since models are never true in real data analysis. Here we characterize the asymptotic behavior of Bayesian model selection and show that when the competing models are equally wrong, Bayesian model selection exhibits surprising and polarized behaviors in large datasets, supporting one model with full force while rejecting the others. If one model is slightly less wrong than the other, the less wrong model will eventually win when the amount of data increases, but the method may become overconfident before it becomes reliable. We suggest that this extreme behavior may be a major factor for the spuriously high posterior probabilities for evolutionary trees. The philosophical implications of our results to the application of Bayesian model selection to evaluate opposing scientific hypotheses are yet to be explored, as are the behaviors of non-Bayesian methods in similar situations.

Genetic Load and Potential Mutational Meltdown in Cancer Cell Populations.

Large genomes with elevated mutation rates are prone to accumulating deleterious mutations more rapidly than natural selection can purge (Muller's ratchet). As a consequence, it may lead to the extinction of small populations. Relative to most unicellular organisms, cancer cells, with large and nonrecombining genome and high mutation rate, could be particularly susceptible to such "mutational meltdown." However, the most common type of mutation in organismal evolution, namely, deleterious mutation, has received relatively little attention in the cancer biology literature. Here, by monitoring single-cell clones from HeLa cell lines, we characterize deleterious mutations that retard the rate of cell proliferation. The main mutation events are copy number variations (CNVs), which, estimated from fitness data, happen at a rate of 0.29 event per cell division on average. The mean fitness reduction, estimated reaching 18% per mutation, is very high. HeLa cell populations therefore have very substantial genetic load and, at this level, natural population would likely face mutational meltdown. We suspect that HeLa cell populations may avoid extinction only after the population size becomes large enough. Because CNVs are common in most cell lines and tumor tissues, the observations hint at cancer cells' vulnerability, which could be exploited by therapeutic strategies.

The Spectre of Too Many Species.

Recent simulation studies examining the performance of Bayesian species delimitation as implemented in the bpp program have suggested that bpp may detect population splits but not species divergences and that it tends to over-split when data of many loci are analyzed. Here, we confirm these results and provide the mathematical justifications. We point out that the distinction between population and species splits made in the protracted speciation model (PSM) has no influence on the generation of gene trees and sequence data, which explains why no method can use such data to distinguish between population splits and speciation. We suggest that the PSM is unrealistic as its mechanism for assigning species status assumes instantaneous speciation, contradicting prevailing taxonomic practice. We confirm the suggestion, based on simulation, that in the case of speciation with gene flow, Bayesian model selection as implemented in bpp tends to detect population splits when the amount of data (the number of loci) increases. We discuss the use of a recently proposed empirical genealogical divergence index (gdi) for species delimitation and illustrate that parameter estimates produced by a full likelihood analysis as implemented in bpp provide much more reliable inference under the gdi than the approximate method phrapl. We distinguish between Bayesian model selection and parameter estimation and suggest that the model selection approach is useful for identifying sympatric cryptic species, while the parameter estimation approach may be used to implement empirical criteria for determining species status among allopatric populations.

Long-term outcomes of single-incision laparoscopic technique in Soave procedure compared with heart-shaped anastomosis for Hirschsprung disease.

PURPOSE: This retrospective study compared the long-term outcomes of single-incision laparoscopy-assisted Soave procedure (SILSP) with single-incision laparoscopy-assisted heart-shaped anastomosis (SILHSA) in patients with Hirschsprung disease (HSCR). METHODS: Patients diagnosed with HSCR that underwent SILSP or SILHSA between January 2009 and January 2015 at our institute were enrolled in this retrospective study. Data on the clinical characteristics, perioperative complications, and postoperative quality of life were retrospectively collected and analyzed. RESULTS: There were 109 patients in the SILSP group and 95 patients in the SILHSA group. No differences in clinical characteristics, including age, weight, hospitalization length, blood loss volume, and operation time, were noted between the two groups. The incidence rates of constipation and soiling were lower in the SILHSA group than those in the SILSP group. The SILHSA group showed lower scores in constipation and soiling compared with the SILSP group, indicating a better surgical outcome for patients receiving SILHSA procedure. CONCLUSION: SILHSA is a feasible and reliable minimally invasive surgical procedure for patients with HSCR. Patients who underwent SILHSA had lower incidence rates of constipation and soiling than patients who underwent SILSP, suggesting that SILHSA could be a better choice for patients with HSCR.

Continuing epidural analgesia during the second stage and ACOG definition of arrest of labor on maternal-fetal outcomes.

BACKGROUND: Despite an increase in the rates of epidural labor analgesia, continuation of epidural labor analgesia in the second stage of labor (CEADSSOL) was interrupted by care providers due to fears of increased risk of operative delivery and adverse neonatal outcomes. Therefore, we evaluated the effect of CEADSSOL and the newer American College of Obstetricians and Gynecologists (ACOG) definition of arrest of labor on the length of secondary stage of labor, newborn outcomes, and mode of delivery. METHODS: This is a retrospective cohort study. Data collection began during March 2014 and ended in May 2015, 1 year after implementation of both interventions. The primary outcome was the length of secondary stage of labor, mode of delivery and neonatal outcome (Apgar < 7, at 5 minutes). The implementation of continuing epidural analgesia during the second stage of labor was performed with 0.08%-0.15% ropivacaine and 0.1-0.2 µg/mL sufentanil. RESULTS: There were a total 10 414 deliveries during the study period. The length of the second stage of labor has no significant differences among groups. The cesarean delivery rate decreased 4.1% (36% vs 40.1%, P = .0038). Moreover, no significant difference was found in neonatal Apgar scores less than 7 at 5 minutes between two phases. Maternal outcomes remained unchanged. Post-intervention neonatal parameters including NICU admissions (P < .001), incidences of antibiotics usage (P < .0001), intubation (P = .0003), and 7 days mortality (P = .0020) were remarkably reduced compared to pre-interventions. CONCLUSION: The important finding of this study was the improvement in neonatal outcomes by implementing two simultaneous interventions without a cost of increased operative delivery.

The role of NMDA receptors in rat propofol self-administration.

BACKGROUND: Propofol is among the most frequently used anesthetic agents, and it has the potential for abuse. The N-methyl-D-aspartate (NMDA) receptors are key mediators neural plasticity, neuronal development, addiction, and neurodegeneration. In the present study, we explored the role of these receptors in the context of rat propofol self-administration. METHODS: Sprague-Dawley Rats were trained to self-administer propofol (1.7 mg/kg/infusion) using a fixed-ratio (FR) schedule over the course of 14 sessions (3 h/day). After training, rats were intraperitoneally administered the non-competitive NDMA receptor antagonist MK-801, followed 10 min later by a propofol self-administration session. RESULTS: After training, rats successfully underwent acquisition of propofol self-administration, as evidenced by a significant and stable rise in the number of active nose-pokes resulting in propofol administration relative to the number of control inactive nose-pokes (P < 0.01). As compared to control rats, rats that had been injected with 0.2 mg/kg MK-801 exhibited a significantly greater number of propofol infusions (F (3, 28) = 4.372, P < 0.01), whereas infusions were comparable in the groups administered 0.1 mg/kg and 0.4 mg/kg of this compound. In addition, MK-801 failed to alter the numbers of active (F (3, 28) = 1.353, P > 0.05) or inactive (F (3, 28) = 0.047, P > 0.05) responses in these study groups. Animals administered 0.4 mg/kg MK-801 exhibited significantly fewer infusions than animals administered 0.2 mg/kg MK-801 (P = 0.006, P < 0.01). In contrast, however, animals in the 0.4 mg/kg MK-801 group displayed a significant reduction in the number of active nose-poke responses (F (3, 20) = 20.8673, P < 0.01) and the number of sucrose pellets (F (3, 20) = 23.77, P < 0.01), while their locomotor activity was increased (F (3, 20) = 22.812, P < 0.01). CONCLUSION: These findings indicate that NMDA receptors may play a role in regulating rat self-administration of propofol.