DYRK1A-haploinsufficiency in mice causes autistic-like features and febrile seizures

Accepted Manuscript DYRK1A-haploinsufficiency in mice causes autistic-like features and febrile seizures

Matthieu Raveau, Atsushi Shimohata, Kenji Amano, Hiroyuki Miyamoto, Kazuhiro Yamakawa PII: DOI: Reference:

S0969-9961(17)30281-4 doi:10.1016/j.nbd.2017.12.003 YNBDI 4075

To appear in:

Neurobiology of Disease

Received date: Revised date: Accepted date:

22 June 2017 28 November 2017 5 December 2017

Please cite this article as: Matthieu Raveau, Atsushi Shimohata, Kenji Amano, Hiroyuki Miyamoto, Kazuhiro Yamakawa , DYRK1A-haploinsufficiency in mice causes autisticlike features and febrile seizures. The address for the corresponding author was captured as affiliation for all authors. Please check if appropriate. Ynbdi(2017), doi:10.1016/ j.nbd.2017.12.003

This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

ACCEPTED MANUSCRIPT

DYRK1A-haploinsufficiency in mice causes autistic-like features and febrile seizures Matthieu Raveau1, Atsushi Shimohata1, Kenji Amano1, Hiroyuki Miyamoto1,

1

PT

Kazuhiro Yamakawa1*

Laboratory for Neurogenetics, RIKEN Brain Science Institute, Saitama 351-0198,

RI

Japan.

AC CE P

TE

D

MA

NU

SC

*To whom correspondence should be addressed. Email: [email protected]

1

ACCEPTED MANUSCRIPT

Abstract Mutations and copy number variants affecting DYRK1A gene encoding the dualspecificity tyrosine phosphorylation-regulated kinase 1A are among the most frequent genetic causes of neurodevelopmental disorders including autism spectrum disorder (ASD) associated with microcephaly, febrile seizures and severe speech acquisition

PT

delay. Here we developed a mouse model harboring a frame-shift mutation in Dyrk1a resulting in a protein truncation and elimination of its kinase activity site. Dyrk1a+/−

RI

mice showed significant impairments in cognition and cognitive flexibility, communicative ultrasonic vocalizations, and social contacts. Susceptibility to

SC

hyperthermia-induced seizures was also significantly increased in these mice. The truncation leading to haploinsufficiency of DYRK1A in mice thus recapitulates the

NU

syndromic phenotypes observed in human patients and constitutes a useful model for further investigations of the mechanisms leading to ASD, speech delay and febrile

MA

seizures.

D

Keywords: DYRK1A; cognitive flexibility; autism; sociability deficit;

AC CE P

TE

communication deficit; febrile seizure

2

ACCEPTED MANUSCRIPT

Introduction Over the past few years, advances in sequencing technologies uncovered genes responsible for neurodevelopmental disorders such as autism spectrum disorder (ASD)(DDD Study, 2015; DDD Study, 2017; Krumm et al., 2014). Among the recurrent high-score candidates lies the dual-specificity tyrosine phosphorylation

PT

regulated kinase 1A gene DYRK1A (DDD Study, 2015; DDD Study, 2017; Krumm et al., 2014). Deletions (van Bon et al., 2011) and mutations in DYRK1A (Courcet et al.,

RI

2012; Møller et al., 2008; O’Roak et al., 2012) were first described as leading to a severe intellectual disability (ID) associated with microcephaly (OMIM 614104). More

SC

recent reports describe cases of partial deletions, missense, nonsense mutations and frameshifts resulting in syndromic phenotypes (Bronicki et al., 2015; Ji et al., 2015; van

NU

Bon et al., 2016). Regardless of the type and location of these genetic variations, all patients show microcephaly, developmental delay, ID and severe speech delay, ASD,

MA

febrile seizures and abnormal gait/clumsiness.

DYRK1A belongs to a family of dual-specificity protein kinases highly conserved across mammalian species (Becker et al., 1998). In rodents, Dyrk1a is widely

D

expressed in the central nervous system and the DYRK1A protein cellular localization

TE

is controlled by a combination of phosphorylation sites (Kaczmarki et al., 2014; Martí et al., 2003). In neuronal cells, DYRK1A localizes in the nucleus as well as in the

AC CE P

cytoplasm associated with proteins of the cytoskeleton (Hämmerle et al., 2008). Interactions with- and phosphorylation by DYRK1A have been identified for a number of proteins involved in various signaling pathways and were extensively reviewed recently (Duchon and Herault, 2016). Among the main roles identified so far, Dyrk1a regulates the cell cycle, differentiation and arborization of neuronal cells, vesicle trafficking and neurotransmission (Dang et al., 2017; Grau et al., 2014; Hämmerle et al., 2011; Ori-McKenney et al., 2016; Park et al., 2012; Soppa et al., 2014). All these studies, mostly realized in-vitro point to a key function of Dyrk1a in both neuronal circuit development and function. Models of gain-of-function and loss-of-function in mouse so far support these observations and suggest a tight regulation of Dyrk1a expression is necessary. Decreases in dendritic branching, synaptic spine formation and synaptic transmission

3

ACCEPTED MANUSCRIPT

have indeed been observed in a model of Dyrk1a overexpression (Ahn et al., 2006; Martinez De Lagran et al., 2012). In these transgenic mice, the overexpression of Dyrk1a leads to a delayed acquisition of neurodevelopmental landmarks and affects learning and memory skills (Ahn et al., 2006; Altafaj et al., 2001; Souchet et al., 2014). Interestingly, the loss-of-function leads to comparable deficits in the development of

PT

dendritic arborization and synaptic spines (Benavides-Piccione et al., 2005). Even if the knowledge on Dyrk1a knockout remains limited, the first reports available so far

RI

describe a global growth delay, microcephaly and decreased grip strength but conserved motor control performance (Fotaki et al., 2002; Fotaki et al., 2004). Additionally, parts

SC

of the hippocampus-dependent memory functions are deficient in Dyrk1a knockout mice: the working memory appears to be deficient (Arqué et al., 2008) whereas the

NU

spatial learning and memory is rather mildly or non-affected depending on the testing paradigm (Arqué et al., 2008; Arqué et al., 2009). Above all, the major syndromic

MA

phenotypes observed in human patients with DYRK1A mutations have not been reported in rodent models so far.

In the present work, we generated a mouse model carrying a deleterious point

D

mutation in Dyrk1a and focused on characterizing its social behavior, communication

TE

and susceptibility to epileptic seizures using multiple experimental approaches.

AC CE P

Materials and methods

