Genome-wide meta-analysis for rheumatoid arthritis

Meta-analysis is being increasingly used as a tool for integrating data from different studies of complex phenotypes, because the power of any one study to identify causal loci is limited. We applied a novel meta-analytical approach (Loesgen et al. in Genet Epidemiol 21(Suppl 1):S142–S147, 2001) in compiling results from four studies of rheumatoid arthritis in Caucasians including two studies from NARAC (Jawaheer et al. in Am J Hum Genet 68:927–936, 2001; Jawaheer et al. in Arthritis Rheum 48:906–916, 2003), one study from the UK (MacKay et al. in Arthritis Rheum 46:632–639, 2001) and one from France (Cornelis et al. in Proc Natl Acad Sci USA 95:10746–10750, 1998). For each study, we obtained NPL scores by performing interval mapping (2 cM intervals) using GeneHunter2 (Kruglyak et al. in Am J Hum Genet 58:1347–1363, 1996; Markianos et al. in Am J Hum Genet 68:963–977, 2001). The marker maps differed among the three consortium groups, therefore, the marker maps were aligned after the interval mapping was completed and the NPL scores that were within 1 cM of each other were combined using the method of Loesgen et al. (Genet Epidemiol 21(Suppl 1):S142–S147, 2001) by calculating the weighted average of the NPL score. This approach avoids some problems in analysis encountered by using GeneHunter2 when some markers in the sample are not genotyped. This procedure provided marginal evidence (P<0.05) of linkage on chromosome 1, 2, 5 and 18, strong evidence (P<0.01) on chromosomes 8 and 16, and overwhelming evidence in the HLA region of chromosome 6.     

The evolutionary origins of tensed language and belief

I outline the debate in metaphysics between those who believe time is tensed and those who believe it is tenseless. I describe the terms in which this debate has been carried out, and the significance to it of ordinary tensed language and widespread common sense beliefs that time is tensed. I then outline a case for thinking that our intuitive beliefs about tense constitute an Adaptive Imaginary Representation (Wilson, in Biol Philos 5:37–62, 1990; Wilson, in Biol Philos 10:77–97, 1995). I also outline a case for thinking that our ordinary tensed beliefs and tensed language owe their tensed nature to its being adaptive to adopt a temporally self-locating perspective on reality. If these conclusions are right, then common sense intuitions and temporal language will be utterly misleading guides to the nature of temporal reality.     

mGluR7 Genetics and Alcohol: Intersection Yields Clues for Addiction

