Single-cell, whole-embryo phenotyping of mammalian developmental disorders.

Mouse models are a critical tool for studying human diseases, particularly developmental disorders1. However, conventional approaches for phenotyping may fail to detect subtle defects throughout the developing mouse2. Here we set out to establish single-cell RNA sequencing of the whole embryo as a scalable platform for the systematic phenotyping of mouse genetic models. We applied combinatorial indexing-based single-cell RNA sequencing3 to profile 101 embryos of 22 mutant and 4 wild-type genotypes at embryonic day 13.5, altogether profiling more than 1.6 million nuclei. The 22 mutants represent a range of anticipated phenotypic severities, from established multisystem disorders to deletions of individual regulatory regions4,5. We developed and applied several analytical frameworks for detecting differences in composition and/or gene expression across 52 cell types or trajectories. Some mutants exhibit changes in dozens of trajectories whereas others exhibit changes in only a few cell types. We also identify differences between widely used wild-type strains, compare phenotyping of gain- versus loss-of-function mutants and characterize deletions of topological associating domain boundaries. Notably, some changes are shared among mutants, suggesting that developmental pleiotropy might be 'decomposable' through further scaling of this approach. Overall, our findings show how single-cell profiling of whole embryos can enable the systematic molecular and cellular phenotypic characterization of mouse mutants with unprecedented breadth and resolution.

TRPS1 maintains luminal progenitors in the mammary gland by repressing SRF/MRTF activity.

The transcription factor TRPS1 is a context-dependent oncogene in breast cancer. In the mammary gland, TRPS1 activity is restricted to the luminal population and is critical during puberty and pregnancy. Its function in the resting state remains however unclear. To evaluate whether it could be a target for cancer therapy, we investigated TRPS1 function in the healthy adult mammary gland using a conditional ubiquitous depletion mouse model where long-term depletion does not affect fitness. Using transcriptomic approaches, flow cytometry and functional assays, we show that TRPS1 activity is essential to maintain a functional luminal progenitor compartment. This requires the repression of both YAP/TAZ and SRF/MRTF activities. TRPS1 represses SRF/MRTF activity indirectly by modulating RhoA activity. Our work uncovers a hitherto undisclosed function of TRPS1 in luminal progenitors intrinsically linked to mechanotransduction in the mammary gland. It may also provide new insights into the oncogenic functions of TRPS1 as luminal progenitors are likely the cells of origin of many breast cancers.

A rationale for considering heart/brain axis control in neuropsychiatric disease.

Neuropsychiatric diseases (NPD) represent a significant global disease burden necessitating innovative approaches to pathogenic understanding, biomarker identification and therapeutic strategy. Emerging evidence implicates heart/brain axis malfunction in NPD etiology, particularly via the autonomic nervous system (ANS) and brain central autonomic network (CAN) interaction. This heart/brain inter-relationship harbors potentially novel NPD diagnosis and treatment avenues. Nevertheless, the lack of multidisciplinary clinical approaches as well as a limited appreciation of molecular underpinnings has stymied progress. Large-scale preclinical multi-systemic functional data can therefore provide supplementary insight into CAN and ANS interaction. We here present an overview of the heart/brain axis in NPD and establish a unique rationale for utilizing a preclinical cardiovascular disease risk gene set to glean insights into heart/brain axis control in NPD. With a top-down approach focusing on genes influencing electrocardiogram ANS function, we combined hierarchical clustering of corresponding regional CAN expression data and functional enrichment analysis to reveal known and novel molecular insights into CAN and NPD. Through 'support vector machine' inquiries for classification and literature validation, we further pinpointed the top 32 genes highly expressed in CAN brain structures altering both heart rate/heart rate variability (HRV) and behavior. Our observations underscore the potential of HRV/hyperactivity behavior as endophenotypes for multimodal disease biomarker identification to index aberrant executive brain functioning with relevance for NPD. This work heralds the potential of large-scale preclinical functional genetic data for understanding CAN/ANS control and introduces a stepwise design leveraging preclinical data to unearth novel heart/brain axis control genes in NPD.

