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.

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 behavioral 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., 2019; 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 behavioral 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 genes 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 is 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. There are C. elegans orthologs for a high number of neurodevelopmental risk genes [e.g. >80% of high-confidence ASD risk genes (McDiarmid et al., 2020a,b)], 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., 2018a,b; Post et al., 2020). C. elegans has 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 enables 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økjær-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 behaviors in real time (Husson et al., 2012; McDiarmid et al., 2018a,b; Swierczek et al., 2011).

A behavior 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 (Blok et al., 2022; 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 behavioral 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 protein products of three genes – EBF3•unc-3, BRN3A•unc-86 and DYNC1H1•dhc-1 (‘•’ denotes orthologous relationship between human and C. elegans genes) – at multiple developmental time points in hundreds 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 if re-expression began 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.

Fig. 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 indole-3-acetic acid (IAA), 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 IAA to inactivate or restore gene function at specific time points in development or adulthood. (C,D) The effects of protein degradation and re-expression across 30 morphological, locomotor, and sensory and learning phenotypes were objectively quantified in hundreds of animals simultaneously using a machine vision tracking system (C) throughout a short-term mechanosensory habituation paradigm (D).

Fig. 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 indole-3-acetic acid (IAA), 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 IAA to inactivate or restore gene function at specific time points in development or adulthood. (C,D) The effects of protein degradation and re-expression across 30 morphological, locomotor, and sensory and learning phenotypes were objectively quantified in hundreds of animals simultaneously using a machine vision tracking system (C) throughout a short-term mechanosensory habituation paradigm (D).

We selected three neurodevelopmental disorder risk genes [as identified by Simons Foundation Autism Research Initiative (Abrahams et al., 2013); Satterstorm et al. (2020)] for reversibility analysis based on the availability of C. elegans strains that contain a neurodevelopmental disorder risk gene ortholog tagged with an auxin-inducible degron at the endogenous locus using CRISPR-Cas9 (see Materials and 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 auxin class plant hormone indole-3-acetic acid (IAA), 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 IAA (McDiarmid et al., 2020a,b; Zhang et al., 2015). Since the AID degron is tagged to the endogenous locus, protein expression is restored using the native regulatory machinery, therefore bypassing the biological confounds associated with conventional approaches that rely on overexpression. Importantly, our laboratory and others have shown that IAA exposure does not cause any overt effects on C. elegans morphology, locomotion, short-term learning or mechanosensory processing phenotypes (McDiarmid et al., 2020a,b; 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 IAA 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 30 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, 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 re-expression 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 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., 2008, 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).

Fig. 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 of UNC-3 (blue) impaired the animals' ability to respond to mechanosensory stimuli compared to the no IAA control (black). Starting IAA exposure at L4 partially induced this impairment. (C) Ending IAA exposure after L4 (48 h post-hatch) partially rescued the impairment in response probability. (B,C) Data shown as mean±s.e.m. using plates as n (n=4 plates each with 40-100 worms per condition). Significant differences between individual components of the response curves across conditions are indicated in D. (D) Full phenotypic profile of UNC-3, indicating all phenotypes induced by continuous degradation (row 1), induced by adult-specific degradation (row 2), or observed following developmental degradation and adult re-expression (row 3). Cells represent directional t-statistics from comparisons to wild-type controls. Directional t-statistics are shown such that phenotypes higher than the no IAA control are progressively more yellow, and phenotypes lower than the no IAA control are progressively more blue. Only significant differences [false discovery rate (FDR)<0.01] are shown and t-statistics are capped at ±10 to aid visualization. The fourth row indicates the degree of phenotypic rescue (difference from continuous degradation vs adult re-expression), and the fifth row indicates whether that difference was significant (FDR<0.01). For the fourth row, progressively darker green indicates a higher amount of rescue and white indicates no rescue occurred or there was no effect observed with continuous degradation to rescue. See Materials and Methods and Table S1 for a complete description of all 30 quantitative phenotypic features. (E-G) Ending IAA exposure at L2 (24 h post-hatch, gray) almost completely rescued the impairment (E), with the level of phenotypic rescue decreasing with later onset of UNC-3 re-expression (F,G). Starting IAA exposure at L2 (yellow) induced an impairment level similar to the continuous degradation group (E), with the degree of phenotypic impairment decreasing with later onset of UNC-3 degradation (F,G). (E-G) Data shown as mean±s.e.m using plates as n (n=4 plates per condition each with 40-100 worms).

Fig. 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 of UNC-3 (blue) impaired the animals' ability to respond to mechanosensory stimuli compared to the no IAA control (black). Starting IAA exposure at L4 partially induced this impairment. (C) Ending IAA exposure after L4 (48 h post-hatch) partially rescued the impairment in response probability. (B,C) Data shown as mean±s.e.m. using plates as n (n=4 plates each with 40-100 worms per condition). Significant differences between individual components of the response curves across conditions are indicated in D. (D) Full phenotypic profile of UNC-3, indicating all phenotypes induced by continuous degradation (row 1), induced by adult-specific degradation (row 2), or observed following developmental degradation and adult re-expression (row 3). Cells represent directional t-statistics from comparisons to wild-type controls. Directional t-statistics are shown such that phenotypes higher than the no IAA control are progressively more yellow, and phenotypes lower than the no IAA control are progressively more blue. Only significant differences [false discovery rate (FDR)<0.01] are shown and t-statistics are capped at ±10 to aid visualization. The fourth row indicates the degree of phenotypic rescue (difference from continuous degradation vs adult re-expression), and the fifth row indicates whether that difference was significant (FDR<0.01). For the fourth row, progressively darker green indicates a higher amount of rescue and white indicates no rescue occurred or there was no effect observed with continuous degradation to rescue. See Materials and Methods and Table S1 for a complete description of all 30 quantitative phenotypic features. (E-G) Ending IAA exposure at L2 (24 h post-hatch, gray) almost completely rescued the impairment (E), with the level of phenotypic rescue decreasing with later onset of UNC-3 re-expression (F,G). Starting IAA exposure at L2 (yellow) induced an impairment level similar to the continuous degradation group (E), with the degree of phenotypic impairment decreasing with later onset of UNC-3 degradation (F,G). (E-G) Data shown as mean±s.e.m using plates as n (n=4 plates per condition each with 40-100 worms).

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 IAA exposure at the final larval stage, L4) resulted in effective protein degradation (Fig. S1) and 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 IAA immediately after L4) (Fig. 2C).

Fig. 3.

EBF3•unc-3 shows diverse temporal patterns of phenotypic reversibility across morphological, locomotor and mechanosensory response phenotypes. The no IAA control group is depicted in black and continuous degradation group is depicted in blue for all panels. (A) Altered animal length could be partially rescued with early post-embryonic re-expression (starting at L2/24 h post-hatch) or fully induced with early post-embryonic degradation. (B,C) The degree of rescue and impairment of animal length increasingly diminished if UNC-3 was re-expression or degraded at L4 (48 h post-hatch) (B) or in adulthood (72 h post-hatch) (C). (D-F) Similarly, impairments in kink could be fully rescued with early post embryonic re-expression (D), but degree of rescue diminished with later re-expression (E,F). Degrading UNC-3 starting at L2 resulted in altered kink to a level similar to the continuous degradation control group (D); the degree of impairment lessened with later onset of degradation (E,F). (A-F) Data shown as mean±s.e.m. using plates as n (n=4 plates per condition each with 40-100 worms). Small points represent individual plate replicates and large points represent the mean±s.e.m. of plate replicates. n.s., not significant. (G-I) Impairments in response duration could not be rescued with UNC-3 re-expression across any of the tested temporal conditions. Degrading UNC-3 in early post embryonic development or development (L4) induced impairments similar to the continuous degradation condition (G,H), yet duration impairments were not strongly induced with UNC-3 degradation starting at 72 h post-hatch (I). (G-I) Data shown as mean±s.e.m. using plates as n (n=4 plates per condition each with 40-100 worms).

Fig. 3.

