Skip to main content

Advertisement

Loss-of-function mutation in inositol monophosphatase 1 (IMPA1) results in abnormal synchrony in resting-state EEG

Article metrics

Abstract

Background

Dysregulation of the inositol cycle is implicated in a wide variety of human diseases, including developmental defects and neurological diseases. A homozygous frameshift mutation in IMPA1, coding for the enzyme inositol monophosphatase 1 (IMPase), has recently been associated with severe intellectual disability (ID) in a geographically isolated consanguineous family in Northeastern Brazil (Figueredo et al., 2016). However, the neurophysiologic mechanisms that mediate the IMPA1 mutation and associated ID phenotype have not been characterized. To this end, resting EEG (eyes-open and eyes-closed) was collected from the Figueredo et al. pedigree. Quantitative EEG measures, including mean power, dominant frequency and dominant frequency variability, were investigated for allelic associations using multivariate family-based association test using generalized estimating equations.

Results

We found that the IMPA1 mutation was associated with relative decreases in frontal theta band power as well as altered alpha-band variability with no regional specificity during the eyes-open condition. For the eyes-closed condition, there was altered dominant theta frequency variability in the central and parietal regions.

Conclusions

These findings represent the first human in vivo phenotypic assessment of brain function disturbances associated with a loss-of-function IMPA1 mutation, and thus an important first step towards an understanding the pathophysiologic mechanisms of intellectual disability associated with the mutation that affects this critical metabolic pathway.

Background

Dysregulation of the inositol cycle is implicated in a wide variety of human diseases including developmental defects, cancer, diabetes, and neurological diseases. A homozygous frameshift mutation in the gene coding for the enzyme inositol monophosphatase 1 (IMPA1) has recently been associated with severe intellectual disability (ID) in a geographically isolated consanguineous family in Northeastern Brazil [1]. Similar findings of ID have been found in another consanguineous cohort in Pakistan associated with an IMPA1 mutation [2]. Preliminary case study MRI findings in a single subject from the Brazilian cohort revealed no structural abnormalities and no reduction of myo-inositol in the basal ganglia (consistent with similar findings in rodents) suggesting that neural disturbances may be more present and detectable at the circuit/systems level [1]. Therefore, we conducted a resting state electroencephalography (EEG) study with a subset of the Figueiredo et al. [1] cohort to test the hypothesis that the observed IMPA1 loss of function mutation impairs neural circuits underlying normal brain functioning.

EEG allows for the non-invasive recording of coordinated activity across large populations of neurons with high temporal resolution. Scalp potentials measured by EEG reflect the summation of post-synaptic potentials along the apical dendrites of pyramidal neurons oriented perpendicular to the cortical surface [3]. EEG and other electrophysiology recording techniques yield measurements of cognitive and sensory networks in the form of quasi-stationary neural oscillations that offer an in vivo index of circuitry-level neurophysiologic function. More significantly, modern EEG equipment is highly mobile and can be easily transported to conduct novel field research in isolated regions.

Unfortunately, there have been no human or animal electrophysiologic studies to date investigating the effects of IMPA1 or IP3 accumulation, a downstream effect of IMPA1 inhibition. In mice, Impa1 knockout is lethal during embryonic development if inositol is not supplemented in the mother’s diet [4]. Postnatally, clear behavioral effects are observed in homozygous Impa1 mutant mice rescued by dietary myo-inositol, including hyperlocomotion and altered circadian rhythms [5]. This developmental lethality was notably absent in the Figueiredo et al. (2016) [1] cohort which, in the presence of the ID phenotype, suggests the presence of imperfect redundancies or compensatory mechanisms in the inositol cycle pathways that alter neurophysiologic function.

Due to the absence of prior human and animal neurophysiology studies of IMPA1/Impa1 mutation, we could not use such literature to develp an a priori hypothesis of an electrophyiological biomarker of IMPA1 mutation. Therefore, we drew upon the closest pharmacological model of IMPA1 inhibition with a rich neurophysiological literature: lithium. Prior rodent work suggests lithium acts through a combination of Impa1 inhibition and reduced Smit1 mRNA expression [6]. In particular, lithium inhibition of inositol monophosphatase was shown in cell culture to increase IP1 (inositol phosphate-1) concentrations which in turn mediated the rate of phosphoinositide synthesis [7]. We reasoned that a loss-of-function mutation of IMPA1 may result in a similar disruption of the inositol metabolic cycle as lithium administration, and as such, we may use human EEG studies of lithium administration to generate testable hypotheses. In human EEG studies, lithium has been found to enhance early sensory potentials, increase low frequency activity in resting state EEG, and increase event-related beta oscillations [8, 9]. While it is unclear the degree to which IMPA1 related mechanisms specifically influence lithium-induced EEG effects, we believe previous research into the effects of lithium on human EEG may offer a useful framework for guiding our predictions.

To that end, we hypothesized that homozygous carriers of a loss of function IMPA1 mutation would resemble healthy controls under chronic lithium administration and show elevated low frequency activity (i.e., delta, theta, and alpha oscillations) compared to carriers of the wild-type form of IMPA1. To test this hypothesis, we collected resting-state EEG recordings under eyes open and eyes closed conditions in the Figueiredo et al. (2016) [1] cohort and computed standard quantitative EEG measures of frequency band power and variability. We identified an electrophysiological phenotype that did not follow the predictions of increased low-frequency power, but rather identified increased low-frequency variability as a potential characteristic of the loss of function mutation identified in this cohort.

