International Mouse Phenotyping Consortium Research Articles

High-Throughput Screening of Mouse Gene Knockouts Identifies Established and Novel High Body Fat Phenotypes.

PurposeObesity is a major public health problem. Understanding which genes contribute to obesity may better predict individual risk and allow development of new therapies. Because obesity of a mouse gene knockout (KO) line predicts an association of the orthologous human gene with obesity, we reviewed data from the Lexicon Genome5000TM high throughput phenotypic screen (HTS) of mouse gene KOs to identify KO lines with high body fat.Materials and MethodsKO lines were generated using homologous recombination or gene trapping technologies. HTS body composition analyses were performed on adult wild-type and homozygous KO littermate mice from 3758 druggable mouse genes having a human ortholog. Body composition was measured by either DXA or QMR on chow-fed cohorts from all 3758 KO lines and was measured by QMR on independent high fat diet-fed cohorts from 2488 of these KO lines. Where possible, comparisons were made to HTS data from the International Mouse Phenotyping Consortium (IMPC).ResultsBody fat data are presented for 75 KO lines. Of 46 KO lines where independent external published and/or IMPC KO lines are reported as obese, 43 had increased body fat. For the remaining 29 novel high body fat KO lines, Ksr2 and G2e3 are supported by data from additional independent KO cohorts, 6 (Asnsd1, Srpk2, Dpp8, Cxxc4, Tenm3 and Kiss1) are supported by data from additional internal cohorts, and the remaining 21 including Tle4, Ak5, Ntm, Tusc3, Ankk1, Mfap3l, Prok2 and Prokr2 were studied with HTS cohorts only.ConclusionThese data support the finding of high body fat in 43 independent external published and/or IMPC KO lines. A novel obese phenotype was identified in 29 additional KO lines, with 27 still lacking the external confirmation now provided for Ksr2 and G2e3 KO mice. Undoubtedly, many mammalian obesity genes remain to be identified and characterized.

An Axin2\xa0mutation and perinatal risk factors contribute to sagittal craniosynostosis: evidence from a Chinese female monochorionic diamniotic twin family