Knockout of the Complex III subunit Uqcrh causes bioenergetic impairment and cardiac contractile dysfunction.

Ubiquinol cytochrome c reductase hinge protein (UQCRH) is required for the electron transfer between cytochrome c1 and c of the mitochondrial cytochrome bc1 Complex (CIII). A two-exon deletion in the human UQCRH gene has recently been identified as the cause for a rare familial mitochondrial disorder. Deletion of the corresponding gene in the mouse (Uqcrh-KO) resulted in striking biochemical and clinical similarities including impairment of CIII, failure to thrive, elevated blood glucose levels, and early death. Here, we set out to test how global ablation of the murine Uqcrh affects cardiac morphology and contractility, and bioenergetics. Hearts from Uqcrh-KO mutant mice appeared macroscopically considerably smaller compared to wildtype littermate controls despite similar geometries as confirmed by transthoracic echocardiography (TTE). Relating TTE-assessed heart to body mass revealed the development of subtle cardiac enlargement, but histopathological analysis showed no excess collagen deposition. Nonetheless, Uqcrh-KO hearts developed pronounced contractile dysfunction. To assess mitochondrial functions, we used the high-resolution respirometer NextGen-O2k allowing measurement of mitochondrial respiratory capacity through the electron transfer system (ETS) simultaneously with the redox state of ETS-reactive coenzyme Q (Q), or production of reactive oxygen species (ROS). Compared to wildtype littermate controls, we found decreased mitochondrial respiratory capacity and more reduced Q in Uqcrh-KO, indicative for an impaired ETS. Yet, mitochondrial ROS production was not generally increased. Taken together, our data suggest that Uqcrh-KO leads to cardiac contractile dysfunction at 9 weeks of age, which is associated with impaired bioenergetics but not with mitochondrial ROS production. Global ablation of the Uqcrh gene results in functional impairment of CIII associated with metabolic dysfunction and postnatal developmental arrest immediately after weaning from the mother. Uqcrh-KO mice show dramatically elevated blood glucose levels and decreased ability of isolated cardiac mitochondria to consume oxygen (O2). Impaired development (failure to thrive) after weaning manifests as a deficiency in the gain of body mass and growth of internal organ including the heart. The relative heart mass seemingly increases when organ mass calculated from transthoracic echocardiography (TTE) is normalized to body mass. Notably, the heart shows no signs of collagen deposition, yet does develop a contractile dysfunction reflected by a decrease in ejection fraction and fractional shortening.

Echo2Pheno: a deep-learning application to uncover echocardiographic phenotypes in conscious mice.

Echocardiography, a rapid and cost-effective imaging technique, assesses cardiac function and structure. Despite its popularity in cardiovascular medicine and clinical research, image-derived phenotypic measurements are manually performed, requiring expert knowledge and training. Notwithstanding great progress in deep-learning applications in small animal echocardiography, the focus has so far only been on images of anesthetized rodents. We present here a new algorithm specifically designed for echocardiograms acquired in conscious mice called Echo2Pheno, an automatic statistical learning workflow for analyzing and interpreting high-throughput non-anesthetized transthoracic murine echocardiographic images in the presence of genetic knockouts. Echo2Pheno comprises a neural network module for echocardiographic image analysis and phenotypic measurements, including a statistical hypothesis-testing framework for assessing phenotypic differences between populations. Using 2159 images of 16 different knockout mouse strains of the German Mouse Clinic, Echo2Pheno accurately confirms known cardiovascular genotype-phenotype relationships (e.g., Dystrophin) and discovers novel genes (e.g., CCR4-NOT transcription complex subunit 6-like, Cnot6l, and synaptotagmin-like protein 4, Sytl4), which cause altered cardiovascular phenotypes, as verified by H&E-stained histological images. Echo2Pheno provides an important step toward automatic end-to-end learning for linking echocardiographic readouts to cardiovascular phenotypes of interest in conscious mice.

