Emergent global patterns of ecosystem structure and function from a mechanistic general ecosystem model

Anthropogenic activities are causing widespread degradation of ecosystems worldwide, threatening the ecosystem services upon which all human life depends. Improved understanding of this degradation is urgently needed to improve avoidance and mitigation measures. One tool to assist these efforts is predictive models of ecosystem structure and function that are mechanistic: based on fundamental ecological principles. Here we present the first mechanistic General Ecosystem Model (GEM) of ecosystem structure and function that is both global and applies in all terrestrial and marine environments. Functional forms and parameter values were derived from the theoretical and empirical literature where possible. Simulations of the fate of all organisms with body masses between 10 µg and 150,000 kg (a range of 14 orders of magnitude) across the globe led to emergent properties at individual (e.g., growth rate), community (e.g., biomass turnover rates), ecosystem (e.g., trophic pyramids), and macroecological scales (e.g., global patterns of trophic structure) that are in general agreement with current data and theory. These properties emerged from our encoding of the biology of, and interactions among, individual organisms without any direct constraints on the properties themselves. Our results indicate that ecologists have gathered sufficient information to begin to build realistic, global, and mechanistic models of ecosystems, capable of predicting a diverse range of ecosystem properties and their response to human pressures

Social determinants of common metabolic risk factors (high blood pressure, high blood sugar, high body mass index and high waist-hip ratio) of major non-communicable diseases in South Asia region:a systematic review protocol

BackgroundPrevalence of non-communicable diseases has been increasing at a greater pace in developing countries and, in particular, the South Asia region. Various behavioral, social and environmental factors present in this region perpetuate common metabolic risk factors of non-communicable diseases. This study will identify social determinants of common metabolic risk factors of major non-communicable diseases in the context of the South Asian region and map their causal pathway.MethodsA systematic review of selected articles will be carried out following Cochrane guidelines. Review will be guided by Social Determinants of Health Framework developed by the World Health Organization to extract social determinants of metabolic risk factors of non-communicable diseases from studies. A distinct search strategy will be applied using key words to screen relevant studies from online databases. Primary and grey literature published from the year 2000 to 2016 and studies with discussion on proximal and distal determinants of non-communicable risk factors among adults of the South Asia region will be selected. They will be further checked for quality, and a matrix illustrating contents of selected articles will be developed. Thematic content analysis will be done to trace social determinants and their interaction with metabolic risk factors. Findings will be illustrated in causal loop diagrams with social determinants of risk factors along with their interaction (feedback mechanism).DiscussionThe review will describe the interplay of social determinants of common NCD metabolic risk factors in the form of causal loop diagram. Findings will be structured in two parts: the first part will explain the linkage between proximal determinants with the metabolic risk factors and the second part will describe the linkage among the risk factors, proximal determinants and distal determinants. Evidences across different regions will be discussed to compare and validate and/or contrast the findings. Possible bias and limitations of this study will also be discussed.Systematic review registrationPROSPERO CRD42017067212 <br/

Nicotiana benthamiana as a Production Platform for Artemisinin Precursors

Background Production of pharmaceuticals in plants provides an alternative for chemical synthesis, fermentation or natural sources. Nicotiana benthamiana is deployed at commercial scale for production of therapeutic proteins. Here the potential of this plant is explored for rapid production of precursors of artemisinin, a sesquiterpenoid compound that is used for malaria treatment. Methodology/Principal Findings Biosynthetic genes leading to artemisinic acid, a precursor of artemisinin, were combined and expressed in N. benthamiana by agro-infiltration. The first committed precursor of artemisinin, amorpha-4,11-diene, was produced upon infiltration of a construct containing amorpha-4,11-diene synthase, accompanied by 3-hydroxy-3-methylglutaryl-CoA reductase and farnesyl diphosphate synthase. Amorpha-4,11-diene was detected both in extracts and in the headspace of the N. benthamiana leaves. When the amorphadiene oxidase CYP71AV1 was co-infiltrated with the amorphadiene-synthesizing construct, the amorpha-4,11-diene levels strongly decreased, suggesting it was oxidized. Surprisingly, no anticipated oxidation products, such as artemisinic acid, were detected upon GC-MS analysis. However, analysis of leaf extracts with a non-targeted metabolomics approach, using LC-QTOF-MS, revealed the presence of another compound, which was identified as artemisinic acid-12-ß-diglucoside. This compound accumulated to 39.5 mg.kg-1 fwt. Apparently the product of the heterologous pathway that was introduced, artemisinic acid, is further metabolized efficiently by glycosyl transferases that are endogenous to N. benthamiana. Conclusion/Significance This work shows that agroinfiltration of N. bentamiana can be used as a model to study the production of sesquiterpenoid pharmaceutical compounds. The interaction between the ectopically introduced pathway and the endogenous metabolism of the plant is discussed

