Skip to main content

Suppressed prefrontal cortex oscillations associate with clinical pain in fibrodysplasia ossificans progressiva



Pain is a highly prevalent symptom experienced by patients across numerous rare musculoskeletal conditions. Much remains unknown regarding the central, neurobiological processes associated with clinical pain in musculoskeletal disease states. Fibrodysplasia ossificans progressiva (FOP) is an inherited condition characterized by substantial physical disability and pain. FOP arises from mutations of the bone morphogenetic protein (BMP) receptor Activin A receptor type 1 (ACVR1) causing patients to undergo painful flare-ups as well as heterotopic ossification (HO) of skeletal muscles, tendons, ligaments, and fascia. To date, the neurobiological processes that underlie pain in FOP have rarely been investigated. We examined pain and central pain mechanism in FOP as a model primary musculoskeletal condition. Central nervous system (CNS) functional properties were investigated in FOP patients (N = 17) stratified into low (0–3; 0–10 Scale) and high (≥ 4) pain cohorts using functional near-infrared spectroscopy (fNIRS). Associations among clinical pain, mental health, and physical health were also quantified using responses derived from a battery of clinical questionnaires.


Resting-state fNIRS revealed suppressed power of hemodynamic activity within the slow-5 frequency sub-band (0.01–0.027 Hz) in the prefrontal cortex in high pain FOP patients, where reduced power of slow-5, prefrontal cortex oscillations exhibited robust negative correlations with pain levels. Higher clinical pain intensities were also associated with higher magnitudes of depressive symptoms.


Our findings not only demonstrate a robust coupling among prefrontal cortex functionality and clinical pain in FOP but lays the groundwork for utilizing fNIRS to objectively monitor and central pain mechanisms in FOP and other musculoskeletal disorders.


Pain is a highly common driver of reduced quality of life across rare, musculoskeletal diseases [4, 21, 32]. Prominent and more accessible features of musculoskeletal disease such as skeletal deformities and lesions, or muscle atrophy and weakness garner extensive clinical and research attention, whereas the neurobiological cause(s) of pain have been less well-studied in this clinical domain. While there have been advances in understanding central abnormalities in musculoskeletal diseases [26, 50], neurological processes underlying pain in a chronic musculoskeletal state remain largely undefined. Demonstrating that musculoskeletal pathology yields or associates with central, pain-related alteration can serve as a neurobiological basis for implementing specific pharmacological or non-pharmacological interventions. A deeper understanding of the central manifestations of pain in rare musculoskeletal diseases, an area that has received minimal attention to date, may help uncover new targetable pain mechanisms present in musculoskeletal illnesses impacting broader populations.