CRISPR/Cas9-mediated generation of a frameshift in Dyrk1a exon 3. Dyrk1a heterozygous knockout mice (Dyrk1a+/−) were engineered using CRISPR/Cas9-mediated mutagenesis as described by Wang and colleagues (Wang et al., 2013). Briefly, two single-guide RNA (sgRNA) were designed to target the 3’ end of Dyrk1a exon 3 (Figure 1A). Forward and reverse oligonucleotides corresponding to these sgRNA (Supplementary table 1) were annealed and cloned into pX330-U6Chimeric_BB-CBh-hSpCas9 plasmid digested with BbsI (Cong et al., 2013). T7 promoter was added to the Cas9 coding fragment by PCR amplification using primers T7_Cas9_Fw and T7_Cas9_rev (Supplementary table 1), purified and transcribed invitro using mMESSAGE mMACHINE T7 ULTRA kit (Ambion #AM1345). T7 promoter was added to the sgRNA fragments by PCR amplification using primers

4

ACCEPTED MANUSCRIPT

T7_guide#1_Fw, T7_guide#2_Fw and T7_guide_Rev (Supplementary table 1), purified and transcribed in-vitro using MEGAshortscript T7 kit (Ambion #1354). Cas9 mRNA and the sgRNAs were purified using MEGAclear Transcription Clean-Up Kit (Ambion #1908). Cas9 mRNA (200 ng/μl) and sgRNAs (50 ng/μl) were microinjected into the

PT

cytoplasm of C57BL/6J fertilized eggs. Two-cells embryos were then implanted into the oviduct of pseudopregnant ICR females. Two pups carrying a frameshift caused by a

RI

single nucleotide insertion were used as founder mice (c.134insC from sgRNA #2 and c.152insC+A>C from sgRNA #1, Figure 1B). These lines were backcrossed for four

SC

generations on C57BL/6J background before the start of the phenotype screening. Genotype was assessed by sequencing a PCR product amplified using the following

NU

primers: Dyrk1a_Fs_Fw: CGTATGTCAGGCTGTGAGCAC and Dyrk1a_Fs_Rev: ATGGTAAGGCAGACACCTGC.

MA

Putative offtargets were estimated using the online “CRISPR Design” tool designed by Dr. Zhang’s laboratory at MIT (http://crispr.mit.edu). For offtarget screening, genomic DNA was prepared from the tail of two WT and two mutant mice

D

for each line. The top five offtarget loci (in order of hit score) were checked by

Western-blotting

TE

sequencing PCR fragments encompassing the potential offtarget sites.

AC CE P

Brains were sampled from 1.5 months old mice (N=3 animals per genotype for each line, total N=6 per genotype) and flash frozen in liquid nitrogen. One hemisphere per sample was homogenized in ice-cold deoxycholate buffer (1% deoxycholate, 50 mM Tris-HCl pH8.5, 150 mM NaCl, 10% glycerol) supplemented with protease inhibitors (Complete Protease Inhibitor Cocktail; Roche Applied Science, Mannheim, Germany) and phosphatases inhibitors (PhosSTOP; Roche Applied Science). Homogenates were centrifuged twice at 23,500 g for 60 min. Supernatants were diluted in sample buffer (2.3% SDS, 62.5 mM Tris-HCl pH6.8, 5% 2-mercaptoethanol) to reach a final total protein concentration of 4.2 g/l. 15 l of each sample was loaded on 520% polyacrylamide SuperSep HG gels (Wako pure reagents, Osaka, Japan), and proteins were transferred to nitro-cellulose membrane (Bio-Rad Laboratories Inc., Hercules, CA). The membrane was blocked using 3.5% skim milk (Wako pure

5

ACCEPTED MANUSCRIPT

reagents) in TBS-T (0.1% Tween-20/TBS pH 7.4) for 1 h and incubated overnight with a rabbit anti DYRK1A C-terminus antibody (1:1,000; Sigma D1819) or a goat antiDYRK1A N-terminus antibody (1:200; Santa Cruz sc-12568). The membrane was then washed in 0.1% TBS-T and incubated for 90 min with a HRP-conjugated anti-rabbit or anti-goat IgG antibody diluted in blocking buffer (1:10 000, Jackson Immuno Research

PT

Laboratories, Inc., West Grove, PA). Labeled proteins were revealed using enhanced chemiluminescence (ECL) detection (Perkin-Elmer, Waltham, MA). Membranes were

RI

then washed using Restore Plus Western Blot stripping buffer (Pierce, Rockford, IL), re-probed with rabbit anti-GAPDH antibody (1:2000; Santa Cruz Biotechnology, Dallas,

SC

TX) and HRP-conjugated goat anti-rabbit antibody (1:10 000; BA-1000, Vector Laboratories) and revealed as described above. Band intensities were quantified using

NU

NIH ImageJ software (developed at U.S. National Institute of Health and available at http://rsb.info.nih.gov/ij/).

MA

Brain morphology

For both frameshift lines, brains from four adult (4 months old) Dyrk1a+/- and four WT littermates were sampled, weighted and post-fixed in 4% PFA. Macroscopic

D

pictures were acquired for gross morphology comparison.

TE

Two pairs of Dyrk1a+/- and WT brains per frameshift line were embedded in paraffin. Equivalent parasagittal sections (6μm thickness) were prepared according to

AC CE P

common landmarks (Paxinos et al. 2001) and hematoxylin-eosin stained for morphometric measurements. Images were acquired using a BZ-X710 light microscope (Keyence) and processed using NIH ImageJ software for area measurement. A total of three sections were measured for each animal in a blinded manner.

Animals and experimental conditions All animal breeding and experimental procedures were performed in accordance with the ARRIVE guidelines and the guidelines of the Animal Experiments Committee of RIKEN Brain Science Institute. Animals were maintained on 12 hours light/dark cycle with ad libitum access to food and water. On experimental days, mice were moved to the testing rooms and given a one-hour habituation before the start of the test. Devices were cleaned and wiped with 70% ethanol and water. In order to comply with the “3R” rule, a total of thirteen

6

ACCEPTED MANUSCRIPT

Dyrk1a+/− and thirteen WT littermates were generated and underwent the screening process as follows: open field, elevated plus maze, 3-chambers social behavior, tube test, tail suspension test and rotarod. Three Dyrk1a+/− and three WT mice were then used for electrocorticogram recording whereas the remaining ten Dyrk1a+/− and ten WT animals were submitted to Barnes maze task and pentylenetetrazol injections.

PT

Independent groups of mice were generated for ultrasonic vocalization (N=13 WT and 11 Dyrk1a+/−), kainic acid injection (N=5 WT and 5 Dyrk1a+/−), febrile seizure

RI

experiments (N=10 WT and 10 Dyrk1a+/−) and buried food test (N=8 WT and 8 Dyrk1a+/−). All behavioral tasks were assessed using exclusively male mice, except for

SC

ultrasonic vocalization (P6 pups; WT: 6 males + 7 females; Dyrk1a+/-: 5 males + 6 females). The detailed ages and mouse number are summarized in the Supplementary

NU

table 2.

Open-field

MA

Mice were placed in a 60 x 60 cm square automated open-field homogeneously illuminated at 70 lux and allowed to freely explore for 30 min. Data was acquired and analyzed using manufacturer’s tracking software (TimeOFCR4; O’Hara & Co, Ltd.,

AC CE P

Elevated Plus maze

TE

the open field.

D

Tokyo, Japan). The center area was defined as a 36 x 36 cm square area at the center of

The elevated plus maze consists of two open arms (25 × 5 cm) crossing two enclosed arms of the same size with 15 cm high transparent walls and was placed 50 cm above the floor. Luminosity was homogeneously set at 70 lux. Mice were placed at the center of the maze, facing one of the open arms and allowed to freely explore for 10 min. Data was acquired and analyzed using manufacturer’s tracking software (TimeEP2; O’Hara & Co, Ltd., Tokyo, Japan).

Barnes maze The Barnes maze apparatus consists of a white circular board (100 cm diameter) with twelve holes evenly spaced at the periphery. Four large visual cues were placed around the apparatus to be used as spatial landmarks. The board was brightly illuminated at 700 lux and elevated 100 cm from the floor. Mice were placed in a box at

7

ACCEPTED MANUSCRIPT

the center of the board and automatically released at the start of each session. Mice were first given a 5 min habituation to the board and 1 min habituation to a black Plexiglas escape box. For four consecutive days, animals were submitted to three daily trials in which they were given a maximum of 5 min to explore the board, find and enter the escape box placed under one of the peripheral holes. The location of the box was

PT

randomized across the animals but remained identical for a given mouse. A probe test (3 min free exploration without escape box) was conducted 24 hours after the last training

RI

session. Reverse training was conducted in similar way (four consecutive days, three trials per day), with the escape box placed 180° from its original location. A “reverse

SC

probe test” was conducted 24 hours after the last reverse training trial. Data was acquired and analyzed using manufacturer’s tracking software (TimeBCM; O’Hara &

NU

Co, Ltd., Tokyo, Japan).

Rotarod

MA

The rotarod apparatus consists of a 30 mm diameter, polyvinyl chloride-coated rotating rod (MK-610A, Muromachi Kikai Co., Ltd., Tokyo, Japan). Mice were first given a 1 min habituation to the rod set at 4 rpm rotation speed. They were then tested

D

four trials a day for four consecutive days on the rod linearly accelerating from 4 rpm to

TE

40 rpm within 5 min. Latency to fall off the rod and maximal rotation speed were automatically recorded. A maximum of 300 sec / 40 rpm was attributed whenever

AC CE P

animals did not fall off the rod. Results of 4 sessions were averaged for each day.

Social behavior: 3-chambers task The 3-chamber apparatus consists in a 43 x 63 cm transparent Plexiglas box, separated into 3 equivalent-sized (21 x 43 cm) chambers by transparent Plexiglas plates with square openings to allow mice to move freely from a chamber to another. One wire quarter-cylinder shaped cage was placed in a corner of the side chambers and used to enclose 10 weeks old C57BL/6J stranger mice. Tested mice were first given a 10 minutes habituation period to the apparatus. A first stranger mouse was then enclosed randomly in one of the side-chambers and the tested mice were given 10 minutes exploration. In a final step a second stranger was enclosed in the opposite side-chamber and the tested mice were given another 10 minutes exploration. Data was acquired and

8

ACCEPTED MANUSCRIPT

analyzed using manufacturer’s tracking software (TimeCSI; O’Hara & Co, Ltd., Tokyo, Japan).

Social dominance: tube test Mice were first allowed to run twice through the tube to habituate to the apparatus (cylindrical transparent tube, diameter: 3 cm; length 30 cm; elevation from

PT

the table: 2 cm). Light intensity was set at 5 lux and videos recorded using an infrared camera. For every testing round, one WT and one Dyrk1a+/− mouse were introduced in

RI

opposite ends of the tube facing each other and released simultaneously. An animal was

SC

considered to have lost the round when both of its posterior feet touched the table, either from retreating or being pushed out of the tube. Each mouse was submitted to four rounds with a different opponent every time (i.e. every Dyrk1a+/- animal was paired

NU

with four different WT mice and every WT animal paired with four different Dyrk1a+/mice). Mice combinations and starting side were randomized to prevent any bias.

MA

Buried food test

After a 16 hours fasting period, every mouse was placed in a cage containing

D

about 3 cm of clean white paper-chip bedding (Alpha-Dri, Shepherd Specialty Products,

TE

Kalamazoo, MI) and allowed to explore freely for 5 minutes for habituation. The mouse was then removed, a small piece of cookie was buried 1 cm beneath the surface in a

AC CE P

corner of the cage, and the mouse was then placed back in the opposite corner, facing the wall. The latency to find the hidden food was recorded. The cutoff time was set at 15 minutes, but all mice were able to find the buried food within the allocated time. We defined “food finding” when the mouse started eating and/or holding the pellet in the forepaws.

Tail suspension test

Mice were suspended 30 cm above the floor of a while plastic box (31 x 41 x 41 cm; O’Hara & Co, Ltd., Tokyo, Japan) using adhesive tape placed 1 cm from the base of the tail. Luminosity was set at 70 lux in the closed box. Behavior was recorded for 10 min and data analyzed automatically using manufacturer’s tracking software (TimeFZ2; O’Hara & Co, Ltd., Tokyo, Japan).

9

ACCEPTED MANUSCRIPT

Ultrasonic vocalizations (USV) Isolation-induced USV were recorded from pups separated from their mother for 5 min at post-natal day 6. Pups were individually removed from the nest and placed in a plastic cylindrical container on top of fresh nesting material. The container was placed in an acoustic attenuation box and a microphone adjusted 10 cm above the pup.

PT

Vocalizations were recorded for 5 minutes at a 250,000 Hz sampling rate in 16-bit format using Avisoft Recorder (version 4.2.8; Avisoft Bioacoustics, Berlin, Germany).

RI

One Dyrk1a+/− pup did not produce any vocalization and was excluded from the calculation to avoid a bias.

SC

Recorded files were processed using Avisoft SAS Lab Pro (version 5.1; Avisoft Bioacoustics, Berlin, Germany). A fast Fourier transform was applied (512 FFT length,

NU

100% frame size, Hamming window and 75% time window overlap) and spectrograms generated with a 0.512 ms time resolution and 488 Hz frequency resolution. A high-

MA

pass filter set at 30 kHz was applied to reduce background noise outside of the USV relevant frequency band. USV calls were automatically detected using a double threshold based algorithm (amplitude: -40dB, hold time: 10 ms). An experimenter blind

D

to the genotype of the pups verified the accuracy of the detected calls and, when

TE

necessary, manually added or corrected missed and incorrectly detected calls. Call duration, peak amplitude, peak frequency, minimum frequency and frequency

AC CE P

modulation (Figure 4A) were averaged for every mouse. Calls were classified into ten categories (Figure 4H) commonly observed in C57BL/6 pups (Roy et al., 2012; Scattoni et al., 2008).

Seizure susceptibility

For hyperthermia-induced seizures, the baseline body temperature and controlled temperature increase were monitored using a rectal probe (IT-18; Physitemp Instruments, Clifton, NJ, USA). Mice were placed in a hermetic plexiglas box (30×30×30 cm) heated from below through a perforated floor plate (81 evenly spaced holes of 2.5 mm diameter) using an Air-Therm heater (World Precision Instruments, Sarasota, FL, USA). Temperature was first maintained at 37°C until animals body temperature reached 39°C. The air temperature was then raised by 0.5°C every minute, animals body temperature monitored and their behavior video-recorded. At the onset of

10

ACCEPTED MANUSCRIPT

a generalized seizure mice were promptly removed from the apparatus and placed in an ice-cold box until their body temperature came back to the baseline value. Kainic acid monohydrate (KA; K0250, SIGMA Aldrich, Saint-Louis, MO, USA) was dissolved in a drop of 1N NaOH, diluted in 0.9% saline and administered intraperitoneously (30 mg/kg; total volume: 300 μl). Animal behavior was then

PT

carefully watched and video-recorded for a maximum of 60 min. Latency to generalized seizures and death were recorded.

RI

Pentylenetetrazole (PTZ; P6500, SIGMA Aldrich, Saint-Louis, MO, USA) dissolved in 1X phosphate-buffered saline (PBS) was administered intraperitoneously

SC

(50 mg/kg; total volume: 300 μl). Animal behavior was then carefully watched and video-recorded for a maximum of 10 min. Severity score were given on a 1 to 5 scale as

NU

follows: 1) myoclonic, 2) tonic, 3) generalized seizure, 4) full body and limbs extension, 5) death. Latency to generalized seizures and death were recorded and a maximum

