Multivariate phenotype analysis enables genome-wide inference of mammalian gene function.

The function of the majority of genes in the human and mouse genomes is unknown. Investigating and illuminating this dark genome is a major challenge for the biomedical sciences. The International Mouse Phenotyping Consortium (IMPC) is addressing this through the generation and broad-based phenotyping of a knockout (KO) mouse line for every protein-coding gene, producing a multidimensional data set that underlies a genome-wide annotation map from genes to phenotypes. Here, we develop a multivariate (MV) statistical approach and apply it to IMPC data comprising 148 phenotypes measured across 4,548 KO lines. There are 4,256 (1.4% of 302,997 observed data measurements) hits called by the univariate (UV) model analysing each phenotype separately, compared to 31,843 (10.5%) hits in the observed data results of the MV model, corresponding to an estimated 7.5-fold increase in power of the MV model relative to the UV model. One key property of the data set is its 55.0% rate of missingness, resulting from quality control filters and incomplete measurement of some KO lines. This raises the question of whether it is possible to infer perturbations at phenotype-gene pairs at which data are not available, i.e., to infer some in vivo effects using statistical analysis rather than experimentation. We demonstrate that, even at missing phenotypes, the MV model can detect perturbations with power comparable to the single-phenotype analysis, thereby filling in the complete gene-phenotype map with good sensitivity. A factor analysis of the MV model's fitted covariance structure identifies 20 clusters of phenotypes, with each cluster tending to be perturbed collectively. These factors cumulatively explain 75% of the KO-induced variation in the data and facilitate biological interpretation of perturbations. We also demonstrate that the MV approach strengthens the correspondence between IMPC phenotypes and existing gene annotation databases. Analysis of a subset of KO lines measured in replicate across multiple laboratories confirms that the MV model increases power with high replicability.

A bioimage informatics platform for high-throughput embryo phenotyping.

High-throughput phenotyping is a cornerstone of numerous functional genomics projects. In recent years, imaging screens have become increasingly important in understanding gene-phenotype relationships in studies of cells, tissues and whole organisms. Three-dimensional (3D) imaging has risen to prominence in the field of developmental biology for its ability to capture whole embryo morphology and gene expression, as exemplified by the International Mouse Phenotyping Consortium (IMPC). Large volumes of image data are being acquired by multiple institutions around the world that encompass a range of modalities, proprietary software and metadata. To facilitate robust downstream analysis, images and metadata must be standardized to account for these differences. As an open scientific enterprise, making the data readily accessible is essential so that members of biomedical and clinical research communities can study the images for themselves without the need for highly specialized software or technical expertise. In this article, we present a platform of software tools that facilitate the upload, analysis and dissemination of 3D images for the IMPC. Over 750 reconstructions from 80 embryonic lethal and subviable lines have been captured to date, all of which are openly accessible at mousephenotype.org. Although designed for the IMPC, all software is available under an open-source licence for others to use and develop further. Ongoing developments aim to increase throughput and improve the analysis and dissemination of image data. Furthermore, we aim to ensure that images are searchable so that users can locate relevant images associated with genes, phenotypes or human diseases of interest.

Redox regulation of pyruvate kinase M2 by cysteine oxidation and S-nitrosation.

We show here that the M2 isoform of human pyruvate kinase (M2PYK) is susceptible to nitrosation and oxidation, and that these modifications regulate enzyme activity by preventing the formation of the active tetrameric form. The biotin-switch assay carried out on M1 and M2 isoforms showed that M2PYK is sensitive to nitrosation and that Cys326 is highly susceptible to redox modification. Structural and enzymatic studies have been carried out on point mutants for three cysteine residues (Cys424, Cys358, and Cys326) to characterise their potential roles in redox regulation. Nine cysteines are conserved between M2PYK and M1PYK. Cys424 is the only cysteine unique to M2PYK. C424S, C424A, and C424L showed a moderate effect on enzyme activity with 80, 100, and 140% activity, respectively, compared with M2PYK. C358 had been previously identified from in vivo studies to be the favoured target for oxidation. Our characterised mutant showed that this mutation stabilises tetrameric M2PYK, suggesting that the in vivo resistance to oxidation for the Cys358Ser mutation is due to stabilisation of the tetrameric form of the enzyme. In contrast, the Cys326Ser mutant exists predominantly in monomeric form. A biotin-switch assay using this mutant also showed a significant reduction in biotinylation of M2PYK, confirming that this is a major target for nitrosation and probably oxidation. Our results show that the sensitivity of M2PYK to oxidation and nitrosation is regulated by its monomer-tetramer equilibrium. In the monomer state, residues (in particular C326) are exposed to oxidative modifications that prevent reformation of the active tetrameric form.