A review of standardized high-throughput cardiovascular phenotyping with a link to metabolism in mice.

Cardiovascular diseases cause a high mortality rate worldwide and represent a major burden for health care systems. Experimental rodent models play a central role in cardiovascular disease research by effectively simulating human cardiovascular diseases. Using mice, the International Mouse Phenotyping Consortium (IMPC) aims to target each protein-coding gene and phenotype multiple organ systems in single-gene knockout models by a global network of mouse clinics. In this review, we summarize the current advances of the IMPC in cardiac research and describe in detail the diagnostic requirements of high-throughput electrocardiography and transthoracic echocardiography capable of detecting cardiac arrhythmias and cardiomyopathies in mice. Beyond that, we are linking metabolism to the heart and describing phenotypes that emerge in a set of known genes, when knocked out in mice, such as the leptin receptor (Lepr), leptin (Lep), and Bardet-Biedl syndrome 5 (Bbs5). Furthermore, we are presenting not yet associated loss-of-function genes affecting both, metabolism and the cardiovascular system, such as the RING finger protein 10 (Rfn10), F-box protein 38 (Fbxo38), and Dipeptidyl peptidase 8 (Dpp8). These extensive high-throughput data from IMPC mice provide a promising opportunity to explore genetics causing metabolic heart disease with an important translational approach.

Knockout mouse models as a resource for the study of rare diseases.

Rare diseases (RDs) are a challenge for medicine due to their heterogeneous clinical manifestations and low prevalence. There is a lack of specific treatments and only a few hundred of the approximately 7,000 RDs have an approved regime. Rapid technological development in genome sequencing enables the mass identification of potential candidates that in their mutated form could trigger diseases but are often not confirmed to be causal. Knockout (KO) mouse models are essential to understand the causality of genes by allowing highly standardized research into the pathogenesis of diseases. The German Mouse Clinic (GMC) is one of the pioneers in mouse research and successfully uses (preclinical) data obtained from single-gene KO mutants for research into monogenic RDs. As part of the International Mouse Phenotyping Consortium (IMPC) and INFRAFRONTIER, the pan-European consortium for modeling human diseases, the GMC expands these preclinical data toward global collaborative approaches with researchers, clinicians, and patient groups.Here, we highlight proprietary genes that when deleted mimic clinical phenotypes associated with known RD targets (Nacc1, Bach2, Klotho alpha). We focus on recognized RD genes with no pre-existing KO mouse models (Kansl1l, Acsf3, Pcdhgb2, Rabgap1, Cox7a2) which highlight novel phenotypes capable of optimizing clinical diagnosis. In addition, we present genes with intriguing phenotypic data (Zdhhc5, Wsb2) that are not presently associated with known human RDs.This report provides comprehensive evidence for genes that when deleted cause differences in the KO mouse across multiple organs, providing a huge translational potential for further understanding monogenic RDs and their clinical spectrum. Genetic KO studies in mice are valuable to further explore the underlying physiological mechanisms and their overall therapeutic potential.

Mouse mutant phenotyping at scale reveals novel genes controlling bone mineral density.

The genetic landscape of diseases associated with changes in bone mineral density (BMD), such as osteoporosis, is only partially understood. Here, we explored data from 3,823 mutant mouse strains for BMD, a measure that is frequently altered in a range of bone pathologies, including osteoporosis. A total of 200 genes were found to significantly affect BMD. This pool of BMD genes comprised 141 genes with previously unknown functions in bone biology and was complementary to pools derived from recent human studies. Nineteen of the 141 genes also caused skeletal abnormalities. Examination of the BMD genes in osteoclasts and osteoblasts underscored BMD pathways, including vesicle transport, in these cells and together with in silico bone turnover studies resulted in the prioritization of candidate genes for further investigation. Overall, the results add novel pathophysiological and molecular insight into bone health and disease.