MA

value of 600 s was given to mice that did not reach these stages.

Somatosensory ECoG (electrocorticogram) and EMG (electromyogram) recordings

D

Adult male mice (20 weeks old) were used in this study. Stainless steel screws

TE

(1.1 mm diameter) serving as ECoG electrodes were inserted through the skull in contact with the somatosensory cortex (±1.5 mm lateral and 1.0 mm posterior from

AC CE P

bregma) under 1–1.5% isoflurane anesthesia. A screw electrode in contact with the cerebellum was used as reference. Stainless wires bipolar electrodes (100 m diameter) were inserted in the cervical region of the trapezius muscle to record EMG. After a minimum of one week of recovery, local field potentials (LFPs) were recorded for 3 continuous days (sampled at 256 Hz) and analyzed off-line (SleepSign, Kissay, Japan). A 3 hours window was chosen within the dark part of the daily light/dark cycle on the third day recorded and was used for further analysis. This time window provides sufficient data to characterize the main behavioral states (waking (W), non-REM sleep (NREM), REM sleep (REM)) that were then assigned manually by an experimenter blinded to sample’s genotype. Power spectrum densities for each behavioral state were automatically generated by a fast Fourier transform.

11

ACCEPTED MANUSCRIPT

Statistical analysis Statistical significance was assessed using one-way ANOVA for parameters with a single value per individual and two-way repeated measures ANOVA for parameters involving repeated measures. Statistical significance was calculated using KyPlot v2.0 software (Kyens Lab, Tokyo, Japan). Experiments were conducted in a

PT

blinded manner.

RI

Results

Creation of a mouse model of Dyrk1a truncation mutation using the

SC

CRISPR/Cas9 system

In order to create a model of DYRK1A truncation in mouse, we designed a

NU

strategy to disrupt the genomic sequence in Dyrk1a exon 3, upstream from the enzymatic and nuclear localization signal sites (Figure 1A). Using the CRISPR/Cas9

MA

system, we generated two mouse lines carrying single nucleotide insertions leading to frameshifts and truncations of the protein products (Figure 1A, B). For both lines the five highest scoring predictive offtargets were located in either non-coding regions or

D

intronic segments (Supplementary Figure 1A). The absence of mutations at these loci

TE

was confirmed by sequencing (Supplementary Figure 1B). Dyrk1a+/− animals from both lines displayed significant brain morphological

AC CE P

abnormalities. As described by Fotaki and colleagues (Fotaki et al. 2002) the mesencephalic tectum was significantly smaller in both of our frameshift lines (Supplementary Figure 2A). The brain weight was significantly lighter and whereas the cerebellum size did not appear to be affected, the cerebrum size was significantly smaller in Dyrk1a+/- mice (Supplementary Figure 2B-F). Bodyweight monitoring showed a short and transient, yet significant, decrease from 2 months old to 4 months old (Supplementary table 3). In early postnatal development (P6) as well as in late adulthood (8 months old) bodyweight was not significantly different between Dyrk1a+/− and WT littermates. These observations reproducing the phenotypes of previously reported knockout mice thus indicate that our frameshift models display morphological changes characteristic of a heterozygous loss of Dyrk1a (Fotaki et al. 2002). Western blot analysis with an antibody targeting the C-terminus of DYRK1A (D1819 in Figure 1A) showed a drastic decrease in DYRK1A amount in the brain of

12

ACCEPTED MANUSCRIPT

mice from both Dyrk1a mutant lines (Figure 1C). Normalized DYRK1A amount was significantly decreased in both lines (Figure 1D; One-way ANOVA, p=4.5E-03 and p=7.4E-03 respectively). No significant differences were observed between the two lines (two-way ANOVA: “genotype” effect p=4.3E-04, “line” effect p=0.491). A second antibody targeting the N-terminus (sc-12568 in Figure 1A) failed to detect any

PT

protein in either WT or Dyrk1a+/− samples (data not shown). Taken together, these results confirm a disruption in Dyrk1a caused by a

RI

frameshift in two equivalent mutated alleles resulting in a similar drop in DYRK1A protein level and abnormal brain morphology identical to Dyrk1a knockout mice

SC

reported in the literature. Mice derived from both lines were thus considered as carrying an equivalent haploinsufficiency and cohorts used for further experiments were

NU

prepared from a mix of mice derived from both lines in equivalent proportion.

MA

Spontaneous activity and exploratory behavior are conserved in Dyrk1a+/− mice We first assessed spontaneous exploratory behavior in Dyrk1a+/− mice using the open field and elevated plus maze tasks. Travelled distance and moving speed were not

D

significantly different between Dyrk1a+/− and WT littermates in both tests

TE

(Supplementary table 3). Spontaneous activity was thus conserved in the Dyrk1a+/− mice. The time spent in the central area in the open field or in the open arms in the

AC CE P

elevated plus maze are commonly used as markers of anxiety in rodents. These parameters were not significantly different between Dyrk1a+/− mice and their WT littermates (Supplementary table 3), suggesting that anxiety level is not significantly affected in Dyrk1a+/− mice. Though the number of rearing events in the open field was significantly decreased in Dyrk1a+/− mice (p=0.0096), the number of circling and the overall behavior in the elevated plus maze did not hint to any obvious sign of persistence-like behavior (Supplementary table 3). Taken together, these results indicate that spontaneous activity and anxiety level are conserved in Dyrk1a+/− mice and that these mice do not show significant signs of ASD-like repetitive behavior.

13

ACCEPTED MANUSCRIPT

Dyrk1a+/− mice display a mild deficit in reverse spatial learning and memory We next assessed spatial learning and memory performance of Dyrk1a+/− mice using the Barnes maze task. In the habituation phase, Dyrk1a+/− mice travelled a mildly longer distance and their average speed (over the total session, including idle events) was mildly higher than their WT littermates but this did not reach the statistical

PT

significance (Supplementary Figure 3A, B). The moving speed (including only periods in which mice move faster than 3 cm/s, i.e. excluding idle events) was comparable to

RI

that of their WT littermates (Supplementary Figure 3C).

We then challenged mice to learn the position of an escape box over a four-days

SC

training course. WT and Dyrk1a+/− mice similarly showed a significant improvement in the evolution in number of errors and travelled distance to reach the target (Figure 2A-

NU

B; Two-way repeated measures ANOVA, Errors to target: “Day” effect: F(3,17)=12.74; p=9.0E-7 – “Genotype” effect: F(1,19)=0.55; p=0.459; Distance to target: “Day” effect:

MA

F(3,17)=11.99; p=2.9E-4 – “Genotype” effect: F(1,19)=2.92; p=0.092). Moving speed was significantly lower in the Dyrk1a+/− group from the second through the fourth day (Figure 2C; Two-way repeated measures ANOVA, “Genotype” effect: F(1,19)=27.94;

D

p=2.8E-4) whereas their average speed was not significantly different from their WT

TE

counterparts (Supplementary Figure 3E; Two-way repeated measures ANOVA, “Genotype” effect: F(1,19)=4.93; p=0.092).

