Rapid assessment of the temporal function and phenotypic reversibility of neurodevelopmental disorder risk genes in Caenorhabditis elegans

ABSTRACT Recent studies have indicated that some phenotypes caused by decreased function of select neurodevelopmental disorder (NDD) risk genes can be reversed by restoring gene function in adulthood. However, few of the hundreds of risk genes have been assessed for adult phenotypic reversibility. We developed a strategy to rapidly assess the temporal requirements and phenotypic reversibility of NDD risk gene orthologs using a conditional protein degradation system and machine-vision phenotypic profiling in Caenorhabditis elegans. We measured how degrading and re-expressing orthologs of EBF3, BRN3A and DYNC1H1 at multiple periods throughout development affect 30 morphological, locomotor, sensory and learning phenotypes. We found that phenotypic reversibility was possible for each gene studied. However, the temporal requirements of gene function and degree of rescue varied by gene and phenotype. This work highlights the critical need to assess multiple windows of degradation and re-expression and a large number of phenotypes to understand the many roles a gene can have across the lifespan. This work also demonstrates the benefits of using a high-throughput model system to prioritize NDD risk genes for re-expression studies in other organisms.


INTRODUCTION
Neurodevelopmental disorders such as Autism Spectrum Disorder (ASD) and Intellectual Disability (ID) are highly genetically heterogeneous and are accompanied by a range of cognitive and behavioural phenotypes including sensory processing and learning impairments (American Psychiatric Association, 2013; Boyle et al., 2011;de la Torre-Ubieta et al., 2016;Iakoucheva et al., 2019;Sanders et al., 2019).More severe cases of these disorders can cause significant challenges for affected individuals and their families (Boyle et al., 2011;de la Torre-Ubieta et al., 2016;Iakoucheva et al., 2019;Sanders et al., 2019).Recently, there has been remarkable progress in identifying genetic risk factors that contribute to diverse neurodevelopmental disorders, with hundreds of genes now implicated in ASD and ID alone (de la Torre-Ubieta et al., 2016;De Rubeis et al., 2014;Deciphering Developmental Disorders Study, 2015;Iakoucheva et al., 2019;Iossifov et al., 2014;Sanders et al., 2019Sanders et al., , 2015;;Satterstrom et al., 2020;Vissers et al., 2016).Variants in the majority of these genes (e.g.89/102 ASD-associated genes; 87%) are thought to confer risk through haploinsufficiency as the individual carries one loss-of-function allele with insufficient residual function from the remaining copy (Satterstrom et al., 2020).The identification of how these variants contribute to disorder pathology suggests re-expression therapies, where a second functional allele is introduced to restore protein levels to compensate for the decreased function of the faulty allele, could be a viable treatment option.
Historically, it was assumed that any treatment targeting neurodevelopmental disorder risk genes would need to be administered very early in development to be effective.This long-held assumption was challenged by reports that re-expression of several risk gene orthologs could reverse multiple altered neurophysiological and/or behavioural phenotypes in adult mice (Creson et al., 2019;Ehninger et al., 2008;Gao et al., 2020;Guy et al., 2007;Mei et al., 2016;Vogel-Ciernia et al., 2013;Zeier et al., 2009).In addition, inactivating orthologs of some of these risk genes in adult mice could also induce the phenotypic impairments previously associated only with altered gene function during development (Ehninger et al., 2008;Guy et al., 2007).Together, these findings suggest there may be a degree of temporal flexibility in the neurodevelopmental processes these gene contribute to, and that some genes typically associated with neurodevelopment may continue to have important functions well into adulthood.
The handful of reports that show the possibility of phenotypic reversibility with gene reactivation in adult mice offer critical insights into which genes are promising candidates for future re-expression-based therapies.However, because of limitations including cost, developmental rate, and technical difficulties (e.g.injection of viral vectors for large numbers of animals, etc.) very few risk genes have been tested for adult phenotypic reversibility.Further, the rapidly growing number of risk genes identified in recent years has exacerbated this problem and drastically increased the need for candidate prioritization to better direct research efforts.While most neurodevelopmental disorder risk genes are highly expressed early in pre-natal development (Jin et al., 2020;Parikshak et al., 2013;Satterstrom et al., 2020;Willsey et al., 2013), many continue to be expressed well into adulthood, and we currently do not know if the relationship between temporal expression patterns and inferred temporal functional windows are a significant predictor of whether a gene will be suitable for re-expression therapies.Assessing the phenotypic reversibility of neurodevelopmental disorder risk genes in more high-throughput model organisms offers the ability to rapidly screen a large number of genes to aid in prioritizing risk genes for further study in mammalian models.
The nematode Caenorhabditis elegans offers multiple advantages to systematically assess the temporal requirements and phenotypic reversibility of neurodevelopmental disorder risk genes.C. elegans have orthologs for a high number of neurodevelopmental risk genes (e.g.>80% of high-confidence ASD risk genes (McDiarmid et al., 2020)) and these genes have repeatedly been shown to be so well-conserved that in many cases expression of the human risk gene can compensate for loss of the C. elegans ortholog (Kaletta and Hengartner, 2006;Levitan et al., 1996;McDiarmid et al., 2018;Post et al., 2020).C. elegans have rapid development, growing from egg through well-characterized larval stages (L1, L2, L3, and L4) to egg-laying adults within 3 days.The hermaphroditic reproduction of C. elegans enable large colonies of genetically identical animals to be rapidly and cheaply cultivated.In addition, multiple genomic tools are available to precisely control the spatial and/or temporal activity of genes in vivo (Ashley et al., 2021;Au et al., 2019;Dickinson and Goldstein, 2016;Nance and Frøkjaer-Jensen, 2019;Wang et al., 2017;Zhang et al., 2015).Lastly, the phenotypic profiles of hundreds of animals can be simultaneously assessed using automated tracking systems which capture and analyze the impact of genetic perturbations on morphological, sensory, and learning behaviours in real time (Husson, Steuer Costa, Schmitt, and Gottschalk, 2012;McDiarmid et al., 2018;Swierczek et al., 2011).
A behaviour that has been increasingly used in high-throughput model organisms to investigate the biological function of neurodevelopmental disorder risk genes and functional impact of disorder-associated variants is habituation (Kepler et al., 2020).Habituation is a highly conserved form of non-associative learning observed as a decrement in responding to a repetitive stimulus (Rankin et al., 2009;Thompson and Spencer, 1966).Habituation is thought to be an important building block for higher cognitive functions and enable ongoing shifts in behavioral strategy (McDiarmid et al., 2019;Schmid et al., 2015).Alterations in the ability to habituate have been reported in ASD, ID, and Schizophrenia (McDiarmid et al., 2017) and are hypothesized to contribute to more complex behavioural symptoms (Green et al., 2015;Kavšek, 2004;Kleinhans et al., 2009;Massa and O'Desky, 2012;Williams et al., 2013).
Here, we developed a strategy to assess the phenotypic reversibility and temporal functional windows of three neurodevelopmental disorder risk genes in vivo using CRISPR-Cas9 Auxin Inducible degradation (AID) and high-throughput machine vision phenotyping (Fig. 1).We took advantage of the genetic tractability and rapid, stereotyped development of C. elegans to precisely investigate the effects of degrading and re-expressing the proteins of interest at multiple developmental time points in thousands of age-synchronized genetically identical animals.We found that some level of phenotypic reversibility was possible for each risk gene if re-expression occurred early in post-embryonic development, but only re-expressing EBF3•unc-3 and DYNC1H1•dhc-1 could reverse multiple phenotypic alterations later in life.More broadly, we provide an adaptable strategy and important examples/criteria that illustrate a path towards prioritizing neurodevelopmental disorder risk genes for further study and therapeutic development.

RESULTS
We selected three neurodevelopmental disorder risk genes (as identified by Simons Foundation Autism Research Initiative (Abrahams et al., 2013), Satterstorm et al. (2020), and/or literature search) for reversibility analysis based on the availability of C. elegans strains that contain a neurodevelopmental disorder risk gene ortholog tagged with the Auxin-Inducible Degron at the endogenous locus using CRISPR-Cas9 (see methods).The AID system relies on tagging the gene of interest with a short degron peptide tag as well as transgenic expression of TIR1 which is an E3 ubiquitin ligase typically found only in plants (Nishimura et al., 2009;Zhang et al., 2015).In the presence of the plant hormone Auxin, TIR1 can associate with the AID degron, adding a poly-ubiquitin chain to the protein of interest causing it to be degraded by the proteasome (Nishimura et al., 2009;Zhang et al., 2015) (Fig. 1A).We chose to use the AID system as it enables temporal control of protein degradation, which can be reversed by transferring populations of worms to culture plates without Auxin (McDiarmid et al., 2020;Zhang et al., 2015).
Since the AID degron is tagged to the endogenous loci, protein expression is restored using the native regulatory machinery, therefore bypassing the biological confounds associated with conventional approaches that rely on overexpression.Importantly, our lab and others have shown that Auxin exposure does not cause any overt effects on C. elegans morphology, locomotion, short-term learning, or mechanosensory processing phenotypes (McDiarmid et al., 2020;Zhang et al., 2015).We systematically assessed the functional consequence of multiple developmental degradation time windows in vivo by transferring animals on or off plates containing Auxin at precise time points in C. elegans development (Fig. 1B).We used our high-throughput machine vision tracking system, the Multi-Worm Tracker (MWT) (Swierczek et al., 2011), to quantify 26 phenotypes spanning morphology, baseline locomotion, sensory responding, and learning while animals were subjected to a short-term mechanosensory habituation behavioral paradigm (Fig. 1C & D).Our phenotypic features included multiple measures of mechanosensory responding and habituation learning, as both are disrupted across neurodevelopmental disorders (Green et al., 2019(Green et al., , 2015;;McDiarmid et al., 2017;Williams et al., 2013), and because we have previously shown that different components of a single behavioral response can be mediated by genetically dissociable underlying mechanisms (Ardiel et al., 2018;Kindt et al., 2007;McDiarmid et al., 2019;Randlett et al., 2019).Inclusion of a range of phenotypes not only aids in characterizing gene function across development, but also enables any unexpected phenotypes caused by protein reexpression to be captured.

The transcription factor EBF3•unc-3 displays a reciprocal pattern of phenotypic induction and reversibility across development
The first neurodevelopmental disorder risk gene we assessed was EBF3•unc-3, which is a highly conserved transcription factor.In C. elegans, unc-3 that acts to specify the identity of different neuronal classes by initiating and maintaining the expression of class-specific effector genes (Fig. 2A) (Kratsios et al., 2017;Prasad et al., 2008Prasad et al., , 1998)).Variants in EBF3 have been associated with multiple neurodevelopmental disorders including ASD and ID, and are thought to confer risk through haploinsufficiency or by interfering with DNA binding (Chao et al., 2017;Lopes et al., 2017;Sleven et al., 2017;Tanaka et al., 2017).In C. elegans, unc-3 loss-of-function results in severe locomotion and coordination defects caused by undifferentiated/abnormal identity of cholinergic motor neurons in the ventral nerve cord (Brenner, 1974;Feng et al., 2020;Kratsios et al., 2017;Prasad et al., 1998).
In our paradigm, continuous degradation of UNC-3 from egg through to adulthood produced uncoordinated locomotion, an inability to respond to mechanosensory stimuli, and severe alterations in several other morphological and behavioral phenotypes (Fig. 2B & D and     Fig 3).Early adult inactivation of unc-3 (achieved by beginning Auxin exposure at the final larval stage, L4) induced impairments for the majority of affected phenotypes, supporting previous findings that unc-3 function is continuously required throughout development (Feng et al., 2020;Kratsios et al., 2017;Li et al., 2020) (Fig. 2B).Importantly, impairments in mechanosensory response probability and several other altered phenotypes were partially rescued in animals when UNC-3 was degraded during development but was re-expressed from the endogenous locus starting in early adulthood (animals taken off Auxin immediately after L4) (Fig. 2C).Phenotypic change

D
Our initial findings of partial rescue of some phenotypes following protein re-expression starting at early adulthood motivated us to explore whether earlier restoration of unc-3 function would produce more effective rescue.Re-expression of UNC-3 during late post-embryonic development (ending Auxin exposure at L2) resulted in an almost complete rescue of impairments in response probability (Fig. 2E).Further, starting UNC-3 degradation during late post-embryonic development produced more severe impairments which were similar to the continuous degradation control, suggesting unc-3 plays a critical role between L2 and L4 for reversal probability (Fig. 2E).Lastly, we explored the phenotypic consequences of exposing or removing 3-day old animals (young adult-72 hrs post-hatch) from Auxin.Re-expression of UNC-3 in adulthood did not rescue impairments in response probability, reaffirming that the crucial functional period of unc-3 occurs during development.In line with this, starting UNC-3 degradation at 72 hrs post-hatch did not induce impairments, suggesting that the role of unc-3 in maintaining the expression of terminal identity genes throughout the lifespan may not be required for normal behaviour once the nervous system has fully developed.
Assessment across the three temporal conditions revealed a striking reciprocal pattern of unc-3 temporal function for response probability (Fig. 2E-G).The degree of phenotypic impairment for response probability induced by inactivation of unc-3 at a given developmental time point almost perfectly mirrored the degree of reversibility possible with re-expression (Fig.

2E-G).
While some other phenotypes followed this reciprocal pattern others did not (Fig. 3A-F).
For example, reversal duration seems to be mediated by a mechanism that acts earlier in development, with reversibility only possible with early re-expression (Fig. 3G-I).Taken together, these results clearly demonstrate that the temporal windows of gene function can drastically vary by phenotype, and that both the degree of reversibility and number of reversible phenotypes can be influenced by how early re-expression occurs in development.Importantly, these findings reveal that unc-3 expression can rescue multiple phenotypic alterations relatively late in life, prioritizing this gene as a candidate for further study.•unc-86 is a POU-type transcription factor that plays conserved roles in the initiation and maintenance neuronal identity across species (Badea et al., 2009;Serrano-Saiz et al., 2018;Xiang et al., 1995;Zou et al., 2012).Variants in POUF4/BRN3A have been associated with abnormal development of sensory structures, including auditory (Huang et al., 2001) and visual cells (Badea et al., 2009).In C. elegans, unc-86 is thought to be required across the lifespan to maintain the expression of terminal identity genes in multiple neuronal subtypes (Serrano-Saiz et al., 2018;Sze and Ruvkun, 2003) (Fig. 4A).
Analysis of 30 quantitative phenotypes revealed that continuous degradation of UNC-86 specifically impaired mechanosensory response probability (Fig. 4B).Animals were hyporesponsive to mechanosensory stimuli, with a lower likelihood of responding throughout the tracking session, while other parameters of the reversal response (e.g.duration and speed) were not affected.Re-expression of UNC-86 in early adulthood did not rescue these impairments (Fig. 3C) and exposing animals to Auxin from L4 onwards did not induce the hyporesponsive phenotype seen with continuous degradation (Fig. 4B).Together, these findings suggest a primarily early developmental role for unc-86 in regulating mechanosensory responses.While we found that impairing unc-86 from L4 onwards did not induce any significant impairments, a previous study found that inactivating unc-86 at L4 using a temperature-sensitive allele resulted in impaired chemotaxis to multiple odorants (Sze and Ruvkun, 2003).These differences may suggest that unc-86 has distinct temporal functional windows in different neuronal classes (i.e.

Phenotypic change
continuous degradation (blue), induced by adult-specific degradation (starting at L4, yellow), and rescued with adult re-expression (gray).E) Exposing animals to Auxin starting at L2 (yellow) induces impairments in response probability.F) Ending Auxin exposure at L2 (24 hrs post-hatch, gray) enabled phenotypic rescue.
We next investigated the phenotypic consequences of degrading and re-expressing UNC-86 starting at an earlier developmental time point (specifically 24 hours after age-synchronization during late post-embryonic development; beginning after L2).Restoration of unc-86 expression beginning at this earlier time point was sufficient to completely rescue impairments in mechanosensory hyporesponsivity (Fig. 4E).In contrast, we found that degrading UNC-86 from this same time point onwards also impaired response probability (Fig. 4F).These results suggest the window of phenotypic reversibility for unc-86 extends into early post-embryonic development.
Moreover, these results indicate that as long as unc-86 has played its role in neurodevelopment in this early critical window, it is no longer required for normal mechanosensory responding in adulthood.

