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

Investigating the effects of a cryptic splice site in the En2 splice acceptor sequence used in the IKMC knockout-first alleles

Targeted mouse mutants are a common tool used to investigate gene function. The International Knockout Mouse Consortium undertook a large-scale screen of mouse mutants, making use of the knockout-first allele design that contains the En2 splice acceptor sequence coupled to the lacZ reporter gene. Although the knockout-first allele was designed to interfere with splicing and thus disrupt gene function, the En2 sequence has been reported to be transcribed within the host gene mRNA due to a cryptic splice site within the En2 sequence which allows splicing to the next exon of the host gene. In some circumstances, this has the potential to permit translation of a mutant protein. Here, we describe our computational analysis of all the mouse protein-coding genes with established knockout-first embryonic stem cell lines, and our predictions of their transcription outcome should the En2 sequence be included. As part of the large-scale mutagenesis program, mutant mice underwent a broad phenotyping screen, and their phenotypes are available. No wide-scale effects on mouse phenotypes reported were found as a result of the predicted En2 insertion. However, the En2 insertion was found experimentally in the transcripts of 24 of 35 mutant alleles examined, including the five already described, two with evidence of readthrough. Splicing from the cryptic splice site also has the potential to disrupt expression of the lacZ reporter gene. It is recommended that mutant transcripts be checked for this insertion as well as for leaky transcription in studies involving knockout-first alleles.

Computational identification of disease models through cross-species phenotype comparison

The use of standardised phenotyping screens to identify abnormal phenotypes in mouse knockouts, together with the use of ontologies to describe such phenotypic features, allows the implementation of an automated and unbiased pipeline to identify new models of disease by performing phenotype comparisons across species. Using data from the International Mouse Phenotyping Consortium (IMPC), approximately half of mouse mutants are able to mimic, at least partially, the human ortholog disease phenotypes as computed by the PhenoDigm algorithm. We found the number of phenotypic abnormalities in the mouse and the corresponding Mendelian disorder, the pleiotropy and severity of the disease, and the viability and zygosity status of the mouse knockout to be associated with the ability of mouse models to recapitulate the human disorder. An analysis of the IMPC impact on disease gene discovery through a publication-tracking system revealed that the resource has been implicated in at least 109 validated rare disease-gene associations over the last decade.

A consensus score to combine inferences from multiple centres

Experiments in which data are collected by multiple independent resources, including multicentre data, different laboratories within the same centre or with different operators, are challenging in design, data collection and interpretation. Indeed, inconsistent results across the resources are possible. In this paper, we propose a statistical solution for the problem of multi-resource consensus inferences when statistical results from different resources show variation in magnitude, directionality, and significance. Our proposed method allows combining the corrected p-values, effect sizes and the total number of centres into a global consensus score. We apply this method to obtain a consensus score for data collected by the International Mouse Phenotyping Consortium (IMPC) across 11 centres. We show the application of this method to detect sexual dimorphism in haematological data and discuss the suitability of the methodology.

KOMPUTE: imputing summary statistics of missing phenotypes in high-throughput model organism data

Motivation

The International Mouse Phenotyping Consortium (IMPC) is striving to build a comprehensive functional catalog of mammalian protein-coding genes by systematically producing and phenotyping gene-knockout mice for almost every protein-coding gene in the mouse genome and by testing associations between gene loss-of-function and phenotype. To date, the IMPC has identified over 90 000 gene-phenotype associations, but many phenotypes have not yet been measured for each gene, resulting in largely incomplete data; ∼75.6% of association summary statistics are still missing in the latest IMPC summary statistics dataset (IMPC release version 16).

Results

To overcome these challenges, we propose KOMPUTE, a novel method for imputing missing summary statistics in the IMPC dataset. Using conditional distribution properties of multivariate normal, KOMPUTE estimates the association Z-scores of unmeasured phenotypes for a particular gene as a conditional expectation given the Z-scores of measured phenotypes. Our evaluation of the method using simulated and real-world datasets demonstrates its superiority over the singular value decomposition matrix completion method in various scenarios.

Availability and implementation

An R package for KOMPUTE is publicly available at https://github.com/statsleelab/kompute, along with usage examples and results for different phenotype domains at https://statsleelab.github.io/komputeExamples.

Knockout mice are an important tool for human monogenic heart disease studies

Mouse models are relevant to studying the functionality of genes involved in human diseases; however, translation of phenotypes can be challenging. Here, we investigated genes related to monogenic forms of cardiovascular disease based on the Genomics England PanelApp and aligned them to International Mouse Phenotyping Consortium (IMPC) data. We found 153 genes associated with cardiomyopathy, cardiac arrhythmias or congenital heart disease in humans, of which 151 have one-to-one mouse orthologues. For 37.7% (57/151), viability and heart data captured by electrocardiography, transthoracic echocardiography, morphology and pathology from embryos and young adult mice are available. In knockout mice, 75.4% (43/57) of these genes showed non-viable phenotypes, whereas records of prenatal, neonatal or infant death in humans were found for 35.1% (20/57). Multisystem phenotypes are common, with 58.8% (20/34) of heterozygous (homozygous lethal) and 78.6% (11/14) of homozygous (viable) mice showing cardiovascular, metabolic/homeostasis, musculoskeletal, hematopoietic, nervous system and/or growth abnormalities mimicking the clinical manifestations observed in patients. These IMPC data are critical beyond cardiac diagnostics given their multisystemic nature, allowing detection of abnormalities across physiological systems and providing a valuable resource to understand pleiotropic effects.

Genome-wide screening reveals the genetic basis of mammalian embryonic eye development

BACKGROUND

