A major challenge facing the genetics of autism spectrum disorders (ASDs) is the large and growing number of candidate risk genes and gene variants of unknown functional significance. Here, we usedCaenorhabditis elegansto systematically functionally characterize ASD-associated genes in vivo. Using our custom machine vision system, we quantified 26 phenotypes spanning morphology, locomotion, tactile sensitivity, and habituation learning in 135 strains each carrying a mutation in an ortholog of an ASD-associated gene. We identified hundreds of genotype–phenotype relationships ranging from severe developmental delays and uncoordinated movement to subtle deficits in sensory and learning behaviors. We clustered genes by similarity in phenomic profiles and used epistasis analysis to discover parallel networks centered onCHD8•chd-7andNLGN3•nlg-1that underlie mechanosensory hyperresponsivity and impaired habituation learning. We then leveraged our data for in vivo functional assays to gauge missense variant effect. Expression of wild-type NLG-1 innlg-1mutantC. elegansrescued their sensory and learning impairments. Testing the rescuing ability of conserved ASD-associated neuroligin variants revealed varied partial loss of function despite proper subcellular localization. Finally, we used CRISPR-Cas9 auxin-inducible degradation to determine that phenotypic abnormalities caused by developmental loss of NLG-1 can be reversed by adult expression. This work charts the phenotypic landscape of ASD-associated genes, offers in vivo variant functional assays, and potential therapeutic targets for ASD.

Functional variomics provides the foundation for personalized medicine by linking genetic variation to disease expression, outcome and treatment, yet its utility is dependent on appropriate assays to evaluate mutation impact on protein function. To fully assess the effects of 106 missense and nonsense variants of PTEN associated with autism spectrum disorder, somatic cancer and PTEN hamartoma syndrome (PHTS), we take a deep phenotypic profiling approach using 18 assays in 5 model systems spanning diverse cellular environments ranging from molecular function to neuronal morphogenesis and behavior. Variants inducing instability occur across the protein, resulting in partial-to-complete loss-of-function (LoF), which is well correlated across models. However, assays are selectively sensitive to variants located in substrate binding and catalytic domains, which exhibit complete LoF or dominant negativity independent of effects on stability. Our results indicate that full characterization of variant impact requires assays sensitive to instability and a range of protein functions.

With rare exception, ciliated cells entering mitosis lose their cilia, thereby freeing basal bodies to serve as centrosomes in the formation of high-fidelity mitotic spindles. Cilia can be lost by shedding or disassembly, but either way, it appears that the final release may be via a coordinated severing of the nine axonemal outer doublet microtubules linking the basal body to the ciliary transition zone. Little is known about the mechanism or regulation of this important process. The stress-induced deflagellation response of Chlamydomonas provides a basis to identifying key players in axonemal severing. In an earlier screen we uncovered multiple alleles for each of three deflagellation genes, ADF1, FA1, and FA2. Products of the two FA genes localize to the site of axonemal severing and encode a scaffolding protein and a member of the NIMA-related family of ciliary-cell cycle kinases. The identity of the ADF1 gene remained elusive. Here, we report a new screen using a mutagenesis that yields point mutations in Chlamydomonas, an enhanced screening methodology, and whole genome sequencing. We isolated numerous new alleles of the three known genes, and one or two alleles each of at least four new genes. We identify ADF1 as a TRP ion channel, which we suggest may reside at the flagellar transition zone.

As sequencing becomes more economical, we are identifying sequence variations in the population faster than ever. For diseaseassociated genes, it is imperative that we differentiate a sequence variant as either benign or pathogenic, such that the appropriate therapeutic interventions or surveillance can be implemented. PTEN is a frequently mutated tumor suppressor that has been linked to the PTEN hamartoma tumor syndrome. Although the domain structure of PTEN and the functional impact of a number of its most common tumor-linked mutations have been characterized, there is a lack of information about many recently identified clinical variants. To address this challenge, we developed a cell-based assay that utilizes a premalignant phenotype of normal mammary epithelial cells lacking PTEN. We measured the ability of PTEN variants to rescue the spheroid formation phenotype of PTEN À/À MCF10A cells maintained in suspension. As proof of concept, we functionalized 47 missense variants using this assay, only 19 of which have clear classifications in ClinVar. We utilized a machine learning model trained with annotated genotypic data to classify variants as benign or pathogenic based on our functional scores. Our model predicted with high accuracy that loss of PTEN function was indicative of pathogenicity. We also determined that the pathogenicity of certain variants may have arisen from reduced stability of the protein product. Overall, this assay outperformed computational predictions, was scalable, and had a short run time, serving as an ideal alternative for annotating the clinical significance of cancer-associated PTEN variants.