Development of addiction to alcohol or other substances can be attributed in part to exposure-dependent modifications at synaptic efficacy leading to an organism which functions at an altered homeostatic setpoint. Genetic factors may also influence setpoints and the stability of the homeostatic system of an organism. Quantitative genetic analysis of voluntary alcohol drinking, and mapping of the involved genes in the quasi-congenic Recombinant QTL Introgression strain system, identified Eac2 as a Quantitative Trait Locus (QTL) on mouse chromosome 6 which explained 18% of the variance with an effect size of 2.09 g/kg/day alcohol consumption, and Grm7 as a quantitative trait gene underlying Eac2 [Vadasz et al. in Neurochem Res 32:1099–1112, 100, Genomics 90:690–702, 102]. In earlier studies, the product of Grm7 mGluR7, a G protein-coupled receptor, has been implicated in stress systems [Mitsukawa et al. in Proc Natl Acad Sci USA 102:18712–18717, 63], anxiety-like behaviors [Cryan et al. in Eur J Neurosci 17:2409–2417, 14], memory [Holscher et al. in Learn Mem 12:450–455, 26], and psychiatric disorders (e.g., [Mick et al. in Am J Med Genet B Neuropsychiatr Genet 147B:1412–1418, 61; Ohtsuki et al. in Schizophr Res 101:9–16, 72; Pergadia et al. in Paper presented at the 38th Annual Meeting of the Behavior Genetics Association, Louisville, Kentucky, USA, 76]. Here, in experiments with mice, we show that (1) Grm7 knockout mice express increased alcohol consumption, (2) sub-congenic, and congenic mice carrying a Grm7 variant characterized by higher Grm7 mRNA drink less alcohol, and show a tendency for higher circadian dark phase motor activity in a wheel running paradigm, respectively, and (3) there are significant genetic differences in Grm7 mRNA abundance in the mouse brain between congenic and background mice identifying brain areas whose function is implicated in addiction related processes. We hypothesize that metabotropic glutamate receptors may function as regulators of homeostasis, and Grm7 (mGluR7) is involved in multiple processes (including stress, circadian activity, reward control, memory, etc.) which interact with substance use and the development of addiction. In conclusion, we suggest that mGluR7 is a significant new therapeutic target in addiction and related neurobehavioral disorders.     

Periodicity and synchronization in blood-stage malaria infection

Malaria fever is highly periodic and is associated with the parasite replication cycles in red blood cells. The existence of periodicity in malaria infection demonstrates that parasite replication in different red blood cells is synchronized. In this article, rigorous mathematical analysis of an age-structured human malaria model of infected red blood cells (Rouzine and McKenzie, Proc Natl Acad Sci USA 100:3473–3478, 2003) is provided and the synchronization of Plasmodium falciparum erythrocytic stages is investigated. By using the replication rate as the bifurcation parameter, the existence of Hopf bifurcation in the age-structured malaria infection model is obtained. Numerical simulations indicate that synchronization with regular periodic oscillations (of period 48 h) occurs when the replication rate increases. Therefore, Kwiatkowski and Nowak’s observation (Proc Natl Acad Sci USA 88:5111–5113, 1991) that synchronization could be generated at modest replication rates is confirmed.     

ROS-Mediated ABA Signaling

In plants, reactive oxygen species (ROS) are short-lived molecules produced through various cellular mechanisms in response to biotic and abiotic stimuli. ROS function as second messengers for hormone signaling, development, oxygen deprivation, programmed cell death, and plant–pathogen interactions. Recent research on ROS-mediated responses has produced stimulating findings such as the specific sources of ROS production, molecular elements that work in ROS-mediated signaling and homeostasis, and a ROS-regulated gene network (Neill et al., Curr Opin Plant Biol 5:388–395, 2002a; Apel and Hirt, Annu Rev Plant Biol 55:373–399, 2004; Mittler et al., Trends Plant Sci 9:490–498, 2004; Mori and Schroeder, Plant Physiol 135:702–708, 2004; Kwak et al., Plant Physiol 141:323–329, 2006; Torres et al., Plant Physiol 141:373–378, 2006; Miller et al., Physiol Plant 133:481–489, 2008). In this review, we highlight new discoveries in ROS-mediated abscisic acid (ABA) signaling. Drs. Daeshik Cho and June M. Kwak are the corresponding authors for this paper.     

1H, 13C and 15N NMR assignments of the C1A and C1B subdomains of PKC-delta

The Protein Kinase C family of enzymes is a group of serine/threonine kinases that play central roles in cell-cycle regulation, development and cancer. A key step in the activation of PKC is translocation to membranes and binding of membrane-associated activators including diacylglycerol (DAG). Interaction of novel and conventional isotypes of PKC with DAG and phorbol esters occurs through the two C1 regulatory domains (C1A and C1B), which exhibit distinct ligand binding selectivity that likely controls enzyme activation by different co-activators. PKC has also been implicated in physiological responses to alcohol consumption and it has been proposed that PKCα (Slater et al. J Biol Chem 272(10):6167–6173, 1997; Slater et al. Biochemistry 43(23):7601–7609, 2004), PKCε (Das et al. Biochem J 421(3):405–413, 2009) and PKCδ (Das et al. J Biol Chem 279(36):37964–37972, 2004; Das et al. Protein Sci 15(9):2107–2119, 2006) contain specific alcohol-binding sites in their C1 domains. We are interested in understanding how ethanol affects signal transduction processes through its affects on the structure and function of the C1 domains of PKC. Here we present the ¹H, ¹⁵N and ¹³C NMR chemical shift assignments for the Rattus norvegicus PKCδ C1A and C1B proteins.     

Identification of <Emphasis Type="Italic">Porphyra yezoensis</Emphasis> (Rhodophyta) meiosis by DNA quantification using confocal laser scanning microscopy

Conchospore germlings of Porphyra yezoensis were stained with a fluorescent dye for DNA and observed with confocal laser scanning microscopy (CLSM). Relative DNA values of the germling nuclei were obtained by measuring fluorescence intensities of nuclear regions of the optically sliced specimens, using the mean value of the smallest blade cells as a reference of the genomic n value. Such quantification revealed that the nuclear DNA amounts of the one-cell, two-cell, and four-cell-stage germlings are approximately 4 × n, 2 × n, and n ∼2 × n values respectively; these values agreed well with the expected ones from the hypothesis that meiosis corresponds to the first successive cell divisions after the conchospore germination. These results are consistent with a previous study on cytogenetic analysis of the chimaera blade formation (Ohme and Miura 1988, Plant Sci 57:135–140) and not consistent with a recent microscopic study (Wang et al. 2006, Phycol Res 54:201–207) which proposed that the first meiotic division occurs at the conchospore formation and the second division at the germination.     

A Computational Analysis of Localized Ca<Superscript>2+</Superscript>-Dynamics Generated by Heterogeneous Release Sites

We investigate the role of heterogeneous expression of IP₃R and RyR in generating diverse elementary Ca²⁺ signals. It has been shown empirically (Wojcikiewicz and Luo in Mol. Pharmacol. 53(4):656–662, 1998; Newton et al. in J. Biol. Chem. 269(46):28613–28619, 1994; Smedt et al. in Biochem. J. 322(Pt. 2):575–583, 1997) that tissues express various proportions of IP₃ and RyR isoforms and this expression is dynamically regulated (Parrington et al. in Dev. Biol. 203(2):451–461, 1998; Fissore et al. in Biol. Reprod. 60(1):49–57, 1999; Tovey et al. in J. Cell Sci. 114(Pt. 22):3979–3989, 2001). Although many previous theoretical studies have investigated the dynamics of localized calcium release sites (Swillens et al. in Proc. Natl. Acad. Sci. U.S.A. 96(24):13750–13755, 1999; Shuai and Jung in Proc. Natl. Acad. Sci. U.S.A. 100(2):506–510, 2003a; Shuai and Jung in Phys. Rev. E, Stat. Nonlinear Soft Matter Phys. 67(3 Pt. 1):031905, 2003b; Thul and Falcke in Biophys. J. 86(5):2660–2673, 2004; DeRemigio and Smith in Cell Calcium 38(2):73–86, 2005; Nguyen et al. in Bull. Math. Biol. 67(3):393–432, 2005), so far all such studies focused on release sites consisting of identical channel types. We have extended an existing mathematical model (Nguyen et al. in Bull. Math. Biol. 67(3):393–432, 2005) to release sites with two (or more) receptor types, each with its distinct channel kinetics. Mathematically, the release site is represented by a transition probability matrix for a collection of nonidentical stochastically gating channels coupled through a shared Ca²⁺ domain. We demonstrate that under certain conditions a previously defined mean-field approximation of the coupling strength does not accurately reproduce the release site dynamics. We develop a novel approximation and establish that its performance in these instances is superior. We use this mathematical framework to study the effect of heterogeneity in the Ca²⁺-regulation of two colocalized channel types on the release site dynamics. We consider release sites consisting of channels with both Ca²⁺-activation and inactivation (“four-state channels”) and channels with Ca²⁺-activation only (“two-state channels”) and show that for the appropriate parameter values, synchronous channel openings within a release site with any proportion of two-state to four-state channels are possible, however, the larger the proportion of two-state channels, the more sensitive the dynamics are to the exact spatial positioning of the channels and the distance between channels. Specifically, the clustering of even a small number of two-state channels interferes with puff/spark termination and increases puff durations or leads to a tonic response.     

Shaping bursting by electrical coupling and noise

Gap-junctional coupling is an important way of communication between neurons and other excitable cells. Strong electrical coupling synchronizes activity across cell ensembles. Surprisingly, in the presence of noise synchronous oscillations generated by an electrically coupled network may differ qualitatively from the oscillations produced by uncoupled individual cells forming the network. A prominent example of such behavior is the synchronized bursting in islets of Langerhans formed by pancreatic β-cells, which in isolation are known to exhibit irregular spiking (Sherman and Rinzel, Biophys J 54:411–425, 1988; Sherman and Rinzel, Biophys J 59:547–559, 1991). At the heart of this intriguing phenomenon lies denoising, a remarkable ability of electrical coupling to diminish the effects of noise acting on individual cells. In this paper, building on an earlier analysis of denoising in networks of integrate-and-fire neurons (Medvedev, Neural Comput 21 (11):3057–3078, 2009) and our recent study of spontaneous activity in a closely related model of the Locus Coeruleus network (Medvedev and Zhuravytska, The geometry of spontaneous spiking in neuronal networks, submitted, 2012), we derive quantitative estimates characterizing denoising in electrically coupled networks of conductance-based models of square wave bursting cells. Our analysis reveals the interplay of the intrinsic properties of the individual cells and network topology and their respective contributions to this important effect. In particular, we show that networks on graphs with large algebraic connectivity (Fiedler, Czech Math J 23(98):298–305, 1973) or small total effective resistance (Bollobas, Modern graph theory, Graduate Texts in Mathematics, vol. 184, Springer, New York, 1998) are better equipped for implementing denoising. As a by-product of the analysis of denoising, we analytically estimate the rate with which trajectories converge to the synchronization subspace and the stability of the latter to random perturbations. These estimates reveal the role of the network topology in synchronization. The analysis is complemented by numerical simulations of electrically coupled conductance-based networks. Taken together, these results explain the mechanisms underlying synchronization and denoising in an important class of biological models.     

Loss of <Emphasis Type="Italic">Tc-arrow</Emphasis> and canonical Wnt signaling alters posterior morphology and pair-rule gene expression in the short-germ insect, <Emphasis Type="Italic">Tribolium castaneum</Emphasis>

Wnt signaling has been implicated in posterior patterning in short-germ insects, including the red flour beetle Tribolium castaneum (Bolognesi et al. Curr Biol 18:1624–1629, 2008b; Angelini and Kaufman Dev Biol 283:409–423, 2005; Miyawaki et al. Mech Dev 121:119–130, 2004). Specifically, depletion of Wnt ligands Tc-Wnt1 and Tc-WntD/8 produces Tribolium embryos lacking abdominal segments. Similar phenotypes are produced by depletion of Tc-porcupine (Tc-porc) or Tc-pangolin (Tc-pan), indicating that the signal is transmitted through the canonical Wnt pathway (Bolognesi et al. Curr Biol 18:1624–1629, 2008b). Here we show that RNAi for the receptor Tc-arrow produced similar truncated phenotypes, providing additional evidence supporting canonical signal transduction. Furthermore, since in Tribolium segments are defined sequentially by a pair-rule gene circuit that, when interrupted, produces truncated phenotypes (Choe et al. Proc Natl Acad Sci U S A 103:6560–6564, 2006), we investigated the relationship between loss of Wnt signaling and this pair-rule gene circuit. After depletion of the receptor Tc-arrow, expression of Tc-Wnt1 was noticeably absent from the growth zone, while Tc-WntD/8 was restricted to a single spot of expression in what remained of the posterior growth zone. The primary pair-rule genes Tc-runt (Tc-run) and Tc-even-skipped (Tc-eve) were expressed normally in the anterior segments, but were reduced to a single spot in the remnants of the posterior growth zone. Thus, expression of pair-rule genes and Tc-WntD/8 are similarly affected by depletion of Wnt signal and disruption of the posterior growth zone.     

Existence and stability of steady states of a reaction convection diffusion equation modeling microtubule formation

We generalize the Dogterom–Leibler model for microtubule dynamics (Dogterom and Leibler in Phys Rev Lett 70(9):1347–1350, 1993) to the case where the rates of elongation as well as the lifetimes of the elongating shortening phases are a function of GTP-tubulin concentration. We analyze also the effect of nucleation rate in the form of a damping term which leads to new steady-states. For this model, we study existence and stability of steady states satisfying the boundary conditions at x = 0. Our stability analysis introduces numerical and analytical Evans function computations as a new mathematical tool in the study of microtubule dynamics.     

Pleiotropic effect, clock genes, and reproductive isolation

 The mechanism by which a clock gene pleiotropically controlling life history and behavioral traits causes reproductive isolation is explained using a model species, the melon fly, Bactrocera cucurbitae (Coquillett) (Diptera: Tephritidae). Melon flies mate once a day, at dusk. The population selected for life history traits exhibits correlated responses in the time of mating during the day. For example, the fly populations selected for faster (slower) development have an earlier (later) time of mating. A circadian rhythm controls the time of mating. The circadian periods in constant darkness were about 22 h in lines selected for a short developmental period and about 31 h in lines selected for a long developmental period. The data on crosses between the selected lines indicated that the developmental period is controlled by a polygene, whereas the circadian period may be controlled by a single clock gene. These results suggest a clock gene pleiotropically controls developmental and circadian periods in the melon fly. Reproductive isolation may often evolve as an indirect (pleiotropic) consequence of adaptation to different environments or habitats. For example, niches that are temporally or seasonally offset can select organisms with different developmental characteristics. These developmental differences can inadvertently cause reproductive isolation by a variety of means including shifts in mating activity patterns. The difference in time of mating between populations selected for developmental period translated into significant prezygotic isolation, as measured by mate choice tests. If the mating time between populations differed more than 1 h, the isolation index was significantly higher than zero. These findings indicate that premating isolation can be established by a pleiotropic effect of a clock gene. There are many examples in which the difference in timing of reproduction prevents gene flow between populations, such as the egg spawning time in marine organisms, the flowering time in angiosperms, and the time of mating in insects. In such organisms, if genetic correlations between circadian rhythm and reproductive traits exist, multifarious divergent selection for life history traits would often accelerate the evolution of reproductive isolation through clock genes. Natural populations may diverge in reproduction time through drift, direct natural selection for time of reproduction, or as a by-product effect of genetic correlations. In any case, clock genes are keys in reproductive isolation. Received: January 31, 2002 / Accepted: July 29, 2002 Acknowledgments I am grateful to Tetsuo Arai, Akira Matsumoto, Takashi Matsuyama, Toru Shimizu, Aya Takahashi, Teiichi Tanimura, Tetsuya Toyosato, and Yasuhiko Watari for useful discussion, and to the responsible editor and two anonymous reviewers for helpful suggestions. I also thank Yoshihiko Chiba, Norio Ishida, Emi Koyama, Kazuhiko Sakai, and Takaomi Sakai for useful information. My work on speciation has been supported by a Grant-in-Aid for Scientific Research (KAKENHI 14340244) from the Ministry of Education, Culture, Sports, Science and Technology of Japan.     

Adding Adhesion to a Chemical Signaling Model for Somite Formation

Somites are condensations of mesodermal cells that form along the two sides of the neural tube during early vertebrate development. They are one of the first instances of a periodic pattern, and give rise to repeated structures such as the vertebrae. A number of theories for the mechanisms underpinning somite formation have been proposed. For example, in the “clock and wavefront” model (Cooke and Zeeman in J. Theor. Biol. 58:455–476, 1976), a cellular oscillator coupled to a determination wave progressing along the anterior-posterior axis serves to group cells into a presumptive somite. More recently, a chemical signaling model has been developed and analyzed by Maini and coworkers (Collier et al. in J. Theor. Biol. 207:305–316, 2000; Schnell et al. in C. R. Biol. 325:179–189, 2002; McInerney et al. in Math. Med. Biol. 21:85–113, 2004), with equations for two chemical regulators with entrained dynamics. One of the chemicals is identified as a somitic factor, which is assumed to translate into a pattern of cellular aggregations via its effect on cell–cell adhesion. Here, the authors propose an extension to this model that includes an explicit equation for an adhesive cell population. They represent cell adhesion via an integral over the sensing region of the cell, based on a model developed previously for adhesion driven cell sorting (Armstrong et al. in J. Theor. Biol. 243:98–113, 2006). The expanded model is able to reproduce the observed pattern of cellular aggregates, but only under certain parameter restrictions. This provides a fuller understanding of the conditions required for the chemical model to be applicable. Moreover, a further extension of the model to include separate subpopulations of cells is able to reproduce the observed differentiation of the somite into separate anterior and posterior halves. N.J. Armstrong was supported by a Doctoral Training Account Studentship from EPSRC. K.J. Painter and J.A. Sherratt were supported in part by Integrative Cancer Biology Program Grant CA113004 from the US National Institute of Health and in part by BBSRC grant BB/D019621/1 for the Centre for Systems Biology at Edinburgh.     

Comparative mapping of the oat <Emphasis Type="Italic">Dw6/dw6</Emphasis> dwarfing locus using NILs and association with vacuolar proton ATPase subunit H

Seven pairs of oat near-isogenic lines (NILs) (Kibite in Crop Sci 41:277–278, 2001) contrasting for the Dw6 dwarfing gene were used to test for correlation between tall/dwarf phenotype and polymorphic genotype using restriction fragment length polymorphism (RFLP) and other molecular markers selected from the Kanota × Ogle (K×O) (Wight et al. in Genome 46:28–47, 2003) and Terra × Marion (De Koeyer et al. in Theor Appl Genet 108:1285–1298, 2004) recombination maps. This strategy located the Dw6/dw6 locus to a small chromosomal region on K×O linkage group (LG) KO33, near or at a putative RFLP locus aco245z. Aco245z and other tightly linked flanking markers have potential for use in marker-assisted selection (MAS), and PCR-based markers were developed from several of these. RFLP genotyping of the Dw6 NILs indicated that 13 of the 14 individual lines were homogeneously maternal or paternal for a large genomic region near Dw6/dw6, an unexpected result for NILs. The cDNA clone aco245 codes for a vacuolar proton ATPase subunit H, a potential candidate gene for Dw6. Vacuolar proton ATPase enzymes have a central role in plant growth and development and a mutation in subunit C is responsible for the det3 dwarfing mutation in Arabidopsis thaliana (Schumacher et al. in Genes Dev 13:3259–3270, 1999). Aco245 affords the potential of designing highly precise diagnostic markers for MAS for Dw6. The Dw6 NILs have potential utility to investigate the role of vacuolar proton ATPases in growth and development in plants.     

Progress in erigonine spider phylogeny—the <Emphasis Type="Italic">Savignia</Emphasis>-group is not monophyletic (Araneae: Linyphiidae)

We present the most inclusive study on the higher-level phylogeny of erigonine spiders, including about 30% of all erigonine genera. By expanding the previously most comprehensive analysis (Miller and Hormiga Cladistics 20:385–442, 2004) we tested the robustness of its results to the addition of closely related taxa, and also the monophyly of the Savignia-group defined by Millidge (Bulletin of the British Arachnological Society 4:1–60, 1977). The character matrix was expanded by adding 18 newly scored species in 15 genera, and also includes all species scored by other authors. This adds up to 98 species in 91 erigonine genera plus 13 outgroup taxa. The parsimony analysis led to eight fully resolved most parsimonious trees (L=1084). The topology concerning the taxa basal to the ‘distal erigonines’ remained unchanged, and the latter clade still shares 67% of all nodes with the original analysis. The Savignia-group is not monophyletic at genus level with respect to Saloca diceros and Alioranus pastoralis, and the same applies at species level in Diplocephalus and Erigonella. From the Savignia-group, Glyphesis servulus, Diplocephalus cristatus, Savignia frontata, and two representatives each of Erigonella, Dicymbium and Araeoncus combine to form a monophyletic clade.     

Competing with oneself: introducing self-interaction in a model of competitive learning

A competitive learning model was introduced in Mehta and Luck (Phys Rev E 60, 5:5218–5230, 1999), in which the learning is outcome-related. Every individual chooses between a pair of existing strategies or types, guided by a combination of two factors: tendency to conform to the local majority, and a preference for the type with higher perceived success among its neighbors, based on their relative outcomes. Here, an extension of the interfacial model of Mehta and Luck (Phys Rev E 60, 5:5218–5230, 1999) is proposed, in which individuals additionally take into account their own outcomes in arriving at their outcome-based choices. Three possible update rules for handling bulk sites are considered. The corresponding phase diagrams, obtained at coexistence, show systematic departures from the original interfacial model. Possible relationships of these variants with the cooperative model of Mehta and Luck (Phys Rev E 60, 5:5218–5230, 1999) are also touched upon.     

Mutational analysis of wheat (<Emphasis Type="Italic">Triticum aestivum</Emphasis> L.) nucleotide pyrophosphatase/phosphodiesterase shows the role of six amino acids in the catalytic mechanism

Nucleotide pyrophosphatases/phosphodiesterases (NPPs, PF01663) release nucleoside 5′-monophosphates from a wide range of nucleotide substrates. Only very recently, the first plant members of the NPP family were characterised (Joye et al. J Cereal Sci 51: 326–336, 2010), and little is known about their substrate-specifying residues. We elucidated the role of six amino acid residues of the recently identified and characterised Triticum aestivum L. NPP (Joye et al. J Cereal Sci 51: 326–336, 2010). Substitution of the highly conserved catalytic Thr132 into Ser or Ala completely abolished enzyme activity. Mutation of a highly conserved His255 residue into an apolar Ala suprisingly increased enzyme activity against most phosphodiester substrates. Four other residues moderately to highly conserved over NPPs of different organisms were studied as well. Mutation of the Asn153, Asn165 and Glu199 into an Arg, Ser and Asp residue, respectively, increased the relative enzyme activity against p-nitrophenyl phosphate. Furthermore, mutation of Phe194 into Ser increased the relative enzyme activity against adenosine 5′-monophosphate-containing substrates, although the overall enzyme activity of this mutant enzyme decreased. We conclude that the structural requirements and the conservation of the amino acids of the catalytic site of TaNPPr and, by extension, probably of all NPPs, are very stringent.