EBF3•unc-3 shows diverse temporal patterns of phenotypic reversibility across morphological, locomotor and mechanosensory response phenotypes. The no IAA control group is depicted in black and continuous degradation group is depicted in blue for all panels. (A) Altered animal length could be partially rescued with early post-embryonic re-expression (starting at L2/24 h post-hatch) or fully induced with early post-embryonic degradation. (B,C) The degree of rescue and impairment of animal length increasingly diminished if UNC-3 was re-expression or degraded at L4 (48 h post-hatch) (B) or in adulthood (72 h post-hatch) (C). (D-F) Similarly, impairments in kink could be fully rescued with early post embryonic re-expression (D), but degree of rescue diminished with later re-expression (E,F). Degrading UNC-3 starting at L2 resulted in altered kink to a level similar to the continuous degradation control group (D); the degree of impairment lessened with later onset of degradation (E,F). (A-F) Data shown as mean±s.e.m. using plates as n (n=4 plates per condition each with 40-100 worms). Small points represent individual plate replicates and large points represent the mean±s.e.m. of plate replicates. n.s., not significant. (G-I) Impairments in response duration could not be rescued with UNC-3 re-expression across any of the tested temporal conditions. Degrading UNC-3 in early post embryonic development or development (L4) induced impairments similar to the continuous degradation condition (G,H), yet duration impairments were not strongly induced with UNC-3 degradation starting at 72 h post-hatch (I). (G-I) Data shown as mean±s.e.m. using plates as n (n=4 plates per condition each with 40-100 worms).

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 IAA 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 those present with continuous degradation, 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 h post-hatch) from IAA. 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 h 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 behavior 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.

Degradation of BRN3A•UNC-86 specifically impairs mechanosensory response probability and displays a reversibility window restricted to early post-embryonic development

BRN3A•unc-86 (also known as POU4F1) 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 BRN3A have been associated with abnormal development of sensory structures, including auditory (Huang et al., 2001) and visual cells (Badea et al., 2009). BRN3A may also be affected by a copy number variant identified in a family with ASD that harbor a deletion 153.5 kp upstream of the BRN3A locus (Salyakina et al., 2011). While the gene sequence is not affected, upstream regulatory components of BRN3A may be altered, leading to abnormal development (Salyakina et al., 2011). 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).

Fig. 4.

Degrading 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 to 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 IAA (black). Staring IAA exposure at L4 (48 h post-hatch, yellow) did not significantly induce phenotypic impairments. (C) Ending IAA exposure at L4 (gray) did not rescue impairments in response probability. (B,C) Data shown as mean±s.e.m. using plates as n (n=4 plates each with 40-100 worms per condition). Significant differences between individual components of the response curves across conditions are indicated in D. (D) Full phenotypic profile of UNC-86, indicating all phenotypes induced by continuous degradation (row 1), induced by adult-specific degradation (row 2), or observed following developmental degradation and adult re-expression (row 3). Cells represent directional t-statistics from comparisons to wild-type controls. Directional t-statistics are shown such that phenotypes higher than the no auxin control are progressively more yellow, and phenotypes lower than the no IAA control are progressively more blue. Only significant differences (FDR<0.01) are shown and t-statistics are capped at ±10 to aid visualization. The fourth row indicates the degree of phenotypic rescue (difference from continuous degradation vs adult re-expression), and the fifth row indicates whether that difference was significant (FDR<0.01). For the fourth row, progressively darker green indicates a higher amount of rescue and white indicates no rescue occurred or there was no effect observed with continuous degradation to rescue. See Materials and Methods and Table S1 for a complete description of all 30 quantitative phenotypic features. (E) Exposing animals to IAA starting at L2 (yellow) induces impairments in response probability. (F) Ending IAA exposure at L2 (24 h post-hatch, gray) enabled phenotypic rescue. (E,F) Data shown as mean±s.e.m. using plates as n (n=4-6 plates per condition each with 40-100 worms).

Fig. 4.

Degrading 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 to 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 IAA (black). Staring IAA exposure at L4 (48 h post-hatch, yellow) did not significantly induce phenotypic impairments. (C) Ending IAA exposure at L4 (gray) did not rescue impairments in response probability. (B,C) Data shown as mean±s.e.m. using plates as n (n=4 plates each with 40-100 worms per condition). Significant differences between individual components of the response curves across conditions are indicated in D. (D) Full phenotypic profile of UNC-86, indicating all phenotypes induced by continuous degradation (row 1), induced by adult-specific degradation (row 2), or observed following developmental degradation and adult re-expression (row 3). Cells represent directional t-statistics from comparisons to wild-type controls. Directional t-statistics are shown such that phenotypes higher than the no auxin control are progressively more yellow, and phenotypes lower than the no IAA control are progressively more blue. Only significant differences (FDR<0.01) are shown and t-statistics are capped at ±10 to aid visualization. The fourth row indicates the degree of phenotypic rescue (difference from continuous degradation vs adult re-expression), and the fifth row indicates whether that difference was significant (FDR<0.01). For the fourth row, progressively darker green indicates a higher amount of rescue and white indicates no rescue occurred or there was no effect observed with continuous degradation to rescue. See Materials and Methods and Table S1 for a complete description of all 30 quantitative phenotypic features. (E) Exposing animals to IAA starting at L2 (yellow) induces impairments in response probability. (F) Ending IAA exposure at L2 (24 h post-hatch, gray) enabled phenotypic rescue. (E,F) Data shown as mean±s.e.m. using plates as n (n=4-6 plates per condition each with 40-100 worms).

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. 4C) and exposing animals to IAA from L4 onwards did not induce the hyporesponsive phenotype seen with continuous degradation (Fig. 4B) even though effective degradation occurred (Fig. S2). 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 (Fig. 4D), 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. unc-86 is continuously required for the function of chemotaxis neurons but is only required in early development for function of mechanosensory neurons).

We next investigated the phenotypic consequences of degrading and re-expressing UNC-86 starting at an earlier developmental time point (specifically 24 h 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 behavior remains relatively uncharacterized.

Fig. 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-D) Starting IAA exposure in early adulthood (72 h post-hatch) did not impair response probability (B) but did deepen habituation of response duration (C) and decreased response speed (D) compared to the no IAA 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. (B-E) Data shown as mean±s.e.m. using plates as n (n=4 plates per condition each with 40-100 worms). (F) Full phenotypic profiles of DHC-1, indicating all phenotypes induced by adult-specific degradation starting at 3 days post-synchronization. Directional t-statistics are shown such that phenotypes higher than the no IAA control are progressively more yellow, and phenotypes lower than the no IAA control are progressively more blue. Only significant differences (FDR<0.01) are shown and t-statistics are capped at ±10 to aid visualization.

Fig. 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-D) Starting IAA exposure in early adulthood (72 h post-hatch) did not impair response probability (B) but did deepen habituation of response duration (C) and decreased response speed (D) compared to the no IAA 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. (B-E) Data shown as mean±s.e.m. using plates as n (n=4 plates per condition each with 40-100 worms). (F) Full phenotypic profiles of DHC-1, indicating all phenotypes induced by adult-specific degradation starting at 3 days post-synchronization. Directional t-statistics are shown such that phenotypes higher than the no IAA control are progressively more yellow, and phenotypes lower than the no IAA control are progressively more blue. Only significant differences (FDR<0.01) are shown and t-statistics are capped at ±10 to aid visualization.