SYNGAP1 is a Ras and Rap GTPase with important roles in regulating excitatory synaptic plasticity. While many SYNGAP1 missense and nonsense mutations have been associated with intellectual disability, epilepsy, schizophrenia and autism spectrum disorder (ASD), there are many variants of unknown significance (VUS). In this report, we characterize 58 variants in nine assays that examine multiple aspects of SYNGAP1 function. Specifically, we used multiplex phospho-flow cytometry to measure the impact of variants on pERK, pGSK3β and pCREB and high-content imaging to examine their subcellular localization. We find variants ranging from complete loss-of-function (LoF) to wildtype (WT)-like in their ability to regulate pERK and pGSK3β, while all variants retain at least partial ability to regulate pCREB. Interestingly, our assays reveal that a high percentage of variants located within the disordered domain of unknown function that makes up the C-terminal half of SYNGAP1 exhibited LoF, compared to the more well studied catalytic domain. Moreover, we find protein instability to be a major contributor to dysfunction only for two missense variants both located within the catalytic domain. Using highcontent imaging, we find variants with nuclear enrichment/exclusion and aberrant nuclear speckle localization. These variants are primarily located within the C2 domain known to mediate membrane lipid interactions. We find that mislocalization is distinct from altered catalytic activity, highlighting multiple independent molecular mechanisms underlying variant dysfunction. Our multidimensional dataset allows clustering of variants based on functional phenotypes and provides high-confidence pathogenicity classification..

ABSTRACTFunctional variomics provides the foundation for personalized medicine by linking genetic variation to disease expression, outcome and treatment, yet its utility is dependent on appropriate assays to evaluate mutation impact on protein function. To fully assess the effects of 106 missense and nonsense variants of PTEN associated with autism spectrum disorder, somatic cancer and PHTS, we take a deep phenotypic profiling approach using 18 assays in 5 model systems spanning diverse cellular environments ranging from molecular function to neuronal morphogenesis and behavior. Variants inducing instability occurred across the protein, resulting in partial to complete loss of function (LoF), which was well correlated across models. However, assays were selectively sensitive to variants located in substrate binding and catalytic domains, which exhibited complete LoF or dominant negativity independent of effects on stability. Our results indicate that full characterization of variant impact requires assays sensitive to instability and a range of protein functions.

Advances in sequencing technology have led to an explosion in the number of known genetic variants of human genes. A major challenge is to now determine which of these variants contribute to diseases as a result of their effect on gene function. Here, we describe a generic approach using the yeast Saccharomyces cerevisiae to quickly develop gene-specific in vivo assays that can be used to quantify the level of function of a genetic variant. Using synthetic dosage lethality screening, ‘sentinel’ yeast strains are identified that are sensitive to overexpression of a human disease gene. Variants of the gene can then be functionalized in a high-throughput fashion through simple growth assays using solid or liquid media. Sentinel interaction mapping (SIM) has the potential to create functional assays for the large majority of human disease genes that do not have a yeast orthologue. Using the tumour suppressor gene PTEN as an example, we show that SIM assays can provide a fast and economical means to screen a large number of genetic variants.

A major challenge facing the genetics of Autism Spectrum Disorders (ASD) is the large and growing number of candidate risk genes and gene variants of unknown functional significance.Here, we used Caenorhabditis elegans to systematically functionally characterize ASDassociated genes in vivo. Using our custom machine vision system we quantified 26 phenotypes spanning morphology, locomotion, tactile sensitivity, and habituation learning in 87 strains each carrying a mutation in an ortholog of an ASD-associated gene. We identified hundreds of novel genotype-phenotype relationships ranging from severe developmental delays and uncoordinated movement to subtle deficits in sensory and learning behaviors. We clustered genes by similarity in phenomic profiles and used epistasis analysis to discover parallel networks centered on CHD8•chd-7 and NLGN3•nlg-1 that underlie mechanosensory hyper-responsivity and impaired habituation learning. We then leveraged our data for in vivo functional assays to gauge missense variant effect. Expression of wild-type NLG-1 in nlg-1 mutant C. elegans rescued their sensory and learning impairments. Testing the rescuing ability of all conserved ASD-associated neuroligin variants revealed varied partial loss-of-function despite proper subcellular localization. Finally, we used CRISPR-Cas9 auxin inducible degradation to determine that phenotypic abnormalities caused by developmental loss of NLG-1 can be reversed by adult expression. This work charts the phenotypic landscape of ASDassociated genes, offers novel in vivo variant functional assays, and potential therapeutic targets for ASD.