Creld1 regulates myocardial development and function.

CRELD1 (Cysteine-Rich with EGF-Like Domains 1) is a risk gene for non-syndromic atrioventricular septal defects in human patients. In a mouse model, Creld1 has been shown to be essential for heart development, particularly in septum and valve formation. However, due to the embryonic lethality of global Creld1 knockout (KO) mice, its cell type-specific function during peri- and postnatal stages remains unknown. Here, we generated conditional Creld1 KO mice lacking Creld1 either in the endocardium (KOTie2) or the myocardium (KOMyHC). Using a combination of cardiac phenotyping, histology, immunohistochemistry, RNA-sequencing, and flow cytometry, we demonstrate that Creld1 function in the endocardium is dispensable for heart development. Lack of myocardial Creld1 causes extracellular matrix remodeling and trabeculation defects by modulation of the Notch1 signaling pathway. Hence, KOMyHC mice die early postnatally due to myocardial hypoplasia. Our results reveal that Creld1 not only controls the formation of septa and valves at an early stage during heart development, but also cardiac maturation and function at a later stage. These findings underline the central role of Creld1 in mammalian heart development and function.

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

MOTIVATION: High-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. RESULTS: Here 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 IMPLEMENTATION: The 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.

A comprehensive and comparative phenotypic analysis of the collaborative founder strains identifies new and known phenotypes.

The collaborative cross (CC) is a large panel of mouse-inbred lines derived from eight founder strains (NOD/ShiLtJ, NZO/HILtJ, A/J, C57BL/6J, 129S1/SvImJ, CAST/EiJ, PWK/PhJ, and WSB/EiJ). Here, we performed a comprehensive and comparative phenotyping screening to identify phenotypic differences and similarities between the eight founder strains. In total, more than 300 parameters including allergy, behavior, cardiovascular, clinical blood chemistry, dysmorphology, bone and cartilage, energy metabolism, eye and vision, immunology, lung function, neurology, nociception, and pathology were analyzed; in most traits from sixteen females and sixteen males. We identified over 270 parameters that were significantly different between strains. This study highlights the value of the founder and CC strains for phenotype-genotype associations of many genetic traits that are highly relevant to human diseases. All data described here are publicly available from the mouse phenome database for analyses and downloads.

Role of pancreatic cancer-derived exosomes in salivary biomarker development.

Recent studies have demonstrated that discriminatory salivary biomarkers can be readily detected upon the development of systemic diseases such as pancreatic cancer, breast cancer, lung cancer, and ovarian cancer. However, the utility of salivary biomarkers for the detection of systemic diseases has been undermined due to the absence of the biological and mechanistic rationale as to why distal diseases from the oral cavity would lead to the development of discriminatory biomarkers in saliva. Here, we examine the hypothesis that pancreatic tumor-derived exosomes are mechanistically involved in the development of pancreatic cancer-discriminatory salivary transcriptomic biomarkers. We first developed a pancreatic cancer mouse model that yielded discriminatory salivary biomarkers by implanting the mouse pancreatic cancer cell line Panc02 into the pancreas of the syngeneic host C57BL/6. The role of pancreatic cancer-derived exosomes in the development of discriminatory salivary biomarkers was then tested by engineering a Panc02 cell line that is suppressed for exosome biogenesis, implanting into the C56BL/6 mouse, and examining whether the discriminatory salivary biomarker profile was ablated or disrupted. Suppression of exosome biogenesis results in the ablation of discriminatory salivary biomarker development. This study supports that tumor-derived exosomes provide a mechanism in the development of discriminatory biomarkers in saliva and distal systemic diseases.

Canonical and Noncanonical Contribution of Thyroid Hormone Receptor Isoforms Alpha and Beta to Cardiac Hypertrophy and Heart Rate in Male Mice.