Methods

Participants

Thirty participants in the family identified in Figueiredo et al. (2016) [1] were recruited to participate in the current study. Of the 30, four participants were homozygous for the mutant IMPA1 allele (c.489_493dupGGGCT) (hereafter called HOM, 4 female, 0 male). Nine participants carried only one copy (HET, 2 female, 7 male), and the remaining 17 carried the normal variant (WT, 11 female, 6 male). All HOM patients showed the intellectual disability phenotype. Several attempts were made to recruit additional HOM patients; however, several affected patients were notably irritable and unable to sit for the EEG. All participants were evaluated for neuropsychiatric co-morbidities through the Mini-International Neuropsychiatric Interview (MINI-6.0, Portuguese version). Eight participants were identified as exhibiting symptoms of psychosis (2 HOM, 3 HET, & 3 WT); however, we observed no statistical relationship between genotype and the presence of psychosis symptoms (Fisher’s Exact Tests, all ps > 0.50). Most participants were taking medications (antipsychotics: 1 HOM, 2 HET, 1 WT; benzodiazepines: 1 HOM, 1 HET, 3 WT; SSRIs: 1 WT; antihypertensives: 1 HOM, 4 HET, 5 WT; oral hypoglycemic: 1 HET; allopurinol: 1 WT). Additional demographic information is reported in Table 1.

Table 1 Subject demographics

Electroencephalographic recordings

Resting EEG data were collected in an air-conditioned room under eyes open and eyes closed conditions (2 min each). EEG data were collected from a 32 channel ActiCAP with active Ag/AgCl electrodes using a BrainAmpMR amplifier (Brain Products, Munich, Germany). Scalp locations were based on the International 10/20 system. Data were digitized at 5000 Hz with a 0.1 to 1000 Hz hardware band-pass filter. Ground and reference electrodes were placed at AFz and FCz respectively. Electrode impedances were maintained at or below 20 kΩ. For the final 11 participants (9 WT, 3 HET), the ActiCAP electrode impedance measurement failed, preventing accurate recording of impedances. Therefore, online EEG data were evaluated by two experts for overall quality. As an additional index of data integrity, an on-line running average of visual evoked potentials from a separate EEG task were reviewed during the recordings to determine that the EEG data collected were valid.

EEG preprocessing and analysis

Data preprocessing was done offline using custom Matlab scripts (Mathworks, Natick, MA). First, Continuous EEG data were notch filtered from 59 to 61 Hz followed by a 0.2 to 150 Hz band-pass filter. Data were then segmented into 2500 ms epochs for artifact identification. Individual trials and channels were rejected based on statistical distance from channel and trial means and variances. Data which were deemed to be contaminated by artifact were removed from the data, the remaining clean data were submitted to independent components analysis (ICA) [10]. Briefly, ICA is a blind signal-source separation technique that decomposes mixed signals (e.g., scalp EEG) into unmixed ‘components’ based on spatial patterns of activity. Each ICA component acts as a spatial filter defining unique sources of activity by a weighted sum of EEG channels which can be subtracted from the raw data to remove modeled artifacts. Stereotyped artifacts such as blinks, eye movements, heartbeats, and muscle artifacts were visually inspected and removed by expert reviewers [11, 12]. Artifact free data were submitted for a second round of trial and channel evaluation. Individual channels/trials exceeding a 3 standard deviations of the sample mean of the data were removed from the final analysis.

To evaluate the influence of the IMPA1 mutation on the spectral profile of resting EEG, two quantitative EEG (qEEG) measures were derived from the power spectrum densities (PSD) of clean EEG data. First, mean band power (MBP) was calculated using Welch’s method of windowed averaging over a time series [13]. Window sizes were set to half of the data segment length (i.e., 1250 ms) with no window overlap to avoid discontinuities across potentially non-adjacent epochs (final frequency resolution = 0.63 Hz). PSDs were converted to relative PSD by normalizing each subject’s estimate by the sum power over all frequencies up to 100 Hz. Band power was then defined as the mean power across 6 canonical frequency band: (1) Delta, δ < 4 Hz; (2) theta, θ = 4–8 Hz; (3) alpha, α = 9–14 Hz; (4) beta, β = 14–30 Hz; (5) low gamma, low γ = 31–55 Hz; (6) high gamma, high γ = 65–100 Hz. Lastly, we estimated oscillatory network stability by computing dominant frequency variability (DFV) over trials. In each 1250 ms window, the frequency exhibiting the maximum power was identified within frequency bands. DFV was defined as the mean of absolute deviation scores (i.e., \( DFV=\frac{\sum \left(\left|X-\mu X\right|\right)}{n} \)). Mean absolute deviation was used to reduce the potential influence of outlier trials on the variability estimation (i.e., compared to standard deviation). MBP and DFV scores were summarized by averaging estimates within 8 scalp regions (i.e., left frontal, middle frontal, right frontal, left temporal, central, right temporal, parietal, occipital). Matlab scripts used for this analysis can be found at http://github.com/cholab/IMPA1-EEG.