AC CE P

In the probe test, mice from both WT and Dyrk1a+/− groups spent significantly more time investigating the target hole than the remaining holes around the board (Figure 2D; Two-way repeated measure ANOVA, “Position” effect: F(11, 9)=14.64; p=5.1E-21) without significant genotype effect (F(1,19)=1.83E-5; p=0.996). Dyrk1a+/− mice tended to make a larger number of mistakes to reach the target but this did not reach the significance level (Figure 2E; One-way ANOVA, p=0.097). In the reverse learning phase mice were trained to learn the new location of the escape box (shifted 180°) for four consecutive days. Dyrk1a+/− mice tended to show a delay in acquiring the task, making more mistakes and travelling longer distances than their WT littermates (Figure 2F-G; Two-way repeated measures ANOVA, Errors to target: “Day” effect: F(3,17)=14.03; p=2.7E-7 – “Genotype” effect: F(1,19)=4.20; p=0.044; Distance to target: “Day” effect: F(3,17)=14.29; p=2.1E-7 – “Genotype” effect: F(1,19)=9.82; p=0.002). The moving speed of Dyrk1a+/− mice was significantly

14

ACCEPTED MANUSCRIPT

slower than their WT littermates (Figure 2H; Two-way repeated measures ANOVA, “Genotype” effect: F(1,19)=49.70; p=8.9E-10) whereas their average speed and latency to reach the target were not significantly affected (Supplementary Figure 3K-L; Twoway repeated measures ANOVA, Average speed - “Genotype” effect: F(1,19)=0.33; p=0.563; Latency to reach target - “Genotype” effect: F(1,19)=1.58; p=0.213). Taken

PT

together, these results suggest a delay in the acquisition of a new paradigm in Dyrk1a+/− mice.

RI

In the reverse probe test, mice from both WT and Dyrk1a+/− groups spent significantly more time investigating the target hole than the remaining holes around the

SC

board (Figure 2I; Two-way repeated measure ANOVA, “Position” effect: F(11, 9)=28.59; p=2.1E-36) without significant genotype effect (F(1,19)=3.7E-8; p=0.999).

(Figure 2E; One-way ANOVA, p=0.042).

NU

Dyrk1a+/− mice however committed significantly more mistakes to reach the target

MA

Though changes in the neuronal cell cycle have been described in the retina of Dyrk1atm1Mla mice, the exact impact on visual acuity in these animals has not been reported (Laguna et al. 2008). In the various behavioral tasks we did not observe

D

behavior that would raise suspicions of visual acuity deficiency. Notably, in the Barnes

TE

maze task, Dyrk1a+/- mice were able to learn the position of the escape box in the first learning phase, and we noted numerous occurrences where mutant mice moved directly

AC CE P

toward the target after a short orientation phase easily identifiable on tracking data, even being able to turn around when they were facing an opposite direction at the start of the test. Visual acuity is thus unlikely to be significantly affected in our Dyrk1a+/- mice. We then used the rotarod task to investigate the locomotor coordination skills of Dyrk1a+/− mice. Mice from both the WT and Dyrk1a+/− groups showed a significant improvement over the four days training course (Figure 2K-L, Two-way repeated measures ANNOVA, “Day” effect: F(3,17)=7.90; p=9.2E-5). The performance and improvement score of Dyrk1a+/− mice was not significantly different from their WT counterparts (Figure 2L, One-way ANOVA p=0.281).

15

ACCEPTED MANUSCRIPT

Dyrk1a+/− mice show impaired sociability without signs of depression-like behavior We next investigated the social behavior of Dyrk1a+/− mice using the 3chambers and the tube “social dominance” tasks. During the habituation phase, no bias or preference for a side of the box was observed (Supplementary Figure 4A-D). During

PT

the “sociability” phase of the 3-chambers test mice from both WT and Dyrk1a+/− groups showed a significant preference for the cage with a stranger mouse (Figure 3A, B; One-

RI

way ANOVA - WT: p=5.7E-06; Dyrk1a+/−: p=1.3E-04). Dyrk1a+/− mice however spent significantly less time investigating the stranger mouse than their WT littermates

SC

(Figure 3B; One-way ANOVA p=0.012). Dyrk1a+/− showed a significant increase in total travelled distance and average speed, whereas their moving speed was comparable

NU

to that of WT mice (Figure 3C-E; One-way ANOVA p=0.025; p=0.025 and p=0.199 respectively). In the “preference for social novelty” phase of the task, no significant

MA

differences were observed between Dyrk1a+/− and WT mice (Supplementary Figure 4EJ). Even though the preference for the stranger mouse over the familiar one did not reach significance in the WT group, the significance level was reached in the Dyrk1a+/−

D

group (Supplementary Figure 4J; One-way ANOVA p=0.571 and p=2.9E-05

TE

respectively). The social memory is thus unlikely to be affected in Dyrk1a+/− mice. We investigated further the social behavior of Dyrk1a+/− mice using the tube test

AC CE P

for social dominance. Dyrk1a+/− mice lost 75% of their 52 confrontations with WT counterparts (Figure 3F; Chi-square test p=3.1E-4). We noticed that Dyrk1a+/− mice tended to retreat from the tube after nose-to-nose contact and short confrontations. The ratio of loss by retreat out of the total number of losses was indeed significantly higher in the Dyrk1a+/− group (Figure 3G, Supplementary Video 1; One-way ANOVA p=0.015).

In order to clarify whether Dyrk1a+/− mice had a stronger tendency to give up, we performed the tail suspension test, commonly used to detect depression-like behavior. No significant difference was observed between Dyrk1a+/− and WT mice (Figure 3H; Two-way repeated measures ANOVA, “Genotype” effect: F(1,19)=1.05; p=0.307). In addition, as social behavior in rodents mainly relies on the olfactory system, we used the buried food task to challenge olfaction in the Dyrk1a+/- mice. The latency to

16

ACCEPTED MANUSCRIPT

find the buried pellet was not significantly different between Dyrk1a+/- animals and their WT littermates, suggesting their olfactory performance was not significantly affected (Supplementary Figure 5). Taken together, these data indicate that although the preference for social over non-social interaction is observed in Dyrk1a+/- mice, their sociability and in particular

PT

their social approach behavior is significantly lower than that of their WT littermates. Moreover, the absence of signs of depression or easiness to give up in non-social

RI

conditions suggests that the confrontation avoidance in the tube test is mainly driven by

SC

a deficit in social behavior.

Dyrk1a+/− pups produce less ultrasonic calls and a shift in call patterns

NU

We then assessed communication skills of Dyrk1a+/− pups by measuring ultrasonic vocalizations (USV) at post-natal day 6 where isolation-induced ultrasonic

MA

vocalizations tend to peak in mice (Mosienko et al., 2015). Data was acquired from male and female pups. As males and females were evenly distributed between the groups (WT: 7 females and 6 males, Dyrk1a+/−: 6 females and 5 males) and no

TE

were pulled together.

D

significant differences were observed between males and females, data from both sexes

Automatic call detection allowed the measurement of duration, peak frequency,

AC CE P

minimum frequency, frequency modulation (i.e. the difference between the minimum and peak frequencies within each call) and amplitude for every single call (Figure 4A). The default-preset references were used in Avisoft SASLabPro software, resulting in negative values for call amplitudes. Dyrk1a+/− pups produced significantly fewer calls than their WT littermates (Figure 4B; One-way ANOVA p=0.038). In a broad quantitative analysis, calls produced by Dyrk1a+/− mice were significantly shorter, of lower amplitude and lower frequency modulation (Figure 4C-E; One-way ANOVA p=0.002, p=0.005 and p=0.016 respectively). The minimum frequency was significantly higher in the Dyrk1a+/− group whereas the peak frequency was not significantly different between Dyrk1a+/− and WT pups (Figure 4F-G; One-way ANOVA p=0.007 and p=0.707 respectively).

17

ACCEPTED MANUSCRIPT

We next classified the calls into the ten main types commonly observed in C57BL/6J pups (Figure 4H)(Scattoni et al., 2008). Dyrk1a+/− pups produced calls of the “upward” and “short” types at a significantly higher rate than their WT littermates (Figure 4H-I; One-way ANOVA p=0.047 and p=0.008 respectively). In contrast, they produced complex calls such as “complex” and “harmonic” types at significantly lower

PT

ratio than their WT littermates (Figure 4H-I; One-way ANOVA p=0.031 and p=0.003 respectively). They also produced fewer “flat” calls and tended to produce fewer “two-

RI

syllable” calls than their WT littermates (Figure 4H-I; One-way ANOVA p=0.011 and p=0.069 respectively).

SC

Dyrk1a+/− pups thus showed a significant deficit in communication by USV,

NU

overall producing fewer calls and of lower complexity.

Dyrk1a+/− mice show a mild susceptibility to seizures induced by hyperthermia

MA

but not by pentylenetetrazole or kainic acid injections As a part of human patients with mutations in DYRK1A develop febrile seizures, evolving in some cases to other forms of epilepsy (Bronicki et al., 2015; Ji et al., 2015;

D

van Bon et al., 2016), we next investigated seizures susceptibility of Dyrk1a+/− mice.

TE

Before testing hyperthermia-induced seizure susceptibility, the baseline body temperature was not significantly different between Dyrk1a+/− and WT mice (Figure

AC CE P

5A; One-way ANOVA p=0.944). During the test, the rate of body temperature elevation did not reveal significant differences between Dyrk1a+/− and WT mice (Figure 5B; Twoway repeated measures ANOVA – “Genotype” effect: F(1,19)=0.12; p=0.726). Generalized seizures however appeared at a slightly yet significantly lower body temperature in Dyrk1a+/− mice (Figure 5C; One-way ANOVA p=0.005). Kainic acid, an agonist of the ionotropic glutamate receptors induces seizures and is commonly used to model temporal lobe epilepsies (Ben-Ari et al., 1979; Lévesque et al., 2016). A 30mg/kg intraperitoneal injection of kainic acid induced generalized seizures and death within the 60 min monitoring period of all Dyrk1a+/− and WT mice. The latencies to generalized seizure and death were not significantly different between Dyrk1a+/− and their WT animals (Figure 5D, E; One-way ANOVA – Latency to GS: p=0.411; Latency to death: p=0.962).