Ubiquitous degradation of the essential protein DYNC1H1•dhc-1 in adult animals reveals
specific roles in mechanosensory responding and habituation DYNC1H1•dhc-1 is an essential motor protein implicated in diverse processes including cell division and cargo transport along microtubules (e.g.retrograde axonal transport in neurons) (Cianfrocco et al., 2015) (Fig. 5A).Variants in DYNC1H1 have been implicated in several neurodevelopmental disorders including ID and ASD (Satterstrom et al., 2020;Willemsen et al., 2013).Determining the biological functions of DYNC1H1•dhc-1 throughout development has been challenging as dynein loss-of-function results in early embryonic lethality in multiple model organisms (Hamill et al., 2002;Harada et al., 1998;Howell and Rose, 1990;Mains et al., 1990;Robinson et al., 1999).As a result, the role of DYNC1H1•dhc-1 in behaviour remains relatively uncharacterized.
Here, we used the AID system to ubiquitously degrade DHC-1.As expected from loss-offunction alleles, continuous degradation of dynein was lethal.Degrading dynein in early adulthood (starting Auxin exposure immediately after L4) also resulted in lethality, indicating that dynein function remains essential throughout the late stages of C. elegans development.To determine whether dynein function is essential in adulthood, we ubiquitously degraded dynein in 3-day old adults (beginning 72 hrs post-hatch).We found that degrading dynein in adulthood was not lethal, allowing us an opportunity to investigate the biological functions of dynein in adult animals.
Despite its broad expression and essential role in early development, degradation of dynein later in life revealed surprisingly specific roles in adult behaviour (Fig. 5A-F).Adult-specific degradation of dynein did not cause severe alterations in morphology or baseline locomotion.
Instead, only select components of the mechanosensory reversal response were affected while others were left intact.Response probability was unaffected, as the proportion of worms that responded to each mechanosensory stimulus was similar to the no Auxin control (Fig. 5B).
However, degrading dynein at 72 hrs post-hatch caused animals to display deeper habituation of response duration (Fig. 5C) and a slower response speed compared to animals not exposed to Auxin (Fig. 5D).Interestingly, the absolute speed trace across the entire experiment shows animals have the ability to respond as fast as control animals but decrement their response speed faster with repeated stimulation (Fig. 5E).These findings suggest that, in adulthood, dynein may function to promote normal habituation of response duration and speed, but not response probability.Taken together, these results support the hypothesis that different components of To further investigate the role of dynein across development, we took advantage of the ability to activate the AID system cell specifically and obtained a line of C. elegans that allowed for specific and reversible degradation of dynein only in neurons by driving TIR1 expression under the rab-3 promoter (Fig. 6A).Continuous pan-neuronal degradation of dynein did not cause lethality, offering us an unprecedented opportunity to determine the phenotypic consequences of decreased dynein function in the nervous system throughout development and whether the resulting impairments were reversible.Starting degradation of neuronal dynein at L4 (yellow) did not induce significant impairments in response duration but did induce impairments in response speed.E) Full phenotypic profile of dhc-1, indicating all phenotypes induced by continuous degradation (blue), induced by adultspecific degradation (starting at L4, yellow), and rescued with adult re-expression (gray).F) Exposing animals to Auxin at L2 (24 hours post-hatch) impaired response probability (yellow) and ending Auxin exposure at L2 rescued impairments in response probability (gray).G and H) Impairments in response duration and speed were rescued with re-expression of dynein starting at L2 (gray).Starting Auxin exposure at L2 (yellow) did not affect response duration, but did impair response speed.
Continuous pan-neuronal degradation of dynein caused a broad range of sensory responding impairments, including a low probability of responding to mechanosensory stimuli, rapid habituation of response duration, and lower response speed (Fig. 6B-E and Fig. 7A-D).In this case, multiple phenotypes showed distinct temporal functional windows and reversibility patterns.Re-expression of dynein in early adulthood (ending Auxin exposure immediately after L4) did not restore normal mechanosensory responding, but instead resulted in severely impaired habituation of response probability (animals were hyperresponsive and did not learn to decrease their likelihood of responding to repeated stimuli) (Fig. 6B).Adult pan-neuronal degradation of dynein (beginning Auxin exposure at L4) did not alter initial response probability but did decrease habituation of response probability (Fig. 6B).For response duration, adult re-expression of dynein was sufficient to fully rescue the impairment seen with continuous degradation, but adult degradation did not induce the response duration impairment (Fig. 6C).Consistent with our findings for ubiquitous dynein degradation in 72 hr old adult animals, continuous degradation of pan neuronal dynein also caused animals to exhibit slower reversal speed.The speed impairment could be rescued with dynein re-expression in early adulthood and was also induced by adult degradation (Fig. 6D) suggesting that dynein is continuously required in neurons to mediate response speed.probability while starting degradation at L2 produced a similar level of impairment in response probability as the continuous degradation condition (Fig. 6F).Importantly, the lack of the novel hyperresponsive reversal probability phenotype that occurred when dynein was re-expressed at L4 suggests that for certain genes there will be crucial windows in development when reexpression must occur to avoid inducing alternative impairments.For reversal duration and speed, impairments in both phenotypes were fully rescued with re-expression starting at L2 (Fig.

6G-H).
Earlier degradation did not induce impairments in reversal duration, suggesting the crucial functional window of DHC-1 for reversal duration is occurring prior to L2 but is not required throughout the lifespan (Fig. 6G).Earlier degradation did lower reversal speed (Fig. 6H), providing more evidence that DHC-1 is continuously required for this phenotype.Together, these results reveal many new roles for dynein in the developing and adult nervous system and illustrate the diversity of temporal functional windows that can be observed for a single gene (Fig. 6E).