Statistical analysis

To test for associations between the IMPA1 genotype and the observed qEEG phenotypes, we employed a multivariate family-based associate test using generalized estimating equations (FBAT-GEE; FBAT-Toolkit v204, http://sites.google.com/view/fbat-web-page) [14]. Briefly, FBAT-GEE is an extension of a traditional FBAT that tests for associations between allelic presence and a given phenotype while conditioning the null hypothesis of the test statistic on the observed phenotype distribution. FBAT-GEE expands on the traditional FBAT statistic given by

$$ {\chi}^2=\frac{{\left(S-E(S)\right)}^2}{V_S} $$
(1)

where \( S=\sum \limits_{i=1}^n{t}_i{x}_i \), the expected value, \( E(S)=\sum \limits_{i=1}^n{t}_iE\left({x}_i|{p}_{i1},{p}_{i2}\right) \) …, and the variance \( {V}_s=\sum \limits_{i=1}^n{t}_i^2 Var\left({x}_i|{p}_{i1},{p}_{i2}\right) \). By this metric, both quantitative (i.e., continuous) or categorical (i.e., bivariate) phenotype may be represented as ti for the ith individual. The value of xi codes for the hypothesized phenotype expression based on allele frequency for the marker of interest (i.e., dominant, additive, or recessive). Both E(S) and Vs set the normalization parameters to be conditional on the genotype expression in the parents.

The multivariate extension replaces S and E(S) with the m-dimensional vector \( \overset{\sim }{S} \) defined as

$$ \overset{\sim }{S}=\sum \limits_{i=1}^n{t}_i\left({x}_i-E\left({x}_i|{p}_{i1},{p}_{i2}\right)\right) $$
(2)

where m is number of phenotypes being tested, and the variance is given by the m x m matrix defined as

$$ {V}_{\overset{\sim }{S}}= Var\left(\overset{\sim }{S}\right)=\sum \limits_{i=1}^n{t}_i{t}_i^t Var\left({x}_i|{p}_{i1},{p}_{i2}\right) $$
(3)

where t signifies a vector transpose. The final FBAT-GEE statistic is thus given by

$$ {\chi}_{FBAT- GEE}^2={\overset{\sim }{S}}^t\ {V}_{\overset{\sim }{S}}^{-1}\ \overset{\sim }{S} $$
(4)

which is asymptotically χ2-distributed with degrees of freedom equal to \( k=\operatorname{rank}\left({V}_{\overset{\sim }{s}}\right) \). The statistical inference follows that a significant test statistic indicates the association between an allele loading and the observed phenotype is unlikely to occur by chance given the observed distribution of alleles in the parents. Since the intellectual disability associated with the IMPA1 mutation has previously demonstrated a strongly recessive pattern, we elected to use a recessive FBAT model where aa carriers are coded as 1, and AA and Aa carriers are coded as 0.

The recessive FBAT-GEE model was employed based on the a priori assumption that the inheritance pattern would follow that of the more global phenotype of intellectual disability. However, since the narrower neurophysiologic phenotypes investigated here may show expression even with partial allelic loading, as an additional exploratory step, we repeated our FBAT-GEE analyses with the additive model (i.e., AA = 0, Aa = 1, and aa = 2) to identify potential phenotypes that are more sensitive to such partial loading.

Statistical analyses were conducted in the FBAT-GEE framework for each frequency band and qEEG measure with scalp region acting as the multivariate dimension. Univariate FBATs were run post-hoc within each region to determine if effects were global or regional. Owing to the small sample size, the uniqueness of the mutation identified in the study population, and the generally exploratory nature of our investigation, all tests were evaluated at α = 0.05, uncorrected. Therefore, we present the following analysis under the framework of hypothesis generation, and stress the need for confirmatory follow-up experiments.

Results

Quantitative EEG measures including mean power and dominant frequency variability were extracted from resting EEG data of 30 participants (17 WT, 9 HET, 4 HOM). From this sample, 27 subjects (14 WT, 9 HET, 4 HOM) were from the eight families informative for the IMPA1 mutation and therefore chosen for our analysis.

Both quantitative EEG measures were natural log transformed prior to FBAT testing to better approximate a normal distribution. We performed FBAT analysis to evaluate the association between the IMPA1 mutation and our quantitative phenotypes (i.e., MBP, and DFV). See Tables 2 and 3 for a summary of results. First, a multivariate FBAT-GEE was performed for each of the 6 frequency bands across all scalp regions; significant associations were found with θ-band power (χ2 = 18.451, p = .018) and dominant α-band variability (χ2 = 19.771, p = .011) for the eyes-open condition, and dominant θ-band variability (χ2 = 15.848, p = .045) for the eyes-closed condition. With these identified frequency bands, the 8 individual scalp regions were subsequently tested with univariate analysis. For the eyes-open condition, lower than expected θ power over the left frontal scalp region was significantly associated (Z = − 2.211, p = .027) with the mutated allele, while higher than expected θ power over the right frontal scalp region was significantly associated (Z = 2.248, p = .025) with the wild-type allele (see Fig. 1). For the eyes-closed condition, dominant θ variability over the central scalp region was significantly associated (Z = 2.411, p = .016) with the mutated allele, while dominant θ variability over the parietal scalp region was significantly associated (Z = − 2.329, p = .020) with the wild-type allele (see Fig. 2). Univariate FBATs showed no significant associations for alpha, indicating a general effect across the scalp for dominant α-band variability.

Table 2 FBAT-GEE for eyes-open condition
Table 3 FBAT-GEE for eyes-closed condition
Fig. 1
figure1

Mean (± SEM) theta-band power by scalp region and genotype for eyes open and eyes closed conditions. (BLUE = WT, TEAL = HET, YELLOW=HOM). Arrows indicate the direction of significant associations between theta-band power and specific IMPA1 alleles under the recessive FBAT model. Arrows on the left of a graph indicate an association with the wild-type allele within a specific region. Arrows on the right of the graph indicate an association with the mutant allele within a specific region. a In the eyes open condition, significantly lower left frontal theta power was associated with the mutant allele, and significantly greater right frontal theta power was associated with the wild-type allele. The right frontal pattern was also present under the additive FBAT model. b In the eyes closed condition, these differences were not observed. This lateral shift in theta-band power representations on the scalp suggests subtle anatomical or network-level differences may arise during development in the presence of the IMPA1 mutation

Fig. 2
figure2

Mean (± SEM) dominant frequency variability by scalp region and phenotype. (BLUE = WT, TEAL = HET, YELLOW=HOM). Arrows indicate the direction of significant associations between theta-band power and specific IMPA1 alleles under the recessive FBAT model. Arrows on the left of a graph indicate an association with the wild-type allele within a specific region. Arrows on the right of the graph indicate an association with the mutant allele within a specific region. a No significant associations were observed for theta-band variability in the eyes open condition. b In the eyes closed condition, the mutant allele was associated with a significantly higher central theta variability, and the wild-type allele was associated with significantly lower parietal theta variability. These patterns supplement the power results suggesting a core theta-band oscillatory disturbance resulting from the IMPA1 mutation. c Eyes open and d eyes closed alpha-band variability did not show any regional associations with specific alleles; however, multivariate-FBAT analyses found significant associations between the wild-type allele and alpha variability in the eyes open condition suggesting a more global phenomenon

As noted above, we conducted FBAT-GEE analyses employing the additive model as an exploratory analysis. This showed a significant association (χ2 = 16.388, p = .037) between θ-band power and the mutated allele for the eyes-open condition. Univariate analysis showed lower than expected θ power over the right frontal scalp region was found to be significantly associated (Z = − 2.542, p = .011).

Discussion

In the current study, we conducted a resting EEG study to investigate the electrophysiologic phenotype associated with the IMPA1 mutation in a geographically isolated, consanguineous cohort in Brazil [1]. Our primary findings with the IMPA1 mutation being associated with relative decreases in left and right frontal θ power and altered α-band variability with no identifiable regional specificity during the eyes-open condition. For the eyes-closed condition, there was altered dominant θ frequency variability in the central and parietal regions. These findings differ from the elevated low-frequency activity that would be hypothesized based on such findings in healthy individuals under lithium administration. However, it is not overly surprising that our findings may depart from these simple predictions based on acute lithium administration, given the complexity of the inositol metabolic pathway, the comparative timeframes involved, and specific population being investigated. Regardless, the patterns identified in this analysis will need to be replicated in future studies to establish the generalizability of these effect to other carriers of the IMPA1 mutation and IMPA1 function more generally.

The IMPA1 mutation in the cohort of the current study was associated with severe ID and disruptive behavior. Most homozygous individuals found it very challenging or were unable to complete simple computer or neuropsychological tests of cognitive and intellectual functioning (and a number of identified individuals could not participate at all due to the severity of their impairment and debilitation). Resting-state EEG placed minimal burden on participants and as such was an ideal approach for probing neural phenotypic expression of the IMPA1 mutation in this population. While no cognitive process is actively engaged through explicit task performance, such intrinsic resting state activity is nevertheless thought to reflect the activity and functional integrity of networks that support various cognitive and perceptual processes [15, 16] and segregates into distinct, separable frequency bands [17].

In the current study, two prominent rhythms in resting EEG, θ-band and α-band activity, showed variation with the IMPA1 mutation. Theta oscillations have been associated with cognitive control functions such as error monitoring [18] and working memory [19, 20] including the modulation of local gamma activity [21] and mediating cortical interactions with the hippocampus [22]. In our cohort, we found evidence of reduced frontal theta power, in association with intellectual disability. While formal assessment of working memory and cognitive control was not possible due to the level of impairment in individuals homozygous for IMPA1 mutation, it seems highly likely that the disturbances in theta were associated with severe impairments in both these cognitive processes. Interestingly, both these cognitive processes are highly correlated with IQ [23,24,25,26] and thus impairments in these fundamental cognitive functions may underlie the significant intellectual disability present in our cohort. Furthermore, altered θ-band and α-band oscillations have been observed in other disorders associated with intellectual disability including fragile X syndrome (FXS) [27, 28] and attention-deficit hyperactivity disorder (ADHD) [29, 30]. The precise alterations differ between specific disorders, but suggest θ-band and α-band oscillations represent a potential convergent mechanism underlying intellectual disability across etiologies. Further research will be necessary to evaluate the usefulness of θ/α oscillations as a generalizable biomarker of ID.

Putative mechanisms of alpha and theta oscillations

Alpha-band oscillations are thought to underlie a different putative function. By and large, alpha oscillations are considered to reflect cortical inhibitory processes which can serve to suppress or gate information transfer in the brain (i.e., the inhibition-timing hypothesis) [31, 32]. In the context of cognitive tasks, alpha oscillations are considered to reflect a top-down control mechanism to inhibit task-irrelevant or task-distracting information [33]. By contrast, decreases in alpha activity are often observed in task relevant brain regions [34,35,36,37,38]. In the sensory domain, the phase of ongoing alpha oscillations have been shown to influence bottom-up perceptual fidelity [39, 40] which similar to theta oscillations may help bind information carried and maintained in high frequency, gamma-band oscillations [41]. This balance between suppression and gating functions are thought to be driven by thalamocortical neurons via muscarinic acetylcholine (mAChR) and metabotropic glutamate type-1 receptors (mGluR1) [42, 43]. Notably, mAChRs are Gq-protein coupled receptors known to upregulate inositol triphosphate (IP3) through activation of phospholipase C [44]. Given the role of IMPA1 in modulating the production of IP3, our finding of a significant multivariate FBAT association between the wild-type allele and alpha-band variability suggests that intact IMPA1 functioning may support the flexible coordination of thalamocortical alpha rhythms, which in the absence of a sufficient supply of IP3, results in a decrease in alpha-band variability for homozygous carriers of the IMPA1 mutation.

The theta disturbances in the current study manifested in both the power and dominant frequency variability of theta in association with the IMPA1 mutation. The power of any oscillatory rhythm is dependent on a number of factors including the numbers of pyramidal neurons and synapses giving rise to the rhythm and how synchronous their activations are. It is not yet known the degree to which IMPA1 mutation may be associated with decreased neuronal number or disturbances in synaptic morphology or function. However, Impa1 knockout has been found to increase autophagy [6] which, in a non-physiologic context, may be deleterious to synaptic function or morphology and thus diminish EEG theta power. How the IMPA1 mutation may affect the dominant frequency is also unclear, but potential mechanisms are suggested by studies of calbindin (CB), which can activate IMPase [45,46,47]. CB suppression in hippocampal excitatory neurons has been associated with memory impairments and is thought to be mediated by downstream effects on IMPase [47]. The firing of hippocampal CB positive pyramidal cells is strongly locked to theta rhythms in the context of spatial navigation and memory [48, 49]. It is not clear that such findings in the hippocampus would have relevance to the cortical theta rhythms that are detectable by EEG as in the current study. However, cortical theta rhythms arise from the coordination of CB positive basket interneurons and pyramidal cells which, in turn, may have monosynaptic input from the ventral hippocampus and thus be modulated by the strong theta rhythms generated there [50]. Thus, if the kinds of CB suppression effects as observed in the hippocampus are mediated by impact on IMPase function, the IMPA1 mutation could lead to disturbances in the coordination of network activity that manifest as instabilities in the dominant frequency of cortical theta rhythms.

Our findings of impaired cortical theta rhythms run counter to the predictions that would be made taking lithium as a model of IMPA1 mutation. Administration of lithium leads to reductions in IMPA1 activity and enhancements in low-frequency EEG rhythms [8]. However, lithium also results in reduced Smit1 mRNA expression [6] and thus it is not clear what may lead to the observed enhancements in theta. Further, even if findings of increased theta with lithium were attributable to decreased IMPA1 function, the impact of an IMPA1 mutation may be fundamentally different from lithium administration. In the case of an IMPA1 mutation, the protracted course of effects could start as early as in utero and impact a neurodevelopment course that may also involve compensatory mechanisms in inositol metabolism and associated pathways that together result in neurophysiologic outcomes altogether different from much shorter time scale lithium administration in adults.

Strengths and limitations

Our study had a number of strengths, notably including the uniqueness of the study sample and the successful collection of electrophysiologic measures in a remote, rural setting. The mobility of our EEG equipment facilitated access to a geographically isolated sample which allowed the study team to bring the lab to the participants. We found that the homozygous carriers of the IMPA1 allele were too impaired to accurately assess IQ or perform cognitive behavioral tasks during EEG. As such, the data was collected under passive resting state conditions, avoiding reliance on subject motivation and attention. Thus, the present findings are not confounded by differences in attentional engagement as may be the case in cognitive or sensory tasks. The present data supports the feasibility of data collection and meaningful analysis using a resting state paradigm in this population.

There are some limitations that should be considered when interpreting the results. While it is the first study of its kind to investigate the IMPA1 mutation employing human electrophysiology, the sample size was limited. This was to large extent unavoidable due to a number of factors including the geographically isolated nature of the population, the inability of a number of homozygous IMPA1 mutation candidates to participate in the study because of the severity of their impairment, and the lack of on-site lab facilities and expertise to conduct longer-term data collection. We conducted a high-powered multivariate statistical analysis of the data to capitalize on the familial pedigree structure of our sample, but were still limited by the sample size which was underpowered relative to the conservative corrections for our post-hoc multiple comparisons. However, given the rarity of the sample, we elected to present the findings of the current study with the understanding that they require replication in additional cohorts. On a related note, despite the much higher proportion of psychosis exhibited in individuals who were homozygous or heterozygous for the mutant allele, the sample size limited a more definitive examination of the phenotypic relationship between intellectual disability and psychosis in the context of the IMPA1 mutation. Such associations between psychiatric symptoms and cognitive impairments are common in neuropsychiatric disorders such as schizophrenia, but the precise nature of the relationship will require much larger samples, perhaps primarily comparing heterozygous to wildtype groups given the severe cognitive impairment manifesting in individuals homozygous for the mutant allele. Finally, future investigations could also attempt higher-density EEG recordings to facilitate anatomic source analyses of the resting state networks.

Conclusion

We conducted the first electrophysiologic phenotyping human study of the IMPA1 mutation in a geographically isolated, consanguineous cohort. Our findings of disturbances in frontal theta and more global alpha band disturbances raises the possibility of perturbations in certain cellular subtypes and aspects of the inositol metabolic pathways. More mechanistically definitive studies will require preclinical models with controlled manipulation of the IMPA1 expression and electrophysiological recordings. Based on these important preliminary findings, now we are expanding the studies to patient-derived neuronal cell lines in order to clearly elucidate the mechanisms by which impairment of IMPA1 can alter important signaling pathways that may lead to development of intellectual disability and to correlate with in vivo measures of brain activity and cognitive capacity to better understand the pathway from genetic variants to behavior. To our knowledge, however, this current study is the first human in vivo phenotypic assessment of brain function disturbances associated with the IMPA1 mutation, thus representing an important first step towards understanding the pathophysiology of intellectual disability associated with the mutation that affects this critical metabolic pathway.

Abbreviations

CB:

Calbindin

DFV:

Dominant frequency variability

EEG:

Electroencephalography

FBAT-GEE:

Family-based association test with generalized estimating equations

HET:

Heterozygous

HOM:

Homozygous

ICA:

Independent components analysis

ID:

Intellectual disability

IMPA1:

Inositol monophosphatase 1

IP3:

Inositol triphosphate

IQ:

Intelligence quotient

mAChR:

Muscarinic acetylcholine receptor

MBP:

Mean band power

mGluR1:

Metabotropic glutamate type-1 receptor

mRNA:

Messenger ribonucleic acid

SEM:

Standard error of the mean

SMIT1:

Sodium-myoinositol cotransporter 1

WT:

Wild-type

References

  1. 1.

    Figueiredo T, Melo US, Pessoa AL, Nobrega PR, Kitajima JP, Rusch H, et al. A homozygous loss-of-function mutation in inositol monophosphatase 1 (IMPA1) causes severe intellectual disability. Mol Psychiatry. 2016;21:1125–9.

  2. 2.

    Harripaul R, Vasli N, Mikhailov A, Rafiq MA, Mittal K, Windpassinger C, et al. Mapping autosomal recessive intellectual disability: combined microarray and exome sequencing identifies 26 novel candidate genes in 192 consanguineous families. Mol Psychiatry. 2018;23:973–84.

  3. 3.

    Luck S. An introduction to the event-related potential technique. Cognitive neuroscience. Cambridge (Massachusetts): MIT Press; 2005.

  4. 4.

    Cryns K, Shamir A, Van Acker N, Levi I, Daneels G, Goris I, et al. IMPA1 is essential for embryonic development and lithium-like pilocarpine sensitivity. Neuropsychopharmacology. 2008;33:674–84.

  5. 5.

    Ohnishi T, Murata T, Watanabe A, Hida A, Ohba H, Iwayama Y, et al. Defective craniofacial development and brain function in a mouse model for depletion of intracellular inositol synthesis. J Biol Chem. 2014;289:10785–96.

  6. 6.

    Sade Y, Toke L, Kara NZ, Einat H, Rapoport S, Moechars D, et al. IP3 accumulation and/or inositol depletion: two downstream lithium’s effects that may mediate its behavioral and cellular changes. Transl Psychiatry. 2016;6(12):e968.

  7. 7.

    Saiardi A, Mudge AW. Lithium and fluoxetine regulate the rate of phosphoinositide synthesis in neurons: a new view of their mechanisms of action in bipolar disorder. Transl Psychiatry. 2018;8:175.

  8. 8.

    Atagün Mİ. Brain oscillations in bipolar disorder and lithium-induced changes. Neuropsychiatr Dis and Treat. 2016;12:589–600.

  9. 9.

    Atagün Mİ, Güntekin B, Tan D, Tülay EE, Başar E. Lithium excessively enhances event related beta oscillations in patients with bipolar disorder. J Affect Disord. 2015;170:59–65.

  10. 10.

    Bell AJ, Sejnowski TJ. An Information-Maximization Approach to Blind Separation and Blind Deconvolution. Neural Comput. 1995;7(6):1129–59.

  11. 11.

    Delorme A, Makeig S. EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. J Neurosci Methods. 2004;134(1):9–21.

  12. 12.

    Jung TP, Makeig S, Humphries C, Lee TW, McKeown MJ, Iragui V, Sejnowski TJ. Removing electroencephalographic artifacts by blind source separation. Psychophysiology. 2000;37(2):163–78.

  13. 13.

    Welch P. The use of fast Fourier transform for the estimation of power spectra: A method based on time averaging over short, modified periodograms. IEEE Trans Audio Electroacoust. 1967;15(2):70–73.

  14. 14.

    Lange C, Silverman EK, Xu X, Weiss ST. Laird NM. A multivariate family-based association test using generalized estimating equations: FBAT-GEE. Biostatistics. 2003;4(2):195–206.

  15. 15.

    Fox MD, Greicius M. Clinical applications of resting state functional connectivity. Front Syst Neurosci. 2010;4:19.

  16. 16.

    Mantini D, Perrucci MG, Del Gratta C, Romani GL, Corbetta M. Electrophysiological signatures of resting state networks in the human brain. Proc Natl Acad Sci U S A. 2007;104:13170–5.

  17. 17.

    Groppe DM, Bickel S, Keller CJ, Jain SK, Hwang ST, Harden C, et al. Dominant frequencies of resting human brain activity as measured by the electrocorticogram. Neuroimage 2013. 2013;79:223–33.

  18. 18.

    Debener S, Ullsperger M, Siegel M, Fiehler K, von Cramon DY, Engel AK. Trial-by-trial coupling of concurrent electroencephalogram and functional magnetic resonance imaging identifies the dynamics of performance monitoring. J Neurosci. 2005;25:11730–7.

  19. 19.

    Onton J, Delorme A, Makeig S. Frontal midline EEG dynamics during working memory. Neuroimage. 2005;27:341–56.

  20. 20.

    Raghavachari S, Lisman JE, Tully M, Madsen JR, Bromfield EB, Kahana MJ. Theta oscillations in human cortex during a working-memory task: evidence for local generators. J Neurophysiol. 2006;95:1630–8.

  21. 21.

    Canolty RT, Edwards E, Dalal SS, Soltani M, Nagarajan SS, Kirsch HE, et al. High gamma power is phase-locked to Theta oscillations in human neocortex. Science. 2006;313:1626–8.

  22. 22.

    Mitchell DJ, McNaughton N, Flanagan D, Kirk IJ. Frontal-midline theta from the perspective of hippocampal "theta". Prog Neurobiol. 2008;86:156–85.

  23. 23.

    Blair RJ. The emergence of psychopathy: implications for the neuropsychological approach to developmental disorders. Cognition. 2006;101:414–42.

  24. 24.

    Shamosh NA, Deyoung CG, Green AE, Reis DL, Johnson MR, Conway AR, et al. Individual differences in delay discounting: relation to intelligence, working memory, and anterior prefrontal cortex. Psychol Sci. 2008;19:904–11.

  25. 25.

    Conway ARA, Kane MJ, Engle RW. Working memory capacity and its relation to general intelligence. Trends Cogn Sci. 2003;7:547–52.

  26. 26.

    Ackerman PL, Beier ME, Boyle MO. Working memory and intelligence: the same or different constructs? Psychol Bull. 2005;131:30–60.

  27. 27.

    Van der Molen MJW, Van der Molen MW. Reduced alpha and exaggerated theta power during resting-state EEG in fragile X syndrome. Biol Psychol. 2013;92:216–9.

  28. 28.

    Van der Molen MJW, Stam CJ, Van der Molen MW. Resting-state EEG oscillatory dynamics in fragile X syndrome: abnormal functional connectivity and brain network organization. PLoS One. 2014;9:e88451.

  29. 29.

    Bresnahan SM, Barry RJ. Specificity of quantitative EEG analysis of adults with attention deficit hyperactivity disorder. Psychaitry Res. 2002;112:133–44.

  30. 30.

    Barry RJ, Clarke AR, Johnstone SJ, McCarthy R, Selikowitz M. Electroencephalogram theta/beta ratio and arousal in attention-deficit/hyperactivity disorder: evidence of independent processes. Biol Psychiatry. 2009;66:398–401.

  31. 31.

    Klimesch W, Sauseng P, Hanslmayr S. EEG alpha oscillations: the inhibition-timing hypothesis. Brain Res Rev. 2007;53:63–88.

  32. 32.

    Klimesch W. Alpha-band oscillations, attention, and controlled access to stored information. Trends Cogn Sci. 2012;16:606–17.

  33. 33.

    Bonnefond M, Jensen O. Alpha oscillations serve to protect working memory maintenance against anticipated distracters. Curr Biol. 2012;22:1969–74.

  34. 34.

    Bastiaansen MC, Brunia CH. Anticipatory attention: an event-related desynchronization approach. Int J Psychophysiol. 2001;43:91–107.

  35. 35.

    Capotosto P, Babiloni C, Romani GL, Corbetta M. Fronto-parietal cortex controls spatial attention through modulation of anticipatory alpha rhythms. J Neurosci. 2009;29:5863–72.

  36. 36.

    Kelly SP, Lalor EC, Reilly RB, Foxe JJ. Increases in alpha oscillatory power reflect an active retinotopic mechanism for distracter suppression during sustained visuospatial attention. J Neurophysiol. 2006;95:3844–51.

  37. 37.

    Rajagovindan R, Ding M. From Prestimulus alpha oscillation to visual-evoked response: an inverted-U function and its attentional modulation. J Cogn Neurosci. 2011;23:1379–94.

  38. 38.

    Sauseng P, Klimesch W, Schabus M, Doppelmayr M. Fronto-parietal EEG coherence in theta and upper alpha reflect central executive functions of working memory. Int J Psychophysiol. 2005;57:97–103.

  39. 39.

    Lou B, Li Y, Philiastides MG, Sajda P. Prestimulus alpha power predicts fidelity of sensory encoding in perceptual decision making. NeuroImage. 2014;87:242–51.

  40. 40.

    Mazaheri A, DiQuattro NE, Bengson J, Geng JJ. Pre-stimulus activity predicts the winner of top-down vs. bottom-up attentional selection. PLoS One. 2011;6:e16243.

  41. 41.

    Roux F, Uhlhaas PJ. Working memory and neural oscillations: α-γ versus θ-γ codes for distinct WM information? Trends Cogn Sci. 2014;18:16–25.

  42. 42.

    Lőrincz ML, Crunelli V, Hughes SW. Cellular dynamics of cholinergically-induced alpha (8-13 Hz) rhythms in sensory thalamic nuclei in vitro. J Neurosci. 2008;28:660–71.

  43. 43.

    Vijayan S, Kopell NJ. Thalamic model of awake alpha oscillations and implications for stimulus processing. Proc Natl Acad Sci U S A. 2012;109:18553–8.

  44. 44.

    Burford NT, Nahorski SR. Muscarinic m1 receptor-stimulated adenylate cyclase activity in Chinese hamster ovary cells is mediated by Gsα and is not a consequence of phosphoinositidase C activation. Biochem J. 1996;315:883–8.

  45. 45.

    Berggard T, Szczepankiewicz O, Thulin E, Linse S. Myo-inositol monophosphatase is an activated target of calbindin D28k. J Bio Chem. 2002;277:41954–9.

  46. 46.

    Levi I, Eskira Y, Eisenstein M, Gilon C, Hoffman A, Tal-Gan Y, et al. Inhibition of inositol monophosphatase (IMPase) at the calbindin-D28k binding site: molecular and behavioral aspects. Eur Neuropsychopharmacol. 2013;12:1806–15.

  47. 47.

    Li JT, Xie XM, Yu JY, Sun YX, Liao XM, Wang XX, et al. Supressed Calbindin levels in hippocampal excitatory neurons mediate stress-induced memory loss. Cell Rep. 2017;21:891–900.

  48. 48.

    Tang Q, Burgalossi A, Ebbesen CL, Ray S, Naumann R, Schmidt H, et al. Pyramidal and stellate cell specificity of grid and border representations in layer 2 of medial entorhinal cortex. Neuron. 2014;84:1191–7.

  49. 49.

    Ray S, Naumann R, Burgalossi A, Tang Q, Schmidt H, Brecht M. Grid-layout and theta-modulation of layer 2 pyramidal neurons in medial entorhinal cortex. Science. 2014;343:891–6.

  50. 50.

    Hartwich K, Pollak T, Klausberger T. Distinct firing patterns of identified basket and dendrite-targeting interneurons in the prefrontal cortex during hippocampal theta and local spindle oscillations. J Neurosci. 2009;29:9563–74.

Download references

Acknowledgments

We would like to thank the families who gave their time and energy and without whose support and hospitality this work would not be possible.

Funding

This work was supported by funding from the Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), process n° 2016/09618–5, Fundação de Apoio à Pesquisa do Estado da Paraíba (FAPESQ/CNPq - PPSUS 015/2014), Universidade Estadual da Paraíba (PROPESQ), and Departments of Psychiatry and Behavioral Science at the University of Texas Health Science Center at Houston and Baylor College of Medicine.