18

ACCEPTED MANUSCRIPT

Pentylenetetrazol (PTZ), an inhibitor of the γ-aminobutyric acid type A (GABAA) receptor, is the most widely used model of acute chemically induced seizures (Huang et al., 2001; Löscher, 2011). PTZ injection (50mg/kg) led to a generalized seizure in 6 out of 10 mice in both Dyrk1a+/− and WT groups (Figure 5F). The latency to generalized seizure and seizure severity score were not significantly different between

PT

Dyrk1a+/− and WT mice (Figure 5F, G; One-way ANOVA – Latency to GS: p=0.543; Severity score: p=0.764).

RI

No abnormal epileptic patterns were seen on basal electrocorticogram recordings in Dyrk1a+/− mice and local field potential traces showed similar properties in waking

SC

state as well as REM (Rapid Eye Movement) and non-REM sleep states (Supplementary Figure 6).

NU

Taken together, these results indicate a mildly yet significantly increased susceptibility to hyperthermia induced seizures in Dyrk1a+/− mice without significant

MA

differences in susceptibilities to kainic acid and PTZ.

Discussion

D

The multiplication of case reports over the past two years have led to the

TE

description of a syndromic phenotype consistently observed in human patients carrying mutations of DYRK1A (Bronicki et al., 2015; Ji et al., 2015; van Bon et al., 2016).

AC CE P

Frameshifts or nonsense mutations resulting in DYRK1A truncation have been identified in about two-third of these patients. In the present work we developed a mouse model in which a frameshift by single nucleotide insertion resulted in a truncation of DYRK1A protein upstream from its nuclear localization signal and kinase activity site. This led to a total amount of DYRK1A in the adult brain to be drastically decreased in Dyrk1a+/− mice. This model is thus genetically close to the most deleterious mutations described in human case reports. According to previous reports, Dyrk1a+/− mice show a delayed acquisition of early post-natal landmarks (Fotaki et al., 2002). Despite a decrease in grip strength and significantly shorter strides, locomotor and neuromuscular performance is conserved in Dyrk1a+/− adult mice (Arqué et al., 2008; Arqué et al., 2009; Fotaki et al., 2002; Fotaki et al., 2004). In addition, signs of hypoactivity (decreased rearing and travelled

19

ACCEPTED MANUSCRIPT

distances) have been described (Fotaki et al., 2004). Whereas the working memory in the novel object recognition task is significantly deficient in Dyrk1a+/− mice, their spatial learning memory is only very mildly affected, not reaching the significance level in the Morris water maze nor in the radial water maze (Arqué et al., 2008; Arqué et al., 2009). Consistent with these observations, the performance of our Dyrk1a+/− mice in the

PT

rotarod task was conserved whereas their moving speed in the Barnes maze task as well as their rearing count in the open field were significantly decreased. However, looking

RI

at the travelled distance in the different behavior tasks, we did not observe significant signs of hypoactivity in our Dyrk1a+/- mice and their latency to reach the target in the

SC

Barnes maze was not significantly affected. Moreover, our work provides new insights in Dyrk1a+/− behavior, refining the cognitive phenotypes in a reverse learning strategy

NU

and introducing for the first time sociability and communication deficits. Indeed, in our present study the most remarkable deficit of Dyrk1a+/− mice in the

MA

Barnes maze task occurred in the reverse learning phase, whereas their performance in the primary learning phase was normal. This suggests that cognitive flexibility, rather than hippocampal-dependent spatial memory, is affected in Dyrk1a+/− mice. Reversal

D

learning is an executive function that relies on the prefrontal cortex (McAlonan and

TE

Brown, 2003). Though no functional study of the prefrontal cortex in Dyrk1a knockout mice is available so far, overexpression of Dyrk1a increases excitatory post-synaptic

AC CE P

currents at the local level ex-vivo (Thomazeau et al., 2014) but decreases overall neuronal activity, affecting gamma waves in-vivo (Ruiz-Mejias et al., 2016). The outcome of the loss of function in Dyrk1a+/− mice is difficult to predict, as it does not always lead to a phenotype opposite to the overexpression models. For instance, both overexpression and silencing of Dyrk1a lead to a similar decrease in neuronal dendrites and spines formation (Dang et al., 2017). Evidence however strongly suggests a critical role of Dyrk1a in the prefrontal cortex and further work will help understanding how mutations affect the performance in cognitive flexibility. ASD is characterized by social behavior impairments, often associated with repetitive behavior. In mouse models, various combinations of autistic-like features have been identified so far. For instance, whereas sociability is significantly decreased in Mecp2 knockout (Rett syndrome model), Mef2c knockout and BTBR (inbred autism model) mice, the Fmr1 knockout mice (Fragile X model) shows normal sociability

20

ACCEPTED MANUSCRIPT

(Harrington et al., 2016; McFarlane et al., 2008; Moretti et al., 2005; Spencer et al., 2005). In the tube test for social dominance Fmr1 knockout mice lose significantly more than their WT littermates, whereas Mecp2 knockout mice have an even win/lose ratio and Mef2c knockout mice have a significantly higher win rate. In addition, among these models only Mef2c knockout and BTBR mice display repetitive behavior in the form of

PT

excessive grooming. In the present study we investigated for the first time social behavior in Dyrk1a+/− model. Dyrk1a+/− mice showed a significantly decreased

RI

sociability and escaped confrontations in the tube test. No sign of repetitive behavior was however seen in Dyrk1a+/− mice.

SC

Severe speech acquisition delay is found in all patients with mutations in DYRK1A (Bronicki et al., 2015; Ji et al., 2015; van Bon et al., 2016). Mimicking this

NU

strong developmental delay, we identified for the first time a communication deficit in juvenile Dyrk1a+/− mice. Dyrk1a+/− pups indeed produce fewer vocalizations with a shift

MA

toward more “simple” types of calls. A variety of changes in isolation-induced ultrasonic vocalization properties have also been identified in several other models of ASD. Shank1, Tph2 and Mef2c knockout pups produce fewer calls than their WT

D

littermates whereas Fmr1 deficient pups produce a normal number of vocalizations

TE

(Harrington et al., 2016; Mosienko et al., 2015; Roy et al., 2012; Wöhr et al., 2011). No significant shift in call types has been described in these models, even though minor

AC CE P

variations in peak frequency or frequency modulations are seen in Shank1, Tph2 and Fmr1 mutant pups (Mosienko et al., 2015; Roy et al., 2012; Wöhr et al., 2011). A shift in the types of vocalizations has been reported in BTBR mice producing more calls and of higher complexity than C57BL/6J pups (Scattoni et al., 2008). Dyrk1a+/− pups is thus the first model to our knowledge to display a decrease in call number associated with a shift toward a lower complexity that could reflect a deficit in communication skills. Social behavior is a complex function for which there is so far no consensus in terms of candidate circuitry or neurotransmission system. Oxytocin signaling in the paraventricular thalamus (Peñagarikano et al., 2015), NMDA receptors in the amygdala (Schoch et al., 2017), GABAA in the forebrain (Han et al., 2012) or more recently glutamatergic transmission in the ventral tegmental area (Krishnan et al., 2017) have all been raised as leading causes to ASD-like sociability deficit. Dyrk1a is expressed in a large amount of neurons throughout the brain and the neurons subtypes expressing

21

ACCEPTED MANUSCRIPT

DYRK1A have not been described so far (Hämmerle et al., 2008). Further extensive work is thus required to determine how Dyrk1a mutations affect social behavior, using region specific and/or neuron subtype specific conditional knockouts. Nearly 60% of the reported patients with a mutation in DYRK1A suffer febrile seizures in infancy, sometimes evolving to other types of epileptic seizures later in life

PT

(Bronicki et al., 2015; Ji et al., 2015; van Bon et al., 2016). Dyrk1a+/− mice were slightly yet significantly more susceptible to hyperthermia-induced seizures. A temperature

RI

elevation in the brain negatively affects inhibitory transmission more than excitatory (Qu et al., 2007). The principal mechanisms identified so far involve both presynaptic

SC

(GABA release and reuptake) and postsynaptic components (GABA receptor trafficking, availability and function) resulting in an overall decrease in GABAergic transmission

NU

(Kang et al., 2006; Qu and Leung, 2008). This phenomenon would be amplified by mutations affecting inhibitory neurons. Indeed, a number of cases of febrile seizures in

MA

human are linked to mutations in the voltage activated sodium channels SCN1A and SCN1B, resulting in a decrease in GABAergic transmission (Steinlein, 2014; Yamakawa, 2016). Dyrk1a may generate susceptibilities to seizures by affecting neurotransmission

D

through its target proteins involved in vesicle trafficking and recycling as well as

TE

regulation of receptors' availability at the cell surface (Chen et al., 2014; Grau et al., 2014; Ori-McKenney et al., 2016; Park et al., 2012).

AC CE P

In conclusion we report for the first time ASD-like sociability and communication deficits in Dyrk1a+/− mice. We also provide a refinement of their cognitive skills pointing to a cognitive flexibility impairment that adds up to the already known working memory deficiency. In addition, we identified a mild susceptibility to hyperthermia-induced seizures in Dyrk1a+/− mice. Autism, intellectual disability, severe speech delay and febrile seizures have been described in every patient with DYRK1A mutations reported so far. Our work shows that Dyrk1a+/− mice thus recapitulate these phenotypes and should help understand the pathology of neurodevelopmental disorders caused by DYRK1A mutations and to the development of therapeutic approaches.