Applying the ARRIVE Guidelines to an In Vivo Database.

The Animal Research: Reporting of In Vivo Experiments (ARRIVE) guidelines were developed to address the lack of reproducibility in biomedical animal studies and improve the communication of research findings. While intended to guide the preparation of peer-reviewed manuscripts, the principles of transparent reporting are also fundamental for in vivo databases. Here, we describe the benefits and challenges of applying the guidelines for the International Mouse Phenotyping Consortium (IMPC), whose goal is to produce and phenotype 20,000 knockout mouse strains in a reproducible manner across ten research centres. In addition to ensuring the transparency and reproducibility of the IMPC, the solutions to the challenges of applying the ARRIVE guidelines in the context of IMPC will provide a resource to help guide similar initiatives in the future.

The International Mouse Phenotyping Consortium (IMPC): a functional catalogue of the mammalian genome that informs conservation.

The International Mouse Phenotyping Consortium (IMPC) is building a catalogue of mammalian gene function by producing and phenotyping a knockout mouse line for every protein-coding gene. To date, the IMPC has generated and characterised 5186 mutant lines. One-third of the lines have been found to be non-viable and over 300 new mouse models of human disease have been identified thus far. While current bioinformatics efforts are focused on translating results to better understand human disease processes, IMPC data also aids understanding genetic function and processes in other species. Here we show, using gorilla genomic data, how genes essential to development in mice can be used to help assess the potentially deleterious impact of gene variants in other species. This type of analyses could be used to select optimal breeders in endangered species to maintain or increase fitness and avoid variants associated to impaired-health phenotypes or loss-of-function mutations in genes of critical importance. We also show, using selected examples from various mammal species, how IMPC data can aid in the identification of candidate genes for studying a condition of interest, deliver information about the mechanisms involved, or support predictions for the function of genes that may play a role in adaptation. With genotyping costs decreasing and the continued improvements of bioinformatics tools, the analyses we demonstrate can be routinely applied.

The International Mouse Phenotyping Consortium Web Portal, a unified point of access for knockout mice and related phenotyping data.