Comparison of temporal profiles reveals shared phenotypic disruptions and prioritizing principles of phenotypic reversibility
All neurodevelopmental disorder risk genes assessed here showed post-embryonic reversibility for at least one impaired phenotype, with several showing phenotypic reversibility relatively late in development (e.g. after stopping Auxin exposure at L4; Fig. 8A).Looking across all genes assessed, at least one phenotype within each phenotypic class (morphology, locomotion, mechanosensory responding, and learning) could be reversed later in life, suggesting a degree of flexibility in when multiple cellular processes can occur during development (Fig. 8A).
However, not all phenotypes caused by the inactivation of a single gene could be reversed, even if reversibility was possible for other affected phenotypes (Fig. 8A).These results suggest that multiple phenotypic disruptions stemming from a single affected gene can show distinct windows of reversibility (Fig. 8A).In addition, there were also cases where the same organism-level phenotypic disruption (e.g.impaired mechanosensory responding) could be rescued later in development by re-expression of one of the genes, but not others (Fig. 8A-C).Based on our results, EBF3•unc-3 is a promising candidate for further research as many impaired phenotypes could be reversed even when gene function was restored later stages of development (Fig. 8A), and re-expression did not induce novel phenotypic disruptions (Fig. 2 & 3).Further, identifying which phenotypes are most commonly affected across risk genes may provide insight into points of convergence in the mechanisms that are altered in neurodevelopmental disorders (Fig. 8B).
We found that reversal probability metrics appeared to be both the most affected and reversible phenotypes across all genes assessed (Fig. 8B & C).This finding supports previous work from multiple model organisms which found that inactivating neurodevelopmental disorder risk genes commonly impaired habituation of response probability (Fenckova et al., 2019;McDiarmid et al., 2020) and that impairments in habituation caused by inactivation of the ASD risk gene ortholog neuroligin (NLGN1/2/3/X•nlg-1) could be partially rescued with adult re-expression (McDiarmid et al., 2020).