22

ACCEPTED MANUSCRIPT

Supplementary Information Supplementary

information

includes

six

supplementary

figures,

three

supplementary tables and one supplementary video.

Acknowledgements We are grateful to RIKEN Brain Science Institute (BSI) Research Resource

PT

Center for the support to generate the Dyrk1a+/− model, taking care of the mice and maintaining behavior-testing equipment. We also thank all the members of the

SC

work was supported in part by RIKEN-BSI [to K.Y.].

RI

laboratory for Neurogenetics in RIKEN-BSI for helpful comments and advice. This

Conflict of Interest

AC CE P

TE

D

MA

NU

None declared.

23

ACCEPTED MANUSCRIPT

References 1. Ahn, K.J., Jeong, H.K., Choi, H.S., Ryoo, S.R., Kim, Y.J., Goo, J.S., Choi, S.Y., Han, J.S., Ha, I., Song, W.J., 2006. DYRK1A BAC transgenic mice show altered synaptic plasticity with learning and memory defects. Neurobiol. Dis. 3, 463-472.

PT

2. Altafaj, X., Dierssen, M., Baamonde, C., Martí, E., Visa, J., Guimerà, J., Oset, M., González, J.R., Flórez, J., Fillat, C., Estivill, X., 2001. Neurodevelopmental

RI

delay, motor abnormalities and cognitive deficits in transgenic mice overexpressing Dyrk1A (minibrain), a murine model of Down's syndrome. Hum.

SC

Mol. Genet. 10, 1915-1923.

3. Arqué, G., Fotaki, V., Fernández, D., Martínez de Lagrán, M., Arbonés, M.L.,

NU

Dierssen, M., 2008. Impaired spatial learning strategies and novel object recognition in mice haploinsufficient for the dual specificity tyrosine-regulated

MA

kinase-1A (Dyrk1A). PLoS One 3, e2575.

4. Arqué, G., de Lagrán, M.M., Arbonés, M.L., Dierssen, M., 2009. Age-associated motor and visuo-spatial learning phenotype in Dyrk1A heterozygous mutant

D

mice. Neurobiol. Dis. 36, 312-319.

TE

5. Becker, W., Weber, Y., Wetzel, K., Eirmbter, K., Tejedor, F.J., Joost, H.G., 1998. Sequence characteristics, subcellular localization, and substrate specificity

AC CE P

of DYRK-related kinases, a novel family of dual specificity protein kinases. Biol. Chem. 273, 25893-25902. 6. Ben-Ari, Y., Lagowska, J., Tremblay, E., Le Gal La Salle, G., 1979. A new model of focal status epilepticus: intra-amygdaloid application of kainic acid elicits repetitive secondarily generalized convulsive seizures. Brain Res. 163, 176-179.

7. Benavides-Piccione, R., Dierssen, M., Ballesteros-Yáñez, I., Martínez de Lagrán, M., Arbonés, M.L., Fotaki, V., DeFelipe, J., Elston, G.N., 2005. Alterations in the phenotype of neocortical pyramidal cells in the Dyrk1A+/− mouse. Neurobiol. Dis. 20, 115-122. 8. Bronicki, L.M., Redin, C., Drunat, S., Piton, A., Lyons, M., Passemard, S., Baumann, C., Faivre, L., Thevenon, J., Rivière, J.B., Isidor, B., Gan, G.,

24

ACCEPTED MANUSCRIPT

Francannet, C., Willems, M., Gunel, M., Jones, J.R., Gleeson, J.G., Mandel, J.L., Stevenson, R.E., Friez, M.J., Aylsworth, A.S., 2015. Ten new cases further delineate the syndromic intellectual disability phenotype caused by mutations in DYRK1A. Eur. J. Hum. Genet. 23, 1482-1487. 9. Chen, C.K., Bregere, C., Paluch, J., Lu, J.F., Dickman, D.K., Chang, K.T., 2014.

PT

Activity-dependent facilitation of Synaptojanin and synaptic vesicle recycling by the Minibrain kinase. Nat. Commun. 5, 4246.

RI

10. Cong, L., Ran, F.A., Cox, D., Lin, S., Barretto, R., Habib, N., Hsu, P.D., Wu, X., Jiang, W., Marraffini, L.A., Zhang, F., 2013. Multiplex genome engineering

SC

using CRISPR/Cas systems. Science 339, 819-823.

11. Courcet, J.B., Faivre, L., Malzac, P., Masurel-Paulet, A., Lopez, E., Callier, P.,

NU

Lambert, L., Lemesle, M., Thevenon, J., Gigot, N., Duplomb, L., Ragon, C., Marle, N., Mosca-Boidron, A.L., Huet, F., Philippe, C., Moncla, A., Thauvin-

MA

Robinet, C., 2012. The DYRK1A gene is a cause of syndromic intellectual disability with severe microcephaly and epilepsy. J. Med. Genet. 49, 731-736. 12. Dang, T., Duan, W.Y., Yu, B., Tong, D.L., Cheng, C., Zhang, Y.F., Wu, W., Ye,

D

K., Zhang, W.X., Wu, M., Wu, B.B., An, Y., Qiu, Z.L., Wu, B.L., 2017.

TE

Autism-associated Dyrk1a truncation mutants impair neuronal dendritic and spine growth and interfere with postnatal cortical development. Mol. Psychiatry

AC CE P

Epub ahead of print.

13. Deciphering Developmental Disorders Study, 2015. Large-scale discovery of novel genetic causes of developmental disorders. Nature 519, 223-228. 14. Deciphering Developmental Disorders Study, 2017. Prevalence and architecture of de novo mutations in developmental disorders. Nature 542, 433-438. 15. Duchon, A., Herault, Y., 2016. DYRK1A, a Dosage-Sensitive Gene Involved in Neurodevelopmental Disorders, Is a Target for Drug Development in Down Syndrome. Front. Behav. Neurosci. 10, 104. 16. Fotaki, V., Dierssen, M., Alcántara, S., Martínez, S., Martí, E., Casas, C., Visa, J., Soriano, E., Estivill, X., Arbonés, M.L., 2002. Dyrk1A haploinsufficiency affects viability and causes developmental delay and abnormal brain morphology in mice. Mol. Cell. Biol. 22, 6636-6647.

25

ACCEPTED MANUSCRIPT

17. Fotaki, V., Martínez De Lagrán, M., Estivill, X., Arbonés, M., Dierssen, M., 2004. Haploinsufficiency of Dyrk1A in mice leads to specific alterations in the development and regulation of motor activity. Behav. Neurosci. 118, 815-821. 18. Grau, C., Arató, K., Fernández-Fernández, J.M., Valderrama, A., Sindreu, C., Fillat, C., Ferrer, I., de la Luna, S., Altafaj, X., 2014. DYRK1A-mediated

PT

phosphorylation of GluN2A at Ser(1048) regulates the surface expression and channel activity of GluN1/GluN2A receptors. Front. Cell. Neurosci. 8, 331.

RI

19. Hämmerle, B., Elizalde, C., Tejedor, F.J., 2008. The spatio-temporal and subcellular expression of the candidate Down syndrome gene Mnb/Dyrk1A in

SC

the developing mouse brain suggests distinct sequential roles in neuronal development. Eur. J. Neurosci. 27, 1061-1074.

NU

20. Hämmerle, B., Ulin, E., Guimera, J., Becker, W., Guillemot, F., Tejedor, F.J., 2011. Transient expression of Mnb/Dyrk1a couples cell cycle exit and

MA

differentiation of neuronal precursors by inducing p27KIP1 expression and suppressing NOTCH signaling. Development 138, 2543-2554. 21. Han, S., Tai, C., Westenbroek, R.E., Yu, F.H., Cheah, C.S., Potter, G.B.,

D

Rubenstein, J.L., Scheuer, T., de la Iglesia, H.O., Catterall, W.A., 2012.

TE

Autistic-like behaviour in Scn1a+/− mice and rescue by enhanced GABAmediated neurotransmission. Nature 489, 385-390.

AC CE P

22. Harrington, A.J., Raissi, A., Rajkovich, K., Berto, S., Kumar, J., Molinaro, G., Raduazzo, J., Guo, Y., Loerwald, K., Konopka, G., Huber, K.M., Cowan, C.W., 2016. MEF2C regulates cortical inhibitory and excitatory synapses and behaviors

relevant

to

neurodevelopmental

disorders.

Elife

doi:

10.7554/eLife.20059.

23. Huang, R.Q., Bell-Horner, C.L., Dibas, M.I., Covey, D.F., Drewe, J.A., Dillon, G.H., 2001. Pentylenetetrazole-induced inhibition of recombinant gammaaminobutyric acid type A (GABA(A)) receptors: mechanism and site of action. J. Pharmacol. Exp. Ther. 298, 986-995. 24. Ji, J., Lee, H., Argiropoulos, B., Dorrani, N., Mann, J., Martinez-Agosto, J.A., Gomez-Ospina, N., Gallant, N., Bernstein, J.A., Hudgins, L., Slattery, L., Isidor, B., Le Caignec, C., David, A., Obersztyn, E., Wiśniowiecka-Kowalnik, B., Fox, M., Deignan, J.L., Vilain, E., Hendricks, E., Horton Harr, M., Noon, S.E.,

26

ACCEPTED MANUSCRIPT