The International Mouse Phenotyping Consortium (IMPC) web portal (http://www.mousephenotype.org) provides the biomedical community with a unified point of access to mutant mice and rich collection of related emerging and existing mouse phenotype data. IMPC mouse clinics worldwide follow rigorous highly structured and standardized protocols for the experimentation, collection and dissemination of data. Dedicated 'data wranglers' work with each phenotyping center to collate data and perform quality control of data. An automated statistical analysis pipeline has been developed to identify knockout strains with a significant change in the phenotype parameters. Annotation with biomedical ontologies allows biologists and clinicians to easily find mouse strains with phenotypic traits relevant to their research. Data integration with other resources will provide insights into mammalian gene function and human disease. As phenotype data become available for every gene in the mouse, the IMPC web portal will become an invaluable tool for researchers studying the genetic contributions of genes to human diseases.

M2 pyruvate kinase provides a mechanism for nutrient sensing and regulation of cell proliferation.

We show that the M2 isoform of pyruvate kinase (M2PYK) exists in equilibrium between monomers and tetramers regulated by allosteric binding of naturally occurring small-molecule metabolites. Phenylalanine stabilizes an inactive T-state tetrameric conformer and inhibits M2PYK with an IC50 value of 0.24 mM, whereas thyroid hormone (triiodo-L-thyronine, T3) stabilizes an inactive monomeric form of M2PYK with an IC50 of 78 nM. The allosteric activator fructose-1,6-bisphosphate [F16BP, AC50 (concentration that gives 50% activation) of 7 µM] shifts the equilibrium to the tetrameric active R-state, which has a similar activity to that of the constitutively fully active isoform M1PYK. Proliferation assays using HCT-116 cells showed that addition of inhibitors phenylalanine and T3 both increased cell proliferation, whereas addition of the activator F16BP reduced proliferation. F16BP abrogates the inhibitory effect of both phenylalanine and T3, highlighting a dominant role of M2PYK allosteric activation in the regulation of cancer proliferation. X-ray structures show constitutively fully active M1PYK and F16BP-bound M2PYK in an R-state conformation with a lysine at the dimer-interface acting as a peg in a hole, locking the active tetramer conformation. Binding of phenylalanine in an allosteric pocket induces a 13° rotation of the protomers, destroying the peg-in-hole R-state interface. This distinct T-state tetramer is stabilized by flipped out Trp/Arg side chains that stack across the dimer interface. X-ray structures and biophysical binding data of M2PYK complexes explain how, at a molecular level, fluctuations in concentrations of amino acids, thyroid hormone, and glucose metabolites switch M2PYK on and off to provide the cell with a nutrient sensing and growth signaling mechanism.

A mouse informatics platform for phenotypic and translational discovery.

The International Mouse Phenotyping Consortium (IMPC) is providing the world's first functional catalogue of a mammalian genome by characterising a knockout mouse strain for every gene. A robust and highly structured informatics platform has been developed to systematically collate, analyse and disseminate the data produced by the IMPC. As the first phase of the project, in which 5000 new knockout strains are being broadly phenotyped, nears completion, the informatics platform is extending and adapting to support the increasing volume and complexity of the data produced as well as addressing a large volume of users and emerging user groups. An intuitive interface helps researchers explore IMPC data by giving overviews and the ability to find and visualise data that support a phenotype assertion. Dedicated disease pages allow researchers to find new mouse models of human diseases, and novel viewers provide high-resolution images of embryonic and adult dysmorphologies. With each monthly release, the informatics platform will continue to evolve to support the increased data volume and to maintain its position as the primary route of access to IMPC data and as an invaluable resource for clinical and non-clinical researchers.

Thermorudis pharmacophila sp. nov., a novel member of the class Thermomicrobia isolated from geothermal soil, and emended descriptions of Thermomicrobium roseum, Thermomicrobium carboxidum, Thermorudis peleae and Sphaerobacter thermophilus.

An aerobic, thermophilic and cellulolytic bacterium, designated strain WKT50.2T, was isolated from geothermal soil at Waikite, New Zealand. Strain WKT50.2T grew at 53-76 °C and at pH 5.9-8.2. The DNA G+C content was 58.4 mol%. The major fatty acids were 12-methyl C18 : 0 and C18 : 0. Polar lipids were all linked to long-chain 1,2-diols, and comprised 2-acylalkyldiol-1-O-phosphoinositol (diolPI), 2-acylalkyldiol-1-O-phosphoacylmannoside (diolP-acylMan), 2-acylalkyldiol-1-O-phosphoinositol acylmannoside (diolPI-acylMan) and 2-acylalkyldiol-1-O-phosphoinositol mannoside (diolPI-Man). Strain WKT50.2T utilized a range of cellulosic substrates, alcohols and organic acids for growth, but was unable to utilize monosaccharides. Robust growth of WKT50.2T was observed on protein derivatives. WKT50.2T was sensitive to ampicillin, chloramphenicol, kanamycin, neomycin, polymyxin B, streptomycin and vancomycin. Metronidazole, lasalocid A and trimethoprim stimulated growth. Phylogenetic analysis of 16S rRNA gene sequences showed that WKT50.2T belonged to the class Thermomicrobia within the phylum Chloroflexi, and was most closely related to Thermorudis peleae KI4T (99.6% similarity). DNA-DNA hybridization between WKT50.2T and Thermorudis peleae DSM 27169T was 18.0%. Physiological and biochemical tests confirmed the phenotypic and genotypic differentiation of strain WKT50.2T from Thermorudis peleae KI4T and other members of the Thermomicrobia. On the basis of its phylogenetic position and phenotypic characteristics, we propose that strain WKT50.2T represents a novel species, for which the name Thermorudis pharmacophila sp. nov. is proposed, with the type strain WKT50.2T ( = DSM 26011T = ICMP 20042T). Emended descriptions of Thermomicrobium roseum, Thermomicrobium carboxidum, Thermorudis peleae and Sphaerobacter thermophilus are also proposed, and include the description of a novel respiratory quinone, MK-8 2,3-epoxide (23%), in Thermomicrobium roseum.

Pyruvate kinases have an intrinsic and conserved decarboxylase activity.

The phosphotransfer mechanism of PYKs (pyruvate kinases) has been studied in detail, but the mechanism of the intrinsic decarboxylase reaction catalysed by PYKs is still unknown. 1H NMR was used in the present study to follow OAA (oxaloacetate) decarboxylation by trypanosomatid and human PYKs confirming that the decarboxylase activity is conserved across distantly related species. Crystal structures of TbPYK (Trypanosoma brucei PYK) complexed with the product of the decarboxylase reaction (pyruvate), and a series of substrate analogues (D-malate, 2-oxoglutarate and oxalate) show that the OAA analogues bind to the kinase active site with similar binding modes, confirming that both decarboxylase and kinase activities share a common site for substrate binding and catalysis. Decarboxylation of OAA as monitored by NMR for TbPYK has a relatively low turnover with values of 0.86 s-1 and 1.47 s-1 in the absence and presence of F26BP (fructose 2,6-bisphosphate) respectively. Human M1PYK (M1 isoform of PYK) has a measured turnover value of 0.50 s-1. The X-ray structures explain why the decarboxylation activity is specific for OAA and is not general for α-oxo acid analogues. Conservation of the decarboxylase reaction across divergent species is a consequence of piggybacking on the conserved kinase mechanism which requires a stabilized enol intermediate.

Trypanosomatid phosphoglycerate mutases have multiple conformational and oligomeric states.

Three structurally distinct forms of phosphoglycerate mutase from the trypanosomatid parasite Leishmania mexicana were isolated by standard procedures of bacterial expression and purification. Analytical size-exclusion chromatography coupled to a multi-angle scattering detector detected two monomeric forms of differing hydrodynamic radii, as well as a dimeric form. Structural comparisons of holoenzyme and apoenzyme trypanosomatid cofactor-independent phosphoglycerate mutase (iPGAM) X-ray crystal structures show a large conformational change between the open (apoenzyme) and closed (holoenzyme) forms accounting for the different monomer hydrodynamic radii. Until now iPGAM from trypanosomatids was considered to be only monomeric, but results presented here show the appearance of a dimeric form. Taken together, these observations are important for the choice of screening strategies to identify inhibitors of iPGAM for parasite chemotherapy and highlight the need to select the most biologically or functionally relevant form of the purified enzyme.

Effect of long-term starvation on the survival, recovery, and carbon utilization profiles of a bovine Escherichia coli O157:H7 isolate from New Zealand.

The ability to maintain a dual lifestyle of colonizing the ruminant gut and surviving in nonhost environments once shed is key to the success of Escherichia coli O157:H7 as a zoonotic pathogen. Both physical and biological conditions encountered by the bacteria are likely to change during the transition between host and nonhost environments. In this study, carbon starvation at suboptimal temperatures in nonhost environments was simulated by starving a New Zealand bovine E. coli O157:H7 isolate in phosphate-buffered saline at 4 and 15°C for 84 days. Recovery of starved cells on media with different nutrient availabilities was monitored under aerobic and anaerobic conditions. We found that the New Zealand bovine E. coli O157:H7 isolate was able to maintain membrane integrity and viability over 84 days and that the level of recovery depended on the nutrient level of the recovery medium as well as the starvation temperature. In addition, a significant difference in carbon utilization was observed between starved and nonstarved cells.

How hospital leaders can bend the anesthesia cost curve.

Factors to consider in evaluating a hospital's anesthesia cost curve include: The efficiency of services provided by anesthesiologists and certified registered nurse anesthetists. Cost-effectiveness of anesthesia staffing. Anesthesia billing performance. Patient satisfaction with anesthesia services.

Investigating the lignocellulolytic gut microbiome of huhu grubs (Prionoplus reticularis) using defined diets and dietary switch.

The huhu beetle (Prionoplus reticularis) is the largest endemic beetle found throughout Aotearoa New Zealand, and is characterised by feeding on wood during its larval stage. It has been hypothesised that its gut microbiome plays a fundamental role in the degradation of wood. To explore this idea we examined the fungal and bacterial community composition of huhu grubs' frass, using amplicon sequencing. Grubs were reared on an exclusive diet of either a predominantly cellulose source (cotton) or lignocellulose source (pine) for 4 months; subsequently a diet switch was performed and the grubs were grown for another 4 months. The fungal community of cellulose-reared huhu grubs was abundant in potential cellulose degraders, contrasting with the community of lignocellulose-reared grubs, which showed abundant potential soft rot fungi, yeasts, and hemicellulose and cellulose degraders. Cellulose-reared grubs showed a less diverse fungal community, however, diet switch from cellulose to lignocellulose resulted in a change in community composition that showed grubs were still capable of utilising this substrate. Conversely, diet seemed to have a limited influence on huhu grub gut bacterial communities.

`In crystallo' substrate binding triggers major domain movements and reveals magnesium as a co-activator of Trypanosoma brucei pyruvate kinase.

The active site of pyruvate kinase (PYK) is located between the AC core of the enzyme and a mobile lid corresponding to domain B. Many PYK structures have already been determined, but the first `effector-only' structure and the first with PEP (the true natural substrate) are now reported for the enzyme from Trypanosoma brucei. PEP soaked into crystals of the enzyme with bound allosteric activator fructose 2,6-bisphosphate (F26BP) and Mg(2+) triggers a substantial 23° rotation of the B domain `in crystallo', resulting in a partially closed active site. The interplay of side chains with Mg(2+) and PEP may explain the mechanism of the domain movement. Furthermore, it is apparent that when F26BP is present but PEP is absent Mg(2+) occupies a position that is distinct from the two canonical Mg(2+)-binding sites at the active site. This third site is adjacent to the active site and involves the same amino-acid side chains as in canonical site 1 but in altered orientations. Site 3 acts to sequester Mg(2+) in a `priming' position such that the enzyme is maintained in its R-state conformation. In this way, Mg(2+) cooperates with F26BP to ensure that the enzyme is in a conformation that has a high affinity for the substrate.

A new family of covalent inhibitors block nucleotide binding to the active site of pyruvate kinase.

PYK (pyruvate kinase) plays a central role in the metabolism of many organisms and cell types, but the elucidation of the details of its function in a systems biology context has been hampered by the lack of specific high-affinity small-molecule inhibitors. High-throughput screening has been used to identify a family of saccharin derivatives which inhibit LmPYK (Leishmania mexicana PYK) activity in a time- (and dose-) dependent manner, a characteristic of irreversible inhibition. The crystal structure of DBS {4-[(1,1-dioxo-1,2-benzothiazol-3-yl)sulfanyl]benzoic acid} complexed with LmPYK shows that the saccharin moiety reacts with an active-site lysine residue (Lys335), forming a covalent bond and sterically hindering the binding of ADP/ATP. Mutation of the lysine residue to an arginine residue eliminated the effect of the inhibitor molecule, providing confirmation of the proposed inhibitor mechanism. This lysine residue is conserved in the active sites of the four human PYK isoenzymes, which were also found to be irreversibly inhibited by DBS. X-ray structures of PYK isoforms show structural differences at the DBS-binding pocket, and this covalent inhibitor of PYK provides a chemical scaffold for the design of new families of potentially isoform-specific irreversible inhibitors.

EDULISS: a small-molecule database with data-mining and pharmacophore searching capabilities.

We present the relational database EDULISS (EDinburgh University Ligand Selection System), which stores structural, physicochemical and pharmacophoric properties of small molecules. The database comprises a collection of over 4 million commercially available compounds from 28 different suppliers. A user-friendly web-based interface for EDULISS (available at http://eduliss.bch.ed.ac.uk/) has been established providing a number of data-mining possibilities. For each compound a single 3D conformer is stored along with over 1600 calculated descriptor values (molecular properties). A very efficient method for unique compound recognition, especially for a large scale database, is demonstrated by making use of small subgroups of the descriptors. Many of the shape and distance descriptors are held as pre-calculated bit strings permitting fast and efficient similarity and pharmacophore searches which can be used to identify families of related compounds for biological testing. Two ligand searching applications are given to demonstrate how EDULISS can be used to extract families of molecules with selected structural and biophysical features.