Background: Stimulation of ventricular hypertrophy and heart rate are two major cardiac effects of thyroid hormone (TH). The aim of this study was to determine in vivo which TH receptor (TR)-α or ß-and which mode of TR action-canonical gene expression or DNA-binding independent noncanonical action-mediate these effects. Methods: We compared global TRα and TRß knockout mice (TRαKO; TRßKO) with wild-type (WT) mice to determine the TR isoform responsible for T3 effects. The relevance of TR DNA binding was studied in mice with a mutation in the DNA-binding domain that selectively abrogates DNA binding and canonical TR action (TRαGS; TRßGS). Hearts were studied with echocardiography at baseline and after 7 weeks of T3 treatment. Gene expression was measured with real-time polymerase chain reaction. Heart rate was recorded with radiotelemetry transmitters for 7 weeks in untreated, hypothyroid, and T3-treated mice. Results: T3 induced ventricular hypertrophy in WT and TRßKO mice, but not in TRαKO mice. Hypertrophy was also induced in TRαGS mice. Thus, hypertrophy is mostly mediated by noncanonical TRα action. Similarly, repression of Mhy7 occurred in WT and TRαGS mice. Basal heart rate was largely dependent on canonical TRα action. But responsiveness to hypothyroidism and T3 treatment as well as expression of pacemaker gene Hcn2 were still preserved in TRαKO mice, demonstrating that TRß could compensate for absence of TRα. Conclusions: T3-induced cardiac hypertrophy could be attributed to noncanonical TRα action, whereas heart rate regulation was mediated by canonical TRα action. TRß could substitute for canonical but not noncanonical TRα action.

Comparative Phenotyping of Mice Reveals Canonical and Noncanonical Physiological Functions of TRα and TRß.

Thyroid hormone (TH) effects are mediated through TH receptors (TRs), TRα1, TRß1, and TRß2. The TRs bind to the DNA and regulate expression of TH target genes (canonical signaling). In addition, they mediate activation of signaling pathways (noncanonical signaling). Whether noncanonical TR action contributes to the spectrum of TH effects is largely unknown. The aim of this study was to attribute physiological effects to the TR isoforms and their canonical and noncanonical signaling. We conducted multiparameter phenotyping in male and female TR knockout mice (TRαKO, TRßKO), mice with disrupted canonical signaling due to mutations in the TR DNA binding domain (TRαGS, TRßGS), and their wild-type littermates. Perturbations in senses, especially hearing (mainly TRß with a lesser impact of TRα), visual acuity, retinal thickness (TRα and TRß), and in muscle metabolism (TRα) highlighted the role of canonical TR action. Strikingly, selective abrogation of canonical TR action often had little phenotypic consequence, suggesting that noncanonical TR action sufficed to maintain the wild-type phenotype for specific effects. For instance, macrocytic anemia, reduced retinal vascularization, or increased anxiety-related behavior were only observed in TRαKO but not TRαGS mice. Noncanonical TRα action improved energy utilization and prevented hyperphagia observed in female TRαKO mice. In summary, by examining the phenotypes of TRα and TRß knockout models alongside their DNA binding-deficient mutants and wild-type counterparts, we could establish that the noncanonical actions of TRα and TRß play a crucial role in modulating sensory, behavioral, and metabolic functions and, thus, contribute to the spectrum of physiological TH effects.

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.

AOX delays the onset of the lethal phenotype in a mouse model of Uqcrh (complex III) disease.