Jackson, J.R., Wilkens, A., Mirzaa, G., Salamon, N., Abramson, J., Zackai, E.H., Krantz, I., Innes, A.M., Nelson, S.F., Grody, W.W., 2015. DYRK1A haploinsufficiency causes a new recognizable syndrome with microcephaly, intellectual disability, speech impairment, and distinct facies. Eur. J. Hum. Genet. 23, 1473-1481.

PT

25. Kaczmarski, W., Barua, M., Mazur-Kolecka, B., Frackowiak, J., Dowjat, W., Mehta, P., Bolton, D., Hwang, Y.W., Rabe, A., Albertini, G., Wegiel, J., 2014.

RI

Intracellular distribution of differentially phosphorylated dual-specificity tyrosine phosphorylation-regulated kinase 1A (DYRK1A). J. Neurosci. Res. 92,

SC

162-173.

26. Kang, J.Q., Shen, W., Macdonald, R.L., 2006. Why does fever trigger febrile

NU

seizures? GABAA receptor gamma2 subunit mutations associated with idiopathic generalized epilepsies have temperature-dependent trafficking

MA

deficiencies. J. Neurosci. 26, 2590-2597.

27. Krishnan, V., Stoppel, D.C., Nong, Y., Johnson, M.A., Nadler, M.J., Ozkaynak, E., Teng, B.L., Nagakura, I., Mohammad, F., Silva, M.A., Peterson, S., Cruz,

D

T.J., Kasper, E.M., Arnaout, R., Anderson, M.P., 2017. Autism gene Ube3a and

TE

seizures impair sociability by repressing VTA Cbln1. Nature 543, 507-512. 28. Krumm, N., O'Roak, B.J., Säendure, J., Eichler, E.E., 2014. A de novo

AC CE P

convergence of autism genetics and molecular neuroscience. Trends Neurosci. 37, 95-105.

29. Laguna, A., Aranda, S., Barallobre, M.J., Barhoum, R., Fernández, E., Fotaki, V., Delabar, J.M., de la Luna, S., de la Villa, P., Arbonés, M.L., 2008. The protein kinase DYRK1A regulates caspase-9-mediated apoptosis during retina development. Dev. Cell 15, 841-853. 30. Lévesque, M., Avoli, M., Bernard, C., 2016. Animal models of temporal lobe epilepsy following systemic chemoconvulsant administration. J Neurosci Methods 260, 45-52. 31. Löscher, W., 2011. Critical review of current animal models of seizures and epilepsy used in the discovery and development of new antiepileptic drugs. Seizure 20, 359-368.

27

ACCEPTED MANUSCRIPT

32. Martí, E., Altafaj, X., Dierssen, M., de la Luna, S., Fotaki, V., Alvarez, M., Pérez-Riba, M., Ferrer, I., Estivill, X., 2003. Dyrk1A expression pattern supports specific roles of this kinase in the adult central nervous system. Brain Res. 964, 250-263. 33. Martinez de Lagran, M., Benavides-Piccione, R., Ballesteros-Yañez, I., Calvo,

PT

M., Morales, M., Fillat, C., Defelipe, J., Ramakers, G.J., Dierssen, M., 2012. Dyrk1A influences neuronal morphogenesis through regulation of cytoskeletal

RI

dynamics in mammalian cortical neurons. Cereb. Cortex 22, 2867-2877. 34. McAlonan, K., Brown, V.J., 2003. Orbital prefrontal cortex mediates reversal

SC

learning and not attentional set shifting in the rat. Behav. Brain Res. 146, 97-103. 35. McFarlane, H.G., Kusek, G.K., Yang, M., Phoenix, J.L., Bolivar, V.J., Crawley,

NU

J.N., 2008. Autism-like behavioral phenotypes in BTBR T+tf/J mice. Genes Brain Behav. 7, 152-163.

MA

36. Møller, R.S., Kübart, S., Hoeltzenbein, M., Heye, B., Vogel, I., Hansen, C.P., Menzel, C., Ullmann, R., Tommerup, N., Ropers, H.H., Tümer, Z., Kalscheuer, V.M., 2008. Truncation of the Down syndrome candidate gene DYRK1A in two

D

unrelated patients with microcephaly. Am. J. Hum. Genet. 82, 1165-1170.

TE

37. Moretti, P., Bouwknecht, J.A., Teague, R., Paylor, R., Zoghbi, H.Y., 2005. Abnormalities of social interactions and home-cage behavior in a mouse model

AC CE P

of Rett syndrome. Hum. Mol. Genet. 14, 205-220. 38. Mosienko, V., Beis, D., Alenina, N., Wöhr, M., 2015. Reduced isolationinduced pup ultrasonic communication in mouse pups lacking brain serotonin. Mol. Autism 6, 13.

39. O'Roak, B.J., Vives, L., Fu, W., Egertson, J.D., Stanaway, I.B., Phelps, I.G., Carvill, G., Kumar, A., Lee, C., Ankenman, K., Munson, J., Hiatt, J.B., Turner, E.H., Levy, R., O'Day, D.R., Krumm, N., Coe, B.P., Martin, B.K., Borenstein, E., Nickerson, D.A., Mefford, H.C., Doherty, D., Akey, J.M., Bernier, R., Eichler, E.E., Shendure, J., 2012. Multiplex targeted sequencing identifies recurrently mutated genes in autism spectrum disorders. Science 338, 1619-1622. 40. Ori-McKenney, K.M., McKenney, R.J., Huang, H.H., Li, T., Meltzer, S., Jan, L.Y., Vale, R.D., Wiita, A.P., Jan, Y.N., 2016. Phosphorylation of β-Tubulin by

28

ACCEPTED MANUSCRIPT

the Down Syndrome Kinase, Minibrain/DYRK1a, Regulates Microtubule Dynamics and Dendrite Morphogenesis. Neuron 90, 551-563. 41. Park, J.H., Jung, M.S., Kim, Y.S., Song, W.J., Chung, S.H., 2012. Phosphorylation of Munc18-1 by Dyrk1A regulates its interaction with Syntaxin 1 and X11α. J. Neurochem. 122, 1081-1091.

PT

42. Paxinos, G., Franklin, K., 2001. The mouse brain in stereotaxic coordinates. 2nd ed. San Diego Academic Press.

RI

43. Peñagarikano, O., Lázaro, M.T., Lu, X.H., Gordon, A., Dong, H., Lam, H.A.,

SC

Peles, E., Maidment, N.T., Murphy, N.P., Yang, X.W., Golshani, P., Geschwind, D.H., 2015. Exogenous and evoked oxytocin restores social behavior in the Cntnap2 mouse model of autism. Sci. Transl. Med. 7, 271ra8.

NU

44. Qu, L., Liu, X., Wu, C., Leung, L.S., 2007. Hyperthermia decreases GABAergic synaptic transmission in hippocampal neurons of immature rats. Neurobiol. Dis.

MA

27, 320-327.

45. Qu, L., Leung, L.S., 2008. Mechanisms of hyperthermia-induced depression of GABAergic synaptic transmission in the immature rat hippocampus. J.

D

Neurochem. 106, 2158-2169.

TE

46. Roy, S., Watkins, N., Heck, D., 2012. Comprehensive analysis of ultrasonic vocalizations in a mouse model of fragile X syndrome reveals limited, call type

AC CE P

specific deficits. PLoS One 7, e44816. 47. Ruiz-Mejias, M., Martinez de Lagran, M., Mattia, M., Castano-Prat, P., PerezMendez, L., Ciria-Suarez, L., Gener, T., Sancristobal, B., García-Ojalvo, J., Gruart, A., Delgado-García, J.M., Sanchez-Vives, M.V., Dierssen, M., 2016. Overexpression of Dyrk1A, a Down Syndrome Candidate, Decreases Excitability and Impairs Gamma Oscillations in the Prefrontal Cortex. J. Neurosci. 36, 3648-3659. 48. Scattoni, M.L., Gandhy, S.U., Ricceri, L., Crawley, J.N., 2008. Unusual repertoire of vocalizations in the BTBR T+tf/J mouse model of autism. PLoS One 3, e3067. 49. Schoch, H., Kreibich, A.S., Ferri, S.L., White, R.S., Bohorquez, D., Banerjee, A., Port, R.G., Dow, H.C., Cordero, L., Pallathra, A.A., Kim, H., Li, H., Bilker,

29

ACCEPTED MANUSCRIPT

W.B., Hirano, S., Schultz, R.T., Borgmann-Winter, K., Hahn, C.G., Feldmeyer, D., Carlson, G.C., Abel, T., Brodkin, E.S., 2017. Sociability Deficits and Altered Amygdala Circuits in Mice Lacking Pcdh10, an Autism Associated Gene. Biol. Psychiatry 81, 193-202. 50. Soppa, U., Schumacher, J., Florencio Ortiz, V., Pasqualon, T., Tejedor, F.J.,

PT

Becker, W., 2014. The Down syndrome-related protein kinase DYRK1A phosphorylates p27(Kip1) and Cyclin D1 and induces cell cycle exit and

RI

neuronal differentiation. Cell Cycle 13, 2084-2100.

51. Souchet, B., Guedj, F., Sahún, I., Duchon, A., Daubigney, F., Badel, A.,

SC

Yanagawa, Y., Barallobre, M.J., Dierssen, M., Yu, E., Herault, Y., Arbones, M., Janel, N., Créau, N., Delabar, J.M., 2014. Excitation/inhibition balance and

NU

learning are modified by Dyrk1a gene dosage. Neurobiol. Dis. 69, 65-75. 52. Spencer, C.M., Alekseyenko, O., Serysheva, E., Yuva-Paylor, L.A., Paylor, R.,

MA

2005. Altered anxiety-related and social behaviors in the Fmr1 knockout mouse model of fragile X syndrome. Genes Brain Behav. 4, 420-430. 53. Steinlein, O.K., 2014. Mechanisms underlying epilepsies associated with sodium