Microphthalmia, anophthalmia, and coloboma (MAC) spectrum disease encompasses a group of eye malformations which play a role in childhood visual impairment. Although the predominant cause of eye malformations is known to be heritable in nature, with 80% of cases displaying loss-of-function mutations in the ocular developmental genes OTX2 or SOX2, the genetic abnormalities underlying the remaining cases of MAC are incompletely understood. This study intended to identify the novel genes and pathways required for early eye development. Additionally, pathways involved in eye formation during embryogenesis are also incompletely understood. This study aims to identify the novel genes and pathways required for early eye development through systematic forward screening of the mammalian genome.

RESULTS

Query of the International Mouse Phenotyping Consortium (IMPC) database (data release 17.0, August 01, 2022) identified 74 unique knockout lines (genes) with genetically associated eye defects in mouse embryos. The vast majority of eye abnormalities were small or absent eyes, findings most relevant to MAC spectrum disease in humans. A literature search showed that 27 of the 74 lines had previously published knockout mouse models, of which only 15 had ocular defects identified in the original publications. These 12 previously published gene knockouts with no reported ocular abnormalities and the 47 unpublished knockouts with ocular abnormalities identified by the IMPC represent 59 genes not previously associated with early eye development in mice. Of these 59, we identified 19 genes with a reported human eye phenotype. Overall, mining of the IMPC data yielded 40 previously unimplicated genes linked to mammalian eye development. Bioinformatic analysis showed that several of the IMPC genes colocalized to several protein anabolic and pluripotency pathways in early eye development. Of note, our analysis suggests that the serine-glycine pathway producing glycine, a mitochondrial one-carbon donator to folate one-carbon metabolism (FOCM), is essential for eye formation.

CONCLUSIONS

Using genome-wide phenotype screening of single-gene knockout mouse lines, STRING analysis, and bioinformatic methods, this study identified genes heretofore unassociated with MAC phenotypes providing models to research novel molecular and cellular mechanisms involved in eye development. These findings have the potential to hasten the diagnosis and treatment of this congenital blinding disease.

Identifying genetic determinants of inflammatory pain in mice using a large-scale gene-targeted screen

ABSTRACT

Identifying the genetic determinants of pain is a scientific imperative given the magnitude of the global health burden that pain causes. Here, we report a genetic screen for nociception, performed under the auspices of the International Mouse Phenotyping Consortium. A biased set of 110 single-gene knockout mouse strains was screened for 1 or more nociception and hypersensitivity assays, including chemical nociception (formalin) and mechanical and thermal nociception (von Frey filaments and Hargreaves tests, respectively), with or without an inflammatory agent (complete Freund's adjuvant). We identified 13 single-gene knockout strains with altered nocifensive behavior in 1 or more assays. All these novel mouse models are openly available to the scientific community to study gene function. Two of the 13 genes (Gria1 and Htr3a) have been previously reported with nociception-related phenotypes in genetically engineered mouse strains and represent useful benchmarking standards. One of the 13 genes (Cnrip1) is known from human studies to play a role in pain modulation and the knockout mouse reported herein can be used to explore this function further. The remaining 10 genes (Abhd13, Alg6, BC048562, Cgnl1, Cp, Mmp16, Oxa1l, Tecpr2, Trim14, and Trim2) reveal novel pathways involved in nociception and may provide new knowledge to better understand genetic mechanisms of inflammatory pain and to serve as models for therapeutic target validation and drug development.

Extensive identification of genes involved in congenital and structural heart disorders and cardiomyopathy

Clinical presentation of congenital heart disease is heterogeneous, making identification of the disease-causing genes and their genetic pathways and mechanisms of action challenging. By using in vivo electrocardiography, transthoracic echocardiography and microcomputed tomography imaging to screen 3,894 single-gene-null mouse lines for structural and functional cardiac abnormalities, here we identify 705 lines with cardiac arrhythmia, myocardial hypertrophy and/or ventricular dilation. Among these 705 genes, 486 have not been previously associated with cardiac dysfunction in humans, and some of them represent variants of unknown relevance (VUR). Mice with mutations in Casz1, Dnajc18, Pde4dip, Rnf38 or Tmem161b genes show developmental cardiac structural abnormalities, with their human orthologs being categorized as VUR. Using UK Biobank data, we validate the importance of the DNAJC18 gene for cardiac homeostasis by showing that its loss of function is associated with altered left ventricular systolic function. Our results identify hundreds of previously unappreciated genes with potential function in congenital heart disease and suggest causal function of five VUR in congenital heart disease.

Transcriptome-wide association study and eQTL colocalization identify potentially causal genes responsible for human bone mineral density GWAS associations

Genome-wide association studies (GWASs) for bone mineral density (BMD) in humans have identified over 1100 associations to date. However, identifying causal genes implicated by such studies has been challenging. Recent advances in the development of transcriptome reference datasets and computational approaches such as transcriptome-wide association studies (TWASs) and expression quantitative trait loci (eQTL) colocalization have proven to be informative in identifying putatively causal genes underlying GWAS associations. Here, we used TWAS/eQTL colocalization in conjunction with transcriptomic data from the Genotype-Tissue Expression (GTEx) project to identify potentially causal genes for the largest BMD GWAS performed to date. Using this approach, we identified 512 genes as significant using both TWAS and eQTL colocalization. This set of genes was enriched for regulators of BMD and members of bone relevant biological processes. To investigate the significance of our findings, we selected PPP6R3, the gene with the strongest support from our analysis which was not previously implicated in the regulation of BMD, for further investigation. We observed that Ppp6r3 deletion in mice decreased BMD. In this work, we provide an updated resource of putatively causal BMD genes and demonstrate that PPP6R3 is a putatively causal BMD GWAS gene. These data increase our understanding of the genetics of BMD and provide further evidence for the utility of combined TWAS/colocalization approaches in untangling the genetics of complex traits.