Discussion
We systematically investigated the effect of degrading and re-expressing multiple neurodevelopmental disorder risk gene orthologs across a suite of morphological, locomotor, sensory, and learning phenotypes in thousands of freely behaving animals using our highthroughput machine vision tracking system.Taking advantage of the CRISPR-Cas9 AID system allowed us to test whether restoring protein levels through re-expression from the endogenous locus was sufficient for phenotypic rescue at multiple time points throughout development.We found that each gene displayed unique temporal functional windows and phenotypic profiles (Fig. processing it is no longer required in adulthood for normal short-term habituation (McDiarmid et al., 2020).Together, these results reveal a remarkable diversity in temporal phenotypic profiles across neurodevelopmental disorder risk genes that would be missed by approaches that focus on a single phenotype or developmental timepoint.The approach established in this study can be used to systematically assess the temporal requirements and phenotypic reversibility of neurodevelopmental disorder risk genes at an unprecedented throughput to prioritize risk genes for further assessment.

8). DYNCH1•dhc
Using neuron-specific reversible protein degradation, we provide the first description of the role of dynein in behaviour across development.We found that continuous degradation of dynein in only neurons affected the majority (22/30) of the phenotypes assessed, including multiple morphology phenotypes.Interestingly, we found that both degrading or re-expressing dynein in neurons during early adulthood impaired habituation of mechanosensory response probability.Previous work from our lab has found habituation of response probability is affected by developmental stage such that habituation becomes deeper with age due to circuit rewiring and reduced sensitivity (Beck and Rankin, 1993;Rai and Rankin, 2007;Timbers et al., 2013).
The impaired habituation phenotypes seen with both development and adult-specific dynein degradation could both stem from an immature nervous system, such that impairing protein function during early development temporarily impedes developmental processes from occurring whereas early adult inactivation freezes the nervous system in an immature state.Determining whether this neurodevelopmental freeze mechanism, or alternative, more complex mechanisms (e.g.changes in synaptic physiology) mediate these habituation impairments requires further study.Taken together these results reveal that dynein has several ongoing functions in the nervous system to modulate sensory and learning behaviors after it's essential period in development.
Using high-throughput model systems to rapidly assess the temporal function and phenotypic reversibility of neurodevelopmental risk genes provides critical insight into the emerging principles that should be considered in future re-expression studies.Importantly, we found restoring protein expression in adulthood may induce novel phenotypes that were not observed when proteins were continuously degraded (e.g.altered habituation with neuronspecific re-expression of dynein).This finding reveals the importance of assessing a large number of morphological and behavioral phenotypes to ensure novel adverse phenotypes do not arise when gene function is restored later in development.These studies also offer a reminder to the diversity of functions of a single gene across development.Although a gene may play an important role in a phenotype under study, it also may have other functions that are not immediately obvious.
Only by studying multiple phenotypes over the span of development can we begin to understand the breadth of what a given gene contributes to the organism.Future re-expression studies should not only aim to assess whether restoring gene function rescues the well-documented cell functions of the gene of interest, but also capture multiple organism-level phenotypes such as forms of learning or other behaviours commonly altered in neurodevelopmental disorders.Overall, our results suggest that earlier protein re-expression will almost always be better, yet identification of genes with longer reversibility windows and multiple reversible phenotypes in high-throughput model systems should be a key principle in prioritizing candidates for further study.
In addition, we found that time windows for when a gene is required and when reexpression can reverse impairments did not always align.While EBF3•unc-3 showed reciprocal functional and reversibility windows for many affected phenotypes, we found other genes (e.g. nlg-1 and dhc-1) where certain phenotypic impairments were not induced if the protein was degraded in early adulthood, but phenotypic rescue was possible if protein levels were restored at that same time point in development.In addition, while EBF3•unc-3 and BRN3A•unc-86 have relatively similar functions and previously described temporal functional windows, we identified stark differences in their reversibility profiles.Together, these findings highlight the need for systematic assessment of the reversibility windows of neurodevelopmental disorder risk genes across different developmental timepoints, even if there is prior indication of when the gene normally functions.
As the number of genes assessed increases, we may uncover patterns in the molecular attributes that enable certain genes to rescue impairments more broadly than others.For example, EBF3•unc-3 is a putative pioneer transcription factor that may be able to more efficiently remodel chromatin and rewire transcriptional networks outside of its typical developmental window compared to other transcription factors.The approach developed in this study can be adapted to determine how gene reactivation reverses more complex behaviours in C. elegans and more conserved model systems.Information gained from high-throughput model organisms is increasingly valuable as they enable rapid assessment of the growing list of risk genes to gain insight into principles governing neurodevelopment and how a nervous system adapts to the reintroduction of a previously inactive protein.