Availability of data and materials

The datasets used and/or analyzed during the currents study is available from the corresponding author on request.

Author information

CPW, MR, NM, RYC processed and analyzed the EEG data presented; CPW, AP, TF, USM, PN, SS, RYC participated in EEG data collection; TF, MZ conducted the genetic analysis to identify the genotypes of the affected participants; AP, TF, PN, SS, RYC performed clinical characterization of all participants; CPW, AP, TF, FK, SS, RYC conceived of the study, participated in its design and coordination; CPW, AP, TF, MR, RYC interpreted the results and contributed to the writing of the manuscript; All authors read and approved the final manuscript.

Correspondence to Christopher P. Walker.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the ethics committee of Universidade Estadual da Paraiba (CAAE: 62513316.0.0000.5187). The participants or their legal guardians gave their informed consent to be included.

Consent for publication

Not applicable.

Competing interests

All authors declare that they have no competing interest.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Walker, C.P., Pessoa, A.L.S., Figueiredo, T. et al. Loss-of-function mutation in inositol monophosphatase 1 (IMPA1) results in abnormal synchrony in resting-state EEG. Orphanet J Rare Dis 14, 3 (2019) doi:10.1186/s13023-018-0977-1

Download citation

Keywords

  • Inosotol monophosphatase
  • IMPA1
  • Electroencephalography
  • EEG
  • Oscillations