The alternative oxidase, AOX, provides a by-pass of the cytochrome segment of the mitochondrial respiratory chain when the chain is unavailable. AOX is absent from mammals, but AOX from Ciona intestinalis is benign when expressed in mice. Although non-protonmotive, so does not contribute directly to ATP production, it has been shown to modify and in some cases rescue phenotypes of respiratory-chain disease models. Here we studied the effect of C. intestinalis AOX on mice engineered to express a disease-equivalent mutant of Uqcrh, encoding the hinge subunit of mitochondrial respiratory complex III, which results in a complex metabolic phenotype beginning at 4-5 weeks, rapidly progressing to lethality within a further 6-7 weeks. AOX expression delayed the onset of this phenotype by several weeks, but provided no long-term benefit. We discuss the significance of this finding in light of the known and hypothesized effects of AOX on metabolism, redox homeostasis, oxidative stress and cell signaling. Although not a panacea, the ability of AOX to mitigate disease onset and progression means it could be useful in treatment.

The human salivary RNA transcriptome revealed by massively parallel sequencing.

BACKGROUND: Evaluation of the salivary transcriptome is an emerging diagnostic technology with discriminatory power for disease detection. This study explored massively parallel sequencing for providing nucleotide-level sequence information for each RNA in saliva. METHODS: Transcriptome profiling with the SOLiD™ system was applied to RNA isolated from unstimulated cell-free saliva (CFS) and whole saliva (WS) from healthy human volunteers. Sequenced reads were aligned to human genome build 18 and the Human Oral Microbiome Database (HOMD). RESULTS: Massively parallel sequencing enabled the acquisition of complete sequence information for each nucleotide position of the human salivary transcriptome through alignment to multiple sequence databases. Approximately 20%-25% of the sequenced reads from CFS aligned to the human genome, and approximately 30% of the sequenced reads aligned to the HOMD. We detected the expression of >4000 coding and noncoding genes in CFS and WS. Gene expression at different genome loci showed that the structural integrity of the transcripts for the annotated genes was preserved in saliva. CONCLUSIONS: A single measurement provided RNA sequence information of gene transcript abundance for both coding and noncoding RNAs and identified sequences from >400 different microbial species within a single sample. Contrary to previous data suggesting that salivary RNA is highly fragmented, in our study the structural integrity of RNA was preserved. The high degree of sequence alignment to annotated exons and introns for each of the respective reference genomes, with sequence coverage spanning the full length of the genes, provides strong evidence that the salivary transcriptome is a complex RNA network.

Salivary diagnostics: moving to the next level.

The evolution of salivary diagnostics has reached a new level toward the goal of using saliva as a powerful fluid for early detection and the first line of diagnosis for life-threatening diseases such as cancer, metabolic disorders, infections, and inflammatory diseases. Newly developed tools such as the novel saliva-based POCT and the SKB are helping to realize the goal of making salivary diagnostics available to clinicians worldwide. This is a unique moment where dentistry may be paving a new path for primary healthcare.

Post-synaptic scaffold protein TANC2 in psychiatric and somatic disease risk.

Understanding the shared genetic aetiology of psychiatric and medical comorbidity in neurodevelopmental disorders (NDDs) could improve patient diagnosis, stratification and treatment options. Rare tetratricopeptide repeat, ankyrin repeat and coiled-coil containing 2 (TANC2)-disrupting variants were disease causing in NDD patients. The post-synaptic scaffold protein TANC2 is essential for dendrite formation in synaptic plasticity and plays an unclarified but critical role in development. We here report a novel homozygous-viable Tanc2-disrupted function model in which mutant mice were hyperactive and had impaired sensorimotor gating consistent with NDD patient psychiatric endophenotypes. Yet, a multi-systemic analysis revealed the pleiotropic effects of Tanc2 outside the brain, such as growth failure and hepatocellular damage. This was associated with aberrant liver function including altered hepatocellular metabolism. Integrative analysis indicates that these disrupted Tanc2 systemic effects relate to interaction with Hippo developmental signalling pathway proteins and will increase the risk for comorbid somatic disease. This highlights how NDD gene pleiotropy can augment medical comorbidity susceptibility, underscoring the benefit of holistic NDD patient diagnosis and treatment for which large-scale preclinical functional genomics can provide complementary pleiotropic gene function information.