Species Used for Drug Testing Reveal Different Inhibition Susceptibility for 17beta-Hydroxysteroid Dehydrogenase Type 1

Steroid-related cancers can be treated by inhibitors of steroid metabolism. In searching for new inhibitors of human 17beta-hydroxysteroid dehydrogenase type 1 (17β-HSD 1) for the treatment of breast cancer or endometriosis, novel substances based on 15-substituted estrone were validated. We checked the specificity for different 17β-HSD types and species. Compounds were tested for specificity in vitro not only towards recombinant human 17β-HSD types 1, 2, 4, 5 and 7 but also against 17β-HSD 1 of several other species including marmoset, pig, mouse, and rat. The latter are used in the processes of pharmacophore screening. We present the quantification of inhibitor preferences between human and animal models. Profound differences in the susceptibility to inhibition of steroid conversion among all 17β-HSDs analyzed were observed. Especially, the rodent 17β-HSDs 1 were significantly less sensitive to inhibition compared to the human ortholog, while the most similar inhibition pattern to the human 17β-HSD 1 was obtained with the marmoset enzyme. Molecular docking experiments predicted estrone as the most potent inhibitor. The best performing compound in enzymatic assays was also highly ranked by docking scoring for the human enzyme. However, species-specific prediction of inhibitor performance by molecular docking was not possible. We show that experiments with good candidate compounds would out-select them in the rodent model during preclinical optimization steps. Potentially active human-relevant drugs, therefore, would no longer be further developed. Activity and efficacy screens in heterologous species systems must be evaluated with caution

Enhanced Notch Activation Is Advantageous but Not Essential for T Cell Lymphomagenesis in Id1 Transgenic Mice

T cell lymphoblastic leukemia (T-ALL) is known to be associated with chromosomal abnormalities that lead to aberrant expression of a number of transcription factors such as TAL1, which dimerizes with basic helix-loop-helix (bHLH) E proteins and inhibits their function. Activated Notch receptors also efficiently induce T cell leukemogenesis in mouse models. Interestingly, gain-of-function mutations or cryptic transcription initiation of the Notch1 gene have been frequently found in both human and mouse T-ALL. However, the correlations between these alterations and overall Notch activities or leukemogenesis have not been thoroughly evaluated. Therefore, we made use of our collection of T cell lymphomas developed in transgenic mice expressing Id1, which like TAL1, inhibits E protein function. By comparing expression levels of Notch target genes in Id1-expressing tumors to those in tumors induced by a constitutively active form of Notch1, N1C, we were able to assess the overall activities of Notch pathways and conclude that the majority of Id1-expressing tumors had elevated Notch function to a varying degree. However, 26% of the Id1-expressing tumors had no evidence of enhanced Notch activation, but that did not delay the onset of tumorigenesis. Furthermore, we examined the genetic or epigenetic alterations thought to contribute to ligand-independent activation or protein stabilization of Notch1 and found that some of the Id1-expressing tumors acquired these changes, but they are not uniformly associated with elevated Notch activities in Id1 tumor samples. In contrast, N1C-expressing tumors do not harbor any PEST domain mutations nor exhibit intragenic transcription initiation. Taken together, it appears that Notch activation provides Id1-expressing tumor cells with selective advantages in growth and survival. However, this may not be absolutely essential for lymphomagenesis in Id1 transgenic mice and additional factors could also cooperate with Id1 to induce T cell lymphoma. Therefore, a broad approach is necessary in designing T-ALL therapy

Home on the Range: Factors Explaining Partial Migration of African Buffalo in a Tropical Environment

