Suppressed prefrontal cortex oscillations associate with clinical pain in fibrodysplasia ossificans progressiva
Orphanet Journal of Rare Diseases volume 16, Article number: 54 (2021)
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) . 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 . 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.
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 . 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.
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.
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 . A battery of clinical questionnaires derived from the Patient-Reported Outcomes Measurement Information System (PROMIS; http://www.healthmeasures.net) 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 (https://www.project-redcap.org). 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.
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 (https://homer-fnirs.org), 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 . 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 . 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.
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.
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 , reward mechanisms , and clinical conditions including major depressive disorder . 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 . 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 . 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  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 . 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 , 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 , postherpetic neuralgia , and chronic back pain [28, 51] populations, a loss of slow-5 power was not only noted in prefrontal cortex  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 . 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 , 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 , 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 . 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
Patient-Reported Outcomes Measurement Information System
Quality assurance/quality control
Baliki MN, Baria AT, Apkarian AV. The cortical rhythms of chronic back pain. J Neurosci. 2011;31(39):13981–90.
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.
Buzsaki G, Draguhn A. Neuronal oscillations in cortical networks. Science. 2004;304(5679):1926–9.
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.
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.
Daut RL, Cleeland CS. The prevalence and severity of pain in cancer. Cancer. 1982;50(9):1913–8.
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.
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.
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.
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.
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.
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.
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.
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.
Hofbauer RK, Rainville P, Duncan GH, Bushnell MC. Cortical representation of the sensory dimension of pain. J Neurophysiol. 2001;86(1):402–11.
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.
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.
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.
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.
Kaplan FS, Al Mukaddam M, Pignolo RJ. Acute unilateral hip pain in fibrodysplasia ossificans progressiva (FOP). Bone. 2018;109:115–9.
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.
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.
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.
Latremoliere A, Woolf CJ. Central sensitization: a generator of pain hypersensitivity by central neural plasticity. J Pain. 2009;10(9):895–926.
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.
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.
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.
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.
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.
Morales-Piga A, Kaplan FS. Osteochondral diseases and fibrodysplasia ossificans progressiva. Adv Exp Med Biol. 2010;686:335–48.
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.
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.
Ortiz-Agapito F, Colmenares-Bonilla D. Quality of life of patients with Fibrodysplasia ossificans progressiva. J Child Orthop. 2015;9(6):489–93.
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.
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.
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.
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.
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.
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.
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.
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.
Shore EM, Kaplan FS. Inherited human diseases of heterotopic bone formation. Nat Rev Rheumatol. 2010;6(9):518–27.
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.
Twycross R, Harcourt J, Bergl S. A survey of pain in patients with advanced cancer. J Pain Symptom Manag. 1996;12(5):273–82.
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.
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.
van Hecke O, Torrance N, Smith BH. Chronic pain epidemiology and its clinical relevance. Br J Anaesth. 2013;111(1):13–8.
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.
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.
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.
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.
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.
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 (https://www.ifopa.org) 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.
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
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
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). https://doi.org/10.1186/s13023-021-01709-4