D

channel mutations. Prog. Brain Res. 213, 97-111.

TE

54. Thomazeau, A., Lassalle, O., Iafrati, J., Souchet, B., Guedj, F., Janel, N., Chavis, P., Delabar, J., Manzoni, O.J., 2014. Prefrontal deficits in a murine model

AC CE P

overexpressing the down syndrome candidate gene dyrk1a. J. Neurosci. 34, 1138-1147.

55. van Bon, B.W., Hoischen, A., Hehir-Kwa, J., de Brouwer, A.P., Ruivenkamp, C., Gijsbers, A.C., Marcelis, C.L., de Leeuw, N., Veltman, J.A., Brunner, H.G., de Vries, B.B., 2011. Intragenic deletion in DYRK1A leads to mental retardation and primary microcephaly. Clin. Genet. 79, 296-399. 56. van Bon, B.W., Coe, B.P., Bernier, R., Green, C., Gerdts, J., Witherspoon, K., Kleefstra, T., Willemsen, M.H., Kumar, R., Bosco, P., Fichera, M., Li, D., Amaral, D., Cristofoli, F., Peeters, H., Haan, E., Romano, C., Mefford, H.C., Scheffer, I., Gecz, J., de Vries, B.B., Eichler, E.E., 2016. Disruptive de novo mutations of DYRK1A lead to a syndromic form of autism and ID. Mol. Psychiatry 21, 126-132.

30

ACCEPTED MANUSCRIPT

57. Wang, H., Yang, H., Shivalila, C.S., Dawlaty, M.M., Cheng, A.W., Zhang, F., Jaenisch, R., 2013. One-step generation of mice carrying mutations in multiple genes by CRISPR/Cas-mediated genome engineering. Cell 153, 910-918. 58. Wöhr, M., Roullet, F.I., Hung, A.Y., Sheng, M., Crawley, J.N., 2011. Communication impairments in mice lacking Shank1: reduced levels of

PT

ultrasonic vocalizations and scent marking behavior. PLoS One 6: e20631. 59. Yamakawa, K., 2016. Mutations of Voltage-Gated Sodium Channel Genes

RI

SCN1A and SCN2A in Epilepsy, Intellectual Disability, and Autism. In Neuronal and Synaptic Dysfunction in Autism Spectrum Disorder and

AC CE P

TE

D

MA

NU

SC

Intellectual Disability. Academic Press, San Diego, pp 233-251.

31

ACCEPTED MANUSCRIPT

Legend to figures Figure 1. Generation of mice with nucleotide insertions in Dyrk1a gene causing premature terminations of translation. (A) The mouse Dyrk1a has highly conserved DNA and protein sequences with its human ortholog. The exon 3 was targeted using two guide RNAs in order to disrupt DYRK1A protein upstream from its nuclear

PT

localization signal (NLS) and Serine-Threonine kinase site (STK). Antibodies targeting the C-terminus (D1819) and N-terminus (sc-12568) were used to detect DYRK1A

RI

protein. (B) Two lines were obtained by CRISPR/Cas9 recombination. DNA sequencing confirmed that both lines carried a single nucleotide insertion and one of them also

SC

carried a conversion (c.134insC and c.152insC+A>C). The resulting frameshifts lead to truncations of the protein in both cases. (C) Western blot using brain samples from 1.5

NU

months old mice confirmed a loss of DYRK1A protein in both lines. (D) DYRK1A normalized protein level was significantly decreased in both Dyrk1a+/− lines compared

MA

to their respective WT littermates. Values are expressed as mean ± standard error of the mean. Statistical significance was assessed using one-way ANOVA with significance

D

set at (**) p < 0.01.

TE

Figure 2. Impaired reversal learning and conserved locomotor coordination in Dyrk1a+/− mice. (A-J) Barnes maze test. In the primary learning phase, errors to target

AC CE P

(A) and distance to target (B) are not significantly different between WT and Dyrk1a+/− mice, whereas moving speed (C) is significantly decreased in Dyrk1a+/− mice. In the probe test, the investigation time per hole (D) did not show significant differences between Dyrk1a+/− and WT mice, even though Dyrk1a+/− mice tended to commit more mistakes to reach the target for the first time (E). In the reverse learning phase, errors to target (F) and distance to target (G) tended to be increased in Dyrk1a+/− mice but this did not reach the significance level. The moving speed was significantly slower in Dyrk1a+/− for the whole training phase (H). In the reverse probe test, Dyrk1a+/− did not show significant differences in the time spent around the different holes (I) but committed significantly more mistakes to reach the hole for the first time (J). (K-L) Rotarod test. The performance of Dyrk1a+/− mice was not significantly different from their WT littermates (K) and their improvement over the course of four days was

32

ACCEPTED MANUSCRIPT

comparable to their WT littermates (L). Values are expressed as mean ± standard error of the mean. Statistical significance was assessed using two-way repeated measures ANOVA followed by day-by-day one-way ANOVA post-hoc test (A-D, F-I and K) or one-way ANOVA (E, J and L) with significance set at (*) p<0.05 and (**) p < 0.01.

PT

Figure 3. Decreased sociability without depression-like behavior in Dyrk1a+/− mice. (A-E) 3-chambers task. In the sociability phase, mount graphs show a cumulative

RI

occupancy targeted around the stranger cage in WT mice and more spread in Dyrk1a+/− mice (A). Investigation time spent on the cage with a stranger mouse was significantly

SC

shorter in Dyrk1a+/− mice than in their WT littermates, whereas the time spent on the empty cage was not significantly different (B). Travelled distance (C) and average

NU

speed (D) were significantly increased in Dyrk1a+/− mice, whereas their moving speed (E) was comparable to their WT littermates. (F,G) Tube test. The win rate of Dyrk1a+/−

MA

mice was significantly lower than their WT littermates (F). The ratio of loss by retreat was significantly higher in the Dyrk1a+/− group (G). (H) Tail suspension test. The percentage of immobility of Dyrk1a+/− mice over the 10 minutes of the test was not

D

significantly different from their WT littermates. Values are expressed as mean ±

TE

standard error of the mean. Statistical significance was assessed using one-way ANOVA (B-E, G), Chi-square test (F) and two-way repeated measures ANOVA (H)

AC CE P

with significance set at (*) p<0.05 and (***) p < 0.001.

Figure 4. Reduced ultrasonic vocalizations with shift toward lower complexity in Dyrk1a+/− pups. (A) Isolation-induced ultrasonic vocalizations recorded at postnatal day 6 were analyzed quantitatively and qualitatively using frequency, duration and amplitude parameters. Dyrk1a+/− pups produced significantly fewer calls than their WT littermates (B). These calls were of shorter duration (C) and lower amplitude (D). The frequency modulation (E) was significantly decreased in calls produced by Dyrk1a+/− pups. Peak frequency (F) of ultrasonic vocalizations was not significantly different between Dyrk1a+/− and WT pups, whereas the minimum frequency (G) was significantly higher in the Dyrk1a+/− group. All ultrasonic vocalizations were sorted into 10 categories (H) and the proportion of calls per category (I) showed a significant shift toward less complex calls in Dyrk1a+/− pups. They produced significantly more “upward”

33

ACCEPTED MANUSCRIPT

and “short” calls but significantly fewer “flat”, “complex” and “harmonic” calls than their WT littermates. They also tended to produce fewer “two-syllable” calls but this did not reach the significance level. Values are expressed as mean ± standard error of the mean. Statistical significance was assessed using one-way ANOVA with significance

PT

set at (*) p<0.05, (**) p<0.01 and (***) p < 0.001. Figure 5. Increased susceptibility to hyperthermia induced seizures in Dyrk1a+/−

RI

mice. Basal body temperature (A) was not significantly different between Dyrk1a+/− and WT mice, and their body temperature regulation (B) while environment temperature

SC

was raised was comparable. (C) Generalized seizures (GS) however occurred at a significantly lower body temperature in the Dyrk1a+/− group. After 30mg/kg injection of

NU

kainic acid the latencies to GS (D) and death (E) were not significantly different between Dyrk1a+/− and WT mice. After 50mg/kg injection of pentylenetetrazole (PTZ)

MA

the latency to GS (F) and seizure severity score (G) were not significantly different between Dyrk1a+/− and WT mice. Severity score for PTZ-induced seizures were given on a 1 to 5 scale as follows: 1) myoclonic, 2) tonic, 3) generalized seizure, 4) full body

D

and limbs extension, 5) death. Values are expressed as mean ± standard error of the

TE

mean. Statistical significance was assessed using two-way repeated measures ANOVA for body temperature regulation in (B) and one-way ANOVA for every other parameter,

AC CE P

with significance set at (*) p<0.05 and (**) p<0.01.

34

ACCEPTED MANUSCRIPT

AC CE P

TE

D

MA

NU

SC

RI

PT

Fig. 1

35

ACCEPTED MANUSCRIPT

AC CE P

TE

D

MA

NU

SC

RI

PT

Fig. 2

36

ACCEPTED MANUSCRIPT

AC CE P

TE

D

MA

NU

SC

RI

PT

Fig. 3

37

ACCEPTED MANUSCRIPT

AC CE P

TE

D

MA

NU

SC

RI

PT

Fig. 4

38

ACCEPTED MANUSCRIPT

AC CE P

TE

D

MA

NU

SC

RI

PT

Fig. 5

39

ACCEPTED MANUSCRIPT

Highlights:

TE

D

MA

NU

SC

RI

PT

Dyrk1a haploinsufficient mice show cognitive flexibility deficits Dyrk1a+/− mice show autistic like sociability and communication deficits Dyrk1a+/− mice are sensitive to hyperthermia-induced seizures Phenotypes in mice mimic the human DYRK1A mutation syndrome

AC CE P

   

40