Animal maintenance
Prior to Auxin experiments, all strains were maintained on Petri plates containing Nematode Growth Medium (NGM) that were seeded with Escherichia coli (E.coli) strain OP50 following standard experimental procedures (Brenner, 1974).96hr post-hatch hermaphrodite animals were used for all experiments.

Auxin inducible degradation strain selection and ortholog identification
All human orthologs of all Auxin-inducible degradation strains available at the CGC were identified using the Alliance of Genome Resources ortholog prediction tool and OrthoList 2 (Agapite et al., 2020;Kim et al., 2018;McDiarmid et al., 2020).Auxin-inducible degradation strains for which the human ortholog corresponded to a known neurodevelopmental disorder risk gene (based on lists generated by recent large-scale sequencing studies and manual literature search (Belmadani et al., 2019;De Rubeis et al., 2014;McDiarmid et al., 2020;Satterstrom et al., 2020)) were selected for analysis.Note that throughout the manuscript the "•" symbol is used to denote the relationship between the human gene and C. elegans ortholog under study (e.g.

Auxin plate preparation
Auxin administration was performed by transferring animals to bacteria-seeded NGM plates containing Auxin (McDiarmid et al., 2020;Zhang et al., 2015).To prepare Auxin plates, a 400 mM stock solution of Auxin indole-3-acetic acid (IAA) (Thermo Fisher, Alfa Aesar™ #A1055614) was created by dissolving Auxin in ethanol.Molten NGM was prepared and allowed to cool to approximately 50°C.The Auxin stock was then diluted into separate flasks of molten NGM agar to final concentrations of 0.025mM, 1mM, and 4mM.The NGM agar + Auxin mixture was then poured into Petri plates and allowed to dry in the dark for 72 hrs.Auxin plates were then seeded with 50 µl of E. coli OP50 liquid culture 48 hrs before use.All plates were stored in the dark at room temperature (20°C) in a temperature and humidity-controlled room (McDiarmid et al., 2020;Zhang et al., 2015).

Population age-synchronization and Auxin administration
Age synchronization by egg lay was used to create the experimental groups for phenotypic analysis as previously described (McDiarmid et al., 2020;McDiarmid et al., 2020).For age synchronization, five gravid adults were placed on to either NGM or Auxin plates and allowed to lay eggs for 4 hours before removal (resulting in 50-100 animals per plate).For the development specific-degradation conditions, approximately 240 progeny were manually transferred from auxin plates onto 6 regular NGM plates (~40 animals per plate) either 24 hrs (at L2), 48 hrs (at L4), or 72 hrs (at early adulthood) after synchronization (egg lay).For adult-specific degradation conditions, approximately 240 progeny were manually transferred from regular NGM plates onto 6 Auxin plates (~40 animals per plate) 24 hrs (at L2), 48 hrs (at L4), or 72 hrs (at early adulthood) after synchronization.All plates remained in the dark other than when animals were being manually transferred to preserve the integrity of Auxin.4-6 plates were run for each experimental condition.

Behavioral paradigm and Multi-Worm Tracker phenotypic analysis
The Multi-Worm Tracker (MWT) was used for all behavioural tracking experiments (Swierczek et al., 2011).Each plate was subjected to the same short-term habituation behavioural paradigm (see results and Fig. 1D) that began with a 5 min period to allow worms to acclimate to being placed on the MWT.After acclimation, we collected data for an additional 5 min period to assess baseline locomotion and morphology features (Fig. 1D).Following this baseline period, thirty mechanosensory stimuli were administered to the side of the Petri plate using an automated push-solenoid at a 10 second inter-stimulus interval (Fig. 1D).These non-localized mechanosensory stimuli cause animals to perform a reversal response, where animals briefly crawl backwards before resuming forward locomotion (Rankin et al., 1990).We quantified multiple mechanosensory sensitivity and habituation learning phenotypes from these reversals which we have previously shown are mediated by genetically dissociable underlying mechanisms.After the 30 th stimulus, a 5 min rest period occurred which was followed by the administration of a final stimulus to assess short-term memory retention of habituation (spontaneous recovery, Fig. 1D).
See The Multi-Worm Tracker user guide (https://sourceforge.net/projects/mwt/) or McDiarmid et al. 2020(McDiarmid et al., 2020) for full description of all phenotypes.All testing occurred in a temperature and humidity-controlled room at approximately 20°C.
We used the MWT software (version 1.2.0.2) to delivery stimuli and acquire images (Swierczek et al., 2011), and Choreography software (version 1.3.0_r103552) to quantify phenotypes.Choreography filters "-shadowless", "-minimum-move-body 2", and "-minimumtime 20" were used to restrict analysis to animals that moved more than 2 body lengths and were tracked for 20 secs or longer.The "MeasureReversal" plug-in was used to identify animals that reversed within 1 sec of the mechanosensory stimulus being administered (Swierczek et al., 2011).Choreography output files were organized using custom R scripts which are freely available at [Github link will be inserted here].All phenotypic features were pooled across the 4-6 plate replicates (each plate replicate capturing 40-100 animals) per strain.The mean of each condition was then compared using an unpaired t-test and Benjamini-Hochberg control of false discovery rate at 0.01.Figures were generated using the ggplot2 package in R (Wickham, 2016).

Strains used
The following strains are available through the Caenorhabditis Genetics Center (CGC):

Figure 1 .
Figure 1.A pipeline to assess the temporal requirements and phenotypic reversibility of neurodevelopmental disorder risk gene orthologs using the AID system and machine vision phenotypic profiling in Caenorhabditis elegans.A) The Auxin-Inducible Degradation (AID) system is a powerful approach that enables temporal and spatial control of protein depletion.CRISPR-Cas9 is used to tag the gene of interest (GOI) with the AID degron along with a fluorescent protein (FP) to visualize protein expression in vivo.In the presence of the small molecule Auxin, TIR1 (an E3 ubiquitin ligase) associates with the AID degron, recruiting endogenous proteosomes to degrade the ubiquitinated protein of interest.B) Temporal degradation conditions were created by manually transferring animals on and off Petri plates containing Auxin to inactive or restore gene function at specific timepoints in development or adulthood.C) The effects of protein degradation and re-expression across 26 morphological, locomotor, and sensory and learning phenotypes were objectively quantified in hundreds of animals simultaneously using a machine vision tracking system throughout a short-term mechanosensory habituation paradigm (D).

