- Open Access
Krüppel-like factor 7 deficiency causes autistic-like behavior in mice via regulating Clock gene
Cell & Bioscience volume 12, Article number: 166 (2022)
Krüppel-like factor 7 (klf7), a transcription factor in the nervous system to regulate cell proliferation and differentiation, has been recently identified as a causal gene for autism spectrum disorder (ASD), but the mechanism behind remains unknown.
To uncover this mechanism, in this study we characterized the involvement of klf7 in circadian rhythm by knocking down klf7 in N2A cells and examining the rhythmic expression of circadian genes, especially Clock gene. We constructed klf7−/− mice and then investigated into klf7 regulation on the expression of rhythm genes in vivo as well as the use of melatonin to rescue the autism behavior. Our results illustrated that circadian rhythm was disrupted in klf7 knockdown cells and that klf7−/− mice showed autism-like behavior. Also, we found that Clock gene was downregulated in the brain of these klf7−/− mice and that the downstream rhythm genes of Clock were disturbed. Melatonin, as a circadian regulation drug, could regulate the expression level and amplitude of rhythm genes in klf7 knockout cells and further rescue the autistic behavior of klf7−/− mice.
Klf7 deficiency causes ASD by disrupting circadian rhythm related genes to trigger rhythm oscillations. To treat ASD, maintaining circadian homeostasis is promising with the use of melatonin.
Krüppel-like factor 7 (klf7) is a transcription factor in the nervous system to regulates cell proliferation and differentiation . Klf7 has been proposed as a candidate gene for autism spectrum disorder (ASD) with the results that patients with deletion in 2q33.3q34, where klf7 is located, show autistic traits [2,3,4,5,6]. More recently, Powis et al., reported 4 unrelated individuals with de novo mutations in klf7 shared similar clinical traits, including developmental delay/intellectual disability, hypotonia, feeding/swallowing issues, psychiatric features, and neuromuscular symptoms . In addition, rare de novo non-coding variants in this gene have been observed in ASD probands [8, 9]. These studies suggest a possible link between klf7 mutation and ASD. In our recent study, klf7 was confirmed to be an autism-related gene by constructing Nestin-Cre conditional knockout mouse model, though the mechanism behind remains unknown .
Klf7 was identified as a regulator related to suprachiasmatic nucleus function and could enrich rhythm genes in mouse suprachiasmatic nucleus [11, 12]. In this regard, we also discovered that target genes of both human klf7 and mouse klf7 were significantly enriched in circadian rhythms , indicating that klf7 is involved in circadian rhythm. More evidences have been recently collected to support the strong link between dysregulated circadian rhythm and ASD. Researches in cognitive and developmental psychology have discovered the important role of rhythm in the early development of social interaction  and impaired circadian clock network causes susceptibility to ASD [14, 15]. Circadian rhythm is regulated by clock-controlled genes (CCGs), variations in CCGs were considered to affect physiological functions, which may lead to disease susceptibility . Investigations of circadian rhythm related genes showed that Clock, Per1, Per2, Per3, Npas2, Bmal1, Cry1, Cry2 and Dbp mutated frequently in ASD individuals [17, 18]. Also, researches suggested that abnormal CCGs had a significant impact on the development of ASD, especially for sleep–wake rhythm and communication disorders [19,20,21,22,23]; and that circadian clock genes played a significant role in ASD , suggesting the link between circadian rhythm and ASD. Based on the above discussion, in this study we hypothesized that klf7 could contribute to ASD via disrupting circadian rhythm, and tested by investigating into the klf7 regulation on the expression of rhythm genes in vivo. Furthermore, we demonstrated that the use of melatonin could rescue the autism behavior.
Klf7 is involved in regulating the circadian rhythm-related genes
We analyzed the expression of klf7 in N2A cells and obtain the diurnal rhythm (Fig. 1A), which peaked at CT16 and nadired at CT0 (CT0 is the onset at the lighting period). In our previous study , we found that the target genes of klf7 were significantly enriched in circadian rhythms. In order to confirm the effects of klf7 on circadian rhythm genes, we knocked down klf7 level in N2A cells with a lentiviral short hairpin RNA (shRNA). Meanwhile, we also transiently transfected klf7 gene into N2A cells. qRT-PCR results confirmed that klf7 could regulate the expression of some circadian rhythm genes (including Clock, Nr1d1, Nr1d2, Nfil3 and Cry1) (Fig. 1B, C). Organoids can simulate early brain development in vitro, which is necessary for the investigation of neurodevelopmental disorders at many developmental stages [25,26,27]. We knocked down klf7 with shklf7 in human brain organoid model, which was derived from hiPSCs (cell name: DYR0100, serial no.: SCSP-1301, National Collection of Authenticated Cell Cultures), and found that the expression level of rhythm genes was also disrupted, further supporting that klf7 was involved in the regulation of rhythm genes (Fig. 1D). These data suggest that klf7 may play an important role in the regulation of circadian rhythm genes.
Klf7 binds to Clock gene and forms a feedback loop with Clock gene
To verify the transcriptional regulation of klf7 on circadian rhythm genes, we performed ChIP-seq assays in N2A cells with the results showing that klf7 binded to Clock, Nr1d1, Nr1d2, Nfil3 and Tef (Fig. 2A). Clock gene is a core transcription factor in the circadian rhythm system, which forms a transcription-translation auto-regulatory feedback loop with Bmal1 and its target genes Per, Cry and Nr1d1. Target genes accumulate rhythmically and form a suppressor complex that interacts with Clock and Bmal1 to inhibit transcription . Considering the significant influence of klf7 on Clock gene expression level and the core role of Clock in circadian system, we selected Clock gene to examine the relationship between circadian gene and ASD. The results from the Luciferase reporter gene assay confirmed that klf7 positively regulated the transcription of Clock gene in a dose dependent manner (Fig. 2B). The promoter of Clock gene was truncated into a 300 bp sequence and a 280 bp sequence, with the results demonstrating that specific 300 bp sequence could be positively bound by klf7 (Fig. 2C). PCR primers could amplify the 300 bp sequence of Clock gene promoter in IP group (Fig. 2D). These results suggest that klf7 can regulate Clock gene transcription.
The circadian rhythm system consists of transcriptional activators and transcriptional repressors, which can be assembled into feedback loops . To determine whether there is such a transcriptional feedback loop between klf7 and Clock gene, we transiently transfected Clock gene in N2A cells by a dose dependent manner and found that Clock could inhibit klf7 expression, forming a negative regulatory loop (Fig. 2E). In order to verify the effect of increased Clock level on downstream rhythm genes, we selected 250 ng Clock plasmid for transient transfection of normal cells and klf7 knockout cells; qPCR was used to examine the expression efficiency after 48 h. Our results showed that high expression of Clock gene was observed in both normal cells and klf7 knockout cells (Fig. 2F). However, the expression level of klf7 was significantly decreased in both normal cells and klf7 knocked out cells after clock overexpression (Fig. 2G), and the changes in protein level were consistent with mRNA result (Fig. 2H). Also, we found that the levels of Per2, Per3 and Rorb were significantly downregulated in klf7 knockout cells. Increasing Clock level did not improve these genes expression levels, but caused them significantly decreased (Fig. 2I). We speculated that this was caused by Clock inhabiting klf7 expression, being associated with the regulation of klf7 on rhythm genes.
Klf7 knockdown cells demonstrate dysregulated circadian phenotypes
To obtain insight into the possible role of klf7 on Clock and other circadian genes of mammalian cells, we examined the expression changes of rhythm genes in klf7 knockdown cells. We knocked down klf7 with a lentiviral short hairpin RNA (shRNA) in N2A, and then found that Clock gene mRNA level and protein level were reduced (Fig. 3A, B). Due to the expression of Clock were decreased significantly in klf7 knockdown cells, we examined to see weather klf7 deficiency led to the abnormal circadian phenotypes. We synchronized these N2A cells and checked mRNA levels of circadian rhythm genes every four hours during the diurnal circadian cycle. We observed that Clock gene exhibited a rhythmic expression in normal cells and that this rhythmic expression was disrupted and even exhibited an inversion rhythm at some time points in klf7 knockdown cells. Moreover, the amplitude of Clock gene between the peak and nadir increased in a cycle (Fig. 3C), which is of an important factor in diurnal biology. Notably, abnormal rhythmic expression of Clock, as a core transcriptional regulatory protein in the circadian rhythm system, may lead to oscillations of downstream rhythm genes. To uncover this, we also assayed some other rhythm genes and found that the amplitudes of Nr1d1, Per2, Per3, Cry2 and Dbp gene were changed in a cycle and that the rhythm of Nr1d1 and Per2 were disturbed (Fig. 3C). These results suggested that klf7 deficiency made Clock exhibits an abnormal rhythmic expression and further triggered the oscillations of downstream rhythm genes. It has been shown that SNP in circadian rhythm genes is associated with ASD and attention deficit hyperactivity disorder [18, 30, 31]. The expression level of Clock gene was decreased in ASD individuals  and rs3762836 (p.H542R) located at Clock gene was found in Japanese ASD individuals . In addition, 4q12 copy number variation including Clock gene have been associated with ASD [33, 34]. Interestingly, mutations in Clock direct transcription target Per1 and dimerization partner Npas2 are also associated with ASD [18, 35]. The preponderance of Clock involvement indicates that this transcription factor plays a variety of roles in regulating processes relevant to human health; some are linked to circadian rhythm and the others may be associated with ASD. This further supports our hypothesis that there is a link between circadian rhythm and klf7 and that klf7 can cause ASD by targeting Clock gene.
klf7 deficiency mice show core autistic behavior
To examine if the circadian rhythm disruption can cause ASD in vivo, we constructed Nestin-Cre klf7 conditional knockout mice, in which exon 2 of klf7 is deleted under the control of loxp-Cre system. To examine whether klf7 deficient mice could reproduce autistic behavior as shown in klf7 mutations patients, we performed behavioral test on male WT mice and male klf7−/− mice aged 2–3 months.
Three-chamber social interaction tests were performed to examine whether klf7 deficient mice had a core symptom of ASD—social deficit. WT mice and klf7−/− mice didn’t show obvious preference for the left and right chambers during the adaptation phase (Fig. 4A). When the first strange mouse was introduced to a chamber, WT mice stayed in the mouse chamber significantly longer than that in the empty chamber. There was no significant difference in the time that klf7−/− mice stayed in the mouse chamber and empty chamber, and klf7−/− mice stayed in the mouse chamber significantly shorter than WT mice stayed in the mouse chamber (Fig. 4B, F). WT mice also spent significantly more time interacting with the mouse than that with the empty cage. However, klf7−/− mice did not interact with the mouse compared with WT mice, the interaction time with mouse was significantly lower than that of WT mice (Fig. 4C). Upon introducing the second mouse to another chamber, WT mice preferred to stay in the second-mouse chamber rather than first-mouse chamber. Klf7−/− mice did not show preference, and klf7−/− mice stayed in the second mouse chamber or strange mouse chamber significantly shorter than WT mice (Fig. 4D, F). In terms of the interaction time, WT mice preferred to interact with the second or strange mouse rather than the first or familiar one. In contrast, the interaction time of klf7−/− mice with the strange mouse and familiar mouse were equal, and the interaction time of klf7−/− mice with the strange mouse was significantly lower than that of WT mice with the strange one (Fig. 4E). These results suggested that klf7 deficient mice had abnormal social behavior.
The self-grooming test and the Y-maze spontaneous selection test were performed to assess whether the klf7 deficient mice exhibit repetitive behaviors, which is also a core behavior of ASD. Our results showed that Klf7−/− mice spent significantly more time on grooming than WT mice (Fig. 4G). In the Y maze spontaneous selection experiment, klf7−/− mice preferred the first target arm and did not like to make a change, while WT mice preferred to explore the new arm (Fig. 4H). These results indicated that klf7 deletion mice had obvious repetitive behaviors.
In addition to the core behaviors of ASD, we also examined other symptoms associated with ASD. In the open field experiment, our results showed that the total distance that klf7−/− mice had moved was significantly higher than that of the WT mice (Fig. 4I, L), and that their movement speeds were significantly increased (Fig. 4J) while their rest durations were also significantly reduced (Fig. 4K). In the Morris water maze test, klf7−/− mice showed an obvious decreased learning ability during the five training days (Fig. 4M). During the testing period, WT mice spent more time in the target region while klf7−/− mice didn’t explore the target region at all (Fig. 4N). Moreover, WT mice crossed the platform location more frequently while klf7−/− mice swam freely (Fig. 4O, P). These data suggested that klf7−/− mice had greater motor capacity, but were severely impaired spatial learning and memory function.
Disrupting of Clock function in klf7−/− mice
To determine whether klf7 deficiency can affect circadian rhythms in vivo, we performed the RNA-seq on 7-day-old male WT mice and klf7−/− mice. KEGG enrichment analysis was performed on the differentially expressed genes (DEGs), and the results showed that DEGs were significantly enriched in circadian rhythms (Fig. 5A). Regarding the effect of klf7 on Clock gene, our quantitative PCR results showed the reduced Clock level (Fig. 5B). Clock protein level was also reduced, being consistent with the RNA level (Fig. 5C). To examine whether the decrease of Clock protein could cause the downstream rhythm genes abnormal, we detected the expression levels of other important rhythm genes in the brain of klf7−/− mice, and found that the expression levels of other rhythm genes were also disturbed (Fig. 5D). In our previous study , klf7± mice showed the mild autistic behavior, which could be improved if klf7 is expressed in the brain by the adeno-associated virus. So we also performed the RNA-seq on klf7± mice treated with adeno-associated virus, and analyzed the expression level of rhythm genes in the brain of klf7± mice before and after treatment. Our results showed that increased klf7 level could reverse the expression level of rhythm genes (Fig. 5E) and rescue the autistic behavior . These results suggest that klf7 is an important member involved in the regulation of rhythmic circadian genes in vivo and links circadian rhythms to ASD (Fig. 5F).
Melatonin regulates the expression of rhythm genes and restores the core autistic behavior
Based on the above-discovered relationship among the ASD, circadian rhythms and klf7, we further examined whether restoring circadian rhythms could be used to treat ASD. Melatonin is an endogenous neurohormone synthesized by the pineal gland, and its main role includes regulating circadian rhythm [36, 37]. Previous in vivo studies have shown that melatonin is involved in regulating the expression of several clock genes in central and peripheral system, such as Clock, Per1, Per2, Bmal1, Rev-ERba and Cry1 [38, 39]. As discussed above, klf7 deficiency disrupted the rhythmic expression of Clock and other rhythmic genes. In this study we subsequently examined whether melatonin could restore the circadian rhythm disorder induced by klf7 deletion and rescue the autistic behavior of klf7−/− mice. For this, we first examined its regulatory effect on rhythm disorders in N2A cells. Upon 24 h of subculture, N2A cells were added with shklf7 or control shRNA lentivirus infection for 48 h, then 20 μm/L melatonin were added into medium and incubated for 24 h. Samples were taken every 6 h, and then the effect of melatonin on the expression level of rhythm genes was analyzed (Fig. 6A). Clock gene showed obvious rhythmic expression in the normal cells and this rhythm was disturbed in the shklf7 group. Once melatonin was introduced and incubated, Clock gene was still rhythmically expressed in the normal cells, but with the reduced amplitude. Also, Clock gene was rhythmically expressed in shklf7 group with melatonin but the amplitude was lower than that of the normal cells (Fig. 6B). The regulation of Clock gene expression level and rhythmic oscillation caused by melatonin brought our attention to Cry and Per, both of which are downstream target genes of Clock. Our results in this regard showed that melatonin could restore the rhythmic expression of Cry2 gene (Fig. 6C) and regulate the oscillation amplitude of Per3 (Fig. 6D). Our data demonstrate that melatonin can regulate the expression and oscillation of rhythm genes and maintain the normal rhythm system.
Also, we were interested in whether it could restore abnormal behavior in mice. Mice were intraperitoneally injected with melatonin at 10 mg/kg for 15 days, followed by behavioral analysis. The tested mice were acclimated in the equipment for 10 min, and then a male mouse of the same age, who had not been reared with tested mice in the same cage, was randomly introduced into the left and right chambers. In the PBS injection group, the time of WT mice staying in the interactive mouse chamber was significantly higher than that in the empty cage chamber, while the time of klf7−/− mice staying in the interactive mouse chamber was significantly lower than that in the empty cage chamber, which is also significantly lower than that of WT mice staying in the interactive mouse chamber. In the melatonin injection group, the time of WT mice staying in the interactive mouse chamber was significantly higher than that in the empty cage chamber. The time of melatonin treated klf7−/− mice staying in the interactive mouse chamber was significantly higher than that in empty cage chamber, although it was significantly lower than that of melatonin treated WT mice staying in the interactive mouse chamber, but significantly higher than PBS treated klf7−/− mice staying in the interactive mouse chamber (Fig. 6E). The interaction time of PBS-treated WT mice with interactive mice was significantly higher than that in the empty cage, while the interaction time of PBS-treated klf7−/− mice with interactive mice was significantly higher than that in empty cage, but significantly lower than that of WT mice with the interactive mice. After melatonin treatment, the interaction time between WT mice and the interactive mice was significantly higher than that in the empty cage. Melatonin-treated klf7−/− mice had significantly longer interaction time with the interactive mice than in the empty cage, which wasn’t significantly different from the wild-type mice (Fig. 6F).
When the second strange mouse was introduced, PBS-treated WT mice preferred to stay in this strange mouse chamber rather than familiar mouse chamber, while PBS-treated klf7−/− mice spent significantly less time in the strange mouse chamber than that in the familiar mouse chamber. The time spent in the strange mouse chamber by PBS-treated klf7−/− mice was also significantly lower than the time that PBS-treated WT mice stayed in the strange mouse chamber. Melatonin-treated WT mice spent significantly more time in the strange mouse chamber than in the familiar mouse chamber. The time of melatonin-treated klf7−/− mice in the strange mouse chamber was not different from that in the familiar mouse chamber, which was significantly lower than that melatonin-treated WT mice spent in the strange mouse chamber, but significantly higher than that PBS-treated klf7−/− mice spent in the strange mouse chamber (Fig. 6G). The interaction time spent by PBS-treated WT mice with the strange mouse was significantly higher than that spent with the familiar mouse, while the interaction time spent by PBS-treated klf7−/− mice with the strange mouse and familiar mouse was equal. The interaction time between PBS-treated klf7−/− mice and the strange mouse was significantly lower than that between PBS-treated WT mice and the strange mouse. The melatonin-treated WT mice still preferred to interact with the strange mouse, but melatonin-treated klf7−/− mice still spent equal interaction time between the strange mouse and familiar mouse, and there is no significant difference between interaction time spent by melatonin-treated klf7−/− mice and PBS-treated klf7−/− mice (Fig. 6H). These data suggest that melatonin can provide some degree of improvement from social deficits.
In the Y maze spontaneous selection experiment, the proportion of PBS-treated klf7−/− mice sticking to the original choice was significantly higher than that of WT mice, and the probability of melatonin-treated klf7−/− mice was significantly lower than that after treatment, suggesting that melatonin can help to reduce the repetitive stereotypical behavior of mice (Fig. 6I). These results suggest that circadian rhythm play a significant role in the development of ASD and likely, is a relevant target for the development of therapeutic strategy in ASD.
In our study, we found that klf7 can target rhythmic gene and form a feedback loop with Clock gene. Klf7 deficiency not only caused N2A cells to produce dysregulated circadian phenotypes, but also led to the autistic behavior in mice. Klf7 deficiency leads to the development of ASD and restore circadian rhythms can rescue the autistic behavior, indicating that circadian rhythm is an important cause for ASD. Our study tested and supported the hypotheses that there was a link among klf7, circadian rhythms, and ASD.
Klf7 is a rhythm transcription factor that can enrich rhythm genes in the suprachiasmatic nucleus (SCNs) and regulate the expression of rhythm genes with clock transcription factors [11, 12]. Our previous study also found that klf7 target genes were also significantly enriched in the circadian rhythm pathway . These evidences indicate that klf7 is involved in the regulation of circadian rhythm and may play an important role in circadian activities. Circadian rhythm controls various biological processes in living organisms; impaired circadian clock networks or rhythm function can increase susceptibility to psychiatric disorders, such as ASD . Recent studies on the cognitive and developmental psychology have shown the important role of rhythms in the development of early social interaction . Meanwhile, other studies have also linked ASD to circadian disturbances [41, 42]. This suggests that there is an implication between circadian rhythms and ASD etiology. Klf7 deficiency leads to the circadian dysregulation, implicating that klf7 may be one link between circadian rhythm and ASD.
Clock gene is a core clock transcription factor to drive the daily transcription patterns of some clock-controlled genes (CCGs) in some tissues. Studies have shown that Clock dysregulation can disrupt the gene co-expression network related to neurological diseases and play a role in controlling the transcriptional cascade in human brain evolution . Clock gene mutation can cause a lengthened rest-activity cycle . Research has shown that polymorphism (rs1801260/rs3762836) of Clock gene is associated with attention deficit hyperactivity disorder (ADHA) and ASD [17, 30, 31, 45,46,47], being consistent with the performance of klf7−/− mice in open field experiment as observed in the present study. Moreover, decreased Clock level was coincided with the results reported from ASD individuals . A more recent study has been in support of the proposed association between Clock gene and ASD . In the present study, we tested and proved the relationship between Clock gene and ASD in vivo for the first time and demonstrated that the decreased Clock level also led to abnormalities of its downstream rhythm genes. Abnormalities in rhythm genes have been associated with memory in mice and communication timing in fruit flies [20, 49]. All of evidences further support that circadian rhythm abnormality is a cause of ASD. Clock gene, as a core transcriptional regulator of rhythm genes, is a major factor in ASD development.
Circadian rhythm dysregulation is associated with many diseases, adjusting circadian rhythm is a promising method in the treatment of diseases [50, 51]. In the light of the demonstrated relationship between ASD and circadian rhythm, we also examined the treatment of ASD by adjusting circadian rhythms. Melatonin is a circadian synchronizer  and is involved in regulating circadian expression of Per1, Per2, Bmal1, Clock and Cry1 in central and peripheral tissues [38, 39]. Melatonin can synchronizes circadian oscillations by affecting circadian expression of Per1 and Bmal1 in the cardiovascular system of rat heart , as well as the metabolism and hormonal function by regulating Per2, Clock and Rev-ERba in adipose tissue [54, 55]. Taken together, these studies demonstrate the important role of melatonin in regulating the expression of rhythmic genes and synchronization of central oscillators. In the present study, we illustrate that melatonin can regulate the expression level and amplitude of Clock gene and its downstream target rhythm genes induced by klf7 deficiency. Meanwhile, melatonin can be used to alleviate the core symptoms of ASD in klf7−/− mice, further supporting the link of klf7 to ASD through the regulation of the circadian genes. The role of melatonin in regulating circadian rhythms makes it a good candidate of drugs for ASD treatment.
The limitation of this study is that there are no detailed evidences to explain how the feedback loop between klf7 and Clock is generated, which would be interesting for future research.
Klf7 has been recently identified as a causal gene for autism spectrum disorder, but leaving a lot to be discovered. Our findings in the present study demonstrate that klf7 causes ASD by targeting Clock genes that trigger circadian genes dysregulated and circadian oscillation. In the light of a strong link between dysregulated circadian rhythm and ASD, our study also supports that circadian rhythm may be a relevant target for the development of therapeutic strategy in ASD.
Mice and drug treatment
Homozygous klf7loxp/loxp mice were purchased from Cyagen Biosciences Inc. (Suzhou, China), with loxp sites inserted into exon 2. Klf7−/− mice were constructed by crossing homozygous klf7loxp/loxp mice with mice expressing Cre. Littermates WT mice were used as control for klf7−/− mice. All mice were placed in a temperature-controlled environment under the cycle of 12 h light and 12 h dark. All animal husbandry and laboratory procedures were approved by Harbin Institute of Technology Laboratory Animal Welfare Ethics Committee (IACUC-2020036). In our study, male mice of 2–3 months old were used for behavioral experiments (either with or without drug treatment). Melatonin treatment (10 mg/kg body weight) [56, 57] or PBS was injected via intraperitoneal injection, once a day over 15 days. Mice were randomly grouped into the melatonin group or PBS group and the examination were blinded to genotypes. Once behavioral tests were completed, the mice were used in the subsequent experiments as described below.
ChIP-seq was performed as described previously . Briefly, the cells were crosslinked with 1% formaldehyde for 10 min and neutralized with glycine (150 mM final solution). Chromatin was cleaved into about 300 bp fragments by ultrasonic treatment after cells were lysed. Chromatin fragments were incubated with antibody (Abcam, ab9110; Abcam, ab172730) overnight at 4 °C, followed by incubating with magnetic beads for 2 h at 4 °C. The proteins were eluted and reversed crosslinking overnight at 65 °C. DNA was extracted with Phenol–chloroform for sequencing. The purified DNA were pair-end sequenced on an illumina platform (Illumina, CA, USA).
Clean reads were aligned to the mouse reference genome (Ensemble_GRCm38.94) using BWA mem (v 0.7.12). IP enrichment regions were identified using MACS2 (version 2.1.0) peak calling software. A q-value threshold of 0.05 was used to identify significant peaks of klf7, which was defined as the overlapping peaks between two biological ChIP-seq replications.
Luciferase reporter assay
293 T cells were transfected with PIRM-control or PIRM-klf7 vector and PGL4-Clock vector at a concentration of 1000 ng/per well. Cells were lysed with 100 μL cold 1 × PLB lysis buffer (Promega), which were then shaken for 15 min at room temperature. The lysis was centrifuged at 12,000 rpm for 2 min and 20 μL supernatant was used for luciferase activity assay (Promega).
RNA-seq was performed as described previously . Total RNA was extracted from klf7 knockdown organoids and control organoids, whole brains of 7-day-old male WT mice and 7-day-old male klf7−/− mice, whole brains of klf7± mice injected with AAV-NC virus and AAV-klf7 virus. The whole brains of 4 mice per genotype and 20 organoids per group were considered as one biological replication.
Hisat2 (v2.0.5) was used to align reads to the human genome (homo_sapiens_grch38_p12) and the mouse genome (mus_musculus_grcm38_p6). The featureCounts v1.5.0-p3 was used to count the reading numbers mapped to each gene; and FPKM of each gene was calculated based on the length of the gene and reading number mapped to this gene. The DESeq2 R package (1.20.0) was used to analyze differential expression gene, where padj < 0.01 and |log2(foldchange)|> 1 were used to detect significantly differential expression genes. ClusterProfiler R package (3.8.1) was used to test the statistical enrichment of differential expression genes in KEGG pathways.
N2A cells were used to study klf7 knockdown, 20 μL (1.0 × 105 infectious units of virus) of either validated shklf7 or negative control shRNA lentivirus were added to the cells. After 48 h, qRT-PCR and western blot were performed to examine the knockout efficiency.
Three-chamber social experiment
The social experiment was performed as previously described  with minor modifications. Male mice of 2–3 months old were used and placed in a testing room for 2 h prior to the behavioral test. After that, the mice to be tested were placed in an apparatus for 10 min for habituation. When a strange mouse was introduced to a wire cage in one side of chambers, the tested mice were allowed to explore all three-chambers for 10 min. Then another strange mouse was placed in the empty cage and the tested mice were allowed to explore all three-chambers for 10 min again. The time duration that the tested mice spent in each chamber and the time duration spent on interacting with the first- and second-introduced mice were recorded for examination.
Self-grooming behaviors were tested as previous study  with minor modifications. The mice were first conditioned in the testing device for 10 min and then examined in terms of the total self-grooming time over the next 10 min.
Y maze spontaneous selection experiment
The selection experiment was performed as previous study  with minor modification. The tested mice were placed in the middle arm, allowing to select between the left and right arms to enter. The tested mice were left in the arms selected by themselves for 10 min. Next, the number of mice entered in the originally-selected arm during 10 trials was recorded for examination.
Open field test
The open field test was performed as previous study . The tested mice were allowed to explore apparatus for 10 min and their activities were recorded by a camera for another period of 10 min.
Morris water maze test
Morris water maze test was performed as previously described in the study  with few modifications. Each mouse was allowed to perform 4 trials a day over 5 days. The time to locate the platform was recorded during the training period. In the testing trial, after removing the platform each tested mouse was allow to search for the platform in 60 s. The time spent in each quadrant and the number of crossing the platform location were recorded.
All data were presented as the mean ± SEM and analyzed by unpaired t test, one-way ANOVA test and two-way ANOVA test.
Availability of data and materials
RNA-seq data and ChIP-seq data that supporting the conclusions of this study are available on zenodo web site (https://zenodo.org/). DOI number for RNA-seq data of 1 month old wild type mice and 1 month old klf7± mice is 5236120 (https://zenodo.org/search?page=1&size=20&q=5236120). DOI number for RNA-seq data of klf7±/AAV-control mice and klf7±/AAV-klf7 mice is 5242635 (https://zenodo.org/search?page=1&size=20&q=5242635). DOI number for RNA-seq data of klf7 knockdown human brain organoid model is 5242821 (https://zenodo.org/search?page=1&size=20&q=5242821). DOI number for ChIP-seq data of klf7 in N2A cells is 5243430 (https://zenodo.org/search?page=1&size=20&q=5243430).
Krüppel-like factor 7
Autism spectrum disorder
Short hairpin RNA
Differentially expressed genes
Attention deficit hyperactivity disorder
Laub F, Lei L, Sumiyoshi H, Kajimura D, Dragomir C, Smaldone S, et al. Transcription factor KLF7 is important for neuronal morphogenesis in selected regions of the nervous system. Mol Cell Biol. 2005;25(13):5699–711.
Courtens W, Speleman F, Messiaen L, Bormans J, Van Roy N, Vamos E. Interstitial deletion 2q33.3-q34 in a boy with a phenotype resembling the Seckel syndrome. Am J Med Genet. 1997;71(4):479–85.
Pescucci C, Meloni I, Bruttini M, Ariani F, Longo I, Mari F, et al. Chromosome 2 deletion encompassing the MAP2 gene in a patient with autism and Rett-like features. Clin Genet. 2003;64(6):497–501.
Bisgaard AM, Kirchhoff M, Tumer Z, Jepsen B, Brondum-Nielsen K, Cohen M, et al. Additional chromosomal abnormalities in patients with a previously detected abnormal karyotype, mental retardation, and dysmorphic features. Am J Med Genet A. 2006;140(20):2180–7.
Brandau DT, Lund M, Cooley LD, Sanger WG, Butler MG. Autistic and dysmorphic features associated with a submicroscopic 2q33.3-q34 interstitial deletion detected by array comparative genomic hybridization. Am J Med Genet A. 2008;146A(4):521–4.
Jang DH, Chae H, Kim M. Autistic and Rett-like features associated with 2q33.3-q34 interstitial deletion. Am J Med Genet A. 2015;167A(9):2213–8.
Powis Z, Hart A, Cherny S, Petrik I, Palmaer E, Tang S, et al. Clinical diagnostic exome evaluation for an infant with a lethal disorder: genetic diagnosis of TARP syndrome and expansion of the phenotype in a patient with a newly reported RBM10 alteration. BMC Med Genet. 2017;18(1):60.
Yuen RK, Merico D, Cao H, Pellecchia G, Alipanahi B, Thiruvahindrapuram B, et al. Genome-wide characteristics of de novo mutations in autism. NPJ Genom Med. 2016;1:160271–1602710.
Turner TN, Coe BP, Dickel DE, Hoekzema K, Nelson BJ, Zody MC, et al. Genomic patterns of de novo mutation in simplex autism. Cell. 2017;171(3):710-722 e712.
Tian H, Qiao S, Zhao Y, Jin X, Wang C, Wang R, et al. Kruppel-like transcription factor 7 is a causal gene in autism development. Int J Mol Sci. 2022. https://doi.org/10.3390/ijms23063376.
Wen S, Ma D, Zhao M, Xie L, Wu Q, Gou L, et al. Spatiotemporal single-cell analysis of gene expression in the mouse suprachiasmatic nucleus. Nat Neurosci. 2020;23(3):456–67.
Cheng AH, Bouchard-Cannon P, Hegazi S, Lowden C, Fung SW, Chiang CK, et al. SOX2-dependent transcription in Clock neurons promotes the robustness of the central circadian pacemaker. Cell Rep. 2019;26(12):3191-3202 e3198.
Tordjman S, Anderson GM, Bellissant E, Botbol M, Charbuy H, Camus F, et al. Day and nighttime excretion of 6-sulphatoxymelatonin in adolescents and young adults with autistic disorder. Psychoneuroendocrinology. 2012;37(12):1990–7.
Yenen AS, Cak HT. Melatonin and circadian rhythm in autism spectrum disorders. Turk Psikiyatri Derg. 2020;31(3):201–11.
Charrier A, Olliac B, Roubertoux P, Tordjman S. Clock genes and altered sleep-wake rhythms: their role in the development of psychiatric disorders. Int J Mol Sci. 2017. https://doi.org/10.3390/ijms18050938.
Takahashi JS, Hong HK, Ko CH, McDearmon EL. The genetics of mammalian circadian order and disorder: implications for physiology and disease. Nat Rev Genet. 2008;9(10):764–75.
Yang Z, Matsumoto A, Nakayama K, Jimbo EF, Kojima K, Nagata K, et al. Circadian-relevant genes are highly polymorphic in autism spectrum disorder patients. Brain Dev. 2016;38(1):91–9.
Nicholas B, Rudrasingham V, Nash S, Kirov G, Owen MJ, Wimpory DC. Association of Per1 and Npas2 with autistic disorder: support for the clock genes/social timing hypothesis. Mol Psychiatry. 2007;12(6):581–92.
Toh KL, Jones CR, He Y, Eide EJ, Hinz WA, Virshup DM, et al. An hPer2 phosphorylation site mutation in familial advanced sleep phase syndrome. Science. 2001;291(5506):1040–3.
Sakai T, Tamura T, Kitamoto T, Kidokoro Y. A clock gene, period, plays a key role in long-term memory formation in Drosophila. Proc Natl Acad Sci U S A. 2004;101(45):16058–63.
Janich P, Pascual G, Merlos-Suarez A, Batlle E, Ripperger J, Albrecht U, et al. The circadian molecular clock creates epidermal stem cell heterogeneity. Nature. 2011;480(7376):209–14.
Marcheva B, Ramsey KM, Buhr ED, Kobayashi Y, Su H, Ko CH, et al. Disruption of the clock components CLOCK and BMAL1 leads to hypoinsulinaemia and diabetes. Nature. 2010;466(7306):627–31.
Wimpory D, Nicholas B, Nash S. Social timing, clock genes and autism: a new hypothesis. J Intellect Disabil Res. 2002;46(Pt 4):352–8.
Liu D, Nanclares C, Simbriger K, Fang K, Lorsung E, Le N, et al. Autistic-like behavior and cerebellar dysfunction in Bmal1 mutant mice ameliorated by mTORC1 inhibition. Mol Psychiatry. 2022. https://doi.org/10.1038/s41380-022-01499-6.
Mertens J, Marchetto MC, Bardy C, Gage FH. Evaluating cell reprogramming, differentiation and conversion technologies in neuroscience. Nat Rev Neurosci. 2016;17(7):424–37.
Pasca SP. The rise of three-dimensional human brain cultures. Nature. 2018;553(7689):437–45.
Gordon A, Yoon SJ, Tran SS, Makinson CD, Park JY, Andersen J, et al. Long-term maturation of human cortical organoids matches key early postnatal transitions. Nat Neurosci. 2021;24(3):331–42.
Takahashi JS. Molecular components of the circadian clock in mammals. Diabetes Obes Metab. 2015;17(Suppl 1):6–11.
Reppert SM, Weaver DR. Coordination of circadian timing in mammals. Nature. 2002;418(6901):935–41.
Kissling C, Retz W, Wiemann S, Coogan AN, Clement RM, Hunnerkopf R, et al. A polymorphism at the 3’-untranslated region of the CLOCK gene is associated with adult attention-deficit hyperactivity disorder. Am J Med Genet B Neuropsychiatr Genet. 2008;147(3):333–8.
Xu X, Breen G, Chen CK, Huang YS, Wu YY, Asherson P. Association study between a polymorphism at the 3’-untranslated region of CLOCK gene and attention deficit hyperactivity disorder. Behav Brain Funct. 2010;6:48.
Hu VW, Sarachana T, Kim KS, Nguyen A, Kulkarni S, Steinberg ME, et al. Gene expression profiling differentiates autism case-controls and phenotypic variants of autism spectrum disorders: evidence for circadian rhythm dysfunction in severe autism. Autism Res. 2009;2(2):78–97.
Sarachana T, Zhou R, Chen G, Manji HK, Hu VW. Investigation of post-transcriptional gene regulatory networks associated with autism spectrum disorders by microRNA expression profiling of lymphoblastoid cell lines. Genome Med. 2010;2(4):23.
Griswold AJ, Ma D, Cukier HN, Nations LD, Schmidt MA, Chung RH, et al. Evaluation of copy number variations reveals novel candidate genes in autism spectrum disorder-associated pathways. Hum Mol Genet. 2012;21(15):3513–23.
Neale BM, Kou Y, Liu L, Ma’ayan A, Samocha KE, Sabo A, et al. Patterns and rates of exonic de novo mutations in autism spectrum disorders. Nature. 2012;485(7397):242–5.
Reiter RJ, Calvo JR, Karbownik M, Qi W, Tan DX. Melatonin and its relation to the immune system and inflammation. Ann N Y Acad Sci. 2000;917:376–86.
Pundir M, Papagerakis S, De Rosa MC, Chronis N, Kurabayashi K, Abdulmawjood S, et al. Emerging biotechnologies for evaluating disruption of stress, sleep, and circadian rhythm mechanism using aptamer-based detection of salivary biomarkers. Biotechnol Adv. 2022;59: 107961.
Dardente H, Menet JS, Poirel VJ, Streicher D, Gauer F, Vivien-Roels B, et al. Melatonin induces Cry1 expression in the pars tuberalis of the rat. Brain Res Mol Brain Res. 2003;114(2):101–6.
von Gall C, Weaver DR, Moek J, Jilg A, Stehle JH, Korf HW. Melatonin plays a crucial role in the regulation of rhythmic clock gene expression in the mouse pars tuberalis. Ann N Y Acad Sci. 2005;1040:508–11.
Tordjman S, Najjar I, Bellissant E, Anderson GM, Barburoth M, Cohen D, et al. Advances in the research of melatonin in autism spectrum disorders: literature review and new perspectives. Int J Mol Sci. 2013;14(10):20508–42.
Cortesi F, Giannotti F, Ivanenko A, Johnson K. Sleep in children with autistic spectrum disorder. Sleep Med. 2010;11(7):659–64.
Richdale AL, Schreck KA. Sleep problems in autism spectrum disorders: prevalence, nature, and possible biopsychosocial aetiologies. Sleep Med Rev. 2009;13(6):403–11.
Fontenot MR, Berto S, Liu Y, Werthmann G, Douglas C, Usui N, et al. Novel transcriptional networks regulated by CLOCK in human neurons. Genes Dev. 2017;31(21):2121–35.
Vitaterna MH, King DP, Chang AM, Kornhauser JM, Lowrey PL, McDonald JD, et al. Mutagenesis and mapping of a mouse gene, Clock, essential for circadian behavior. Science. 1994;264(5159):719–25.
Serretti A, Benedetti F, Mandelli L, Lorenzi C, Pirovano A, Colombo C, et al. Genetic dissection of psychopathological symptoms: insomnia in mood disorders and CLOCK gene polymorphism. Am J Med Genet B Neuropsychiatr Genet. 2003;121B(1):35–8.
Ozburn AR, Purohit K, Parekh PK, Kaplan GN, Falcon E, Mukherjee S, et al. Functional implications of the CLOCK 3111T/C single-nucleotide polymorphism. Front Psychiatry. 2016;7:67.
Shi SQ, White MJ, Borsetti HM, Pendergast JS, Hida A, Ciarleglio CM, et al. Molecular analyses of circadian gene variants reveal sex-dependent links between depression and clocks. Transl Psychiatry. 2016;6: e748.
Huang MC, Ho CW, Chen CH, Liu SC, Chen CC, Leu SJ. Reduced expression of circadian clock genes in male alcoholic patients. Alcohol Clin Exp Res. 2010;34(11):1899–904.
Garcia JA, Zhang D, Estill SJ, Michnoff C, Rutter J, Reick M, et al. Impaired cued and contextual memory in NPAS2-deficient mice. Science. 2000;288(5474):2226–30.
Ruan W, Yuan X, Eltzschig HK. Circadian rhythm as a therapeutic target. Nat Rev Drug Discov. 2021;20(4):287–307.
Lee Y. Roles of circadian clocks in cancer pathogenesis and treatment. Exp Mol Med. 2021;53(10):1529–38.
Pevet P, Challet E. Melatonin: both master clock output and internal time-giver in the circadian clocks network. J Physiol Paris. 2011;105(4–6):170–82.
Zeman M, Herichova I. Melatonin and clock genes expression in the cardiovascular system. Front Biosci (Schol Ed). 2013;5(2):743–53.
Kennaway DJ, Owens JA, Voultsios A, Wight N. Adipokines and adipocyte function in Clock mutant mice that retain melatonin rhythmicity. Obesity (Silver Spring). 2012;20(2):295–305.
Delezie J, Dumont S, Dardente H, Oudart H, Grechez-Cassiau A, Klosen P, et al. The nuclear receptor REV-ERBalpha is required for the daily balance of carbohydrate and lipid metabolism. FASEB J. 2012;26(8):3321–35.
Wang HB, Tahara Y, Luk SHC, Kim YS, Hitchcock ON, MacDowell Kaswan ZA, et al. Melatonin treatment of repetitive behavioral deficits in the Cntnap2 mouse model of autism spectrum disorder. Neurobiol Dis. 2020;145: 105064.
Wang Q, Zhu D, Ping S, Li C, Pang K, Zhu S, et al. Melatonin recovers sleep phase delayed by MK-801 through the melatonin MT2 receptor- Ca(2+) -CaMKII-CREB pathway in the ventrolateral preoptic nucleus. J Pineal Res. 2020;69(3): e12674.
Odawara J, Harada A, Yoshimi T, Maehara K, Tachibana T, Okada S, et al. The classification of mRNA expression levels by the phosphorylation state of RNAPII CTD based on a combined genome-wide approach. BMC Genomics. 2011;12:516.
Peca J, Feliciano C, Ting JT, Wang W, Wells MF, Venkatraman TN, et al. Shank3 mutant mice display autistic-like behaviours and striatal dysfunction. Nature. 2011;472(7344):437–42.
Shoji H, Hagihara H, Takao K, Hattori S, Miyakawa T. T-maze forced alternation and left-right discrimination tasks for assessing working and reference memory in mice. J Vis Exp. 2012. https://doi.org/10.3791/3300.
Katayama Y, Nishiyama M, Shoji H, Ohkawa Y, Kawamura A, Sato T, et al. CHD8 haploinsufficiency results in autistic-like phenotypes in mice. Nature. 2016;537(7622):675–9.
We thank Zhiwei Huang for his critical reading and comments on this manuscript.
This research is funded by the National Natural Science Foundation of China (Grants 51773050) and Heilongjiang Touyan Team (HITTY-20190034).
Ethics approval and consent to participate
All institutional and national guidelines for the care and use of laboratory animals were followed. All animal husbandry and experimental procedures were approved by the Experimental Animal Welfare Ethics Committee of Harbin Institute of Technology (IACUC-2020036).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Tian, H., Jiao, Y., Guo, M. et al. Krüppel-like factor 7 deficiency causes autistic-like behavior in mice via regulating Clock gene. Cell Biosci 12, 166 (2022). https://doi.org/10.1186/s13578-022-00903-6
- Clock gene
- Circadian rhythm