Partial migration (when only some individuals in a population undertake seasonal migrations) is common in many species and geographical contexts. Despite the development of modern statistical methods for analyzing partial migration, there have been no studies on what influences partial migration in tropical environments. We present research on factors affecting partial migration in African buffalo (Syncerus caffer) in northeastern Namibia. Our dataset is derived from 32 satellite tracking collars, spans 4 years and contains over 35,000 locations. We used remotely sensed data to quantify various factors that buffalo experience in the dry season when making decisions on whether and how far to migrate, including potential man-made and natural barriers, as well as spatial and temporal heterogeneity in environmental conditions. Using an information-theoretic, non-linear regression approach, our analyses showed that buffalo in this area can be divided into 4 migratory classes: migrants, non-migrants, dispersers, and a new class that we call “expanders”. Multimodel inference from least-squares regressions of wet season movements showed that environmental conditions (rainfall, fires, woodland cover, vegetation biomass), distance to the nearest barrier (river, fence, cultivated area) and social factors (age, size of herd at capture) were all important in explaining variation in migratory behaviour. The relative contributions of these variables to partial migration have not previously been assessed for ungulates in the tropics. Understanding the factors driving migratory decisions of wildlife will lead to better-informed conservation and land-use decisions in this area

Transglutaminase 6: a protein associated with central nervous system development and motor function.

Transglutaminases (TG) form a family of enzymes that catalyse various post-translational modifications of glutamine residues in proteins and peptides including intra- and intermolecular isopeptide bond formation, esterification and deamidation. We have characterized a novel member of the mammalian TG family, TG6, which is expressed in a human carcinoma cell line with neuronal characteristics and in mouse brain. Besides full-length protein, alternative splicing results in a short variant lacking the second β-barrel domain in man and a variant with truncated β-sandwich domain in mouse. Biochemical data show that TG6 is allosterically regulated by Ca(2+) and guanine nucleotides. Molecular modelling indicates that TG6 could have Ca(2+) and GDP-binding sites related to those of TG3 and TG2, respectively. Localization of mRNA and protein in the mouse identified abundant expression of TG6 in the central nervous system. Analysis of its temporal and spatial pattern of induction in mouse development indicates an association with neurogenesis. Neuronal expression of TG6 was confirmed by double-labelling of mouse forebrain cells with cell type-specific markers. Induction of differentiation in mouse Neuro 2a cells with NGF or dibutyryl cAMP is associated with an upregulation of TG6 expression. Familial ataxia has recently been linked to mutations in the TGM6 gene. Autoantibodies to TG6 were identified in immune-mediated ataxia in patients with gluten sensitivity. These findings suggest a critical role for TG6 in cortical and cerebellar neurons

Body mass index and outcome in renal transplant recipients: a systematic review and meta-analysis

BACKGROUND: Whether overweight or obese end stage renal disease (ESRD) patients are suitable for renal transplantation (RT) is often debated. The objective of this review and meta-analysis was to systematically investigate the outcome of low versus high BMI recipients after RT. METHODS: Comprehensive searches were conducted in MEDLINE OvidSP, Web of Science, Google Scholar, Embase, and CENTRAL (the Cochrane Library 2014, issue 8). We reviewed four major guidelines that are available regarding (potential) RT recipients. The methodology was in accordance with the Cochrane Handbook for Systematic Reviews of Interventions and written based on the PRISMA statement. The quality assessment of studies was performed by using the GRADE tool. A meta-analysis was performed using Review Manager 5.3. Random-effects models were used. RESULTS: After identifying 5,526 studies addressing this topic, 56 studies were included. We extracted data for 37 outcome measures (including data of more than 209,000 RT recipients), of which 26 could be meta-analysed. The following outcome measures demonstrated significant differences in favour of low BMI (<30) recipients: mortality (RR = 1.52), delayed graft function (RR = 1.52), acute rejection (RR = 1.17), 1-, 2-, and 3-year graft survival (RR = 0.97, 0.95, and 0.97), 1-, 2-, and 3-year patient survival (RR = 0.99, 0.99, and 0.99), wound infection and dehiscence (RR = 3.13 and 4.85), NODAT (RR = 2.24), length of hospital stay (2.31 days), operation duration (0.77 hours), hypertension (RR = 1.35), and incisional hernia (RR = 2.72). However, patient survival expressed in hazard ratios was in significant favour of high BMI recipients. Differences in other outcome parameters were not significant. CONCLUSIONS: Several of the pooled outcome measurements show significant benefits for ‘low’ BMI (<30) recipients. Therefore, we postulate that ESRD patients with a BMI >30 preferably should lose weight prior to RT. If this cannot be achieved with common measures, in morbidly obese RT candidates, bariatric surgery could be considered. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s12916-015-0340-5) contains supplementary material, which is available to authorized users

