Abstract
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.
Similar content being viewed by others
Introduction
Mus musculus is the most used animal model in scientific research. It has high similarity with humans at the molecular level with 99% of human genes having homologs in the mouse genome1. Mice can model many human diseases, making them suitable to study rare monogenic disorders and complex multigenic diseases such as cancer, diabetes, and even anxiety2,3,4,5,6,7,8,9. Current genome manipulation techniques to knock out or silence a specific gene have allowed many human conditions to be reproduced in mice, enabling the study of disease mechanism and progression10,11,12,13,14.
These studies were largely performed using high throughput analytical methods. Analysis of mouse tissues in the context of health and disease has been done previously using microarray and deep sequencing technologies15,16. Although genes are the original template for proteins, it is the expressed proteins and their differential abundance that principally determine the function of cells and tissues. Hence, parallel to the various sequencing efforts, comprehensive studies at the proteome level have been performed in recent years and provided insight into the proteins that are differentially expressed between cells, tissues, organs, or organ systems, or are related to a specific condition or disease. These studies have revealed functional genomics insights beyond that derived from sequencing alone17,18,19,20. Such efforts have included the analysis of cells and tissues from wild type, transgenic, knockin, and knockout strains, and mice labeled in vivo with isotope using mass spectrometry-based methods21,22,23,24. Mass spectrometry is a versatile technique that allows system-wide study of the proteome25. In a typical bottom-up workflow, proteins are digested into peptides for analysis using liquid chromatography coupled to mass spectrometry (LC-MS/MS)26. Differential expression is inferred from mass spectrum signal intensity and good comparability across groups can be achieved using labeling approaches such as isobaric tagging using Tandem Mass Tag (TMT), Isobaric Tag for Relative and Absolute Quantitation (iTRAQ), or Stable Isotope Labeling with Amino acids in Cell culture (SILAC) 27. Multiple reaction monitoring (MRM) is considered the gold standard in quantitative measurements28,29,30. When combined with heavy labeled internal standards, high precision and accuracy were achieved while multiplexing assays for hundreds of proteins within a single experiment.
We set out to conduct a systematic comparison using large-scale, targeted proteomic analysis of the impacts caused by single-gene disruption. Two hundred and eighteen plasma samples from 90 female and 90 male mice for 30 knockout (KO) strains and 38 corresponding wild-type controls were analyzed. All KO strains and controls were on the C57BL/6N genetic background. The mutant mice were produced and phenotyped through a standardized pipeline of sequential tests by the International Mouse Phenotyping Consortium (IMPC). The KO gene targets were selected on the basis of their known involvement in diverse biological processes, with the goal of evaluating how plasma proteomics can complement clinical in vivo and terminal phenotyping tests (Table 1). We first chose homozygous (HOM) and heterozygous (HET) strains to study the effect of protein ablation (HOM) and reduced protein abundance (HET). Approximately 30% of the KO strains produced by the IMPC are embryo lethal or subviable11, so it was important to test if the proteomic analysis was sensitive enough to detect changes in heterozygous mice. We also included female and male mice to study sexual dimorphism at the plasma proteome level and determine possible interaction with gene KO related protein abundances. Further, we purposely included KO strains with various protein expression profiles including secreted, widely expressed, ubiquitous, and tissue-specific, as well as proteins with no known tissue specificity.
Selection of proteins measured was based on their involvement in various biological pathways and detectability31. The abundances of 375 plasma proteins were measured using MRM assays validated according to the CPTAC guidelines32 (Supplementary Table 1 and Supplementary Fig. 1). The measured plasma protein concentrations provided a molecular phenotype for each KO strain in addition to the clinical in vivo and terminal test phenotype data from the IMPC. To our knowledge, this is the first large-scale analysis of plasma proteins in KO mice.
Results
We proteotyped 30 mouse KO strains and corresponding wild-type controls using quantitative targeted proteomics. We realize that the number of strains analyzed is small compared to other phenotyping test and interpretation studies; e.g. Karp et al.33 analyzed 2186 strains for sexual dimorphism in 238 standard IMPC phenotyping tests. Proteotyping on that scale, i.e. >54,000 samples would require a large coordinated effort in addition to the associated high operational costs. Recognizing this limitation, strain selection was particularly important to address the questions of proteotyping capabilities to detect protein abundance differences between HOM and HET genotypes, female and male mice, and protein expression profiles. Our results identified differences for all three criteria suggesting that proteotyping by current state-of-the-art quantitative methods is possible, biologically relevant, and scalable.
Of the 375 measured proteins, 284 were detectable, and 234 were quantifiable with a minimum of 5% of all measurements above the lower limit of quantification (LLOQ). Two hundred and twenty-six proteins were quantified within the dynamic range of the assays in all three mice of at least one mouse strain and sex (Fig. 1, Supplementary Fig. 2 and Supplementary Table 1); therefore, we used the minimal set of 226 proteins in our subsequent analyses. The determined concentrations of these proteins spanned five orders of magnitude, ranging from 0.27 to 6.2 × 104 fmol/μL plasma, demonstrating the large dynamic range that is quantifiable using LC-MRM/MS (Fig. 1a). Overall, these measurements had very good precision34,35,36 with an average coefficient of variation (CV) of 9.3%, and all were below 23% (Fig. 1b).
Proteins originating from erythrocytes and platelets
Due to daily differences in sample collection and processing, plasma samples routinely contain variable amounts of proteins originating from red blood cells and platelets. Recently, Geyer et al.33 identified the contaminating proteins from erythrocytes and platelets in human plasma, which can be used as indicators of differences in sample processing and should be excluded from inference analysis between biological groups, unless independence of sample handling can be established. We measured 12 erythrocyte- and 10 platelet-specific intracellular proteins that were previously identified by Geyer et al. as common contaminants. Correlation analysis in all samples between all proteins (Supplementary Fig. 3 and Fig. 2 filtered for minimum absolute Pearson coefficient of 0.8) showed the clustering of these proteins in correlated groups, indicating their amounts measured in some of our samples are in fact artifacts of sample processing. In addition to the erythrocyte proteins identified by Geyer et al.37, a strong correlation was also observed with Ubiquitin-like protein ISG15. This intracellular protein is involved in erythroid differentiation38; therefore, we concluded ISG15 also originated from erythrocytes during sample collection. In our further analyses, these 22 reported erythrocyte contaminants plus ISG15 were closely examined. Specifically, if any of these 23 proteins were significantly altered in a comparison between groups, we determined if the other erythrocyte proteins showed a similar trend. This allowed us to determine if the differential expression was an artifact of sample collection or a signature of the gene KO. As our sample collection method produces platelet-rich plasma, we opted to consider platelet proteins as part of our samples. Our results (Table 1) emphasize the importance of carefully considering whether the presence of intracellular proteins in plasma reflects a biological condition, or if they are the result of sample processing.
Sexual dimorphism
Various studies have previously demonstrated sex-specific differences in protein expression39,40,41. Similar to findings in brain, urine, and apheresis platelet supernatant, we observed a clear protein signature associated with sex for mouse plasma, in both KO and control mice. Using principle component analysis (PCA), wild type and KO mice clustered clearly according to sex (with the exception of one outlier from SRA1) in PC1–PC2 plane, which explained 31.4% of the variance in our dataset (Fig. 3a). Proteins with a high significant difference in expression (p value < 0.01 and fold change of 2) were identified (Fig. 3b), and can be used to create a sex-specific protein signature (Fig. 3c). The top discriminating proteins have been previously reported to have differential expression in males and females, including Adiponectin (ADIPOQ), Alpha-1-antitrypsin (SERPINA1E), Alpha-1B-glycoprotein (A1BG), Alpha-2-macroglobulin-P (A2M), and the complement components (C5, C8A, C8B, C8G)37,42,43,44,45,46,47. The concentrations of these proteins in the plasma of male and female mice are shown in Fig. 3d. Similarly, Transcortin (SERPINA6)48, Epidermal growth factor receptor (EGFR)49, Glycosylation-dependent cell adhesion molecule-1 (GLYCAM1)50, Haptoglobin (HP)51, Murinoglobulin-1 (MUG1)52,53, Serum amyloid A proteins (SAA2, SAA4, SAA1)54, and Thyroxine-binding globulin (SERPINA7)55 were reported previously as sexually dimorphic.
We were able to measure good discrimination over a wide dynamic range spanning from a few fmol/μL in glycosylation-dependent cell adhesion molecule (1 GLYCAM1), up to thousands in Alpha-1-antitrypsin 1–5 (SERPINA1E) and corticosteroid-binding globulin (SERPINA6) as shown in Fig. 3d.
When we compared the profile obtained in our work to a recent study on sexual dimorphism in human plasma proteins56, in which 142 proteins were identified to be differential between females and males, only three proteins were shared: adiponectin, a1-antitrypsin, and thyroxine-binding globulin.
Previously it was shown that 56.6% of the phenotypic continuous (non-categorical) measurements performed by the IMPC are associated with sex57. Our results extend these findings to show sexual dimorphism at the molecular level in plasma.
Correlation with standard phenotyping tests
The mice used here were characterized as part of the IMPC program58 using standardized tests to measure biological parameters from the hematological, metabolic, cardiovascular, musculoskeletal, and neurological systems59. Since our study focused on blood plasma, we compared our proteomic data with available clinical chemistry, hematology, and body composition measurements. We obtained several good correlations between the proteomic and traditional phenotyping measurements despite their separation in time, space, and technology, i.e. correlated values were obtained from measurements on frozen samples performed years apart at different locations using different technologies. Figure 4 shows selected correlations found with Spearman correlation of around 0.8. Strong correlations identified were between high-density lipoprotein (HDL) and cholesterol with apolipoproteins A1 and A2, which was expected given the role of these proteins as major structural components of the high-density lipoprotein complex. Aspartate aminotransferase (AST) is an enzyme involved in amino acid metabolism and its level in blood is often used as an indicator of liver function and damage. In our data, we identified a correlation between AST and beta-enolase, both are enzymes essential for glycolysis/gluconeogenesis60. H-2 class I histocompatibility antigen Q10 has been noted to associate with lipids in C57BL/6 mouse plasma in other studies61, and was found to correlate with measured HDL in our data as well. For these correlations, sexual dimorphism was a clear confounding factor as can be seen in Fig. 4. However, despite decreases in strength, correlations persist after adjusting for sex effect and performing the regression on the residuals. Similarly, regression on stratified data showed similar trends.
The goal of phenotyping KO mice is to identify the consequences of gene dysfunction, which in turn can provide insight into gene function, gene pleiotropy (comorbidities), and generate hypotheses for mechanisms of disease. For the strains examined in this study, none of the standard tests considered here (clinical chemistry, hematology, or body composition) discriminated knockouts from their corresponding controls using the IMPC’s standard statistical analyses62. In such cases, molecular level investigation may identify differences, as discussed below, aiding the detailed characterization of a KO mouse strain.
Proteomic phenotyping of gene deficiency in knockout mice
Although the largest variation in plasma protein abundance was linked to sexual dimorphism (Fig. 3), we were able to determine proteomic signatures specific to 28 gene knockouts (Fig. 5). Here we used simple PCA on the proteins selected by Least Absolute Shrinkage and Selection Operator—LASSO63 (Supplementary Table 2) to demonstrate the possible grouping of samples in the PC1 and PC2 plane. For two of the KO strains, G6pd2 and Sra1, no discriminating proteins were found. In this analysis we removed all erythrocyte-specific proteins for simplicity. The discrimination observed highlights how targeted proteomics with simple data analysis can be used for molecular phenotyping. We have also previously shown possible discrimination between co-housed and co-raised littermate wild type and KO mice (thus much less effect of possible environmental variables) using our targeted proteomics assays and data analysis64.
For each KO strain, we next identified proteins significantly affected by the absence of the gene using Mann–Whitney–Wilcoxon test and calculated the fold change of proteins based on the mean values. We continued our analysis with the proteins differentially expressed (twofold difference in abundance between groups with p value < 0.05) which are listed in Table 1. The number of these proteins ranged from zero as seen in Idh1−/− and A2m−/− mice, to a strong effect with up to 10 and more differentiating proteins as seen in Npc2+/− and Iqgap1−/− mice. Our analysis confirmed the expected absence of protein, when measured, in the corresponding gene KO mouse, as in the case of C8A in the C8a−/− strain. Similar confirmation was also demonstrated in a parallel work, in which we confirmed expected absence of proteins in gene knockdown experiments by targeted proteomics64.
We further performed multiple overrepresentation analyses (ORAs) using the differentially expressed proteins obtained by Mann–Whitney–Wilcoxon test in combination with the discriminating proteins selected by LASSO63. ORA allows identification of known functions, processes, and diseases that are associated with a set of genes or proteins of interest65. We performed systematic ORA using multiple knowledgebases including Gene Ontology Terms—GO66, Molecular Signatures—MsigDB67, molecular pathway using Kyoto Encyclopedia of Genes and Genomes—KEGG68 as well as Reactome60, Disease Ontology—DO69, diseases and their gene associations using DisGeNET70, and Medical Subject Headings—MeSH for processes and diseases71. While some of these resources are overlapping in context, they differ in content and curation method, hence reporting different views. For disease-related analyses, the human orthologs were used. When both mouse proteins and human orthologs were available in a resource, as for Reactome and MeSH processes, we performed parallel analyses. In total, we performed 10 ORA for each mouse KO mouse strain. The results are included in Supplementary ORA-report 1 for discriminating proteins from Mann–Whitney–Wilcoxon test, and Supplementary ORA-report 2 for using the combined protein list of the statistical test and LASSO regression.
C8a −/− and Npc2 +/− strains
Here we report in-depth analysis of two knockouts, C8a−/− and Npc2+/−. Figure 6 shows volcano plots with differentially abundant proteins for these two KO strains, while Fig. 7 represents part of the functional analyses performed. The complete results from all KO strains are included in Table 1 and in the Supplementary Materials.
C8a −/− mice
The C8 alpha (C8A) protein combines with C8 beta (C8B) and C8 gamma (C8G) to form the complement component 8 (C8) protein complex, which plays a key role in the immune response by participating in the assembly of the membrane attack complex (MAC)72. In response to infection, MAC forms a pore in the pathogen cell membranes, resulting in cell lysis and death. Full KO of the C8a gene was confirmed by MRM analysis, since C8A was measured in control mice but not detected in the C8a−/− (C8atm1b(EUCOMM)Hmgu) mice (Fig. 6a). The concentration of C8B was decreased in all C8a−/− males, and no C8B was detected in any of the C8a−/− females. Concentration of C8G was also decreased in all C8a KO mice. C8A, C8B, and C8G are encoded by separate genes, indicating that in the absence of the C8A, the C8 complex does not form and C8B and C8G are cleared from the circulation.
The difference between wild type and C8a−/− KO mice on the plasma protein level can clearly be seen in Fig. 6a and in Fig. 5. Two hundred and forty-six phenotyping tests for C8a−/− mice performed originally by the IMPC reported no significance between wild type and KO using the IMPC’s standard statistical analysis (Supplementary Fig. 4 and Supplementary Table 3). The IMPC’s automated statistical analysis uses significance at a threshold of 0.0001 for unadjusted p value obtained by regression analysis. When applying the criteria we used to evaluate differences in protein abundances (Benjamini–Hochberg adjusted p value threshold of 0.05) to the IMCP phenotyping tests, we obtained a single significance corresponding to the difference in neutrophil differential count between wild type and KO. Furthermore, we used a non-parametric test to compare protein abundances, which is more suitable for low sample number but usually also more stringent (a parametric t-test in case of C8a−/−—which corresponds to a regression analysis with varying intercept—produced 13 additional differentiated proteins besides the 10 included in our study).
ORAs of C8a−/− included five proteins, C8A, C8B, C8G, MASP2, and Ig alpha chain C region, which were obtained by hypothesis testing (Table 1) and selected by LASSO regression (Supplementary Table 2) as discriminators (Fig. 7a–d). While ISG15 levels showed a significant change between KO and control samples, it had a high correlation with erythrocyte-originating proteins and thus we concluded it was a sample preparation contaminant although not reported as such by Geyer et al.33 Nonetheless, we performed the ORAs with and without ISG15 to test its effect, which agreed with the rest of the protein set. Processes and functions overrepresented in this set of proteins showed various immune system related entries which reflected the role of the complement 8 complex and MASP2 in the complement system (as well as ISG15 in the innate immune system when included). Both mediated and adaptive immune responses were represented, which was expected given the terminal role of complement 8 in the innate immune response through its classical, alternative, and lectin pathways. MASP2, on the other hand, takes part only in the lectin pathway of the complement. When included, ISG15 also enriched immune system related functions through its antiviral and DNA repair roles. To that end, the proteins differentiated in their abundance in C8a−/− mice fall into two categories: those with direct interaction with C8A within the complement 8 complex and those indirectly affected through the impact of C8A absence on the innate immune system.
Disease-related ORA results could be linked to an impaired innate immune system, including various infections73 and leukemia74. Neisseriaceae infections, for example, are linked specifically to the deficiency of complement 8 which hinders the formation of MAC75,76. Ig alpha chain C region was the only protein upregulated in C8a−/− compared to the control mice. Having this protein in a set of discriminating proteins that enriches for innate immune system functions is noteworthy. Among the 20 measured immunoglobins, Ig alpha chain C region is the only statistically significant discriminator in C8a−/− KO mice. Currently, little is known about this protein77, and ORA results did not link it to any available annotation in the different knowledge bases we used. Recent studies associated Ig alpha chain C region with Duchenne muscular dystrophy (in Mdx4cv mouse model), prion effect on liver (in PrPC KO mice), as well as the glycoproteome of prion infected mice78,79,80. While our experiments were not sufficient to conclude a direct link of the Ig alpha chain C region to the complement system or immune system, our data suggest a possible link with additional studies necessary to confirm this. In conclusion, targeted proteomics analysis using C8a−/− mice was able to detect the effect of immunodeficiency resulting from an impaired complement system.
Npc2 +/− mice
NPC1 and NPC2 are endosomal/lysosomal proteins involved in the transport of cholesterol. In humans, mutations in either NPC1 or NPC2 lead to the development of Niemann–Pick disease type C (NPC disease), a lysosomal storage disorder with a broad spectrum of visceral and neurological symptoms resulting from cellular accumulation of cholesterol and glycolipids. Individual lysosomal storage disorders are rare but collectively affect 1 in 5000 births with NPC disease affecting 1 in 10,000 (ref. 81). In addition to the aggressive cerebral and visceral inflammation which are hallmarks of NPC disease, generalized immune dysfunction and hematological defects, such as thrombocytopenia and anemia, also occur82,83. NPC disease, as with most lysosomal storage diseases, is inherited in an autosomal recessive manner, affecting homozygous individuals only. In our study we included heterozygous KO and wild-type mice for the proteomics analysis, and performed histopathology on tissue sections from wild type, heterozygous, and homozygous animals (Fig. 6b). Heterozygous mice were primarily phenotyped because Npc2−/− animals were emaciated, ataxic, and needed to be euthanized at clinical endpoint at ~10 weeks of age. Consequently, all Npc2−/− histopathology was done on 10-week-old animals. Phenotyping tests for Npc2+/− mice performed by IMPC are included in Supplementary Fig. 5 and Supplementary Table 4. Analysis of the plasma proteome of Npc2+/− (Npc2tm1e.1(EUCOMM)Wtsi) mice revealed dysregulation of proteins involved in hemostasis, particularly in platelet degranulation (actin, ACTG1; P-selectin, SELP; vinculin, VCL). Several proteins associated with exosomes were also upregulated (actin, ACTG1; CD97 antigen, CD97; Elongation factor 1-alpha-1, EEF1S1; vinculin, VCL) indicating trafficking dysregulation. Histopathological examination of spleen, lymph node, bone marrow, brain, and cerebellum sections revealed no difference between wild type and heterozygous mice; in contrast, homozygous mice showed common lysosomal storage disease phenotypes (Fig. 6c). This corroborates previous histopathological observations obtained in a zebrafish model for Npc1 (ref. 84), in which liver tissue sections of wild type and heterozygous larvae were similar, but different from homozygous larvae. While the known effect of NPC2 dysregulation, i.e. NPC disease is autosomal recessive, which can be confirmed by the absence of pathological phenotype in heterozygous mice, changes at the plasma proteome level in Npc2+/− were quantifiable. In an attempt to investigate these changes, we performed multiple ORA.
Initially we included 18 proteins for ORAs, ENO1, CD97, EEF1A1, Ig heavy chain V region MOPC 47A, SERPINF1, PFN1, SELP, TNC, TALDO1, VCL, ORM2, PZP, CP, FETUB, FCN1, SPINT1, CTLA2A, and SAA1. This set was obtained by combining the results from hypothesis testing and LASSO regression (Table 1 and Supplementary Table 2). Although various associations were found, these had high p values. Reducing the analysis to only those proteins found significant by Mann–Whitney–Wilcoxon test (Table 1) improved the ORA adjusted p values; nevertheless, these were still above 0.1 (Fig. 7e–h, Supplementary ORA-report 1 and Supplementary ORA-report 2). Taking into account the investigatory nature of such analysis to drive hypothesis generation and future research, we investigated the overrepresented entries based on ordered p values. Multiple entries were related to neural development specifically and to cell growth in general. While NPC2 is mainly associated with metabolism and NPC disease, disease ORAs resulted in multiple cancer associations. A direct characterization of the relation of serum NPC2 to cancer has been reported previously85, and has linked upregulated NPC2 levels to breast, colon, and lung cancers, and downregulated levels to kidney and liver cancers in humans. Our results extend these findings and suggest that mice with NPC2 deficiency express a cancer-related protein profile in blood plasma. Further validation of these results is needed and may shed light on the less understood role of NPC2 in cancer. Including a targeted proteomics assay for NPC2 in future analyses will be beneficial to assess its level in the heterozygous KO mice. Furthermore, proteomics analysis of brain, spleen, lymph node, liver, and other tissues will advance the characterization of the Npc2+/− KO mice. The identification of heterologous pathways and disease areas affected by gene ablation (homozygous null) or dosage (heterozygous null) demonstrates the potential for proteomic analyses to increase knowledge about gene and protein function. We believe that complementary proteomic analyses may augment current methodologies to assign significance to variants86,87 or disease risk88 by assessing impacts on pathways known to be involved in disease.
Discussion
Proteomic phenotyping of KO mice using MRM mass spectrometry is a promising method for studying and understanding the function of genes beyond what can be determined through clinical in vivo and terminal test phenotyping alone89. Here we presented a broad proteotyping approach that can be incorporated as a complimentary test in large or small-scale phenotyping studies. We characterized the plasma protein profile of single-gene KO strains deficient for 30 genes. Our validated assays successfully quantified 226 proteins covering five orders of magnitude. All protein measurements had excellent precision with an average CV of 9.3%.
A strong sex-specific signature in measured plasma proteins was identified including 19 up- and downregulated proteins between female and male mice with C-statistics of 0.97, hence a sexually dimorphic blood plasma proteome signature. The differentiating proteins spanned a wide dynamic concentration range, from a few fmol/μL in Glycosylation-dependent cell adhesion molecule-1 (GLYCAM1), up to thousands in Alpha-1-antitrypsin 1–5 (SERPINA1E) and corticosteroid-binding globulin (SERPINA6). Alpha-1B-glycoprotein (A1BG) was undetectable in male animals, acting as a clear binary discriminator. We carefully investigated intracellular erythrocyte-originating proteins present in the measured plasma samples using correlation analysis and comparison to previous work33. It was possible to determine whether these erythrocyte-specific proteins are likely an artifact of sample processing, or a true effect of the deficiency of the gene KO.
The effect of gene KO observed in plasma ranged from no measured effect as seen in Idh1−/− and A2m−/− mice to a strong effect as seen in Npc2+/− and Iqgap1−/− mice where multiple proteins differentiated significantly compared to wild-type controls. We were able to detect changes in protein abundances in homozygous as well as heterozygous KO mice. We also carried out ORAs of the plasma protein profile of all knockouts covering protein functions, involvement in biological processes, and association with diseases. We highlighted insights from C8a−/− and Npc2+/− mice, where a clear plasma molecular profile was observed. Absence of C8A in C8a−/− mice was confirmed by our measurement and resulted in a plasma signature associated with (impaired) complement system. The presence of Ig alpha chain C region in C8a knockouts highlighted how proteotyping approaches help to generate hypotheses for less characterized proteins—in this case suggesting a role in the innate immune system. Functional studies using the mice described here, or other models, are needed to test this hypothesis. Mutations in human NPC2 leads to the development of Niemann–Pick disease type C, an autosomal recessive disorder90,91. Histopathological examination of various tissues from the Npc2 KO strain confirmed the presence of disease-related phenotypes in homozygous mice, but not in heterozygous mice. However, we were able to quantify changes in the plasma proteome of Npc2+/− mice. This clearly shows that proteomics is complementary to other standardized phenotyping tests. The proteomic signature detected in blood plasma of the NPC2-deficient mice was associated with cancer. Confirming previous studies that associated NPC2 levels with various cancers as measured directly by ELISA85, our results associated the blood plasma protein signature of NPC2 deficiency to cancer. We expect that measurement of additional tissues will provide a more comprehensive proteomics phenotype of the gene KO. Indeed, these types of studies may be developed to complement standard genetic screens to assess disease predisposition and risk, particularly for polygenic diseases or when assessing variants of unknown significance. We also compared our proteomic measurements to standard phenotyping tests relevant to plasma, including clinical chemistry, hematology, and body composition measurements. Several correlations were identified between plasma protein concentration and these biological parameters. While we focused our discussion on two KO strains, C8a−/− and Npc2+/−, our work includes measured abundances, determined discriminating proteins, and ORAs for all 30 KO mouse strains that we studied. Our data, in conjunction with available IMPC phenotyping results, provide an enriched resource and will help researchers interested in these proteins, or the pathways and functions their absence affects, to better formulate their hypotheses and develop experiments to test them.
Materials
Mouse plasma samples
Plasma samples for 30 KO strains (Table 1) were obtained from The Centre for Phenogenomics, which is part of the International Mouse Phenotyping Consortium (IMPC)58. Samples were collected from three male and three female mice of each KO line, as well as 19 female and 19 male C57BL/6NCrl wild-type mice collected at a similar time. All sample collection was performed in the morning before noon. Whole-blood samples were collected in tubes containing heparin from the retro-orbital sinus under isoflurane anesthetic. Samples were spun at 5000g for 10 min at 8 °C. The plasma layer was removed, aliquoted, stored at −80 °C, shipped on dry ice to the University of Victoria, and stored again at −80 °C until analyzed. All experimental procedures on animals received approval from the Animal Care Committee of The Centre for Phenogenomics and were conducted in accordance with the guidelines of the Canadian Council on Animal Care. The corresponding license numbers are AUPs 153, 275, 277, and 279. All mutant mouse lines used for plasma proteotyping are available from the Canadian Mouse Mutant Repository (CMMR) at The Centre for Phenogenomics.
Pathology
Wild type and homozygous mice were euthanized at 10 weeks of age, heterozygous mice were euthanized at 16 weeks of age, and a complete necropsy and comprehensive tissue collection for histopathology was done. Fresh tissues were immersion fixed in 10% neutral buffered formalin, paraffin-embedded, sectioned at 4–5 μm, and stained with HE. The tissues collected and processed from each mouse for histopathology included lung, thyroid, trachea, esophagus, heart, thymus, brown adipose tissue, mesenteric lymph node, adrenal gland, liver, spleen, kidney, urinary bladder, mammary gland, uterus, and ovary (from females) or testis, epididymiis, prostate, and seminal vesicle (from males), sternum, pancreas, skeletal muscle, salivary glands, stomach, duodenum, ileum, jejunum with Peyer’s patch, cecum, colon, rectum, eye, ear, spinal cord, brain, femur, tibia, knee joint, and skin (snout, pinna, dorsal, ventral, tail base)92. Histopathology evaluation was done by veterinary pathologists (H.A.A., C.M.) and images were captured using a microscope-mounted Olympus DP71 digital camera (Olympus Life Science Imaging Systems Inc., Markham, ON, Canada).
Surrogate peptide internal standards and assays
Proteotypic peptide surrogates were selected for each protein and chemically synthesized31. First, surrogates were selected by in-silico using PeptidePicker93. For synthesis of the heavy labeled peptides, 13C/15N N-Fmoc l-arginine and l-lysine (98% isotopic enrichment, Cambridge Isotope Laboratories, Andover, MA, USA) were coupled to TentaGelTM R TRT resins (RAPP Polymere, Tübingen, Germany). For synthesis of unlabeled peptides, Wang resins preloaded with non-modified N-Fmoc lysine and arginine were purchased from Matrix Innovations (Quebec City, QC, Canada). All peptides were synthesized and purified in house94. Synthesis was performed using dimethylformamide with a 10× or 20× amino acid excess, using 40% piperidine for Fmoc deprotection, and HCTU(1 eq)/NMM (2 eq) as activator/base reagents. After cleaved from the resin, the synthetic peptides were purified by reverse-phase HPLC on an Onyx silica monolithic C18 column (100 × 10 mm id, 2 μm particles; Phenomenex; Torrance, CA, USA). The peptide elution profiles were monitored by UV absorbance at 230 nm (Ultimate 3000; Dionex; Sunnyvale, CA, USA) and the fractions of interest were measured by MALDI‐TOF‐MS using an Ultraflex III TOF/TOF mass spectrometer (Bruker Daltonik; Bremen, Germany). Fractions containing more than 80% of the target peptide were pooled and lyophilized. Each synthetic peptide was characterized by capillary zone electrophoresis (CZE) to assess the purity, and by amino acid analyses (AAA) to determine its absolute concentration. The results of CZE and AAA were later used to estimate the endogenous surrogate peptide concentration by reference to the exact amount of the spiked-in synthetic heavy labeled peptide. Peptide specific instrument parameters were characterized using an Agilent 6495 Triple Quadrupole mass spectrometer. Peptide assays were validated according to the Clinical Proteome Tumor Analysis Consortium (CPTAC) guidelines for assay characterization32 to assess the response curve, repeatability, selectivity, stability, and reproducible detection of endogenous peptide31. In total, assays measuring 375 peptide surrogates covering same number of proteins were established.
Sample preparation and measurement
A brief explanation is included here with additional details provided in the supplementary materials. Mouse plasma samples were processed using the Tecan Evo (Männedorf, Switzerland) liquid handling robot and all 218 samples were randomized over three 96-well plates. A pooled reference plasma sample (BioReclamationIVT; Westbury, NY, USA) was used for quality control and normalization with 9–12 reference samples per plate inserted semi-randomly. Additional eight samples for establishing the standard curve were included on the first plate, and three curve quality control samples were included on each plate. Tryptic digestion and sample measurement were performed in a standardized way as detailed in Supplementary Materials. An 8-point external calibration curve was established for quantification using synthetic light peptides (ranging in concentration from 1 to 1000× assay LLOQ) spiked at known concentration into digested bovine serum albumin (Sigma Aldrich, Oakville, ON, CA) as a simplified background matrix95, while synthetic heavy labeled peptides were added to all samples at 100× assay LLOQ as the normalizer.
Quantification and data analysis
Endogenous analyte concentrations were calculated from the endogenous/heavy ratio using regression analysis of the standard curves (1/x2 weighting)96. Raw data were processed using Skyline97, including inspection and correction of peak integration. This step ensures that the beginning and end of the eluted peptides are included. Normalization was performed within each plate against the pooled control sample, which were measured on each plate multiple time. If the measured concentration of a specific protein was below the assay’s LLOQ for more than half of the pooled control samples within a plate, the original reported value of each sample for that specific protein was considered more trustworthy and kept unchanged. LASSO was used for identifying the minimal set of best discriminators between KO and wild-type mice that allow best discrimination63. Two-sided Mann–Whitney–Wilcoxon test was used to compare protein abundances between KO and wild-type mice and p values were adjusted with the Benjamini–Hochberg method for multiple testing. Protein fold changes were determined by calculating the ratio of mean concentrations of KO to wild-type mice. Volcano plots were used to represent p values and fold change. ORAs and the required hypergeometric test were performed using the quantified proteins as a background. Entries in seven knowledgebases were used for ORA including GO66, MsigDB67, molecular pathway using KEGG68 and Reactome60, DO69, diseases and their gene associations using DisGeNET70, and MeSH for processes and diseases71.
All data analysis and visualization were performed using R and various libraries including ggplot98 for visualization, glmnet99 for regression and statistical analysis, and ClusterProfiler100 for ORA.
Integration with standard phenotyping tests
Raw data for phenotyping tests were retrieved from the IMPC (www.mousephenotype.org) for our individual mice. Tests for integration with proteomic data included clinical chemistry (i.e. levels of sodium (mmol/L), potasium (mmol/L), chloride (mmol/L), BUN (mg/dL), creatinine (mg/dL), protein (g/L), albumin (g/L), bilirubin (mg/dL), calcium (mg/dL), phosphate (mg/dL), AST (U/L), ALT (U/L), ALP (U/L), cholesterol (mg/dL), HDL (mg/dL), triglycerol (mg/dL), glucose (mg/dL), body composition (body weight (g), fat mass (g), lean mass(g), bone mineral density—BMD (g/cm2), bone mineral content—BMC (g), body length (cm), BMC/weight (ratio), lean/weight (ratio), fat/weight (ratio), head weight (g), bone area (cm2) (BMC/BMD)), and hematology (WBC (103/μL), RBC (106/μL), Hgb (g/dL), HCT (%), MCV (fL), MCH (pg), MCHC (g/dL), PLT (103/μL), MPV (fL), RDW %, NE %, LY %, MO %, EO %, BA %, NE (103 μm/L), LY (103 μm/L), MO (103 μm/L), EO (103 μm/L), and BA (103 μm/L)). Relationships between the phenotyping tests and proteomic results were assessed using Pearson correlation coefficients.
Reporting summary
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
Data availability
All protein concentration values are available in the supplementary material file in Dataset 1.
Code availability
Data processing and analysis methods used are described in the “Quantification and data analysis” section and are all based on publically available open source software tools and packages. No custom code or mathematical algorithms were used for the data analysis.
References
Rosenthal, N. & Brown, S. The mouse ascending: perspectives for human-disease models. Nat. Cell Biol. 9, 993–999 (2007).
Cacheiro, P., Haendel, M. A. & Smedley, D., International Mouse Phenotyping C, the Monarch I. New models for human disease from the International Mouse Phenotyping Consortium. Mamm. Genome 30, 143–150 (2019).
Chen, S. P., Tolner, E. A. & Eikermann-Haerter, K. Animal models of monogenic migraine. Cephalalgia 36, 704–721 (2016).
Perlman, R. L. Mouse models of human disease: an evolutionary perspective. Evol. Med. Public Health 2016, 170–176 (2016).
Justice, M. J. & Dhillon, P. Using the mouse to model human disease: increasing validity and reproducibility. Dis. Model. Mech. 9, 101–103 (2016).
King, A. J. The use of animal models in diabetes research. Br. J. Pharmacol. 166, 877–894 (2012).
Steimer, T. Animal models of anxiety disorders in rats and mice: some conceptual issues. Dialogues Clin. Neurosci. 13, 495–506 (2011).
Cheon, D. J. & Orsulic, S. Mouse models of cancer. Annu. Rev. Pathol. 6, 95–119 (2011).
Frese, K. K. & Tuveson, D. A. Maximizing mouse cancer models. Nat. Rev. Cancer 7, 645–658 (2007).
Meehan, T. F. et al. Disease model discovery from 3,328 gene knockouts by The International Mouse Phenotyping Consortium. Nat. Genet. 49, 1231–1238 (2017).
Dickinson, M. E. et al. High-throughput discovery of novel developmental phenotypes. Nature 537, 508–514 (2016).
Rozman, J. et al. Identification of genetic elements in metabolism by high-throughput mouse phenotyping. Nat. Commun. 9, 288 (2018).
Brommage, R., Powell, D. R. & Vogel, P. Predicting human disease mutations and identifying drug targets from mouse gene knockout phenotyping campaigns. Dis. Model Mech. 12, dmm038224 (2019).
Gurumurthy, C. B. & Lloyd, K. C. K. Generating mouse models for biomedical research: technological advances. Dis. Model Mech. 12, dmm029462 (2019).
Ramskold, D., Wang, E. T., Burge, C. B. & Sandberg, R. An abundance of ubiquitously expressed genes revealed by tissue transcriptome sequence data. PLoS Comput. Biol. 5, e1000598 (2009).
Su, A. I. et al. A gene atlas of the mouse and human protein-encoding transcriptomes. Proc. Natl Acad. Sci. USA 101, 6062–6067 (2004).
Becker, K. et al. Quantifying post-transcriptional regulation in the development of Drosophila melanogaster. Nat. Commun. 9, 4970 (2018).
Jagannathan, S., Ogata, Y., Gafken, P. R., Tapscott, S. J. & Bradley, R. K. Quantitative proteomics reveals key roles for post-transcriptional gene regulation in the molecular pathology of facioscapulohumeral muscular dystrophy. Elife 8, e41740 (2019).
Wang, X., Liu, Q. & Zhang, B. Leveraging the complementary nature of RNA-Seq and shotgun proteomics data. Proteomics 14, 2676–2687 (2014).
Gianazza, E. et al. What if? Mouse proteomics after gene inactivation. J. Proteomics 199, 102–122 (2019).
Geiger, T. et al. Initial quantitative proteomic map of 28 mouse tissues using the SILAC mouse. Mol. Cell. Proteomics 12, 1709–1722 (2013).
Walther, D. M. & Mann, M. Accurate quantification of more than 4000 mouse tissue proteins reveals minimal proteome changes during aging. Mol. Cell. Proteomics 10, M110 004523 (2011).
Malmstrom, E. et al. Large-scale inference of protein tissue origin in gram-positive sepsis plasma using quantitative targeted proteomics. Nat. Commun. 7, 10261 (2016).
Davis, R. G. et al. Top-down proteomics enables comparative analysis of brain proteoforms between mouse strains. Anal. Chem. 90, 3802–3810 (2018).
Aebersold, R. & Mann, M. Mass-spectrometric exploration of proteome structure and function. Nature 537, 347–355 (2016).
Steen, H. & Mann, M. The ABC’s (and XYZ’s) of peptide sequencing. Nat. Rev. Mol. Cell. Biol. 5, 699–711 (2004).
Pappireddi, N., Martin, L. & Wuhr, M. A review on quantitative multiplexed proteomics. Chembiochem 20, 1210–1224 (2019).
Method of the Year 2012. Nat. Methods 10, 1 (2013).
Carr, S. A. et al. Targeted peptide measurements in biology and medicine: best practices for mass spectrometry-based assay development using a fit-for-purpose approach. Mol. Cell Proteomics 13, 907–917 (2014).
Addona, T. A. et al. Multi-site assessment of the precision and reproducibility of multiple reaction monitoring-based measurements of proteins in plasma. Nat. Biotechnol. 27, 633–641 (2009).
Michaud, S. A. et al. Molecular phenotyping of laboratory mouse strains using 500 multiple reaction monitoring mass spectrometry plasma assays. Commun. Biol. 1, 78 (2018).
Whiteaker, J. R. et al. CPTAC Assay Portal: a repository of targeted proteomic assays. Nat. Methods 11, 703–704 (2014).
Geyer, P. E. et al. Plasma Proteome Profiling to detect and avoid sample-related biases in biomarker studies. EMBO Mol. Med. 11, e10427 (2019).
Piehowski, P. D. et al. Sources of technical variability in quantitative LC-MS proteomics: human brain tissue sample analysis. J. Proteome Res. 12, 2128–2137 (2013).
Zhang, T. et al. Block design with common reference samples enables robust large-scale label-free quantitative proteome profiling. J. Proteome Res. 19, 2863–2872 (2020).
Collins, B. C. et al. Multi-laboratory assessment of reproducibility, qualitative and quantitative performance of SWATH-mass spectrometry. Nat. Commun. 8, 291 (2017).
Baba, A., Fujita, T. & Tamura, N. Sexual dimorphism of the fifth component of mouse complement. J. Exp. Med. 160, 411–419 (1984).
Maragno, A. L. et al. ISG15 modulates development of the erythroid lineage. PLoS ONE 6, e26068 (2011).
Shao, C. et al. Comprehensive analysis of individual variation in the urinary proteome revealed significant gender differences. Mol. Cell Proteomics 18, 1110–1122 (2019).
Dzieciatkowska, M., D’Alessandro, A., Hill, R. C. & Hansen, K. C. Plasma QconCATs reveal a gender-specific proteomic signature in apheresis platelet plasma supernatants. J. Proteomics 120, 1–6 (2015).
Yang, J. et al. Proteomics reveals intersexual differences in the rat brain hippocampus. Anat. Rec. (Hoboken) 296, 462–469 (2013).
Stehle, J. R. Jr et al. Mass spectrometry identification of circulating alpha-1-B glycoprotein, increased in aged female C57BL/6 mice. Biochim. Biophys. Acta 1770, 79–86 (2007).
Gui, Y., Silha, J. V. & Murphy, L. J. Sexual dimorphism and regulation of resistin, adiponectin, and leptin expression in the mouse. Obes. Res. 12, 1481–1491 (2004).
Tunstall, A. M., Merriman, J. M., Milne, I. & James, K. Normal and pathological serum levels of alpha2-macroglobulins in men and mice. J. Clin. Pathol. 28, 133–139 (1975).
Kotimaa, J. et al. Sex matters: systemic complement activity of female C57BL/6J and BALB/cJ mice is limited by serum terminal pathway components. Mol. Immunol. 76, 13–21 (2016).
Shand, B. I., Scott, R. S., Elder, P. A. & George, P. M. Plasma adiponectin in overweight, nondiabetic individuals with or without insulin resistance. Diabetes Obes. Metab. 5, 349–353 (2003).
Minnich, M., Kueppers, F. & James, H. Alpha-1-antitrypsin from mouse serum isolation and characterization. Comp. Biochem. Physiol. B 78, 413–419 (1984).
Adams, J. M. et al. Somatostatin is essential for the sexual dimorphism of GH secretion, corticosteroid-binding globulin production, and corticosterone levels in mice. Endocrinology 156, 1052–1065 (2015).
Koshibu, K. & Levitt, P. Sex differences in expression of transforming growth factor-alpha and epidermal growth factor receptor mRNA in Waved-1 and C57Bl6 mice. Neuroscience 134, 877–887 (2005).
Dill-Garlow, R., Chen, K. E. & Walker, A. M. Sex differences in mouse popliteal lymph nodes. Sci. Rep. 9, 965 (2019).
Lamason, R. et al. Sexual dimorphism in immune response genes as a function of puberty. BMC Immunol. 7, 2 (2006).
van Nas, A. et al. Elucidating the role of gonadal hormones in sexually dimorphic gene coexpression networks. Endocrinology 150, 1235–1249 (2009).
Clodfelter, K. H. et al. Role of STAT5a in regulation of sex-specific gene expression in female but not male mouse liver revealed by microarray analysis. Physiol. Genomics 31, 63–74 (2007).
Bain, C. C. et al. Rate of replenishment and microenvironment contribute to the sexually dimorphic phenotype and function of peritoneal macrophages. Sci. Immunol. 5, eabc4466 (2020).
Boulard, M. et al. Histone variant macroH2A1 deletion in mice causes female-specific steatosis. Epigenet. Chromatin 3, 8 (2010).
Curran, A. M. et al. Sexual dimorphism, age, and fat mass are key phenotypic drivers of proteomic signatures. J. Proteome Res. 16, 4122–4133 (2017).
Karp, N. A. et al. Prevalence of sexual dimorphism in mammalian phenotypic traits. Nat. Commun. 8, 15475 (2017).
Munoz-Fuentes, V. et al. The International Mouse Phenotyping Consortium (IMPC): a functional catalogue of the mammalian genome that informs conservation. Conserv. Genet. 19, 995–1005 (2018).
de Angelis, M. H. et al. Analysis of mammalian gene function through broad-based phenotypic screens across a consortium of mouse clinics. Nat. Genet. 47, 969–978 (2015).
Jassal, B. et al. The reactome pathway knowledgebase. Nucleic Acids Res. 48, D498–D503 (2020).
Gordon, S. M. et al. A comparison of the mouse and human lipoproteome: suitability of the mouse model for studies of human lipoproteins. J. Proteome Res. 14, 2686–2695 (2015).
Kurbatova, N., Mason, J. C., Morgan, H., Meehan, T. F. & Karp, N. A. PhenStat: a tool kit for standardized analysis of high throughput phenotypic data. PLoS ONE 10, e0131274 (2015).
Tibshirani, R. Regression shrinkage and selection via the Lasso. J. R. Stat. Soc. B 58, 267–288 (1996).
Tilburg, J. et al. Plasma protein signatures of a murine venous thrombosis model and Slc44a2 knockout mice using quantitative-targeted proteomics. Thromb. Haemost. 120, 423–436 (2020).
Boyle, E. I. et al. GO::TermFinder–open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinformatics 20, 3710–3715 (2004).
The Gene Ontology C. The Gene Ontology Resource: 20 years and still GOing strong. Nucleic Acids Res. 47, D330–D338 (2019).
Liberzon, A. et al. The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst. 1, 417–425 (2015).
Kanehisa, M. & Goto, S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30 (2000).
Schriml, L. M. et al. Human Disease Ontology 2018 update: classification, content and workflow expansion. Nucleic Acids Res. 47, D955–D962 (2019).
Pinero, J. et al. The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic Acids Res. 48, D845–D855 (2020).
Dhammi, I. K. & Kumar, S. Medical subject headings (MeSH) terms. Indian J. Orthop. 48, 443–444 (2014).
Bubeck, D. et al. Structure of human complement C8, a precursor to membrane attack. J. Mol. Biol. 405, 325–330 (2011).
Janeway, C. Immunobiology: The Immune System in Health and Disease 6th edn (Garland Science, 2005).
Schlesinger, M., Broman, I. & Lugassy, G. The complement system is defective in chronic lymphatic leukemia patients and in their healthy relatives. Leukemia 10, 1509–1513 (1996).
Mayilyan, K. R. Complement genetics, deficiencies, and disease associations. Protein Cell 3, 487–496 (2012).
Kugelberg, E., Gollan, B. & Tang, C. M. Mechanisms in Neisseria meningitidis for resistance against complement-mediated killing. Vaccine 26, I34–139 (2008).
Robinson, E. A. & Appella, E. Amino acid sequence of a mouse myeloma immunoglobin heavy chain (MOPC 47 A) with a 100-residue deletion. J. Biol. Chem. 254, 11418–11430 (1979).
Arora, A. S. et al. The role of cellular prion protein in lipid metabolism in the liver. Prion 14, 95–108 (2020).
Murphy, S. et al. Proteomic profiling of liver tissue from the mdx-4cv mouse model of Duchenne muscular dystrophy. Clin. Proteomics 15, 34 (2018).
Wei, X., Herbst, A., Ma, D., Aiken, J. & Li, L. A quantitative proteomic approach to prion disease biomarker research: delving into the glycoproteome. J. Proteome Res. 10, 2687–2702 (2011).
Platt, F. M., d’Azzo, A., Davidson, B. L., Neufeld, E. F. & Tifft, C. J. Lysosomal storage diseases. Nat. Rev. Dis. Primers 4, 27 (2018).
Louwette, S. et al. NPC1 defect results in abnormal platelet formation and function: studies in Niemann-Pick disease type C1 patients and zebrafish. Hum. Mol. Genet. 22, 61–73 (2013).
Spiegel, R. et al. The clinical spectrum of fetal Niemann-Pick type C. Am. J. Med. Genet. A 149A, 446–450 (2009).
Tseng, W. C. et al. Modeling Niemann-Pick disease type C1 in zebrafish: a robust platform for in vivo screening of candidate therapeutic compounds. Dis. Model Mech. 11, dmm034165 (2018).
Liao, Y. J. et al. Characterization of Niemann-Pick Type C2 protein expression in multiple cancers using a novel NPC2 monoclonal antibody. PLoS ONE 8, e77586 (2013).
Oulas, A., Minadakis, G., Zachariou, M. & Spyrou, G. M. Selecting variants of unknown significance through network-based gene-association significantly improves risk prediction for disease-control cohorts. Sci. Rep. 9, 3266 (2019).
Schulz, W. L., Tormey, C. A. & Torres, R. Computational approach to annotating variants of unknown significance in clinical next generation sequencing. Lab. Med. 46, 285–289 (2015).
Lewis, C. M. & Vassos, E. Polygenic risk scores: from research tools to clinical instruments. Genome Med. 12, 44 (2020).
White, J. K. et al. Genome-wide generation and systematic phenotyping of knockout mice reveals new roles for many genes. Cell 154, 452–464 (2013).
Verot, L. et al. Niemann-Pick C disease: functional characterization of three NPC2 mutations and clinical and molecular update on patients with NPC2. Clin. Genet. 71, 320–330 (2007).
Vanier, M. T. & Millat, G. Niemann-Pick disease type C. Clin. Genet. 64, 269–281 (2003).
Adissu, H. A. et al. Histopathology reveals correlative and unique phenotypes in a high-throughput mouse phenotyping screen. Dis. Model. Mech. 7, 515–524 (2014).
Mohammed, Y. et al. PeptidePicker: a scientific workflow with web interface for selecting appropriate peptides for targeted proteomics experiments. J. Proteomics 106, 151–161 (2014).
Percy, A. J., Chambers, A. G., Yang, J., Hardie, D. B. & Borchers, C. H. Advances in multiplexed MRM-based protein biomarker quantitation toward clinical utility. Biochim. Biophys. Acta 1844, 917–926 (2014).
LeBlanc, A. et al. Multiplexed MRM-based protein quantitation using two different stable isotope-labeled peptide isotopologues for calibration. J. Proteome Res. 16, 2527–2536 (2017).
Mohammed, Y., Pan, J., Zhang, S., Han, J. & Borchers, C. H. ExSTA: External Standard Addition Method for Accurate High-Throughput Quantitation in Targeted Proteomics Experiments. Proteomics Clin. Appl. 12, 1600180 (2018).
MacLean, B. et al. Skyline: an open source document editor for creating and analyzing targeted proteomics experiments. Bioinformatics 26, 966–968 (2010).
Hadley, W. Ggplot2 (Springer Science+Business Media, LLC, 2016).
Friedman, J., Hastie, T. & Tibshirani, R. Regularization paths for generalized linear models via coordinate descent. J. Stat. Softw. 33, 1–22 (2010).
Yu, G., Wang, L. G., Han, Y. & He, Q. Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS 16, 284–287 (2012).
Acknowledgements
The University of Victoria-Genome BC Proteomics Center is grateful to Genome Canada and Genome British Columbia for financial support (project codes 181-MRM-GAPP, 204PRO, 214PRO, 264PRO, 234DMP, 282PQP). The Centre for Phenogenomics acknowledges funding support from Genome Canada and Ontario Genomics (OGI-051) and from the National Institutes of Health (USA) through subawards from grants U54HG006364, U42 OD011175, and UM1 OD023221. The study was also supported by the MegaGrant of the Ministry of Science and Higher Education of the Russian Federation (Agreement with Skolkovo Institute of Science and Technology No. 075-10-2019-083).
Author information
Authors and Affiliations
Contributions
C.M.K. and C.H.B. conceived the idea of the study; Y.M., S.A.M., and D.S. designed the study; L.M.J.N. produced the mice; M.G. and A.M.F. collected the samples and did the in-life phenoytping; S.A.M., H.P., and J.Y. carried out the proteomic experiments; S.A.M. and J.Y. processed the proteomic data; H.A.A. and C.M.K. did the histopathological analysis and imaging; Y.M. performed the data analysis and visualization and generated the figures and tables; K.C.K.L., C.M.K. and C.H.B. provided financial support; Y.M. and S.A.M. wrote the first draft of the manuscript; all authors read and contributed to the final version of the manuscript.
Corresponding authors
Ethics declarations
Competing interests
C.H.B. is the Chief Scientific Officer of MRM Proteomics, Inc., the co-founder and Chief Technology Officer of Creative Molecules, Inc., and Chief Technology Officer of Molecular You. The remaining authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Mohammed, Y., Michaud, S.A., Pětrošová, H. et al. Proteotyping of knockout mouse strains reveals sex- and strain-specific signatures in blood plasma. npj Syst Biol Appl 7, 25 (2021). https://doi.org/10.1038/s41540-021-00184-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41540-021-00184-8