Here, we used the AID system to ubiquitously degrade DHC-1. As expected from loss-of-function alleles, continuous degradation of dynein was lethal. Degrading dynein in early adulthood (starting IAA 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 h 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 behavior (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 IAA control (Fig. 5B). However, degrading dynein at 72 h post-hatch caused animals to display deeper habituation of response duration (Fig. 5C) and a slower response speed compared to animals not exposed to IAA (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 habituation can be mediated by distinct mechanisms, and reveal novel, adult-specific roles for DYNC1H1•dhc-1 in mechanosensory responding and habituation (Fig. 5F).

Pan-neuronal degradation of DYNCH1•dhc-1 is not lethal and causes multiple habituation impairments with distinct reversibility profiles

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 even though effective degradation was observed (Fig. S3), 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.

Fig. 6.

Pan-neuronal degradation of DYNC1H1•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 IAA control group (black). Re-expressing (gray) or degrading (yellow) neuronal dynein starting at L4 (48 h post-hatch) caused animals to be hyperresponsive to mechanosensory stimuli. (C,D) Impairments in response duration and response speed could be rescued with re-expression of dynein starting at L4 (gray). Starting degradation of neuronal dynein at L4 (yellow) did not induce significant impairments in response duration but did induce impairments in response speed. (B-D) Data shown as mean±s.e.m. using plates as n (n=4-6 plates each with 40-100 worms per condition). (E) Full phenotypic profile of DHC-1, indicating all phenotypes induced by continuous degradation (row 1), induced by adult-specific degradation (row 2), or observed following developmental degradation and adult re-expression (row 3). Cells represent directional t-statistics from comparisons to wild-type controls. Directional t-statistics are shown such that phenotypes higher than the no IAA control are progressively more yellow, and phenotypes lower than the no IAA control are progressively more blue. Only significant differences (FDR<0.01) are shown and t-statistics are capped at ±10 to aid visualization. The fourth row indicates the degree of phenotypic rescue (difference from continuous degradation vs adult re-expression), and the fifth row indicates whether that difference was significant (FDR<0.01). For the fourth row, progressively darker green indicates a higher amount of rescue and white indicates no rescue occurred or there was no effect observed with continuous degradation to rescue. (F) Exposing animals to IAA at L2 (24 h post-hatch) impaired response probability (yellow) and ending IAA exposure at L2 rescued impairments in response probability (gray). (G,H) Impairments in response duration and speed were rescued with re-expression of dynein starting at L2 (gray). Starting IAA exposure at L2 (yellow) did not affect response duration but did impair response speed. (F-H) Data shown as mean±s.e.m. using plates as n (n=6 plates per condition each with 40-100 worms).

Fig. 6.

Pan-neuronal degradation of DYNC1H1•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 IAA control group (black). Re-expressing (gray) or degrading (yellow) neuronal dynein starting at L4 (48 h post-hatch) caused animals to be hyperresponsive to mechanosensory stimuli. (C,D) Impairments in response duration and response speed could be rescued with re-expression of dynein starting at L4 (gray). Starting degradation of neuronal dynein at L4 (yellow) did not induce significant impairments in response duration but did induce impairments in response speed. (B-D) Data shown as mean±s.e.m. using plates as n (n=4-6 plates each with 40-100 worms per condition). (E) Full phenotypic profile of DHC-1, indicating all phenotypes induced by continuous degradation (row 1), induced by adult-specific degradation (row 2), or observed following developmental degradation and adult re-expression (row 3). Cells represent directional t-statistics from comparisons to wild-type controls. Directional t-statistics are shown such that phenotypes higher than the no IAA control are progressively more yellow, and phenotypes lower than the no IAA control are progressively more blue. Only significant differences (FDR<0.01) are shown and t-statistics are capped at ±10 to aid visualization. The fourth row indicates the degree of phenotypic rescue (difference from continuous degradation vs adult re-expression), and the fifth row indicates whether that difference was significant (FDR<0.01). For the fourth row, progressively darker green indicates a higher amount of rescue and white indicates no rescue occurred or there was no effect observed with continuous degradation to rescue. (F) Exposing animals to IAA at L2 (24 h post-hatch) impaired response probability (yellow) and ending IAA exposure at L2 rescued impairments in response probability (gray). (G,H) Impairments in response duration and speed were rescued with re-expression of dynein starting at L2 (gray). Starting IAA exposure at L2 (yellow) did not affect response duration but did impair response speed. (F-H) Data shown as mean±s.e.m. using plates as n (n=6 plates per condition each with 40-100 worms).

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 IAA 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 IAA 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-h-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. In addition, while continuously degrading dynein did not alter animal curvature, both adult-specific degradation and re-expression of DHC-1 led to alterations in curvature (Fig. 7E).

Fig. 7.

Degrading and re-expressing dynein in neurons reveals distinct temporal functional windows for morphological and baseline locomotion features. The no IAA control group is depicted in black and continuous degradation group is depicted in blue for all panels. (A-D) Dynein is continuously required in neurons for normal aspect ratio (A), kinked body posture (B) animal speed (C) and body length (D). Degrading DHC-1 in neurons beginning at L4 (48 h post-hatch) induced impairment in animal aspect, kink and baseline speed that were similar to the continuous degradation control. Pan-neuronal re-expression of DHC-1 at L4 (gray) rescued impairments in all three phenotypes. (D) Dynein is continuously required in neurons for animal length, but re-expression at L4 can partially rescue impairment. (E) Continuously degrading DHC-1 did not affect animal curvature; however, a novel impairment 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 IAA control and animals continuously exposed to IAA, whereas re-expressing pan-neuronal DHC-1 at L4 caused animals to exhibit a lower body curvature than controls. (A-E) Data shown as mean±s.e.m. using plates as n (n=4-6 plates per condition each with 40-100 worms). Small points represent individual plate replicates and large points represent the mean±s.e.m. of plate replicates. n.s., not significant.

Fig. 7.

Degrading and re-expressing dynein in neurons reveals distinct temporal functional windows for morphological and baseline locomotion features. The no IAA control group is depicted in black and continuous degradation group is depicted in blue for all panels. (A-D) Dynein is continuously required in neurons for normal aspect ratio (A), kinked body posture (B) animal speed (C) and body length (D). Degrading DHC-1 in neurons beginning at L4 (48 h post-hatch) induced impairment in animal aspect, kink and baseline speed that were similar to the continuous degradation control. Pan-neuronal re-expression of DHC-1 at L4 (gray) rescued impairments in all three phenotypes. (D) Dynein is continuously required in neurons for animal length, but re-expression at L4 can partially rescue impairment. (E) Continuously degrading DHC-1 did not affect animal curvature; however, a novel impairment 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 IAA control and animals continuously exposed to IAA, whereas re-expressing pan-neuronal DHC-1 at L4 caused animals to exhibit a lower body curvature than controls. (A-E) Data shown as mean±s.e.m. using plates as n (n=4-6 plates per condition each with 40-100 worms). Small points represent individual plate replicates and large points represent the mean±s.e.m. of plate replicates. n.s., not significant.

We next investigated the phenotypic consequence of degrading and re-expressing dynein in neurons at an earlier time point in development. Re-expressing dynein during early post-embryonic development (starting IAA exposure at L2) fully rescued impairments in response 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 re-expression 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 occurs 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 IAA 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 at later stages of development (Fig. 8A), and re-expression did not induce novel phenotypic disruptions (Figs 2 and 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., 2020a,b) 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., 2020a,b).

Fig. 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 24 h (L2 stage) or 48 h (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. Altered phenotypes observed across all three genes are highlighted in pink. (C) Heatmap showing the phenotypes that could be rescued with protein re-expression starting at L2 (24 h post-synchronization). Phenotypes that were reversible across all three genes are highlighted in green. (B,C) Only significant differences (FDR<0.01) are shown.

Fig. 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 24 h (L2 stage) or 48 h (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. Altered phenotypes observed across all three genes are highlighted in pink. (C) Heatmap showing the phenotypes that could be rescued with protein re-expression starting at L2 (24 h post-synchronization). Phenotypes that were reversible across all three genes are highlighted in green. (B,C) Only significant differences (FDR<0.01) are shown.

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 high-throughput 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. 8). DYNCH1•dhc-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 three genes tested in this study, previous findings from our laboratory provide the temporal requirements for another neurodevelopmental disorder risk gene, the synaptic cell adhesion molecule NLGN1/2/3/4/X•nlg-1. Adulthood re-expression 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 processing, it is no longer required in adulthood for normal short-term habituation (McDiarmid et al., 2020a,b). 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 time point. 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.

Using neuron-specific reversible protein degradation, we provide the first description of the role of dynein in behavior across development. We found that continuous degradation of dynein specifically in neurons affected the majority (23/30) of the phenotypes assessed, including multiple morphology phenotypes. Interestingly, we found that both degrading and re-expressing dynein in neurons during early adulthood impaired habituation of mechanosensory response probability. Previous work from our laboratory has found that 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 its 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 neuron-specific 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 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 behaviors commonly altered in neurodevelopmental disorders. Overall, our results suggest that earlier protein re-expression will almost always enable more phenotypic reversibility, 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 re-expression 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 temporal mRNA expression levels (Fig. S4), we identified stark differences in their reversibility profiles, indicating that temporal mRNA levels are not necessarily predictive of when a gene in needed for a given phenotype. Together, these findings highlight the need for systematic assessment of the reversibility windows of neurodevelopmental disorder risk genes across different developmental time points, 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 the temporal requirements and effects of gene reactivation for cellular phenotypes or more complex behaviors in C. elegans and other model systems. It is important to note that the temporal requirements of gene function and phenotypic reversibility observed in more evolutionarily distant model organisms, such as C. elegans, may not always be conserved in mammals. Nonetheless, information gained from high-throughput model organisms is increasingly valuable as it enables rapid assessment and prioritization of the growing list of risk genes to gain insight into principles governing neurodevelopment and how a nervous system adapts to the re-introduction of a previously inactive protein.

Animal maintenance

Prior to IAA experiments, all strains were maintained on Petri plates containing nematode growth medium (NGM) that were seeded with Escherichia coli strain OP50 following standard experimental procedures (Brenner, 1974). 96 h post-hatch hermaphrodite animals were used for all experiments.

AID strain selection and ortholog identification

All human orthologs of all AID strains available at the Caenorhabditis Genetics Center (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., 2020a,b). AID 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., 2020a,b; Satterstrom et al., 2020)] were selected for analysis. Note that throughout the paper the ‘’ symbol is used to denote the relationship between the human gene and C. elegans ortholog under study (e.g. DYNCH1•dhc-1).

IAA plate preparation

Auxin administration was performed by transferring animals to bacteria-seeded NGM plates containing IAA (McDiarmid et al., 2020a,b; Zhang et al., 2015). To prepare IAA plates, a 400 mM stock solution of IAA (Thermo Fisher, Alfa Aesar™ #A1055614) was created by dissolving IAA in ethanol. Molten NGM was prepared and allowed to cool to approximately 50°C. The IAA stock was then diluted into separate flasks of molten NGM agar to final concentrations of 0.025 mM, 1 mM and 4 mM. The NGM agar+IAA mixture was then poured into Petri plates and allowed to dry in the dark for 72 h. IAA plates were then seeded with 50 µl of E. coli OP50 liquid culture 48 h before use. All plates were stored in the dark at room temperature (20°C) in a temperature and humidity-controlled room (McDiarmid et al., 2020a,b; 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., 2020a,b). For age synchronization, five gravid adults were placed onto either NGM or IAA plates and allowed to lay eggs for 4 h before removal (resulting in 50-100 animals per plate). For the development specific-degradation conditions, approximately 240 progeny were manually transferred from IAA plates onto six regular NGM plates (∼40 animals per plate) either 24 h (at L2), 48 h (at L4) or 72 h (at early adulthood) after synchronization (egg lay). For adult-specific degradation conditions, approximately 240 progeny were manually transferred from regular NGM plates onto six IAA plates (∼40 animals per plate) 24 h (at L2), 48 h (at L4) or 72 h (at early adulthood) after synchronization. All plates remained in the dark other than when animals were being manually transferred to preserve the integrity of IAA. Four to six plates were run for each experimental condition.

Behavioral paradigm and MWT phenotypic analysis

The MWT was used for all behavioral tracking experiments (Swierczek et al., 2011). Each plate was subjected to the same short-term habituation behavioral 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, 30 mechanosensory stimuli were administered to the side of the Petri plate using an automated push-solenoid at a 10 s 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 30th 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 Table S1 for full description of all phenotypes. All testing occurred in a temperature and humidity-controlled room at approximately 20°C (McDiarmid et al., 2020a,b; Swierczek et al., 2011).

We used the MWT software (version 1.2.0.2) to deliver mechanosensory 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 ‘–minimum-time 20’ were used to restrict analysis to animals that moved more than two body lengths and were tracked for 20 s or longer. The ‘MeasureReversal’ plug-in was used to identify animals that reversed within 1 s of the mechanosensory stimulus being administered (Swierczek et al., 2011). Choreography output files were organized using custom R scripts which are freely available at https://github.com/troymcdiarmid/reversibility. All phenotypic features were pooled across the four to six plate replicates (each plate replicate capturing 40-100 animals) per strain. The mean of each condition was then compared using a two-tailed unpaired t-test and Benjamini–Hochberg control of false discovery rate at 0.01. Figures were generated using the ggplot2 and complex heatmap packages in R (Gu et al., 2016; Wickham, 2016).

Confocal imaging

Degradation groups were created as described above. All images were taken using 96-h-old animals. Adult animals were anesthetized by being placed in 10 µl of 7.5 mM Levamisole (Sigma-Aldrich) and 0.225 M BDM (2,3-butanedione monoxime) (Sigma-Aldrich) on glass microscope slide containing a 2% agar pad. After animals were paralyzed, a 1.5 coverslip was placed on top the agar pad. A Leica SP8 white light laser confocal microscope and 63× oil immersion lens was used for imaging. Step size was 0.3 µm. GFP was excited using a 488 nm wavelength laser with emitted light collected through a 490-778 nm bandpass filter. mNeonGreen was excited using a 506 nm wavelength laser with emitted light collected through a 512-742 nm bandpass filter. For each condition, four to five animals were scored for the presence of GFP or mNeonGreen. Final figures were generated using ImageJ (National Institutes of Health, Bethesda, MD, USA).

Strains used

The following strains used are available through the CGC: CA1200 ieSi57[eft-3p::TIR1::mRuby::unc-54 3'UTR+cbr-unc-119(+)] II; OH13988 ieSi57[eft-3p::TIR1::mRuby::unc-54 3'UTR+cbr-unc-119(+)] II; unc-3(ot837[unc-3::mNeonGreen::AID]) X; OH15227 unc-86(ot893[unc-86::3xFlag::mNeonGreen::AID]) III; CA1207 dhc-1(ie28[dhc-1::degron::GFP]) I; CA1210 ieSi57[eft-3p::TIR1::mRuby::unc-54 3'UTR+cbr-unc-119(+)] II; dhc-1(ie28[dhc-1::degron::GFP]) I.

The following strains were generated using standard genetic crosses: VG937 ieSi57[eft-3p::TIR1::mRuby::unc-54 3'UTR+cbr-unc-119(+)] II; unc-86(ot893[unc-86::3xFlag::mNeonGreen::AID]) III; VG946 mizSi6[rab-3p::TIR1::unc-54 3′UTR+LoxP pmyo-2::GFP::unc-54 3′UTR prps27::NeoR::unc-54 3′UTR; LoxP] V; dhc-1(ie28[dhc-1::degron::GFP]) I.

Genotype confirmation

Successful crosses were determined through visual confirmation of fluorescent reporters as well as PCR-based genotyping using the following primers: TIR1 sequence: forward, 5′-GACCGTAACTCCGTCTCC-3′; reverse, 5′-CGTTGGTGGTGATGATTTGAC-3′; AID degron sequence: forward, 5′-CCTAAAGATCCAGCCAAACC-3′; reverse, 5′-CTTCACGAACGCCGC-3′; or forward, 5′-GATCCAGCCAAACCTCCGGC-3′; reverse, 5′-CTTCACGAACGCCGCCGC-3′.

The authors thank Dr Kota Mizumoto for the rab-3p::TIR1 strain that was used for neuron-specific degradation and the CGC (funded by National Institutes of Health Office of Research Infrastructure Programs, P40 OD010440) for other strains used within this research.

Author contributions

Conceptualization: L.D.K., T.A.M.; Methodology: T.A.M.; Software: T.A.M.; Validation: L.D.K.; Formal analysis: T.A.M.; Investigation: L.D.K., T.A.M.; Data curation: T.A.M.; Writing - original draft: L.D.K., T.A.M.; Writing - review & editing: L.D.K., T.A.M., C.H.R.; Visualization: L.D.K., T.A.M.; Supervision: C.H.R.; Project administration: L.D.K.; Funding acquisition: C.H.R.

Funding

This project was supported by a Canadian Institutes of Health Research (CIHR) Doctoral Research Award to L.D.K.; a CIHR Doctoral Research Award to T.A.M.; and a CIHR project grant (CIHR MOP PJT-165947) to C.H.R. Open access funding provided by the CIHR project grant (CIHR MOP PJT-165947) to C.H.R. Deposited in PMC for immediate release.

Data availability

All analysis code is freely available at https://github.com/troymcdiarmid/reversibility. All data are freely available at https://dataverse.scholarsportal.info/dataset.xhtml?persistentId=doi:10.5683/SP3/NOZWXZ.

Abrahams
,
B. S.
,
Arking
,
D. E.
,
Campbell
,
D. B.
,
Mefford
,
H. C.
,
Morrow
,
E. M.
,
Weiss
,
L. A.
,
Menashe
,
I.
,
Wadkins
,
T.
,
Banerjee-Basu
,
S.
and
Packer
,
A.
(
2013
).
SFARI Gene 2.0: a community-driven knowledgebase for the autism spectrum disorders (ASDs)
.
Mol. Autism
4
,
36
.
Agapite
,
J.
,
Albou
,
L.-P.
,
Aleksander
,
S.
,
Argasinska
,
J.
,
Arnaboldi
,
V.
,
Attrill
,
H.
,
Bello
,
S. M.
,
Blake
,
J. A.
,
Blodgett
,
O.
,
Bradford
,
Y. M.
et al. 
(
2020
).
Alliance of Genome Resources Portal: unified model organism research platform
.
Nucleic Acids Res.
48
,
D650
-
D658
.
American Psychiatric Association
. (
2013
).
Diagnostic and Statistical Manual of Mental Disorders
.
American Psychiatric Association
.
Ardiel
,
E. L.
,
McDiarmid
,
T. A.
,
Timbers
,
T. A.
,
Lee
,
K. C. Y.
,
Safaei
,
J.
,
Pelech
,
S. L.
and
Rankin
,
C. H.
(
2018
).
Insights into the roles of CMK-1 and OGT-1 in interstimulus interval-dependent habituation in Caenorhabditis elegans
.
Proc. R. Soc. B Biol. Sci.
285
,
20182084
.
Ashley
,
G. E.
,
Duong
,
T.
,
Levenson
,
M. T.
,
Martinez
,
M. A. Q.
,
Johnson
,
L. C.
,
Hibshman
,
J. D.
,
Saeger
,
H. N.
,
Palmisano
,
N. J.
,
Doonan
,
R.
,
Martinez-Mendez
,
R.
et al. 
(
2021
).
An expanded auxin-inducible degron toolkit for Caenorhabditis elegans
.
Genetics
217
,
iyab006
.
Au
,
V.
,
Li-Leger
,
E.
,
Raymant
,
G.
,
Flibotte
,
S.
,
Chen
,
G.
,
Martin
,
K.
,
Fernando
,
L.
,
Doell
,
C.
,
Rosell
,
F. I.
,
Wang
,
S.
et al. 
(
2019
).
CRISPR/Cas9 methodology for the generation of knockout deletions in Caenorhabditis elegans
.
G3
9
,
135
-
144
.
Badea
,
T. C.
,
Cahill
,
H.
,
Ecker
,
J.
,
Hattar
,
S.
and
Nathans
,
J.
(
2009
).
Distinct roles of transcription factors Brn3a and Brn3b in controlling the development, morphology, and function of retinal ganglion cells
.
Neuron.
61
,
852
-
864
.
Beck
,
C. D. O.
and
Rankin
,
C. H.
(
1993
).
Effects of aging on habituation in the nematode Caenorhabditis elegans
.
Behav. Processes.
28
,
145
-
163
.
Belmadani
,
M.
,
Jacobson
,
M.
,
Holmes
,
N.
,
Phan
,
M.
,
Nguyen
,
T.
,
Pavlidis
,
P.
and
Rogic
,
S.
(
2019
).
VariCarta: a comprehensive database of harmonized genomic variants found in autism spectrum disorder sequencing studies
.
Autism Res.
12
,
1728
-
1736
.
Blok
,
L. E. R.
,
Boon
,
M.
,
van Reijmersdal
,
B.
,
Höffler
,
K. D.
,
Fenckova
,
M.
and
Schenck
,
A.
(
2022
).
Genetics, molecular control and clinical relevance of habituation
.
Preprints
2022010186
.
Boyle
,
C. A.
,
Boulet
,
S.
,
Schieve
,
L. A.
,
Cohen
,
R. A.
,
Blumberg
,
S. J.
,
Yeargin-Allsopp
,
M.
,
Visser
,
S.
and
Kogan
,
M. D.
(
2011
).
Trends in the prevalence of developmental disabilities in US children, 1997-2008
.
Pediatrics
127
,
1034
-
1042
.
Brenner
,
S.
(
1974
).
The genetics of Caenorhabditis elegans
.
Genetics.
77
,
71
-
94
.
Chao
,
H.-T.
,
Davids
,
M.
,
Burke
,
E.
,
Pappas
,
J. G.
,
Rosenfeld
,
J. A.
,
McCarty
,
A. J.
,
Davis
,
T.
,
Wolfe
,
L.
,
Toro
,
C.
,
Tifft
,
C.
et al. 
(
2017
).
A syndromic neurodevelopmental disorder caused by de novo variants in EBF3
.
Am. J. Hum. Genet.
100
,
128
-
137
.
Cianfrocco
,
M. A.
,
DeSantis
,
M. E.
,
Leschziner
,
A. E.
and
Reck-Peterson
,
S. L.
(
2015
).
Mechanism and regulation of cytoplasmic dynein
.
Annu. Rev. Cell Dev. Biol.
31
,
83
-
108
.
Creson
,
T. K.
,
Rojas
,
C.
,
Hwaun
,
E.
,
Vaissiere
,
T.
,
Kilinc
,
M.
,
Jimenez-Gomez
,
A.
,
Holder
,
J. L.
, Jr
,
Tang
,
J.
,
Colgin
,
L. L.
,
Miller
,
C. A.
et al. 
(
2019
).
Re-expression of SynGAP protein in adulthood improves translatable measures of brain function and behavior
.
eLife
8
,
e46752
.
de la Torre-Ubieta
,
L.
,
Won
,
H.
,
Stein
,
J. L.
and
Geschwind
,
D. H.
(
2016
).
Advancing the understanding of autism disease mechanisms through genetics
.
Nat. Med.
22
,
345
-
361
.
De Rubeis
,
S.
,
He
,
X.
,
Goldberg
,
A. P.
,
Poultney
,
C. S.
,
Samocha
,
K.
,
Ercument Cicek
,
A.
,
Kou
,
Y.
,
Liu
,
L.
,
Fromer
,
M.
,
Walker
,
S.
et al. 
(
2014
).
Synaptic, transcriptional and chromatin genes disrupted in autism
.
Nature
515
,
209
-
215
.
Deciphering Developmental Disorders Study
. (
2015
).
Large-scale discovery of novel genetic causes of developmental disorders
.
Nature
519
,
223
-
228
.
Dickinson
,
D. J.
and
Goldstein
,
B.
(
2016
).
CRISPR-based methods for Caenorhabditis elegans genome engineering
.
Genetics
202
,
885
-
901
.
Ehninger
,
D.
,
Li
,
W.
,
Fox
,
K.
,
Stryker
,
M. P.
and
Silva
,
A. J.
(
2008
).
Reversing neurodevelopmental disorders in adults
.
Neuron
60
,
950
-
960
.
Fenckova
,
M.
,
Blok
,
L. E. R.
,
Asztalos
,
L.
,
Goodman
,
D. P.
,
Cizek
,
P.
,
Singgih
,
E. L.
,
Glennon
,
J. C.
,
IntHout
,
J.
,
Zweier
,
C.
,
Eichler
,
E. E.
et al. 
(
2019
).
Habituation learning is a widely affected mechanism in Drosophila models of intellectual disability and autism spectrum
.
Disorders. Biol. Psychiatry
86
,
294
-
305
.
Feng
,
W.
,
Li
,
Y.
,
Dao
,
P.
,
Aburas
,
J.
,
Islam
,
P.
,
Elbaz
,
B.
,
Kolarzyk
,
A.
,
Brown
,
A. E.
and
Kratsios
,
P.
(
2020
).
A terminal selector prevents a Hox transcriptional switch to safeguard motor neuron identity throughout life
.
eLife
9
,
e50065
.
Gao
,
Y.
,
Irvine
,
E. E.
,
Eleftheriadou
,
I.
,
Naranjo
,
C. J.
,
Hearn-Yeates
,
F.
,
Bosch
,
L.
,
Glegola
,
J. A.
,
Murdoch
,
L.
,
Czerniak
,
A.
,
Meloni
,
I.
et al. 
(
2020
).
Gene replacement ameliorates deficits in mouse and human models of cyclin-dependent kinase-like 5 disorder
.
Brain
143
,
811
-
832
.
Green
,
S. A.
,
Hernandez
,
L.
,
Tottenham
,
N.
,
Krasileva
,
K.
,
Bookheimer
,
S. Y.
and
Dapretto
,
M.
(
2015
).
Neurobiology of sensory overresponsivity in youth with autism spectrum disorders
.
JAMA Psychiatry
72
,
778
-
786
.
Green
,
S. A.
,
Hernandez
,
L.
,
Lawrence
,
K. E.
,
Liu
,
J.
,
Tsang
,
T.
,
Yeargin
,
J.
,
Cummings
,
K.
,
Laugeson
,
E.
,
Dapretto
,
M.
and
Bookheimer
,
S. Y.
(
2019
).
Distinct patterns of neural habituation and generalization in children and adolescents with autism with low and high sensory overresponsivity
.
Am. J. Psychiatry
176
,
1010
-
1020
.
Gu
,
Z.
,
Eils
,
R.
and
Schlesner
,
M.
(
2016
).
Complex heatmaps reveal patterns and correlations in multidimensional genomic data
.
Bioinformatics
32
,
2847
-
2849
.
Guy
,
J.
,
Gan
,
J.
,
Selfridge
,
J.
,
Cobb
,
S.
and
Bird
,
A.
(
2007
).
Reversal of neurological defects in a mouse model of Rett Syndrome
.
Science
315
,
1143
-
1147
.
Hamill
,
D. R.
,
Severson
,
A. F.
,
Carter
,
J. C.
and
Bowerman
,
B.
(
2002
).
Centrosome maturation and mitotic spindle assembly in C. elegans require SPD-5, a protein with multiple coiled-coil domains
.
Dev. Cell
3
,
673
-
684
.
Harada
,
A.
,
Takei
,
Y.
,
Kanai
,
Y.
,
Tanaka
,
Y.
,
Nonaka
,
S.
and
Hirokawa
,
N.
(
1998
).
Golgi vesiculation and lysosome dispersion in cells lacking cytoplasmic dynein
.
J. Cell Biol.
141
,
51
-
59
.
Howell
,
A. M.
and
Rose
,
A. M.
(
1990
).
Essential genes in the hDf6 region of chromosome I in Caenorhabditis elegans
.
Genetics
126
,
583
-
592
.
Huang
,
E. J.
,
Liu
,
W.
,
Fritzsch
,
B.
,
Bianchi
,
L. M.
,
Reichardt
,
L. F.
and
Xiang
,
M.
(
2001
).
Brn3a is a transcriptional regulator of soma size, target field innervation and axon pathfinding of inner ear sensory neurons
.
Development
128
,
2421
-
2432
.
Husson
,
S. J.
,
Costa
,
W. S.
,
Schmitt
,
C.
and
Gottschalk
,
A.
(
2013
).
Keeping track of worm trackers
.
WormBook: the online review of C. elegans biology
,
1
-
17
.
Iakoucheva
,
L. M.
,
Muotri
,
A. R.
and
Sebat
,
J.
(
2019
).
Getting to the cores of autism
.
Cell
178
,
1287
-
1298
.
Iossifov
,
I.
,
O'Roak
,
B. J.
,
Sanders
,
S. J.
,
Ronemus
,
M.
,
Krumm
,
N.
,
Levy
,
D.
,
Stessman
,
H. A.
,
Witherspoon
,
K. T.
,
Vives
,
L.
,
Patterson
,
K. E.
et al. 
(
2014
).
The contribution of de novo coding mutations to autism spectrum disorder
.
Nature
515
,
216
-
221
.
Jin
,
X.
,
Simmons
,
S. K.
,
Guo
,
A.
,
Shetty
,
A. S.
,
Ko
,
M.
,
Nguyen
,
L.
,
Jokhi
,
V.
,
Robinson
,
E.
,
Oyler
,
P.
,
Curry
,
N.
et al. 
(
2020
).
In vivo Perturb-Seq reveals neuronal and glial abnormalities associated with autism risk genes
.
Science
370
,
eaaz6063
.
Kaletta
,
T.
and
Hengartner
,
M. O.
(
2006
).
Finding function in novel targets: C. elegans as a model organism
.
Nat. Rev. Drug Discov.
5
,
387
-
399
.
Kavšek
,
M.
(
2004
).
Predicting later IQ from infant visual habituation and dishabituation: a meta-analysis
.
J. Appl. Dev. Psychol.
25
,
369
-
393
.
Kepler
,
L. D.
,
McDiarmid
,
T. A.
and
Rankin
,
C. H.
(
2020
).
Habituation in high-throughput genetic model organisms as a tool to investigate the mechanisms of neurodevelopmental disorders
.
Neurobiol. Learn. Mem.
171
,
107208
.
Kim
,
W.
,
Underwood
,
R. S.
,
Greenwald
,
I.
and
Shaye
,
D. D.
(
2018
).
OrthoList 2: a new comparative genomic analysis of human and Caenorhabditis elegans genes
.
Genetics.
210
,
445
-
461
.
Kindt
,
K. S.
,
Quast
,
K. B.
,
Giles
,
A. C.
,
De
,
S.
,
Hendrey
,
D.
,
Nicastro
,
I.
,
Rankin
,
C. H.
and
Schafer
,
W. R.
(
2007
).
Dopamine mediates context-dependent modulation of sensory plasticity in C. elegans
.
Neuron
55
,
662
-
676
.
Kleinhans
,
N. M.
,
Johnson
,
L. C.
,
Richards
,
T.
,
Mahurin
,
R.
,
Greenson
,
J.
,
Dawson
,
G.
and
Aylward
,
E.
(
2009
).
Reduced neural habituation in the amygdala and social impairments in autism spectrum disorders
.
Am. J. Psychiatry
166
,
467
-
475
.
Kratsios
,
P.
,
Kerk
,
S. Y.
,
Catela
,
C.
,
Liang
,
J.
,
Vidal
,
B.
,
Bayer
,
E. A.
,
Feng
,
W.
,
De La Cruz
,
E. D.
,
Croci
,
L.
,
Consalez
,
G. G.
et al. 
(
2017
).
An intersectional gene regulatory strategy defines subclass diversity of C. elegans motor neurons
.
eLife
6
,
e25751
.
Levitan
,
D.
,
Doyle
,
T. G.
,
Brousseau
,
D.
,
Lee
,
M. K.
,
Thinakaran
,
G.
,
Slunt
,
H. H.
,
Sisodia
,
S. S.
and
Greenwald
,
I.
(
1996
).
Assessment of normal and mutant human presenilin function in Caenorhabditis elegans
.
Proc. Natl. Acad. Sci. USA
93
,
14940
-
14944
.
Li
,
Y.
,
Osuma
,
A.
,
Correa
,
E.
,
Okebalama
,
M. A.
,
Dao
,
P.
,
Gaylord
,
O.
,
Aburas
,
J.
,
Islam
,
P.
,
Brown
,
A. E.
and
Kratsios
,
P.
(
2020
).
Establishment and maintenance of motor neuron identity via temporal modularity in terminal selector function
.
eLife
9
,
e59464
.
Lopes
,
F.
,
Soares
,
G.
,
Gonçalves-Rocha
,
M.
,
Pinto-Basto
,
J.
and
Maciel
,
P.
(
2017
).
Whole gene deletion of EBF3 supporting haploinsufficiency of this gene as a mechanism of neurodevelopmental disease
.
Front. Genet.
8
.
Mains
,
P. E.
,
Sulston
,
I. A.
and
Wood
,
W. B.
(
1990
).
Dominant maternal-effect mutations causing embryonic lethality in Caenorhabditis elegans
.
Genetics
125
,
351
-
369
.
Massa
,
J.
and
O'Desky
,
I. H.
(
2012
).
Impaired Visual Habituation in Adults With ADHD
.
J. Atten. Disord.
16
,
553
-
561
.
McDiarmid
,
T. A.
,
Bernardos
,
A. C.
and
Rankin
,
C. H.
(
2017
).
Habituation is altered in neuropsychiatric disorders—A comprehensive review with recommendations for experimental design and analysis
.
Neurosci. Biobehav. Rev.
80
,
286
-
305
.
McDiarmid
,
T. A.
,
Au
,
V.
,
Loewen
,
A. D.
,
Liang
,
J.
,
Mizumoto
,
K.
,
Moerman
,
D. G.
and
Rankin
,
C. H.
(
2018a
).
CRISPR-Cas9 human gene replacement and phenomic characterization in Caenorhabditis elegans to understand the functional conservation of human genes and decipher variants of uncertain significance
.
Dis. Model. Mech.
11
,
dmm036517
.
McDiarmid
,
T. A.
,
Yu
,
A. J.
and
Rankin
,
C. H.
(
2018b
).
Beyond the response-High throughput behavioral analyses to link genome to phenome in Caenorhabditis elegans
.
Genes, Brain Behav.
17
,
e12437
.
McDiarmid
,
T. A.
,
Yu
,
A. J.
and
Rankin
,
C. H.
(
2019
).
Habituation is more than learning to ignore: multiple mechanisms serve to facilitate shifts in behavioral strategy
.
BioEssays
41
,
1900077
.
McDiarmid
,
T. A.
,
Belmadani
,
M.
,
Liang
,
J.
,
Meili
,
F.
,
Mathews
,
E. A.
,
Mullen
,
G. P.
,
Hendi
,
A.
,
Wong
,
W.-R.
,
Rand
,
J. B.
,
Mizumoto
,
K.
et al. 
(
2020a
).
Systematic phenomics analysis of autism-associated genes reveals parallel networks underlying reversible impairments in habituation
.
Proc. Natl. Acad. Sci. USA
117
,
656
-
667
.
McDiarmid
,
T. A.
,
Kepler
,
L. D.
and
Rankin
,
C. H.
(
2020b
).
Auxin does not affect a suite of morphological or behavioral phenotypes in two wild-type C. elegans strains
.
microPublication Biol.
2020
.
Mei
,
Y.
,
Monteiro
,
P.
,
Zhou
,
Y.
,
Kim
,
J.-A.
,
Gao
,
X.
,
Fu
,
Z.
and
Feng
,
G.
(
2016
).
Adult restoration of Shank3 expression rescues selective autistic-like phenotypes
.
Nature
530
,
481
-
484
.
Nance
,
J.
and
Frøkjær-Jensen
,
C.
(
2019
).
The Caenorhabditis elegans transgenic toolbox
.
Genetics
212
,
959
-
990
.
Nishimura
,
K.
,
Fukagawa
,
T.
,
Takisawa
,
H.
,
Kakimoto
,
T.
and
Kanemaki
,
M.
(
2009
).
An auxin-based degron system for the rapid depletion of proteins in nonplant cells
.
Nat. Methods
6
,
917
-
922
.
Parikshak
,
N. N.
,
Luo
,
R.
,
Zhang
,
A.
,
Won
,
H.
,
Lowe
,
J. K.
,
Chandran
,
V.
,
Horvath
,
S.
and
Geschwind
,
D. H.
(
2013
).
Integrative Functional Genomic Analyses Implicate Specific Molecular Pathways and Circuits in Autism
.
Cell
155
,
1008
-
1021
.
Post
,
K. L.
,
Belmadani
,
M.
,
Ganguly
,
P.
,
Meili
,
F.
,
Dingwall
,
R.
,
McDiarmid
,
T. A.
,
Meyers
,
W. M.
,
Herrington
,
C.
,
Young
,
B. P.
,
Callaghan
,
D. B.
et al. 
(
2020
).
Multi-model functionalization of disease-associated PTEN missense mutations identifies multiple molecular mechanisms underlying protein dysfunction
.
Nat. Commun.
11
,
2073
.
Prasad
,
B. C.
,
Ye
,
B.
,
Zackhary
,
R.
,
Schrader
,
K.
,
Seydoux
,
G.
and
Reed
,
R. R.
(
1998
).
unc-3, a gene required for axonal guidance in Caenorhabditis elegans, encodes a member of the O/E family of transcription factors
.
Development.
125
,
1561
-
1568
.
Prasad
,
B.
,
Karakuzu
,
O.
,
Reed
,
R. R.
and
Cameron
,
S.
(
2008
).
unc-3-dependent repression of specific motor neuron fates in Caenorhabditis elegans
.
Dev. Biol.
323
,
207
-
215
.
Rai
,
S.
and
Rankin
,
C. H.
(
2007
).
Critical and sensitive periods for reversing the effects of mechanosensory deprivation on behavior, nervous system, and development in Caenorhabditis elegans
.
Dev. Neurobiol.
67
,
1443
-
1456
.
Randlett
,
O.
,
Haesemeyer
,
M.
,
Forkin
,
G.
,
Shoenhard
,
H.
,
Schier
,
A. F.
,
Engert
,
F.
and
Granato
,
M.
(
2019
).
Distributed plasticity drives visual habituation learning in larval zebrafish
.
Curr. Biol.
29
,
1337
-
1345.e4
.
Rankin
,
C. H.
,
Beck
,
C. D. O.
and
Chiba
,
C. M.
(
1990
).
Caenorhabditis elegans: a new model system for the study of learning and memory
.
Behav. Brain Res.
37
,
89
-
92
.
Rankin
,
C. H.
,
Abrams
,
T.
,
Barry
,
R. J.
,
Bhatnagar
,
S.
,
Clayton
,
D. F.
,
Colombo
,
J.
,
Coppola
,
G.
,
Geyer
,
M. A.
,
Glanzman
,
D. L.
,
Marsland
,
S.
et al. 
(
2009
).
Habituation revisited: an updated and revised description of the behavioral characteristics of habituation
.
Neurobiol. Learn. Mem.
92
,
135
-
138
.
Robinson
,
J. T.
,
Wojcik
,
E. J.
,
Sanders
,
M. A.
,
McGrail
,
M.
and
Hays
,
T. S.
(
1999
).
Cytoplasmic dynein is required for the nuclear attachment and migration of centrosomes during mitosis in Drosophila
.
J. Cell Biol.
146
,
597
-
608
.
Salyakina
,
D.
,
Cukier
,
H. N.
,
Lee
,
J. M.
,
Sacharow
,
S.
,
Nations
,
L. D.
,
Ma
,
D.
,
Jaworski
,
J. M.
,
Konidari
,
I.
,
Whitehead
,
P. L.
,
Wright
,
H. H.
et al. 
(
2011
).
Copy number variants in extended autism spectrum disorder families reveal candidates potentially involved in autism risk
.
PLoS ONE
6
,
e26049
.
Sanders
,
S. J.
,
He
,
X.
,
Willsey
,
A. J.
,
Ercan-Sencicek
,
A. G.
,
Samocha
,
K. E.
,
Cicek
,
A. E.
,
Murtha
,
M. T.
,
Bal
,
V. H.
,
Bishop
,
S. L.
,
Dong
,
S.
et al. 
(
2015
).
Insights into autism spectrum disorder genomic architecture and biology from 71 risk loci
.
Neuron
87
,
1215
-
1233
.
Sanders
,
S. J.
,
Sahin
,
M.
,
Hostyk
,
J.
,
Thurm
,
A.
,
Jacquemont
,
S.
,
Avillach
,
P.
,
Douard
,
E.
,
Martin
,
C. L.
,
Modi
,
M. E.
,
Moreno-De-Luca
,
A.
et al. 
(
2019
).
A framework for the investigation of rare genetic disorders in neuropsychiatry
.
Nat. Med.
25
,
1477
-
1487
.
Satterstrom
,
F. K.
,
Kosmicki
,
J. A.
,
Wang
,
J.
,
Breen
,
M. S.
,
De Rubeis
,
S.
,
An
,
J.-Y.
,
Peng
,
M.
,
Collins
,
R.
,
Grove
,
J.
,
Klei
,
L.
, et al. 
(
2020
).
Large-scale exome sequencing study implicates both developmental and functional changes in the neurobiology of autism
.
Cell
180
,
P568
-
P584.e23
.
Schmid
,
S.
,
Wilson
,
D. A.
and
Rankin
,
C. H.
(
2015
).
Habituation mechanisms and their importance for cognitive function
.
Front. Integr. Neurosci.
8
,
97
.
Serrano-Saiz
,
E.
,
Leyva-Díaz
,
E.
,
De La Cruz
,
E.
and
Hobert
,
O.
(
2018
).
BRN3-type POU homeobox genes maintain the identity of mature postmitotic neurons in nematodes and mice
.
Curr. Biol.
28
,
2813
-
2823.e2
.
Sleven
,
H.
,
Welsh
,
S. J.
,
Yu
,
J.
,
Churchill
,
M. E. A.
,
Wright
,
C. F.
,
Henderson
,
A.
,
Horvath
,
R.
,
Rankin
,
J.
,
Vogt
,
J.
,
Magee
,
A.
et al. 
(
2017
).
De novo mutations in EBF3 cause a neurodevelopmental syndrome
.
Am. J. Hum. Genet.
100
,
138
-
150
.
Swierczek
,
N. A.
,
Giles
,
A. C.
,
Rankin
,
C. H.
and
Kerr
,
R. A.
(
2011
).
High-throughput behavioral analysis in C. elegans
.
Nature
8
,
592
-
598
.
Sze
,
J. Y.
and
Ruvkun
,
G.
(
2003
).
Activity of the Caenorhabditis elegans UNC-86 POU transcription factor modulates olfactory sensitivity
.
Proc. Natl. Acad. Sci. USA
100
,
9560
-
9565
.
Tanaka
,
A. J.
,
Cho
,
M. T.
,
Willaert
,
R.
,
Retterer
,
K.
,
Zarate
,
Y. A.
,
Bosanko
,
K.
,
Stefans
,
V.
,
Oishi
,
K.
,
Williamson
,
A.
,
Wilson
,
G. N.
et al. 
(
2017
).
De novo variants in EBF3 are associated with hypotonia, developmental delay, intellectual disability, and autism
.
Mol. Case Stud.
3
,
a002097
.
Thompson
,
R. F.
and
Spencer
,
W. A.
(
1966
).
Habituation: a model phenomenon for the study of neuronal substrates of behavior
.
Psychol. Rev.
73
,
16
-
43
.
Timbers
,
T. A.
,
Giles
,
A. C.
,
Ardiel
,
E. L.
,
Kerr
,
R. A.
and
Rankin
,
C. H.
(
2013
).
Intensity discrimination deficits cause habituation changes in middle-aged Caenorhabditis elegans
.
Neurobiol. Aging
34
,
621
-
631
.
Vissers
,
L. E. L. M.
,
Gilissen
,
C.
and
Veltman
,
J. A.
(
2016
).
Genetic studies in intellectual disability and related disorders
.
Nat. Rev. Genet.
17
,
9
-
18
.
Vogel-Ciernia
,
A.
,
Matheos
,
D. P.
,
Barrett
,
R. M.
,
Kramár
,
E. A.
,
Azzawi
,
S.
,
Chen
,
Y.
,
Magnan
,
C. N.
,
Zeller
,
M.
,
Sylvain
,
A.
,
Haettig
,
J.
et al. 
(
2013
).
The neuron-specific chromatin regulatory subunit BAF53b is necessary for synaptic plasticity and memory
.
Nat. Neurosci.
16
,
552
-
561
.
Wang
,
S.
,
Tang
,
N. H.
,
Lara-Gonzalez
,
P.
,
Zhao
,
Z.
,
Cheerambathur
,
D. K.
,
Prevo
,
B.
,
Chisholm
,
A. D.
,
Desai
,
A.
and
Oegema
,
K.
(
2017
).
A toolkit for GFP-mediated tissue-specific protein degradation in C. elegans
.
Development
144
,
2694
-
2701
.
Wickham
,
H.
(
2016
).
ggplot2: Elegant Graphics for Data Analysis
.
New York
:
Springer-Verlag
.
Willemsen
,
M. H.
,
Nijhof
,
B.
,
Fenckova
,
M.
,
Nillesen
,
W. M.
,
Bongers
,
E. M. H. F.
,
Castells-Nobau
,
A.
,
Asztalos
,
L.
,
Viragh
,
E.
,
van Bon
,
B. W. M.
,
Tezel
,
E.
et al. 
(
2013
).
GATAD2B loss-of-function mutations cause a recognisable syndrome with intellectual disability and are associated with learning deficits and synaptic undergrowth in Drosophila
.
J. Med. Genet.
50
,
507
-
514
.
Williams
,
L. E.
,
Blackford
,
J. U.
,
Luksik
,
A.
,
Gauthier
,
I.
and
Heckers
,
S.
(
2013
).
Reduced habituation in patients with schizophrenia
.
Schizophr. Res.
151
,
124
-
132
.
Willsey
,
A. J.
,
Sanders
,
S. J.
,
Li
,
M.
,
Dong
,
S.
,
Tebbenkamp
,
A. T.
,
Muhle
,
R. A.
,
Reilly
,
S. K.
,
Lin
,
L.
,
Fertuzinhos
,
S.
,
Miller
,
J. A.
et al. 
(
2013
).
Coexpression networks implicate human midfetal deep cortical projection neurons in the pathogenesis of autism
.
Cell
155
,
997
-
1007
.
Xiang
,
M.
,
Zhou
,
L.
,
Macke
,
J. P.
,
Yoshioka
,
T.
,
Hendry
,
S. H.
,
Eddy
,
R. L.
,
Shows
,
T. B.
and
Nathans
,
J.
(
1995
).
The Brn-3 family of POU-domain factors: primary structure, binding specificity, and expression in subsets of retinal ganglion cells and somatosensory neurons
.
J. Neurosci.
15
,
4762
-
4785
.
Zeier
,
Z.
,
Kumar
,
A.
,
Bodhinathan
,
K.
,
Feller
,
J. A.
,
Foster
,
T. C.
and
Bloom
,
D. C.
(
2009
).
Fragile X mental retardation protein replacement restores hippocampal synaptic function in a mouse model of fragile X syndrome
.
Gene Ther.
16
,
1122
-
1129
.
Zhang
,
L.
,
Ward
,
J. D.
,
Cheng
,
Z.
and
Dernburg
,
A. F.
(
2015
).
The auxin-inducible degradation (AID) system enables versatile conditional protein depletion in C. elegans
.
Development
142
,
4374
-
4384
.
Zou
,
M.
,
Li
,
S.
,
Klein
,
W. H.
and
Xiang
,
M.
(
2012
).
Brn3a/Pou4f1 regulates dorsal root ganglion sensory neuron specification and axonal projection into the spinal cord
.
Dev. Biol.
364
,
114
-
127
.

Competing interests

The authors declare no competing or financial interests.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

Supplementary information