BackgroundCraniosynostosis, defined as premature fusion of one or more cranial sutures, affects approximately 1 in every 2000–2500 live births. Sagittal craniosynostosis (CS), the most prevalent form of isolated craniosynostosis, is caused by interplay between genetic and perinatal environmental insults. However, the underlying details remain largely unknown.MethodsThe proband (a female monochorionic twin diagnosed with CS), her healthy co-twin sister and parents were enrolled. Obstetric history was extracted from medical records. Genetic screening was performed by whole exome sequencing (WES) and confirmed by Sanger sequencing. Functional annotation, conservation and structural analysis were predicted in public database. Phenotype data of Axin2 knockout mice was downloaded from The International Mouse Phenotyping Consortium (IMPC, http://www.mousephenotype.org).ResultsObstetric medical records showed that, except for the shared perinatal risk factors by the twins, the proband suffered additional persistent breech presentation and intrauterine growth restriction. We identified a heterozygous mutation of Axin2 (c.1181G > A, p.R394H, rs200899695) in monochorionic twins and their father, but not in the mother. This mutation is not reported in Asian population and results in replacement of Arg at residue 394 by His (p.R394H). Arg 394 is located at the GSK3β binding domain of Axin2 protein, which is highly conserved across species. The mutation was predicted to be potentially deleterious by in silico analysis. Incomplete penetrance of Axin2 haploinsufficiency was found in female mice.ConclusionsAxin2 (c.1181G > A, p.R394H, rs200899695) mutation confers susceptibility and perinatal risk factors trigger the occurrence of sagittal craniosynostosis. Our findings provide a new evidence for the gene-environment interplay in understanding pathogenesis of craniosynostosis in Chinese population.

Proteotyping of knockout mouse strains reveals sex- and strain-specific signatures in blood plasma

We proteotyped blood plasma from 30 mouse knockout strains and corresponding wild-type mice from the International Mouse Phenotyping Consortium. We used targeted proteomics with internal standards to quantify 375 proteins in 218 samples. Our results provide insights into the manifested effects of each gene knockout at the plasma proteome level. We first investigated possible contamination by erythrocytes during sample preparation and labeled, in one case, up to 11 differential proteins as erythrocyte originated. Second, we showed that differences in baseline protein abundance between female and male mice were evident in all mice, emphasizing the necessity to include both sexes in basic research, target discovery, and preclinical effect and safety studies. Next, we identified the protein signature of each gene knockout and performed functional analyses for all knockout strains. Further, to demonstrate how proteome analysis identifies the effect of gene deficiency beyond traditional phenotyping tests, we provide in-depth analysis of two strains, C8a−/− and Npc2+/−. The proteins encoded by these genes are well-characterized providing good validation of our method in homozygous and heterozygous knockout mice. Ig alpha chain C region, a poorly characterized protein, was among the differentiating proteins in C8a−/−. In Npc2+/− mice, where histopathology and traditional tests failed to differentiate heterozygous from wild-type mice, our data showed significant difference in various lysosomal storage disease-related proteins. Our results demonstrate how to combine absolute quantitative proteomics with mouse gene knockout strategies to systematically study the effect of protein absence. The approach used here for blood plasma is applicable to all tissue protein extracts.

A resource of targeted mutant mouse lines for 5,061 genes.

The International Mouse Phenotyping Consortium reports the generation of new mouse mutant strains for more than 5,000 genes, including 2,850 novel null, 2,987 novel conditional-ready and 4,433 novel reporter alleles.

LAMA: automated image analysis for the developmental phenotyping of mouse embryos.

ABSTRACTAdvanced 3D imaging modalities, such as micro-computed tomography (micro-CT), have been incorporated into the high-throughput embryo pipeline of the International Mouse Phenotyping Consortium (IMPC). This project generates large volumes of raw data that cannot be immediately exploited without significant resources of personnel and expertise. Thus, rapid automated annotation is crucial to ensure that 3D imaging data can be integrated with other multi-dimensional phenotyping data. We present an automated computational mouse embryo phenotyping pipeline that harnesses the large amount of wild-type control data available in the IMPC embryo pipeline in order to address issues of low mutant sample number as well as incomplete penetrance and variable expressivity. We also investigate the effect of developmental substage on automated phenotyping results. Designed primarily for developmental biologists, our software performs image pre-processing, registration, statistical analysis and segmentation of embryo images. We also present a novel anatomical E14.5 embryo atlas average and, using it with LAMA, show that we can uncover known and novel dysmorphology from two IMPC knockout lines.

Online control of the familywise error rate.

Biological research often involves testing a growing number of null hypotheses as new data are accumulated over time. We study the problem of online control of the familywise error rate, that is testing an a priori unbounded sequence of hypotheses (p-values) one by one over time without knowing the future, such that with high probability there are no false discoveries in the entire sequence. This paper unifies algorithmic concepts developed for offline (single batch) familywise error rate control and online false discovery rate control to develop novel online familywise error rate control methods. Though many offline familywise error rate methods (e.g., Bonferroni, fallback procedures and Sidak's method) can trivially be extended to the online setting, our main contribution is the design of new, powerful, adaptive online algorithms that control the familywise error rate when the p-values are independent or locally dependent in time. Our numerical experiments demonstrate substantial gains in power, that are also formally proved in an idealized Gaussian sequence model. A promising application to the International Mouse Phenotyping Consortium is described.

OpenStats: A robust and scalable software package for reproducible analysis of high-throughput phenotypic data.

Reproducibility in the statistical analyses of data from high-throughput phenotyping screens requires a robust and reliable analysis foundation that allows modelling of different possible statistical scenarios. Regular challenges are scalability and extensibility of the analysis software. In this manuscript, we describe OpenStats, a freely available software package that addresses these challenges. We show the performance of the software in a high-throughput phenomic pipeline in the International Mouse Phenotyping Consortium (IMPC) and compare the agreement of the results with the most similar implementation in the literature. OpenStats has significant improvements in speed and scalability compared to existing software packages including a 13-fold improvement in computational time to the current production analysis pipeline in the IMPC. Reduced complexity also promotes FAIR data analysis by providing transparency and benefiting other groups in reproducing and re-usability of the statistical methods and results. OpenStats is freely available under a Creative Commons license at www.bioconductor.org/packages/OpenStats.

High-throughput genotyping of high-homology mutant mouse strains by next-generation sequencing.

Genotyping of knockout alleles in mice is commonly performed by end-point PCR or gene-specific/universal cassette qPCR. Both have advantages and limitations in terms of assay design and interpretation of results. As an alternative method for high-throughput genotyping, we investigated next generation sequencing (NGS) of PCR amplicons, with a focus on CRISPR-mediated exon deletions where antibiotic selection markers are not present. By multiplexing the wild type and mutant-specific PCR reactions, the genotype can be called by the relative sequence counts of each product. The system is highly scalable and can be applied to a variety of different allele types, including those produced by the International Mouse Phenotyping Consortium and associated projects. One potential challenge with any assay design is locating unique areas of the genome, especially when working with gene families or regions of high homology. These can result in misleading or ambiguous genotypes for either qPCR or end-point assays. Here, we show that genotyping by NGS can negate these issues by simple, automated filtering of undesired sequences. Analysis and genotype calls can also be fully automated, using FASTQ or FASTA input files and an in-house Perl script and SQL database.

Role of UDP-Sugar Receptor P2Y14 in Murine Osteoblasts.

The purinergic (P2) receptor P2Y14 is the only P2 receptor that is stimulated by uridine diphosphate (UDP)-sugars and its role in bone formation is unknown. We confirmed P2Y14 expression in primary murine osteoblasts (CB-Ob) and the C2C12-BMP2 osteoblastic cell line (C2-Ob). UDP-glucose (UDPG) had undiscernible effects on cAMP levels, however, induced dose-dependent elevations in the cytosolic free calcium concentration ([Ca2+]i) in CB-Ob, but not C2-Ob cells. To antagonize the P2Y14 function, we used the P2Y14 inhibitor PPTN or generated CRISPR-Cas9-mediated P2Y14 knockout C2-Ob clones (Y14KO). P2Y14 inhibition facilitated calcium signalling and altered basal cAMP levels in both models of osteoblasts. Importantly, P2Y14 inhibition augmented Ca2+ signalling in response to ATP, ADP and mechanical stimulation. P2Y14 knockout or inhibition reduced osteoblast proliferation and decreased ERK1/2 phosphorylation and increased AMPKα phosphorylation. During in vitro osteogenic differentiation, P2Y14 inhibition modulated the timing of osteogenic gene expression, collagen deposition, and mineralization, but did not significantly affect differentiation status by day 28. Of interest, while P2ry14-/- mice from the International Mouse Phenotyping Consortium were similar to wild-type controls in bone mineral density, their tibia length was significantly increased. We conclude that P2Y14 in osteoblasts reduces cell responsiveness to mechanical stimulation and mechanotransductive signalling and modulates osteoblast differentiation.

The effects of microchipping C57BL/6N mice on standard phenotyping tests.

The C57BL/6N inbred lines of mice are widely used in genetic research. They are particularly favoured in large scale studies such as the International Mouse Phenotyping Consortium (IMPC), where C57BL/6N mice are genetically altered to generate a collection of null alleles (currently more than 8500 null alleles have been generated). In this project, mice carrying null alleles are subjected to a pipeline of broad-based phenotyping tests to produce wide ranging phenotyping data on each model. We have previously described the development of a Home Cage Analysis system that automatically tracks the activity of group housed mice from a microchip inserted in the groin. This platform allows assessment of multiple biologically relevant phenotypes over long periods of time without experimenter interference, and therefore is particularly suited for high through-put studies. To investigate the impact of microchips on other tests carried out in the IMPC pipeline, we inserted microchips in 12 male and 12 female C57BL/6Ntac mice at seven weeks of age. Starting at nine weeks of age these mice underwent standard phenotyping tests, concurrently with 20 unchipped C57BL/6Ntac mice (10 females, 10 males). Tissues from a subset of the microchipped mice (six males and six females), chosen at random, were also sent for histopathological examination at the end of the phenotyping pipeline. No significant impact of insertion of microchip was observed in any of the phenotyping tests apart from bone mineral density measurement at DEXA due to the nature of the microchip. We therefore recommend that the microchip be inserted during the DEXA procedure, after the measurement is taken but before the mouse has recovered from the anaesthetic. This would avoid multiple anaesthetic exposures and prevent the potential variability in DEXA analysis output.

High Fidelity of Mouse Models Mimicking Human Genetic Skeletal Disorders.

The 2019 International Skeletal Dysplasia Society nosology update lists 441 genes for which mutations result in rare human skeletal disorders. These genes code for enzymes (33%), scaffolding proteins (18%), signal transduction proteins (16%), transcription factors (14%), cilia proteins (8%), extracellular matrix proteins (5%), and membrane transporters (4%). Skeletal disorders include aggrecanopathies, channelopathies, ciliopathies, cohesinopathies, laminopathies, linkeropathies, lysosomal storage diseases, protein-folding and RNA splicing defects, and ribosomopathies. With the goal of evaluating the ability of mouse models to mimic these human genetic skeletal disorders, a PubMed literature search identified 260 genes for which mutant mice were examined for skeletal phenotypes. These mouse models included spontaneous and ENU-induced mutants, global and conditional gene knockouts, and transgenic mice with gene over-expression or specific base-pair substitutions. The human X-linked gene ARSE and small nuclear RNA U4ATAC, a component of the minor spliceosome, do not have mouse homologs. Mouse skeletal phenotypes mimicking human skeletal disorders were observed in 249 of the 260 genes (96%) for which comparisons are possible. A supplemental table in spreadsheet format provides PubMed weblinks to representative publications of mutant mouse skeletal phenotypes. Mutations in 11 mouse genes (Ccn6, Cyp2r1, Flna, Galns, Gna13, Lemd3, Manba, Mnx1, Nsd1, Plod1, Smarcal1) do not result in similar skeletal phenotypes observed with mutations of the homologous human genes. These discrepancies can result from failure of mouse models to mimic the exact human gene mutations. There are no obvious commonalities among these 11 genes. Body BMD and/or radiologic dysmorphology phenotypes were successfully identified for 28 genes by the International Mouse Phenotyping Consortium (IMPC). Forward genetics using ENU mouse mutagenesis successfully identified 37 nosology gene phenotypes. Since many human genetic disorders involve hypomorphic, gain-of-function, dominant-negative and intronic mutations, future studies will undoubtedly utilize CRISPR/Cas9 technology to examine transgenic mice having genes modified to exactly mimic variant human sequences. Mutant mice will increasingly be employed for drug development studies designed to treat human genetic skeletal disorders.SignificanceGreat progress is being made identifying mutant genes responsible for human rare genetic skeletal disorders and mouse models for genes affecting bone mass, architecture, mineralization and strength. This review organizes data for 441 human genetic bone disorders with regard to heredity, gene function, molecular pathways, and fidelity of relevant mouse models to mimic the human skeletal disorders. PubMed weblinks to citations of 249 successful mouse models are provided.

Human and mouse essentiality screens as a resource for disease gene discovery

The identification of causal variants in sequencing studies remains a considerable challenge that can be partially addressed by new gene-specific knowledge. Here, we integrate measures of how essential a gene is to supporting life, as inferred from viability and phenotyping screens performed on knockout mice by the International Mouse Phenotyping Consortium and essentiality screens carried out on human cell lines. We propose a cross-species gene classification across the Full Spectrum of Intolerance to Loss-of-function (FUSIL) and demonstrate that genes in five mutually exclusive FUSIL categories have differing biological properties. Most notably, Mendelian disease genes, particularly those associated with developmental disorders, are highly overrepresented among genes non-essential for cell survival but required for organism development. After screening developmental disorder cases from three independent disease sequencing consortia, we identify potentially pathogenic variants in genes not previously associated with rare diseases. We therefore propose FUSIL as an efficient approach for disease gene discovery.

Incremental data integration for tracking genotype-disease associations.

Functional annotation of genes remains a challenge in fundamental biology and is a limiting factor for translational medicine. Computational approaches have been developed to process heterogeneous data into meaningful metrics, but often do not address how findings might be updated when new evidence comes to light. To address this challenge, we describe requirements for a framework for incremental data integration and propose an implementation based on phenotype ontologies and Bayesian probability updates. We apply the framework to quantify similarities between gene annotations and disease profiles. Within this scope, we categorize human diseases according to how well they can be recapitulated by animal models and quantify similarities between human diseases and mouse models produced by the International Mouse Phenotyping Consortium. The flexibility of the approach allows us to incorporate negative phenotypic data to better prioritize candidate genes, and to stratify disease mapping using sex-dependent phenotypes. All our association scores can be updated and we exploit this feature to showcase integration with curated annotations from high-precision assays. Incremental integration is thus a suitable framework for tracking functional annotations and linking to complex human pathology.

Fam151b, the mouse homologue of C.elegans menorin gene, is essential for retinal function

Fam151b is a mammalian homologue of the C. elegans menorin gene, which is involved in neuronal branching. The International Mouse Phenotyping Consortium (IMPC) aims to knock out every gene in the mouse and comprehensively phenotype the mutant animals. This project identified Fam151b homozygous knock-out mice as having retinal degeneration. We show they have no photoreceptor function from eye opening, as demonstrated by a lack of electroretinograph (ERG) response. Histological analysis shows that during development of the eye the correct number of cells are produced and that the layers of the retina differentiate normally. However, after eye opening at P14, Fam151b mutant eyes exhibit signs of retinal stress and rapidly lose photoreceptor cells. We have mutated the second mammalian menorin homologue, Fam151a, and homozygous mutant mice have no discernible phenotype. Sequence analysis indicates that the FAM151 proteins are members of the PLC-like phosphodiesterase superfamily. However, the substrates and function of the proteins remains unknown.

The effects of microchipping C57BL/6N mice on standard phenotyping tests

The C57BL/6N inbred lines of mice are widely used in genetic research. They are particularly favoured in large scale studies such as the International Mouse Phenotyping Consortium (IMPC), where C57BL/6N mice are genetically altered to generate a collection of null alleles (currently more than 8500 null alleles have been generated). In this project, mice carrying null alleles are subjected to a pipeline of broad-based phenotyping tests to produce wide ranging phenotyping data on each model. We have previously described the development of a Home Cage Analysis system that automatically tracks the activity of group housed mice from a microchip inserted in the groin. This platform allows assessment of multiple biologically relevant phenotypes over long periods of time without experimenter interference, and therefore is particularly suited for high through-put studies. To investigate the impact of microchips on other tests carried out in the IMPC pipeline, we inserted microchips in 12 male and 12 female C57BL/6Ntac mice at seven weeks of age. Starting at nine weeks of age these mice underwent standard phenotyping tests, concurrently with 20 unchipped C57BL/6Ntac mice (10 females, 10 males). Tissues from a subset of the microchipped mice (six males and six females), chosen at random, were also sent for histopathological examination at the end of the phenotyping pipeline. No significant impact of insertion of microchip was observed in any of the phenotyping tests apart from bone mineral density measurement at DEXA due to the nature of the microchip. We therefore recommend that the microchip be inserted during the DEXA procedure, after the measurement is taken but before the mouse has recovered from the anaesthetic. This would avoid multiple anaesthetic exposures and prevent the potential variability in DEXA analysis output.

High-throughput discovery of genetic determinants of circadian misalignment.

Circadian systems provide a fitness advantage to organisms by allowing them to adapt to daily changes of environmental cues, such as light/dark cycles. The molecular mechanism underlying the circadian clock has been well characterized. However, how internal circadian clocks are entrained with regular daily light/dark cycles remains unclear. By collecting and analyzing indirect calorimetry (IC) data from more than 2000 wild-type mice available from the International Mouse Phenotyping Consortium (IMPC), we show that the onset time and peak phase of activity and food intake rhythms are reliable parameters for screening defects of circadian misalignment. We developed a machine learning algorithm to quantify these two parameters in our misalignment screen (SyncScreener) with existing datasets and used it to screen 750 mutant mouse lines from five IMPC phenotyping centres. Mutants of five genes (Slc7a11, Rhbdl1, Spop, Ctc1 and Oxtr) were found to be associated with altered patterns of activity or food intake. By further studying the Slc7a11tm1a/tm1a mice, we confirmed its advanced activity phase phenotype in response to a simulated jetlag and skeleton photoperiod stimuli. Disruption of Slc7a11 affected the intercellular communication in the suprachiasmatic nucleus, suggesting a defect in synchronization of clock neurons. Our study has established a systematic phenotype analysis approach that can be used to uncover the mechanism of circadian entrainment in mice.

The importance of sex differences in pharmacology research

This article is part of a themed section on The Importance of Sex Differences in Pharmacology Research. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v176.21/issuetoc.

Soft windowing application to improve analysis of high-throughput phenotyping data.

MotivationHigh-throughput phenomic projects generate complex data from small treatment and large control groups that increase the power of the analyses but introduce variation over time. A method is needed to utlize a set of temporally local controls that maximizes analytic power while minimizing noise from unspecified environmental factors.ResultsHere we introduce ‘soft windowing’, a methodological approach that selects a window of time that includes the most appropriate controls for analysis. Using phenotype data from the International Mouse Phenotyping Consortium (IMPC), adaptive windows were applied such that control data collected proximally to mutants were assigned the maximal weight, while data collected earlier or later had less weight. We applied this method to IMPC data and compared the results with those obtained from a standard non-windowed approach. Validation was performed using a resampling approach in which we demonstrate a 10% reduction of false positives from 2.5 million analyses. We applied the method to our production analysis pipeline that establishes genotype–phenotype associations by comparing mutant versus control data. We report an increase of 30% in significant P-values, as well as linkage to 106 versus 99 disease models via phenotype overlap with the soft-windowed and non-windowed approaches, respectively, from a set of 2082 mutant mouse lines. Our method is generalizable and can benefit large-scale human phenomic projects such as the UK Biobank and the All of Us resources.Availability and implementationThe method is freely available in the R package SmoothWin, available on CRAN http://CRAN.R-project.org/package=SmoothWin.Supplementary information Supplementary data are available at Bioinformatics online.

Genome-wide screening of mouse knockouts reveals novel genes required for normal integumentary and oculocutaneous structure and function

Oculocutaneous syndromes are often due to mutations in single genes. In some cases, mouse models for these diseases exist in spontaneously occurring mutations, or in mice resulting from forward mutatagenesis screens. Here we present novel genes that may be causative for oculocutaneous disease in humans, discovered as part of a genome-wide screen of knockout-mice in a targeted single-gene deletion project. The International Mouse Phenotyping Consortium (IMPC) database (data release 10.0) was interrogated for all mouse strains with integument abnormalities, which were then cross-referenced individually to identify knockouts with concomitant ocular abnormalities attributed to the same targeted gene deletion. The search yielded 307 knockout strains from unique genes with integument abnormalities, 226 of which have not been previously associated with oculocutaneous conditions. Of the 307 knockout strains with integument abnormalities, 52 were determined to have ocular changes attributed to the targeted deletion, 35 of which represent novel oculocutaneous genes. Some examples of various integument abnormalities are shown, as well as two examples of knockout strains with oculocutaneous phenotypes. Each of the novel genes provided here are potentially relevant to the pathophysiology of human integumentary, or oculocutaneous conditions, such as albinism, phakomatoses, or other multi-system syndromes. The novel genes reported here may implicate molecular pathways relevant to these human diseases and may contribute to the discovery of novel therapeutic targets.

A Comprehensive Plasma Metabolomics Dataset for a Cohort of Mouse Knockouts within the International Mouse Phenotyping Consortium.

Mouse knockouts facilitate the study ofgene functions. Often, multiple abnormal phenotypes are induced when a gene is inactivated. The International Mouse Phenotyping Consortium (IMPC) has generated thousands of mouse knockouts and catalogued their phenotype data. We have acquired metabolomics data from 220 plasma samples from 30 unique mouse gene knockouts and corresponding wildtype mice from the IMPC. To acquire comprehensive metabolomics data, we have used liquid chromatography (LC) combined with mass spectrometry (MS) for detecting polar and lipophilic compounds in an untargeted approach. We have also used targeted methods to measure bile acids, steroids and oxylipins. In addition, we have used gas chromatography GC-TOFMS for measuring primary metabolites. The metabolomics dataset reports 832 unique structurally identified metabolites from 124 chemical classes as determined by ChemRICH software. The GCMS and LCMS raw data files, intermediate and finalized data matrices, R-Scripts, annotation databases, and extracted ion chromatograms are provided in this data descriptor. The dataset can be used for subsequent studies to link genetic variants with molecular mechanisms and phenotypes.