Here, we aimed to determine whether and how clinical pain is centrally represented in fibrodysplasia ossificans progressiva (FOP; OMIM #135100) patients. FOP is an ultra-rare disease with an estimated 1:800,000–1:3,000,000 prevalence worldwide (285 cases in the United States) [2, 30, 38]. FOP is an inherited disorder that arises from missense mutations of the type I bone morphogenetic protein (BMP) receptor Activin A receptor type 1 (ACVR1) [43]. The mutant form of ACVR1 confers activin-A-dependent osteogenic signaling, which results in the heterotopic ossification (HO) of soft, connective tissue structures such as skeletal muscles, tendons, ligaments, and fascia [7, 13, 14], often leading to progressive immobilization and persistent pain. Prior observations indicate that the sole presence of the ACVR1 mutation is not sufficient for HO formation, but rather, exogenous events involving soft tissue injury, inflammation, or trauma are necessary to trigger the pathological ACVR1 signaling cascade and HO in FOP [42, 49].

The induction as well as the expansion of HO lesions are often preceded by soft tissue swellings termed “flare-ups” [13, 37, 46], which are characterized by local inflammation or edema, redness or warmth of skin, joint stiffness, and pain [8, 20, 36]. In the event of a flare-up, inflammatory responses alongside sustained pressure on hard or soft tissue structures are factors that can evoke high levels of pain. Interestingly, our prior analysis of data collected from a sample of approximately 100 FOP patients demonstrate that while moderate to severe pain was reported during flare-ups, sub-populations of FOP patients in a quiescent period and showing no objective signs of a flare-up also experienced an identical level of pain [34]. This finding points to the concept that pain in FOP is associated with but not entirely driven by peripheral, flare-up-related mechanisms. Additionally, FOP patients have reported experiencing recurrent severe headaches (migraine, cluster, or tension-type headaches), neuropathic pain and somatosensory abnormalities (e.g., allodynia, hyperalgesia, numbness, or tingling), which may impact or be driven by central pain processes [22, 37, 45]. Therefore, we sought to determine if modulated central mechanisms relate to clinical pain in FOP patients using resting-state, functional near-infrared spectroscopy (fNIRS). fNIRS is a non-invasive, flexible, optical neuroimaging method that quantifies the cortical hemodynamics changes based on the neuro-vascular mechanisms. fNIRS revealed for the first time a disruption in the spontaneous oscillations localized to the prefrontal cortex in FOP patients experiencing moderate to severe clinical pain.

Materials and methods

This investigation was approved by the Boston Children’s Hospital Institutional Review Board and met the Helsinki criteria for the study of human subjects. Each participant was given a detailed overview of study procedures as well as read and provided informed written consent prior to study participation. For individuals below 18 years of age, informed consent was also obtained from their parent or legal guardian.

Study participants

A total of 17 consecutive patients diagnosed with FOP between 7 and 61 years of age (11 females and 6 males) gave consent and were enrolled in the current study (Table 1). Of the 17 FOP patients, 13 individuals were confirmed carriers of the classic R206H mutation. It is estimated that 97% of FOP patients harbor a gain-of-function R206H mutation in ACVR1 [43]. One FOP patient (patient 16) had a clinical and genetic variant of FOP (personal communication with Dr. Fredrick S. Kaplan at Center for Research in FOP and Related Disorders at University of Pennsylvania), while 3 individuals had not undergone genetic testing, but possessed other FOP diagnostic features such as malformation of the big toes. Following enrollment, FOP patients provided medical history, analgesic use, demographics, clinical pain intensity rating, and fNIRS data was acquired. fNIRS data was collected in multiple settings (i.e., patients’ home or community gathering) in order to accommodate patients and families. An additional clinical pain intensity rating and completion of online clinical questionnaires occurred after the fNIRS session and within a two- to three-week period. Thus, clinical pain intensity ratings were first collected at the time of fNIRS data acquisition and subsequently, during the administration of the full battery of clinical questionnaires. Clinical questionnaires were completed at a slightly later time as some patients required special equipment or software to perform online study-related task.

Table 1 Patient characteristics and pain profiles

FOP patients harbored HO of varying degree (Table 1). While some individuals presented with known HO localized to a single site (e.g., craniofacial region or middle of the back), others possessed HO lesions in multiple upper and lower body regions. FOP patients also had varying levels of physical disability with some patients being completely ambulatory and required no assistance, while others used a cane or motorized wheelchair. Nine FOP patients utilized a device (i.e., motorized wheel chair or walking cane) to assist with mobilization. One FOP patient (Patient 1) was potentially in a flare-up state during time of fNIRS acquisition and all other patients were in a quiescent state. Patients did not report experiencing recent acute illnesses (e.g., fever or flu) or physical trauma.

Clinical questionnaires

We designated a 0–3 pain rating as low pain, while ratings ≥ 4 were classified as high pain. The rationale for choosing this criterion stems from prior work noting pain interference in patients’ daily life occurring at pain levels of 4 or more [6, 44] or that analgesic treatment is frequently sought when pain intensities are higher than 3 [11]. A battery of clinical questionnaires derived from the Patient-Reported Outcomes Measurement Information System (PROMIS; database were utilized to capture patient reported levels of (1) pain intensity (0–10 numerical rating scale), where 0 is no pain and 10 is the worst pain imaginable, (2) manifestations of pain, (3) quality of pain (i.e., sensory & affective elements), (4) physical stress, (5) physical ability, (6) psychological stress, (7) anxiety, and (8) depressive mood. Parents or guardians of pediatric FOP patients also completed the Adult Response to Children’s Symptoms (ARCS) questionnaire. The ACRS (29 items) was utilized to determine the parents’ or guardians’ behaviors in response to the child’s pain. Study participants also provided demographic data (age and gender). Each ACRS question was scored on a five-point Likert scale (Never: 0 to Always: 4). All questionnaires were administered and completed online using REDCap ( It is noted that two study participants did not fully complete questionnaires and thus the cross-correlation analysis was restricted to an N = 15 FOP patient population.

fNIRS data acquisition

While methods such as functional magnetic resonance imaging (fMRI) remain critical for characterizing CNS pain pathways in healthy subjects and clinical populations, undergoing an MRI-based assessment can be extremely difficult and, in some cases, impossible, due to the physical limitations experienced by FOP patients. Therefore, CNS hemodynamic properties (i.e., concentration changes of oxygenated hemoglobin (HbO), deoxygenated hemoglobin (HbR) and total hemoglobin) in FOP patients were measured using a multichannel continuous wave fNIRS system (CW7 System, TechEn, Milford, MA, USA) operating at 690 nm and 830 nm wavelengths. Each patient was first comfortably positioned either in a chair or in many cases, remained in their own wheelchair. Once study participants were positioned and at rest, pain intensity ratings (0–10 numerical rating scale) were obtained. Subsequently, an fNIRS cap—one for adults and one for children—consisting of 6 light sources, 8 standard-separation detectors, and 6 short-separation detectors was placed on the participant’s head (Fig. 1). Short-separation detectors were utilized to measure and mitigate any signal changes arising from non-cerebral sources (i.e., scalp or skull). Following the completion of fNIRS cap positioning (~ 20 to 30 min), each FOP patient was asked to remain still and visually fixate on a white cross displayed on computer screening. Resting-state fNIRS data were collected over a six-minute period.

Fig. 1

fNIRS channel arrangement. fNIRS optodes were distributed across two distinct cortical regions. A total of 6 light sources (red circles), 8 standard-separation detectors (blue circles), and 6 short-separation detectors (green circles) were utilized, where a source-detection pair (yellow lines) forms a channel. A standard source-detector distance of 30 mm was used, to yield an approximate light penetration depth of 30 mm, while short-separation detectors were positioned 8 mm from the light source to yield a shallow penetration depth. Channels C1-C8 covered primarily the bilateral medial and anterior prefrontal cortex (superior, middle, and inferior frontal gyri) and to a lesser extent, more lateral areas. Channels C9-C16 were positioned along the primary and secondary somatosensory cortex (pre- and postcentral gyri) on the right hemisphere. The sensitivity profiles for prefrontal and somatosensory channel arrangement for detecting CNS hemodynamics is represented on a logarithmic color scale (arbitrary units). Blue: 0 (minimum sensitivity); Red: 2 (maximum sensitivity)

Data analyses

PROMIS questionnaires: For each individual patient, the summed raw scores derived from short-form, PROMIS-based questionnaires were converted to a T-score metric. The ensuing score was represented by a T-score, a standardized score with a mean ± standard deviation (SD) of 50 ± 10. Thus, a t-score of 40 is one SD below the mean. Higher scores for symptom measures (i.e., depressive symptoms) suggest worse symptom presentation, while higher scores for functional measures (i.e., cognitive function) correspond with higher functioning. Spearman's rank correlation coefficients were calculated to explore the association between each PROMIS-based measure. A false discovery rate (FDR) adjusted p-value (α = 0.05) was subsequently calculated. Statistical comparisons between low and high pain FOP study cohorts was performed using a two-tailed t-test (α = 0.05).

fNIRS data: An open-source MATLAB (Mathworks, Natick, MA, USA) toolbox, HomER2 (, was utilized to preprocess all fNIRS datasets. Briefly, the raw optical intensity data were first transferred to changes in optical density by taking the logarithm of the signal. Optical density time courses were then visually inspected for quality assurance/control (QA/QC) across the 6-min acquisition period for any signal drifts or spikes, indicative of head motion (Additional file 1). During QA/QC of prefrontal cortex signals, 4 subject datasets were excluded due to excessive head motion or low signal to noise ratio, yielding 13 subject datasets for power-spectral analysis (Pain Level: 0–3; N = 7 & Pain Level: 4–10; N = 6). Separate QA/QC of somatosensory cortex signal resulted in 10 viable datasets (Pain Level: 0–3; N = 5 & Pain Level: 4–10; N = 5). Filtered (0.01–0.5 Hz) optical density time courses were further converted to concentration changes in HbO, HbR and total hemoglobin using the modified Beer-Lambert Law with a partial pathlength factor of 6. Finally, a linear regression model was employed to regress out superficial noises from the time series of each cortical fNIRS channel (i.e., having a distance of 3 cm) by including the hemoglobin time course of the short separation channel (i.e. having a distance of 8 mm) with the highest correlation as a covariate in the linear model. A constant regressor with the value of 1 was also added to the model.

Power spectral analysis of low-frequency oscillations on the recorded fNIRS hemodynamic signals was performed to investigate the changes in neuronal oscillatory dynamics with pain. Here, the HbO data were utilized rather than HbR data, based on the availability of higher signal-to-noise ratio of HbO compared to HbR time courses [23]. For each patient, the 6-min HbO time courses of each channel were first transformed to the frequency domain using the fast Fourier transform. The power of each frequency component was obtained by taking the square of its absolute amplitude value. Furthermore, prior work has shown that cortical oscillations occur in distinct frequency bands, where each range is associated with unique physiological functions [3]. Thus, we further subdivided the low-frequency oscillations range into five sub-bands: slow-5 (0.01–0.027 Hz), slow-4 (0.027–0.073 Hz), slow-3 (0.073–0.198 Hz), slow-2 (0.198–0.25 Hz), and slow-1 (0.25–0.5 Hz). For each sub-band, we conducted two-sample t-tests on the mean power values of each normal fNIRS channel between the high pain group and low pain group to test the null hypothesis that patients with high pain levels do not have a statistically significant reduction of the sub-band power of a particular channel (and therefore the underlying brain region). Channel-wise FDR correction was then applied to control the number of false positives. FDR adjusted (α = 0.05) p-values were utilized to determine significance. Resting-state functional connectivity was also explored (Additional file 1).


Pain, mental health and physical health

Of the 17 enrolled FOP patients (Table 1), 11 individuals (6 males, 5 females) were categorized in the low pain (Pain Level: 0–3/10) cohort at the time of fNIRS data acquisition, while 6 female patients were allotted to the high pain (Pain Level: ≥ 4) group. In order to compare patients’ gender with low pain vs. high pain, the Fisher’s Exact test was utilized. We observed a significant association between gender and pain (p = 0.043). Pain ratings between participants less than 18 years of age vs. adult patients was moderately significant (p = 0.046). A comparison of pain levels between the first (fNIRS data acquisition) and second (clinical questionnaire completion) evaluation time points demonstrated that while pain intensity remained relatively stable for many FOP patients, others showed increases or decreases in pain (Fig. 2). Moreover, no relationship among gender and pain was present at the second clinical pain intensity assessment (p = 0.14). Similarly, a significant difference between FOP patients less that 18 years compared to individuals over the age of 18 was not observed at the subsequent time point (p = 0.18). Analgesics, including prescription opioids, were used by both study populations with some patients utilizing treatments on an intermittent or as needed basis, while others had more consistent usage. A Fisher’s Exact test was utilized to assess the association between pain (≥ 4 vs < 4) and analgesic usage (Any usage = 1 vs. Otherwise = 0). No association between pain and analgesics usage was observed with the p-value being highly non-significant (p = 1.0). Clinical questionnaire data acquired from FOP patients showed that individuals in the high pain cohort (N = 7) relative to the low pain group (N = 9) had altered emotional functioning in domains such as anxiety in conjunction with physical symptoms, mainly physical stress (Table 2). The interactions between pain intensity, the quality of pain (affective vs. sensory), mental health and physical health in FOP, all of which were self-reported assessments, were also characterized (Fig. 3). Of note, the intensity of pain was significantly correlated magnitude of self-reported depressive symptoms (r = 0.74, p = 0.0016) in addition to the magnitude physical stress experienced (r = 0.65, p = 0.009). Interestingly, significant relationships were further observed between level of anxiety experienced by patients, sensory aspect of pain and physical stress, suggesting an interaction among these three domains in FOP. Finally, across pediatric FOP patients, the mean ± SE ARCS questionnaire response was 2.31 ± 0.26. Low ARCS scores were likely reflective of low pain levels experienced by pediatric FOP patients.

Fig. 2

Dynamics of pain in FOP. Pain levels in FOP patients were quantified at two time points approximately 2–3 weeks apart (N = 15). Patients 4 and 15 were lost to follow-up. In some patients, pain levels were constant between the two assessments, while for others, there were increases or decreases in self-reported pain intensity. Median self-reported pain intensities were 3 and 4 for fNIRS and clinical questionnaire sessions, respectively

Table 2 Group-level PROMIS scores
Fig. 3

Associations among clinical pain, mental health, and physical health in FOP. a The cross-correlation matrix demonstrates the relationships or lack thereof between measures of pain, mental health and physical health as measured by questionnaires derived from the PROMIS database. Dotted white boxes denote significant effects surpassing multiple comparison correction (FDR-pcorr < 0.05). Correlation analysis was performed using data collected from 15 FOP patients. b Example correlation plots depicting significant relationships among measures of pain, mental health, and physical health. Within the top and bottom plots, 2 patients had overlapping data points. Data from 15 FOP patients are shown

Power spectral analysis (Low vs. high pain FOP patients)

As shown in Fig. 4a, relative to the low pain cohort (N = 7), the amplitude of low-frequency oscillations localized to the prefrontal cortex and within the slow-5 (0.01–0.027 Hz) sub-band were significantly lower in high pain FOP patients (N = 6). A robust suppression of slow-5 power was observed in the high pain group across medial (channels C2 (FDR-pcorr = 0.012) and C7 (FDR-pcorr = 0.013)) as well as more lateral (channel C5 (FDR-pcorr = 0.015)) components of the prefrontal cortex (Cohen’s d: C2 = 1.51; C5 = 1.45; C7 = 1.46). Across all study participants, self-reported pain levels were negatively correlated with prefrontal cortex, slow-5 power (Fig. 4b). Pain score obtained from the second pain evaluation were also negatively correlated with the magnitude of prefrontal cortex, slow-5 power (C2: r = − 0.60, p = 0.039; C5: r = − 0.28, p = 0.39; C7: r = − 0.73, p = 0.007). A significant correlation was not detected between slow-5 power and severity of depressive symptoms (r = − 0.36 to − 0.41; p > 0.05). The age of the patient did not correlate with the magnitude of prefrontal cortex, slow-5 power (C2: r = − 0.06, p = 0.84; C5: r = 0.16, p = 0.60; C7: r = − 0.27, p = 0.37). As low-frequency oscillations at the slow-5 range are involved with long-range interaction of brain regions, these finding suggests that FOP patients experiencing higher levels of pain likely possess a greater magnitude of dysfunctional long-range communication involving the prefrontal cortex [3, 5, 40]. (See also resting-state functional connectivity analysis results in Additional file 1). In other words, higher self-reported pain in FOP was associated with greater disruption of large-scale CNS networks encompassing the prefrontal cortex. Figure 5 demonstrates that slow-5 amplitude alterations were not significantly different over the primary and secondary somatosensory cortices (pre- and post-central gyrus) when comparing low (N = 5) vs. high pain (N = 5) FOP patients. Calculation of effect sizes (Cohen’s d (Range): 0.03–0.79) indicated a low to moderate pain-dependent effect on slow-5 power in somatosensory cortices. A correlation between clinical pain intensity ratings and slow-5 power for any somatosensory channel, channels C9-C16, was not observed (r = − 0.04 to 0.23; FDR-adjusted p > 0.05). Significant differences between low and high pain FOP patients were not detected in other frequency sub-bands (i.e., slow-1 to slow-4) and in either prefrontal or somatosensory cortices.

Fig. 4

Decreased prefrontal cortex slow-5 (0.01–0.027 Hz) power in FOP patients with high pain levels. a Frequency analyses reveal decrease power within the slow-5 sub-band across multiple channels distributed over the medial and anterior aspects of the prefrontal cortex. A significant effect after FDR correction was observed in channels C2, C5, and C7. * denotes FDR correct p-values, while light green bars in C2, C5, and C7 plots further define the slow-5 sub-band. Mean power spectral density and standard error (shaded area) across all sub-bands (i.e., slow-1 to slow-5) are plotted for low (Pain Level: 0–3; N = 7) and high (Pain Level: 4–10; N = 6) clinical pain FOP cohorts. Clinical pain intensity levels for each patient were collected from each FOP patient just prior to acquisition of the six-minute resting-state fNIRS scan. b Correlation analyses performed across all FOP patients showed that for channels C2, C5, and C7, a greater extent of power suppression within the slow-5 frequency range corresponded with greater clinical pain intensity

Fig. 5

No change in primary and secondary somatosensory cortex slow-5 (0.01–0.027 Hz) power (low pain vs. high pain FOP patients). Frequency analyses reveal no significant differences in power within the slow-5 sub-band across all channels distributed over the primary and secondary somatosensory cortex. Calculation of effect sizes (Cohen’s d (Range): 0.03–0.79) demonstrated a low to moderate pain-dependent effect on slow-5 power in somatosensory cortices. Mean power spectral density and standard error (shaded area) across all sub-bands (i.e., slow-1 to slow-5) are plotted for low clinical pain (Pain Level: 0–3; N = 5) and high clinical pain (Pain Level: 4–10; N = 5) FOP cohorts. Clinical pain intensity levels from each FOP patient were collected just prior to acquisition of a six-minute resting-state fNIRS scan. Significant correlation was not detected between clinical pain intensity ratings and slow-5 power (Correlation (r): − 0.04 to 0.23). While a trend towards higher (relative to the low pain FOP group) power in the slow-4 sub-band (Channel C10) for the high pain FOP cohort was observed, a significant intergroup difference was not present (p = 0.137)


This work highlights the central representation of a critical element of FOP, namely, pain. We report that the pain experienced by FOP patients was integrated with an individual’s emotional health, which points to central involvement in FOP pain. By utilizing fNIRS, we demonstrated that the severity of clinical pain was negatively associated with the power of low-frequency oscillations localized to the prefrontal cortex, a region implicated in emotional appraisal [27], reward mechanisms [29], and clinical conditions including major depressive disorder [18]. The observed loss of the prefrontal cortex, low-frequency activity within the slow-5 sub-band, strongly points to the disintegration of large-scale CNS networks in a high pain, FOP state. The current findings revealed fNIRS's utility for characterizing central pain mechanisms in FOP and suggests that this approach might be generalized for discovering central features of clinical pain in other musculoskeletal diseases.

Peripheral pathological features (nerve compression) or events (flare-ups) can cause patients to experience pain. However, across rare and also more common musculoskeletal diseases, there is a frequent discrepancy between objective measures of disease burden or peripheral trauma and patient-reported pain levels [21]. Finding such as those centered around prefrontal cortex functionality in FOP pain patients in conjunction with the elucidation of diminished emotional health may help explain this discordance between pain and tissue pathology, and better direct the implementation of pharmacological and non-pharmacological pain treatment interventions. Therefore, while the current set of results are specific to FOP, there are implications that extend beyond this rare disease and into other clinical conditions characterized by pain and musculoskeletal phenotypes.

In FOP, HO of musculoskeletal tissue causes long-term, severe disability. However, alongside the physical limitations that many FOP patients undergo and which can begin during childhood stages of life [10, 33, 36], the core pathological elements of FOP such as HO induction and growth may lead to other downstream manifestations [34]. Pain in FOP, is most often associated with recurrent soft tissue swelling or flare-ups, which in many circumstances precedes new HO lesions or expansion of existing ones [36, 37, 46]. Yet, as previous work [34] as well as the current study have revealed, clinical pain is frequently reported by FOP patients during what are considered ‘quiescent’ or non-flare-up states.

We found that pain in FOP can be variable, both across subjects and time. Monitoring patients across a 2–3-week period demonstrated intra-subject fluctuations in pain intensity with some individuals reporting increased levels of pain during this interval, while others showed little to no change in pain intensity. The site of pain often corresponded to the location of known HO lesions; however, there were cases where pain was altogether absent, or was experienced only across a subset of HO lesion sites. A pain level ≥ 4 was moderately associated with gender. A trend of higher pain levels or higher chronic pain prevalence in female FOP patients is in accord with the large body of data demonstrating higher perception of pain and more frequent occurrence of chronic pain conditions in women [31, 47]. Additionally, slightly higher pain levels were measured in older FOP patients relative to younger patients; an observation that was previously observed in other rare musculoskeletal conditions [21]. From a battery of clinical questionnaires, the perception of pain in individuals reporting a high level of pain consisted of both sensory and affective elements. Cross-correlation analyses further indicated the interaction between severity of pain, pain quality, physical health, and emotional health; however, two key features stemming from this approach were derived. First, there was clear evidence that the physical (i.e., physical stress) and psychological (i.e., depressive mood or anxiety) challenges that FOP patients harbor are tightly integrated. Second, there was a robust association between higher levels of clinical pain and greater depressive symptoms. While higher pain severity often coincides with a greater depressive symptomology in many clinical conditions [16], the observation of this relationship in a FOP population has additional implications. Pain in FOP can of course be driven by peripheral pathological events, for example, the soft tissue edema or inflammation occurring during a flare-up episode. However, superimposed upon peripheral, pain-inducing occurrences, are abnormalities anchored within the CNS, which can result in a top-down dysregulation of pain perception and relatedly, heightened pain sensitivity or maintain a persistent state of pain [17, 24].

Implementation of fNIRS enabled an objective assessment of CNS function in FOP patients and equally important, an understanding of whether and how cortical activity relates to subjective accounts of pain in this disease. Patients with moderate to severe pain were characterized by a reduced amplitude of slow-5 (0.01–0.027 Hz), low-frequency oscillations occurring within the prefrontal cortex, where higher clinical pain severity negatively correlated with more suppressed prefrontal cortex fluctuations. Using resting-state functional magnetic resonance imaging (fMRI), highly similar frequency-dependent changes have been reported in multiple chronic pain conditions. In trigeminal neuralgia [52], postherpetic neuralgia [12], and chronic back pain [28, 51] populations, a loss of slow-5 power was not only noted in prefrontal cortex [52] but also in subcortical regions such as the nucleus accumbens, which share projections with sub-regions of the frontal cortex areas amongst other (meso-)limbic structures. Thus, in general, the suppression of the frontal cortex (potentially along the frontostriatal pathway) low-frequency oscillation in the slow-5 sub-band may represent a maladaptive cortical process that is derived from persistent pain or mechanisms that sustain a pain state. Further investigation is necessary to determine if and how there is a causal link among suppressed frontal cortex oscillations and chronic pain in FOP.

The alterations in low-frequency dynamics have been previously linked with changes in functional connectivity [1]. Slow-5 oscillations are more closely associated with spontaneous neuronal oscillation intrinsic to cortical regions, including the prefrontal cortex, while higher frequency activity (e.g., slow-4 sub-band: 0.027–0.073 Hz) involves cortical and subcortical regions such as the basal ganglia [39, 48]. Slow-5 oscillations enable long-range connectivity, particularly that which is necessary for proper communication amongst CNS hubs (i.e., prefrontal cortex) embedded within large-scale networks (i.e., default-mode or mesolimbic network). A disruption of mechanisms that facilitate long-range CNS communication may underpin the decreased prefrontal, and somatosensory cortex functional connectivity observed in FOP patients with higher self-reported clinical pain levels. The somatosensory cortex encodes the sensory dimension of the pain [15], while the prefrontal cortex with its connections to the pain modulation centers (periaqueductal grey) is known to exhibit endogenous antinociceptive effects. A decrease in the connectivity of the somatosensory cortex and prefrontal cortex may imply an imbalance in the pain encoding process and the ensuing descending control of pain leading to prolonged pain status. Interestingly, relative to short-range connectivity, there is also a significantly stronger coupling between long-range connectivity and CNS metabolic properties (cerebral blood flow, cerebral metabolic rate for oxygen or glucose [25], possibly suggesting that cerebral metabolic alterations may especially perturbate slow-5 or slow-4 low-frequency oscillations and in parallel, connectivity within large-scale networks.

We project that the frequency-dependent and pain-dependent effects observed in our FOP cohort of patients are a collective and cyclical outcome of living with a chronic illness such as FOP, withstanding recurrent and severely painful flare-up episodes as well as enduring persistent pain over several months or even years. Nonetheless, the impact of aberrant activin-A- or ACVR1/ALK2-dependent activity with the CNS must be considered. In FOP patients [19, 41] lesions or focal MRI hyper-intensities along white matter pathways have been noted. Disruption of white matter integrity can indeed modulate CNS functional properties.

This investigation represents the first functional CNS imaging study in FOP and a novel application of fNIRS. fNIRS was implemented as an alternate approach to probe CNS function, as MRI-based assessments are not easily feasible for many FOP patients. This study, in conjunction with other recent efforts, demonstrates fNIRS's utility towards elucidating central pain mechanisms [9, 35]. However, limitations are noted. fNIRS provides an assessment of cortical (dys-)function, where alterations in the brainstem or subcortical structures and their association with clinical pain cannot be elucidated. Compared to fMRI, fNIRS offers flexibility in terms of implementation, allowing bedside measurement to be made. However, a low signal to noise ratio at specific optode locations, possibly due to hair, was a limiting factor. Some FOP patients presented with HO in the neck, shoulder, and upper back region, which anatomically altered their head position and caused difficultly in fNIRS positioning of optodes, particularly over the somatosensory cortices. The use of wireless fNIRS setups or implementation of advanced light source technology may facilitate improved fNIRS data quality throughout the brain. Additionally, in the current cohort of FOP patients, those categorized in the high pain group were primarily women, which is in line with prior FOP studies [22]. Nonetheless, future investigations may aim to enroll more male FOP patients with moderate to severe pain. Lastly, some patients continued analgesic treatment during the study, but the varied medication usage between patients makes this unlikely.


In conclusion, fNIRS was applied to characterize cortical function in FOP patients in low and high pain states. Suppressed power within the slow-5 frequency band of prefrontal cortex in FOP patients was associated with greater clinical pain. We also have demonstrated the interaction amongst pain, mental health and physical health in FOP. This investigation provides further rationale for utilizing fNIRS to further elucidate mechanisms of pain modulation in FOP and other musculoskeletal disorders.

Availability of data and materials

The datasets generated and/or analyzed during the current study are not publicly available due institutional regulations, but are available from the corresponding author on reasonable request.



Activin A receptor type 1


Adult Response to Children’s Symptoms


Bone morphogenetic protein


Central nervous system


Fibrodysplasia ossificans progressiva


Functional near-infrared spectroscopy


Heterotopic ossification


Oxygenated hemoglobin


Deoxygenated hemoglobin


Patient-Reported Outcomes Measurement Information System


Quality assurance/quality control


  1. 1.

    Baliki MN, Baria AT, Apkarian AV. The cortical rhythms of chronic back pain. J Neurosci. 2011;31(39):13981–90.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  2. 2.

    Baujat G, Choquet R, Bouee S, Jeanbat V, Courouve L, Ruel A, Michot C, Le Quan Sang KH, Lapidus D, Messiaen C, Landais P, Cormier-Daire V. Prevalence of fibrodysplasia ossificans progressiva (FOP) in France: an estimate based on a record linkage of two national databases. Orphanet J Rare Dis. 2017;12(1):123.

    PubMed  PubMed Central  Article  Google Scholar 

  3. 3.

    Buzsaki G, Draguhn A. Neuronal oscillations in cortical networks. Science. 2004;304(5679):1926–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  4. 4.

    Castori M, Camerota F, Celletti C, Danese C, Santilli V, Saraceni VM, Grammatico P. Natural history and manifestations of the hypermobility type Ehlers-Danlos syndrome: a pilot study on 21 patients. Am J Med Genet A. 2010;152A(3):556–64.

    PubMed  Article  PubMed Central  Google Scholar 

  5. 5.

    Cohen AD, Tomasi D, Shokri-Kojori E, Nencka AS, Wang Y. Functional connectivity density mapping: comparing multiband and conventional EPI protocols. Brain Imaging Behav. 2018;12(3):848–59.

    PubMed  Article  PubMed Central  Google Scholar 

  6. 6.

    Daut RL, Cleeland CS. The prevalence and severity of pain in cancer. Cancer. 1982;50(9):1913–8.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  7. 7.

    Dey D, Bagarova J, Hatsell SJ, Armstrong KA, Huang L, Ermann J, Vonner AJ, Shen Y, Mohedas AH, Lee A, Eekhoff EM, van Schie A, Demay MB, Keller C, Wagers AJ, Economides AN, Yu PB. Two tissue-resident progenitor lineages drive distinct phenotypes of heterotopic ossification. Sci Transl Med. 2016;8(366):366ra163.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  8. 8.

    Eekhoff EMW, Botman E, Coen Netelenbos J, de Graaf P, Bravenboer N, Micha D, Pals G, de Vries TJ, Schoenmaker T, Hoebink M, Lammertsma AA, Raijmakers P. [18F]NaF PET/CT scan as an early marker of heterotopic ossification in fibrodysplasia ossificans progressiva. Bone. 2018;109:143–6.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  9. 9.

    Fernandez Rojas R, Huang X, Ou KL. A Machine Learning Approach for the Identification of a Biomarker of Human Pain using fNIRS. Sci Rep. 2019;9(1):5645.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  10. 10.

    Gencer-Atalay K, Ozturk EC, Yagci I, Ata P, Delil K, Ozgen Z, Akyuz G. Challenges in the treatment of fibrodysplasia ossificans progressiva. Rheumatol Int. 2019;39(3):569–76.

    PubMed  Article  PubMed Central  Google Scholar 

  11. 11.

    Gilron I, Bailey JM, Tu D, Holden RR, Weaver DF, Houlden RL. Morphine, gabapentin, or their combination for neuropathic pain. N Engl J Med. 2005;352(13):1324–34.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  12. 12.

    Gu L, Hong S, Jiang J, Liu J, Cao X, Huang Q, Zeng X, Zhou F, Zhang D. Bidirectional alterations in ALFF across slow-5 and slow-4 frequencies in the brains of postherpetic neuralgia patients. J Pain Res. 2019;12:39–47.

    PubMed  Article  PubMed Central  Google Scholar 

  13. 13.

    Hatsell SJ, Idone V, Wolken DM, Huang L, Kim HJ, Wang L, Wen X, Nannuru KC, Jimenez J, Xie L, Das N, Makhoul G, Chernomorsky R, D’Ambrosio D, Corpina RA, Schoenherr CJ, Feeley K, Yu PB, Yancopoulos GD, Murphy AJ, Economides AN. ACVR1R206H receptor mutation causes fibrodysplasia ossificans progressiva by imparting responsiveness to activin A. Sci Transl Med. 2015;7(303):303ra137.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  14. 14.

    Hino K, Ikeya M, Horigome K, Matsumoto Y, Ebise H, Nishio M, Sekiguchi K, Shibata M, Nagata S, Matsuda S, Toguchida J. Neofunction of ACVR1 in fibrodysplasia ossificans progressiva. Proc Natl Acad Sci USA. 2015;112(50):15438–43.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  15. 15.

    Hofbauer RK, Rainville P, Duncan GH, Bushnell MC. Cortical representation of the sensory dimension of pain. J Neurophysiol. 2001;86(1):402–11.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  16. 16.

    IsHak WW, Wen RY, Naghdechi L, Vanle B, Dang J, Knosp M, Dascal J, Marcia L, Gohar Y, Eskander L, Yadegar J, Hanna S, Sadek A, Aguilar-Hernandez L, Danovitch I, Louy C. Pain and depression: a systematic review. Harv Rev Psychiatry. 2018;26(6):352–63.

    PubMed  Article  PubMed Central  Google Scholar 

  17. 17.

    Ji RR, Nackley A, Huh Y, Terrando N, Maixner W. Neuroinflammation and central sensitization in chronic and widespread pain. Anesthesiology. 2018;129(2):343–66.

    PubMed  PubMed Central  Article  Google Scholar 

  18. 18.

    Kaiser RH, Whitfield-Gabrieli S, Dillon DG, Goer F, Beltzer M, Minkel J, Smoski M, Dichter G, Pizzagalli DA. Dynamic resting-state functional connectivity in major depression. Neuropsychopharmacology. 2016;41(7):1822–30.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  19. 19.

    Kan L, Kitterman JA, Procissi D, Chakkalakal S, Peng CY, McGuire TL, Goldsby RE, Pignolo RJ, Shore EM, Kaplan FS, Kessler JA. CNS demyelination in fibrodysplasia ossificans progressiva. J Neurol. 2012;259(12):2644–55.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  20. 20.

    Kaplan FS, Al Mukaddam M, Pignolo RJ. Acute unilateral hip pain in fibrodysplasia ossificans progressiva (FOP). Bone. 2018;109:115–9.

    PubMed  Article  PubMed Central  Google Scholar 

  21. 21.

    Kelly MH, Brillante B, Collins MT. Pain in fibrous dysplasia of bone: age-related changes and the anatomical distribution of skeletal lesions. Osteoporos Int. 2008;19(1):57–63.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  22. 22.

    Kitterman JA, Strober JB, Kan L, Rocke DM, Cali A, Peeper J, Snow J, Delai PL, Morhart R, Pignolo RJ, Shore EM, Kaplan FS. Neurological symptoms in individuals with fibrodysplasia ossificans progressiva. J Neurol. 2012;259(12):2636–43.

    PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Lachert P, Janusek D, Pulawski P, Liebert A, Milej D, Blinowska KJ. Coupling of Oxy- and Deoxyhemoglobin concentrations with EEG rhythms during motor task. Sci Rep. 2017;7(1):15414.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  24. 24.

    Latremoliere A, Woolf CJ. Central sensitization: a generator of pain hypersensitivity by central neural plasticity. J Pain. 2009;10(9):895–926.

    PubMed  PubMed Central  Article  Google Scholar 

  25. 25.

    Liang X, Zou Q, He Y, Yang Y. Coupling of functional connectivity and regional cerebral blood flow reveals a physiological basis for network hubs of the human brain. Proc Natl Acad Sci USA. 2013;110(5):1929–34.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  26. 26.

    Loggia ML, Chonde DB, Akeju O, Arabasz G, Catana C, Edwards RR, Hill E, Hsu S, Izquierdo-Garcia D, Ji RR, Riley M, Wasan AD, Zurcher NR, Albrecht DS, Vangel MG, Rosen BR, Napadow V, Hooker JM. Evidence for brain glial activation in chronic pain patients. Brain. 2015;138(Pt 3):604–15.

    PubMed  PubMed Central  Article  Google Scholar 

  27. 27.

    Ma ST, Abelson JL, Okada G, Taylor SF, Liberzon I. Neural circuitry of emotion regulation: Effects of appraisal, attention, and cortisol administration. Cogn Affect Behav Neurosci. 2017;17(2):437–51.

    PubMed  Article  PubMed Central  Google Scholar 

  28. 28.

    Makary MM, Polosecki P, Cecchi GA, DeAraujo IE, Barron DS, Constable TR, Whang PG, Thomas DA, Mowafi H, Small DM, Geha P. Loss of nucleus accumbens low-frequency fluctuations is a signature of chronic pain. Proc Natl Acad Sci USA. 2020;117(18):10015–23.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  29. 29.

    Malvaez M, Shieh C, Murphy MD, Greenfield VY, Wassum KM. Distinct cortical-amygdala projections drive reward value encoding and retrieval. Nat Neurosci. 2019;22(5):762–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  30. 30.

    Morales-Piga A, Kaplan FS. Osteochondral diseases and fibrodysplasia ossificans progressiva. Adv Exp Med Biol. 2010;686:335–48.

    PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    Naylor JC, Ryan Wagner H, Brancu M, Shepherd-Banigan M, Elbogen E, Kelley M, Fecteau T, Goldstein K, Kimbrel NA, Marx CE, Group VAM-AMW, Group VAM-AMWVW, Strauss JL. Self-reported pain in male and female Iraq/Afghanistan-era veterans: associations with psychiatric symptoms and functioning. Pain Med. 2017;18(9):1658–67.

    PubMed  PubMed Central  Google Scholar 

  32. 32.

    Nghiem T, Louli J, Treherne SC, Anderson CE, Tsimicalis A, Lalloo C, Stinson JN, Thorstad K. Pain experiences of children and adolescents with osteogenesis imperfecta: an integrative review. Clin J Pain. 2017;33(3):271–80.

    PubMed  PubMed Central  Google Scholar 

  33. 33.

    Ortiz-Agapito F, Colmenares-Bonilla D. Quality of life of patients with Fibrodysplasia ossificans progressiva. J Child Orthop. 2015;9(6):489–93.

    PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    Peng K, Cheung K, Lee A, Sieberg C, Borsook D, Upadhyay J. Longitudinal evaluation of pain, flare-up, and emotional health in Fibrodysplasia ossificans progressiva: analyses of the international FOP Registry. JBMR Plus. 2019;3(8):e10181.

    PubMed  PubMed Central  Article  Google Scholar 

  35. 35.

    Peng K, Yucel MA, Steele SC, Bittner EA, Aasted CM, Hoeft MA, Lee A, George EE, Boas DA, Becerra L, Borsook D. Morphine attenuates fNIRS signal associated with painful stimuli in the medial Frontopolar cortex (medial BA 10). Front Hum Neurosci. 2018;12:394.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  36. 36.

    Pignolo RJ, Baujat G, Brown MA, De Cunto C, Di Rocco M, Hsiao EC, Keen R, Al Mukaddam M, Sang KLQ, Wilson A, White B, Grogan DR, Kaplan FS. Natural history of fibrodysplasia ossificans progressiva: cross-sectional analysis of annotated baseline phenotypes. Orphanet J Rare Dis. 2019;14(1):98.

    PubMed  PubMed Central  Article  Google Scholar 

  37. 37.

    Pignolo RJ, Bedford-Gay C, Liljesthrom M, Durbin-Johnson BP, Shore EM, Rocke DM, Kaplan FS. The natural history of flare-ups in Fibrodysplasia Ossificans Progressiva (FOP): a comprehensive global assessment. J Bone Miner Res. 2016;31(3):650–6.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  38. 38.

    Pignolo RJ, Cheung K, Kile S, Fitzpatrick MA, De Cunto C, Al Mukaddam M, Hsiao EC, Baujat G, Delai P, Eekhoff EMW, Di Rocco M, Grunwald Z, Haga N, Keen R, Levi B, Morhart R, Scott C, Sherman A, Zhang K, Kaplan FS. Self-reported baseline phenotypes from the International Fibrodysplasia Ossificans Progressiva (FOP) Association Global Registry. Bone. 2020;134:115274.

    PubMed  Article  PubMed Central  Google Scholar 

  39. 39.

    Ren P, Lo RY, Chapman BP, Mapstone M, Porsteinsson A, Lin F. Alzheimer’s disease neuroimaging I. Longitudinal alteration of intrinsic brain activity in the striatum in mild cognitive impairment. J Alzheimers Dis. 2016;54(1):69–78.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  40. 40.

    Salvador R, Suckling J, Schwarzbauer C, Bullmore E. Undirected graphs of frequency-dependent functional connectivity in whole brain networks. Philos Trans R Soc Lond B Biol Sci. 2005;360(1457):937–46.

    PubMed  PubMed Central  Article  Google Scholar 

  41. 41.

    Severino M, Bertamino M, Tortora D, Morana G, Uccella S, Bocciardi R, Ravazzolo R, Rossi A, Di Rocco M. Novel asymptomatic CNS findings in patients with ACVR1/ALK2 mutations causing fibrodysplasia ossificans progressiva. J Med Genet. 2016;53(12):859–64.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  42. 42.

    Shore EM, Kaplan FS. Inherited human diseases of heterotopic bone formation. Nat Rev Rheumatol. 2010;6(9):518–27.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  43. 43.

    Shore EM, Xu M, Feldman GJ, Fenstermacher DA, Cho TJ, Choi IH, Connor JM, Delai P, Glaser DL, LeMerrer M, Morhart R, Rogers JG, Smith R, Triffitt JT, Urtizberea JA, Zasloff M, Brown MA, Kaplan FS. A recurrent mutation in the BMP type I receptor ACVR1 causes inherited and sporadic fibrodysplasia ossificans progressiva. Nat Genet. 2006;38(5):525–7.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  44. 44.

    Twycross R, Harcourt J, Bergl S. A survey of pain in patients with advanced cancer. J Pain Symptom Manag. 1996;12(5):273–82.

    CAS  Article  Google Scholar 

  45. 45.

    Upadhyay J, Geber C, Hargreaves R, Birklein F, Borsook D. A critical evaluation of validity and utility of translational imaging in pain and analgesia: Utilizing functional imaging to enhance the process. Neurosci Biobehav Rev. 2018;84:407–23.

    PubMed  Article  PubMed Central  Google Scholar 

  46. 46.

    Upadhyay J, Xie L, Huang L, Das N, Stewart RC, Lyon MC, Palmer K, Rajamani S, Graul C, Lobo M, Wellman TJ, Soares EJ, Silva MD, Hesterman J, Wang L, Wen X, Qian X, Nannuru K, Idone V, Murphy AJ, Economides AN, Hatsell SJ. The expansion of heterotopic bone in Fibrodysplasia ossificans progressiva is activin A-dependent. J Bone Miner Res. 2017;32(12):2489–99.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  47. 47.

    van Hecke O, Torrance N, Smith BH. Chronic pain epidemiology and its clinical relevance. Br J Anaesth. 2013;111(1):13–8.

    PubMed  Article  PubMed Central  Google Scholar 

  48. 48.

    Wang L, Kong Q, Li K, Su Y, Zeng Y, Zhang Q, Dai W, Xia M, Wang G, Jin Z, Yu X, Si T. Frequency-dependent changes in amplitude of low-frequency oscillations in depression: a resting-state fMRI study. Neurosci Lett. 2016;614:105–11.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  49. 49.

    Yu PB, Deng DY, Lai CS, Hong CC, Cuny GD, Bouxsein ML, Hong DW, McManus PM, Katagiri T, Sachidanandan C, Kamiya N, Fukuda T, Mishina Y, Peterson RT, Bloch KD. BMP type I receptor inhibition reduces heterotopic [corrected] ossification. Nat Med. 2008;14(12):1363–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  50. 50.

    Zhang A, Lee YC. Mechanisms for joint pain in rheumatoid arthritis (RA): from cytokines to central sensitization. Curr Osteoporos Rep. 2018;16(5):603–10.

    PubMed  PubMed Central  Article  Google Scholar 

  51. 51.

    Zhang B, Jung M, Tu Y, Gollub R, Lang C, Ortiz A, Park J, Wilson G, Gerber J, Mawla I, Chan ST, Wasan A, Edwards R, Lee J, Napadow V, Kaptchuk T, Rosen B, Kong J. Identifying brain regions associated with the neuropathology of chronic low back pain: a resting-state amplitude of low-frequency fluctuation study. Br J Anaesth. 2019;123(2):e303–11.

    PubMed  PubMed Central  Article  Google Scholar 

  52. 52.

    Zhang Y, Mao Z, Pan L, Ling Z, Liu X, Zhang J, Yu X. Frequency-specific alterations in cortical rhythms and functional connectivity in trigeminal neuralgia. Brain Imaging Behav. 2019;13(6):1497–509.

    PubMed  Article  PubMed Central  Google Scholar 

Download references


The authors would like to thank all FOP patients and their families for their involvement and patience during the course of the study. The authors also thank the International FOP Association ( for their assistance during patient outreach, Dr. Kin Cheung for input on statistical analyses, and Dr. Seward Rutkove for helpful comments during the preparation of this manuscript.


This study was supported by the Department of Anesthesiology, Critical Care and Pain Medicine, Boston Children’s Hospital. JU receives supported from MAYDAY Fund, Orphan Disease Center Grant Program, and RF1 MH122967-01. PBY receives support from NIH-NIAMS R01-AR057374.

Author information




KP, KDP, RL, AL, PBY and JU designed the study. KP, KDP, RL, and JU acquired the data. KP, KDP, MV, CC, and JU analyzed data. All authors contributed to data interpretation. KP, KDP, MV, PBY, and JU wrote and edited the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jaymin Upadhyay.

Ethics declarations

Ethics approval and consent to participate

This work was approved by the Boston Children’s Hospital, Institutional Review Board. All participants provided informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1

. QA/QC and functional connectivity analyses of fNIRS data.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Peng, K., Karunakaran, K.D., Labadie, R. et al. Suppressed prefrontal cortex oscillations associate with clinical pain in fibrodysplasia ossificans progressiva. Orphanet J Rare Dis 16, 54 (2021).

Download citation


  • Pain
  • Central nervous system
  • Fibrodysplasia ossificans progressiva
  • Functional near-infrared spectroscopy