Figure 2 .
Figure 2. The transcription factor EBF3•unc-3 displays a reciprocal pattern of phenotypic induction and reversibility across development.A) The transcription factor EBF3•unc-3 acts to specify neuronal identity.B) Continuous degradation unc-3 (blue) impaired the animals' ability to respond to mechanosensory stimuli compared to the no-Auxin control (black).Starting Auxin exposure at L4 partially induced this impairment.C) Ending Auxin exposure after L4 (48 hrs posthatch) partially rescued the impairment in response probability.D) Full phenotypic profile of unc-3, indicating all phenotypes induced by continuous degradation (blue), induced by adult-specific degradation (starting at L4, yellow), and rescued with adult re-expression (gray).E) Ending Auxin exposure at L2 (24 hrs post-hatch, gray) almost completely rescued the impairment, F-G) with the level of phenotypic rescue decreasing with later onset of UNC-3 re-expression.E) Starting Auxin exposure at L2 (yellow) induced an impairment level similar to the continuous degradation group, F-G) with the degree of phenotypic impairment decreasing with later onset of UNC-3 degradation.

Figure 3 .F
Figure 3. EBF3•unc-3 shows diverse temporal patterns of phenotypic reversibility across morphological, locomotor, and mechanosensory response phenotypes.A) The no-Auxin control group is depicted in black and continuous degradation group is depicted in blue for all panels.Altered animal length could be partially rescued with early post embryonic re-expression (starting at L2/24 hrs post-hatch) or fully induced with early post embryonic degradation.B) The degree of rescue and impairment of animal length increasingly diminished if UNC-3 was reexpression or degraded at L4 (48 hrs post-hatch) C) or in adulthood (72 hrs post-hatch).D)Similarly, impairments in kink could be fully rescued with early post embryonic re-expression, E-F) but degree of rescue diminished with later re-expression.D) Degrading UNC-3 starting at L2 resulted in impaired kinkiness to a level similar to the continuous degradation control group, E-F) the degree of impairment lessened with later onset of degradation.G-I) Impairments in response unc-86 is continuously required for the function of chemotaxis neurons but is only required in early development for function of mechanosensory neurons).