HAD hydrolase function unveiled by substrate screening: enzymatic characterization of Arabidopsis thaliana subclass I phosphosugar phosphatise AtSgpp

[EN] This work presents the isolation and the biochemical characterization of the Arabidopsis thaliana gene AtSgpp. This gene shows homology with the Arabidopsis low molecular weight phosphatases AtGpp1 and AtGpp2 and the yeast counterpart GPP1 and GPP2, which have a high specificity for dl-glycerol-3-phosphate. In addition, it exhibits homology with DOG1 and DOG2 that dephosphorylate 2-deoxy-d-glucose-6-phosphate. Using a comparative genomic approach, we identified the AtSgpp gene as a conceptual translated haloacid dehalogenase-like hydrolase HAD protein. AtSgpp (locus tag At2g38740), encodes a protein with a predicted Mw of 26.7 kDa and a pI of 4.6. Its sequence motifs and expected structure revealed that AtSgpp belongs to the HAD hydrolases subfamily I, with the C1-type cap domain. In the presence of Mg2+ ions, the enzyme has a phosphatase activity over a wide range of phosphosugars substrates (pH optima at 7.0 and K (m) in the range of 3.6-7.7 mM). AtSgpp promiscuity is preferentially detectable on d-ribose-5-phosphate, 2-deoxy-d-ribose-5-phosphate, 2-deoxy-d-glucose-6-phosphate, d-mannose-6-phosphate, d-fructose-1-phosphate, d-glucose-6-phosphate, dl-glycerol-3-phosphate, and d-fructose-6-phosphate, as substrates. AtSgpp is ubiquitously expressed throughout development in most plant organs, mainly in sepal and guard cell. Interestingly, expression is affected by abiotic and biotic stresses, being the greatest under Pi starvation and cyclopentenone oxylipins induction. Based on both, substrate lax specificity and gene expression, the physiological function of AtSgpp in housekeeping detoxification, modulation of sugar-phosphate balance and Pi homeostasis, is provisionally assigned.We acknowledge Professors Montserrat Pages (CSIC Barcelona, Spain), Thomas Kupke (University of Heidelberg, Germany) and Manuel Hernandez (University Polytechnic of Valencia, Spain) for their warm support. We also thank the advice and provision of plasmid pSBETa by Dr. Florence Vignols and Yves Meyer (University of Perpignan, France); the computer software helps by Ramon Nogales-Rangel and Alexis Gonzalez-Policarpo; Eugenio Grau-Ferrando for kind advice and help for sequencing. This work was funded by the 10 month research contract MEC-FEDER to J.A.C.-M., 10 month research contract JAE-DOC to I.M.-S. and by the research project BIO2006-10138 from the MEC-FEDER of Spain to F.A.C.-M. In memoriam of Dr. Mari Cruz Cutanda-Perez.Caparrós Martín, JA.; Mccarthy Suarez, I.; Culiañez Macia, FA. (2013). HAD hydrolase function unveiled by substrate screening: enzymatic characterization of Arabidopsis thaliana subclass I phosphosugar phosphatise AtSgpp. Planta. 237(4):943-954. https://doi.org/10.1007/s00425-012-1809-5S9439542374Allen KN, Dunaway-Mariano D (2004) Phosphoryl group transfer: evolution of a catalytic scaffold. Trends Biochem Sci 29:495–503Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ (1990) Basic local alignment search tool. J Mol Biol 215:403–410Ames BN (1966) Assay of inorganic phosphate, total phosphate, and phosphatases. Methods Enzymol 8:115–118Böhmer M, Schroeder JI (2011) Quantitative transcriptomic analysis of abscisic acid-induced and reactive oxygen species-dependent expression changes and proteomic profiling in Arabidopsis suspension cells. Plant J 67:105–118Bradford MM (1976) A rapid and sensitive method for the quantization of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem 72:248–254Burroughs AM, Allen KN, Dunaway-Mariano D, Aravind L (2006) Evolutionary genomics of the HAD superfamily: understanding the structural adaptations and catalytic diversity in a superfamily of phosphoesterases and allied enzymes. J Mol Biol 361:1003–1034Caparrós-Martín JA, Reiland S, Köchert K, Cutanda MC, Culiáñez-Macia FA (2007) Arabidopsis thaliana AtGpp 1 and AtGpp2: two novel low molecular weight phosphatases involved in plant glycerol metabolism. Plant Mol Biol 63:505–517Collet JF, Stroobant V, Pirard M, Delpierre G, Van Schaftingen E (1998) A new class of phosphotransferases phosphorylated on an aspartate residue in an amino-terminal DXDX(T/V) motif. J Biol Chem 273:14107–14112Corpet F, Servantm F, Gouzy J, Kahn D (2000) ProDom and ProDom-CG: tools for protein domain analysis and whole genome comparisons. Nucleic Acids Res 28:267–269Cutanda MC (2003) Effect of altering levels of hexoses phosphate in carbohydrate metabolism and glucose signalling in yeast and plants. PhD thesis, Polytechnic University of Valencia, Valencia, SpainHiggins D, Thompson J, Gibson T, Thompson JD, Higgins DG, Gibson TJ (1994) CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22:4673–4680Koonin EV, Tatusov RL (1994) Computer analysis of bacterial haloacid dehalogenases defines a large superfamily of hydrolases with diverse specificity. Application of an iterative approach to database search. J Mol Biol 244:125–132Kupke T, Caparrós-Martín JA, Malquichagua Salazar KJ, Culiàñez-Macià FA (2009) Biochemical and physiological characterization of Arabidopsis thaliana AtCoAse: a Nudix CoA hydrolyzing protein that improves plant development. Physiol Plant 135:365–378Kuznetsova E, Proudfoot M, Sanders SA, Reinking J, Savchenko A, Arrowsmith CH, Edwards AM, Yakunin AF (2005) Enzyme genomics: application of general enzymatic screens to discover new enzymes. FEMS Microbiol Rev 29:263–279Kuznetsova E, Proudfoo M, Gonzalez CF, Brown G, Omelchenko MV, Borozan I, Carmel L, Wolf YI, Mori H, Savchenko AV, Arrowsmith CH, Koonin EV, Edwards AM, Yakunin AF (2006) Genome-wide analysis of substrate specificities of the Escherichia coli haloacid dehalogenase-like phosphatase family. J Biol Chem 281:36149–36161Lahiri SD, Zhang G, Dai J, Dunaway-Mariano D, Allen KN (2004) Analysis of the substrate specificity loop of the HAD superfamily cap domain. Biochemistry 43:2812–2820Lahiri SD, Zhang G, Dunaway-Mariano D, Allen KN (2006) Diversification of function in the haloacid dehalogenase enzyme superfamily: the role of the cap domain in hydrolytic phosphorus—carbon bond cleavage. Bioorganic Chem 34:394–409Lambert C, Leonard N, De Bolle X, Depiereux E (2002) ESyPred3D: prediction of proteins 3D structures. Bioinformatics 18:1250–1256Lu Z, Dunaway-Mariano D, Allen KN (2005) HAD superfamily phosphotransferase substrate diversification: structure and function analysis of HAD subclass IIB sugar phosphatase BT4131. Biochemistry 44:8684–8696Lu Z, Dunaway-Mariano D, Allen KN (2008) The catalytic scaffold of the haloalkanoic acid dehalogenase enzyme superfamily acts as a mold for the trigonal bipyramidal transition state. Proc Natl Acad Sci USA 105:5687–5692Maniatis T, Fritsch EF, Sambrook J (1982) Molecular cloning: a laboratory manual. Cold Spring Harbor Laboratory Press, Cold Spring HarborMorais MC, Zhang W, Baker AS, Zhang G, Dunaway-Mariano D, Allen KN (2000) The crystal structure of Bacillus cereus phosphonoacetaldehyde hydrolase: insight into catalysis of phosphorus bond cleavage and catalytic diversification within the HAD enzyme superfamily. Biochemistry 39:10385–10396Morais MC, Zhang G, Zhang W, Olsen DB, Dunaway-Mariano D, Allen KN (2004) X-ray crystallographic and site-directed mutagenesis analysis of the mechanism of Schiff-base formation in phosphonoacetaldehyde hydrolase catalysis. J Biol Chem 279:9353–9361Mueller WS, Hilbert B, Dueckershoff K, Roitsch T, Krischke M, Mueller MJ, Berger S (2008) General detoxification and stress responses are mediated by oxidized lipids through TGA transcription factors in Arabidopsis. Plant Cell 20:768–785Murashige T, Skoog F (1962) A revised medium for rapid growth and bioassays with tobacco cultures. Physiol Plant 15:473–497Norbeck J, Pahlman AK, Akhtar N, Blomberg A, Adler L (1996) Purification and characterization of two isoenzymes of dl-glycerol-3-phosphatase from Saccharomyces cerevisiae. Identification of the corresponding GPP1 and GPP2 genes and evidence for osmotic regulation of Gpp 2p expression by the osmosensing mitogen-activated protein kinase signal transduction pathway. J Biol Chem 271:13875–13881Rández-Gil F, Blasco A, Prieto JA, Sanz P (1995) DOGR1 and DOGR2: two genes from Saccharomyces cerevisiae that confer 2-deoxyglucose resistance when overexpressed. Yeast 11:1233–1240Rao KN, Kumaran D, Seetharaman J, Bonanno JB, Burley SK, Swaminathan S (2006) Crystal structure of trehalose-6-phosphate phosphatase-related protein: biochemical and biological implications. Protein Sci 15:1735–1744Schagger H, von Jagow G (1987) Tricine-sodium dodecyl sulfatepolyacrylamide gel electrophoresis for the separation of proteins in the range from 1 to 100 kDa. Anal Biochem 166:368–379Schenk PM, Baumann S, Mattes R, Steinbiss HH (1995) Improved high-level expression system for eukaryotic genes in Escherichia coli using T7 RNA polymerase and rare ArgtRNAs. Biotechniques 19:196–200Selengut JD (2001) MDP-1 is a new and distinct member of the haloacid dehalogenase family of aspartate-dependent phosphohydrolases. Biochemistry 40:12704–12711Selengut JD, Levine RL (2000) MDP-1: a novel eukaryotic magnesium-dependent phosphatase. Biochemistry 39:8315–8324Shin DH, Roberts A, Jancarik J, Yocota H, Kim R, Wemmer DE, Kim S-H (2003) Crystal structure of a phosphatase with a unique substrate binding domain from Thermotoga maritime. Protein Sci 12:1464–1472Sussman I, Avron M (1981) Characterization and partial puri-fication of dl-glycerol-1-phosphatase from Dunaliella salina. Biochim Biophys Acta 661:199–204The Arabidopsis Genome Initiative (2000) Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature 408:796–815Tremblay LW, Dunaway-Mariano D, Allen KN (2006) Structure and activity analyses of Escherichia coli K-12 NagD provide insight into the evolution of biochemical function in the haloalkanoic acid dehalogenase superfamily. Biochemistry 45:1183–1193Vicient CM, Delseny M (1999) Isolation of total RNA from Arabidopsis thaliana seeds. Anal Biochem 268:412–413Wang W, Cho HS, Kim R, Jancarik J, Yokota H, Nguyen HH, Grigoriev IV, Wemmer DE, Kim S-H (2002) Structural characterization of the reaction pathway in phosphoserine phosphatase: crystallographic “snapshots” of intermediate states. J Mol Biol 319:421–431Zhang G, Mazurkie AS, Dunaway-Mariano D, Allen KN (2002) Kinetic evidence for a substrate-induced fit in phosphonoacetaldehyde hydrolase catalysis. Biochemistry 41:13370–13377Zhang G, Morais MC, Dai J, Zhang W, Dunaway-Mariano D, Allen KN (2004) Investigation of metal Ion binding in phosphonoacetaldehyde hydrolase identifies sequence markers for metal-activated enzymes of the HAD enzyme superfamily. Biochemistry 43:4990–4997Zimmermann P, Hirsch-Hoffmann M, Hennig L, Gruissem W (2004) GENEVESTIGATOR: Arabidopsis microarray database and analysis toolbox. Plant Physiol 136:2621–263