Figure 4 .
Figure 4. BRN3A•unc-86 specifically impairs mechanosensory response probability and displays a reversibility window restricted to early post-embryonic development.A) The transcription factor BRNA3•unc-86 acts maintain the expression of terminal identity genes in multiple neuron types.B) Continuous degradation of UNC-86 (blue) specifically impaired response probability to mechanosensory stimuli compared to animals that were not exposed to Auxin (black).Staring Auxin exposure at L4 (48 hrs post-hatch, yellow) did not significantly induce phenotypic impairments.C) Ending Auxin exposure at L4 (gray) did not rescue impairments in response probability.D) Full phenotypic profile of unc-86, indicating all phenotypes induced by

Figure 5 .
Figure 5. Ubiquitous degradation of the essential protein DYNC1H1•dhc-1 in adult animals reveals specific roles in mechanosensory responding and habituation.A) The essential gene DYNC1H1•dhc-1 acts in cargo transport and stabilization of microtubule dynamics.B) Starting Auxin exposure in early adulthood (72 hrs post-hatch) did not impair response probability but C) did deepen habituation of response duration and D) decreased response speed compared to the no-Auxin control animals (black).E) Degrading dynein in adult animals (yellow) decreased average speed during the acclimation period and caused deeper habituation of response speed across the mechanosensory stimuli resulting in a lower average speed during the rest period post mechanosensory stimulation.F) Full phenotypic profiles of dhc-1, indicating all phenotypes induced by adult-specific degradation starting at 3 days post synchronization (yellow).

Figure 6 .
Figure 6.Pan-neuronal degradation of DYNCH1•dhc-1 is not lethal and causes multiple habituation impairments with distinct reversibility profiles.A) Pan neuronal degradation of dynein was achieved by crossing the dhc-1(ie28[dhc-1::degron::GFP]) strain with a strain where TIR1 expression is driven by the rab-3 promoter.B) Continuous degradation of neuronal dynein (blue) impaired response probability compared to the no-Auxin control group (black).Reexpressing (gray) or degrading (yellow) neuronal dynein starting at L4 (48 hrs post-hatch) caused

Figure 7 .
Figure 7. Pan-neuronally degrading and re-expressing dynein reveals distinct temporal functional windows for morphological and baseline locomotion features.The no-Auxin control group is depicted in black and continuous degradation group is depicted in blue for all panels.(A-C) Dynein is continuously required in neurons for normal kinked body posture, aspectratio, and movement direction bias.Degrading DHC-1 in neurons beginning at L4 (48 hrs posthatch) induced impairment levels across these phenotypes similar to the continuous degradation control.Pan neuronal re-expression of DHC-1 at L4 (grey) rescued impairments in all three phenotypes.D) Re-expression of DHC-1 at L4 rescued speed before the onset mechanosensorystimuli, however only partially rescued speed deficits in the 5 min rest period.Degrading dynein at L4 (yellow) caused animals to exhibit lower speed throughout the behavioural paradigm that was similar to the continuous degradation control.E) Dynein is continuously required in neurons for animal length, but re-expression at L4 can partially rescue impairments F) Novel impairments in animal curvature occurred when dynein was re-expressed at L4. Beginning protein degradation at L4 caused animals to have a higher body curvature than the no-Auxin control, whereas reexpressing pan neuronal DHC-1 at L4 caused animals to exhibit a lower body curvature than controls.

Figure 8 .
Figure 8.Comparison of temporal profiles reveals shared phenotypic disruptions and prioritizing principles of phenotypic reversibility.A) Number and kind of phenotypes that could be induced by continuously degrading each gene or reversed by re-expressing the gene 24h (L2 stage) or 48h (L4 stage) after synchronization (i.e.early or late in post-embryonic development).B) Heatmap showing the phenotypes affected by continuous degradation of each gene.Phenotypes observed in all 3 genes are highlighted in red.C) Heatmap showing thephenotypes that could be rescued with protein re-expression starting at L2 (24 hours post synchronization).Phenotypes that were reversible across all 3 genes are highlighted in green.
-1 function is continuously essential for development and plays a specific role in adult mechanosensory behavior.The transcription factor BRN3A•unc-86 can only reverse phenotypic disruptions early in post-embryonic development and is not required for adult function, suggesting a primarily developmental role in mechanosensory responding.The transcription factor EBF3•unc-3 displayed a range of temporal requirements throughout the lifespan and can reverse multiple phenotypic disruptions later in life.In addition to the 3 genes tested in this study, previous findings from our lab provide the temporal requirements for another neurodevelopmental disorder risk gene, the synaptic cell adhesion molecule NLGN1/2/3/4/X•nlg-1. Adulthood reexpression of NLGN1/2/3/4/X•nlg-1 can partially reverse impairments in habituation of response probability, however, once neuroligin has functioned to build a circuit capable of normal sensory