Exercise effect on pain is associated with negative and positive affective components: A large-scale internet-based cross-sectional study in Japan.

Pain is a global health problem that leads to sedentary behavior and tends to cause negative emotion. In contrast, exercise is widely recommended for a health promotion, while pain often worsens with physical activity. Although exercise therapy is often prescribed to people with pain, the mechanisms of exercise effect on pain remains unclear. In this study, we tried to identify a universal association factor between regular exercise and pain intensity utilizing a cross-sectional web-based survey involving 52,353 adult participants from a large national study conducted in Japan. Using principal component analysis, we uncovered a mediation model of exercise effect on pain through psychological components. Analyses were performed in half of the population with pain (n = 20,330) and validated in the other half (n = 20,330), and showed that high-frequency exercise had a significant association with reduction in pain intensity. We also found Negative Affect and Vigor, two psychological components, are fully associating the exercise effect on pain (indirect effect = - 0.032, p < 0.001; association proportion = 0.99) with a dose-dependent response corresponding to the frequency of exercise. These findings were successfully validated (indirect effect of high-frequency exercise = - 0.028, p < 0.001; association proportion = 0.85). Moreover, these findings were also identified in subpopulation analyses of people with low back, neck, knee pain, and the tendency of the exercise effect on pain was increased with older people. In conclusion, the effect of exercise on pain is associated with psychological components and these association effects increased in parallel with the frequency of exercise habit regardless pain location.

Association between smoking and central sensitization pain: a web-based cross-sectional study.

PURPOSE: This study aimed to investigate whether smoking is an independent risk factor for central sensitization syndrome (CSS) in individuals with pain as measured by the Central Sensitization Inventory (CSI). METHODS: In 2020, we conducted an Internet survey targeting 2000 ordinary residents of Japan (aged 20-69 years) who had pain symptoms from October to November 2020. A multiple regression analysis was performed on the association between smoking status (nonsmokers and current smokers; Brinkman index) and CSI values. Moreover, compared to nonsmokers, the relative risk (RR) of the CSI cut-off score of 40 points or higher among current smokers was calculated using a modified Poisson regression model. Covariates included age, sex, body mass index, marital status, equivalized income, exercise habits, history of hypertension, history of hyperlipidemia, history of diabetes, pain chronicity, and Pain Catastrophizing Scale score. RESULTS: This study analyzed 1,822 individuals (1,041 men and 781 women). Among those experiencing pain, current smoking was associated with the increase in CSI values (ß = 0.07). The Brinkman index was also significantly associated with the increase in CSI values (ß = 0.06). Current smoking also increased the risk of being over the CSI cut-off score, with a relative risk (RR) of 1.29 (95% confidence intervals, 1.04-1.60). Younger age, being women, experiencing chronic pain, and higher pain catastrophizing thinking were also significantly associated with increased CSS severity, independent of smoking status. CONCLUSION: Smoking is an independent risk factor for CSS. This indicates that smoking may be an important factor in the management of central pain disorders.

Temporal Relations Between Pain Catastrophizing and Adverse Health and Mental Health Outcomes After Whiplash Injury.

OBJECTIVES: Pain catastrophizing has been shown to be a prognostic indicator for pain severity and the co-occurrence of mental health conditions such as depression and post-traumatic stress disorder after whiplash injury. However, the pattern of available findings is limited in its implications for the possible "antecedent" or "causal" role of pain catastrophizing. The purpose of the present study was to examine the temporal relations between pain catastrophizing, pain severity, depressive symptoms, and post-traumatic stress symptoms (PTSS) in individuals receiving treatment for whiplash injury. MATERIALS AND METHODS: The sample consisted of 388 individuals enrolled in a multidisciplinary program for whiplash injury. Participants completed self-report measures of pain catastrophizing, pain severity, depressive symptoms, and PTSS at the time of admission, mid-treatment (4 week), and treatment completion (7 week). A cross-lagged panel analysis was used to examine the temporal relations between pain catastrophizing, pain severity, depressive symptoms, and PTSS across all 3 timepoints. RESULTS: Model fit was acceptable after the inclusion of modification indices. Pain catastrophizing at the time of admission predicted all other variables at 4 weeks. Pain catastrophizing at 4 weeks also predicted all other variables at 7 weeks. In addition, some bidirectional relations were present, particularly for variables assessed at week 4 and week 7. DISCUSSION: Findings support the view that pain catastrophizing might play a transdiagnostic role in the onset and maintenance of health and mental health conditions. The findings call for greater emphasis on the development of treatment techniques that target pain catastrophizing in intervention programs for whiplash injury.

Surgical strategy for patients with papillary thyroid carcinoma invading the trachea: a comparison of tracheal sleeve resection with end-to-end anastomosis and window resection with tracheocutaneous fistula.

Background: Sleeve resection with end-to-end anastomosis (Procedure A) and window resection with a tracheocutaneous fistula (Procedure B) are the major surgical procedures for patients with papillary thyroid carcinoma (PTC) exhibiting transluminal tracheal invasion. For each procedure, the indications, postoperative course, and treatment results were examined retrospectively. Methods: Of 1,456 patients with PTC (maximum tumor diameter >1 cm) who received initial treatment between 1993 and 2013, we reviewed 51 patients. Of these 51 cases, 45 showed full-layer tracheal invasion, and 6 did not reach the tracheal mucosa, but required full-layer tracheal resection. Twenty-four patients underwent Procedure A, and 27 patients underwent Procedure B. Results: Regarding surgical procedure selection, Procedure B was selected significantly more frequently than Procedure A for cases with preoperative recurrent laryngeal nerve (RLN) palsy, tumor invasion of the esophagus, clinical lymph node metastasis, or a large number of resected tracheal rings. Postoperative airway-related complications were not significantly different between the procedures, but decreased with the use of intraoperative neuromonitoring (IONM). The postoperative hospital stay was significantly longer for Procedure B than for Procedure A. In addition, the rate of a permanent postoperative tracheostoma was higher with Procedure B than with Procedure A. Local recurrence-free survival (LRFS) and cause-specific survival (CSS) did not differ significantly between the two procedures. Conclusions: Certain patients may benefit from Procedure A with IONM in terms of a shorter hospital stay and avoiding the need for a permanent tracheostoma. Although Procedure B was indicated for patients with more advanced disease than Procedure A, treatment outcomes were similar.

Time-varying overdispersion of SARS-CoV-2 transmission during the periods when different variants of concern were circulating in Japan.

Japan has implemented a cluster-based approach for coronavirus disease 2019 (COVID-19) from the pandemic's beginning based on the transmission heterogeneity (overdispersion) of severe acute respiratory coronavirus 2 (SARS-CoV-2). However, studies analyzing overdispersion of transmission among new variants of concerns (VOCs), especially for Omicron, were limited. Thus, we aimed to clarify how the transmission heterogeneity has changed with the emergence of VOCs (Alpha, Delta, and Omicron) using detailed contact tracing data in Yamagata Prefecture, Japan. We estimated the time-varying dispersion parameter ([Formula: see text]) by fitting a negative binomial distribution for each transmission generation. Our results showed that even after the emergence of VOCs, there was transmission heterogeneity of SARS-CoV-2, with changes in [Formula: see text] during each wave. Continuous monitoring of transmission dynamics is vital for implementing appropriate measures. However, a feasible and sustainable epidemiological analysis system should be established to make this possible.

Clinicopathological study of surgically treated non-neoplastic diseases of the pancreas with special reference to autoimmune pancreatitis.

PURPOSE: After the popularization of serum immunoglobulin G4 (IgG4) measurement and endoscopic ultrasound-guided fine-needle aspiration (EUS-FNA) in our institute, surgical resection for non-neoplastic diseases of the pancreas became less common. Although the incidence of such false-positive cases was clarified in the 10-year period after the introduction of these measures (2009-2018), these data were not compared with the 30 years before 2009 (1979-2008). This study was performed to determine the percentage of autoimmune pancreatitis (AIP) that was included during the latter period and how the numbers of false-positive cases differed between the two periods. METHODS: From 1979 to 2008, 51 patients had clinical suspicion of pancreatic carcinoma (false-positive disease). Among these 51 patients, 32 non-alcoholic patients who had tumor-forming chronic pancreatitis (TFCP) were clinically, histologically, and immunohistochemically compared with 11 patients who had TFCP during the latter 10-year period. RESULTS: Retrospective IgG4 immunostaining of false-positive TFCP revealed 14 (35.0%) cases of AIP in the former 30 years versus 5 (45.5%) in the latter 10 years. There were 40 (5.9%) cases of TFCP among 675 patients in the former 30 years and 11 (0.9%) among 1289 patients in the latter 10 years. CONCLUSIONS: When the TFCP ratio of pancreatic resections and the AIP ratio of false-positive TFCPs were compared between the two periods, the TFCP ratio was 5.9% versus 0.9% and the AIP ratio was 35.0% versus 45.5%, respectively. It can thus be speculated that IgG4 measurement and EUS-FNA are absolutely imperative for the diagnosis of TFCP.

Longitudinal trajectories of smoking status using 25 year annually-updated data and all-cause mortality followed over 30 years: A community-based prospective cohort study.

We aimed to evaluate the long-term risk of smoking for all-cause mortality according to smoking status trajectories using 25-year annually-repeated input, traced by group-based trajectory modeling with an extension to account for non-random participant attrition or truncation due to death. We examined 2682 men and 4317 women aged 40 to 59 years who participated in annual health checks for the community-based prospective cohort study, 1975-1984 enrollment in Japan. The main outcome measure was all-cause mortality (follow-up period: median 30.2 years in men and 32.2 years in women). We traced annual smoking trajectories, stratified by sex and smoking status at baseline. For smokers at baseline, we identified five trajectories in both sexes, with different patterns of smoking cessation (e.g., early quitters and lifelong smokers). We calculated HRs and 95% CI of all-cause mortality using Cox proportional hazards regression modeling adjusted for age, body mass index, alcohol intake, blood pressure category, dyslipidemia and glucose category. Compared with one-time-point-based smokers, trajectory-based lifelong smokers had an increased risk of all-cause mortality; HRs were 1.31 (95% CI, 1.18-1.46) in men and 1.26 (95% CI, 0.91-1.73) in women. Among community residents aged 40 to 59 years, 25-year-trajectory-based lifelong smokers had an approximately 30% increased risk for all-cause mortality compared to one-time-point-based smokers. Risk of all-cause mortality among smokers with earlier cessation varied materially. It is necessary to consider the trajectories of smoking status to clarify the long-term excess risk of smoking.

Risk of herpes zoster according to past history in the general population: The Japanese Shozu herpes zoster study.

Immunity is known to persist after vaccination for varicella zoster virus, but the duration of immunity in patients who develop herpes zoster (HZ) remains unknown. To investigate the association between a past history of HZ and its occurrence in the general population. The Shozu HZ (SHEZ) cohort study included data for 12 299 individuals aged ≥50 years with information on their HZ history. Cross-sectional and 3-year follow-up studies were carried out to analyze the associations between a history of HZ (yes <10 years, yes ≥10 years, no) and the proportion of positive varicella zoster virus skin test results (erythema diameter ≥5 mm) and the risk of HZ after adjusting for potential confounding factors including age, sex, body mass index, smoking status, sleep duration, and mental stress. The incidences of positive skin test results were 87.7% (470/536) for individuals with a history of HZ <10 years ago, 82.2% (396/482) for those with a history of HZ ≥10 years, and 80.2% (3614/4509) for those with no history of HZ. The multivariable odds ratios (95% confidence intervals) of erythema diameter ≥5 mm were 2.07 (1.57-2.73) and 1. 39 (1.08-1.80) for individuals with a history <10 years and ≥10 years ago, respectively, compared with no history. The corresponding multivariable hazard ratios of HZ were 0.54 (0.34-0.85) and 1.16 (0.83-1.61), respectively. A past history of HZ <10 years ago may reduce the occurrence of HZ.

Predictors of herpes zoster severity and immune responses according to pain trajectories: A community-based prospective cohort study.

The authors aimed to identify determinants of the clinical course of herpes zoster and immunological responses, focusing on pain trajectories. This prospective community-based cohort study involved the analysis of responses to a valid pain survey provided by 375 patients diagnosed with herpes zoster based on clinical symptoms and virus identification by polymerase chain reaction. The authors analyzed most patients for humoral/cell-mediated immune response against varicella-zoster virus at the onset and 3 months post-onset. Six months post-initial visit, patients self-reported pain on a scale of 0 (no pain) to 5 (extremely strong pain) at up to 18 time points. Moreover, the pain trajectories were traced using group-based trajectory modeling. Subsequently, the authors used analysis of covariance to explore predictors and the humoral/cell-mediated immune response according to the pain trajectories. In addition, humoral/cell-mediated immune responses were assessed among each trajectory using paired t tests. Amon the five identified trajectories, two were isolated that particularly developed postherpetic neuralgia, with or without severe acute pain. Cancer therapy and corticosteroid use before herpes zoster onset specifically predicted postherpetic neuralgia without severe acute pain. In contrast, prescription of nonsteroidal anti-inflammatory drugs was uniquely associated with postherpetic neuralgia accompanied by severe acute pain. The aforementioned trajectories with postherpetic neuralgia showed increased antibodies and decreased cell-mediated immunity compared with those without postherpetic neuralgia. The authors could successfully distinguish between postherpetic neuralgia trajectories with and without severe acute pain. The identified key predictors and immunological responses against varicella-herpes zoster contribute further evidence to our understanding of the clinical features of herpes zoster and postherpetic neuralgia.

Maternal autistic traits and antenatal pain by cross-sectional analysis of the Japan Environment and Children's Study.

The aim of cross-sectional study was to investigate whether the presence of autistic traits in pregnant women was positively associated with the prevalence and severity of antenatal pain. We analyzed 89,068 pregnant women from a Japanese national birth cohort cross-sectionally. Autistic traits were assessed using the Japanese version of the Autism-Spectrum Quotient short form (AQ-10-J). Antenatal pain was measured using the SF-8 bodily pain item (SF-8-Pain). Antenatal pain in the second to third trimester during pregnancy was categorized into three groups: without pain, mild pain, and moderate-to-severe pain. Participants were divided into eight groups by AQ-10-J score: seven consecutive scoring groups (scores 0-6), and those above the cut-off (≥ 7) for probable autistic spectrum disorders. Odds ratios (OR) for the prevalence of mild and moderate-to-severe pain were calculated for each AQ-10-J scoring group (reference: without pain group) using multinominal logistic regression analysis. Autistic traits were positively associated with mild and moderate-to-severe pain in a dose-response manner, but the association with moderate-to-severe pain was strongest. Fully-adjusted ORs (95% confidence intervals) for moderate-to-severe pain were: 1.01 (0.91-1.13) for 1 point, 1.13 (1.02-1.25) for 2 points, 1.16 (1.04-1.29) for 3 points, 1.20 (1.07-1.34) for 4 points, 1.23 (1.09-1.40) for 5 points, 1.27 (1.10-1.47) for 6 points, and 1.24 (1.05-1.46) for ≥ 7 points (AQ-10-J cut-off). We identified an association between maternal autistic traits and antenatal pain. Maternal autistic traits may need to be considered when addressing antenatal pain during healthcare for expectant mothers.

Aprosamine Derivatives Active against Multidrug-Resistant Gram-Negative Bacteria.

Novel aprosamine derivatives were synthesized for the development of aminoglycoside antibiotics active against multidrug-resistant Gram-negative bacteria. The synthesis of aprosamine derivatives involved glycosylation at the C-8' position and subsequent modification (epimerization and deoxygenation at the C-5 position and 1-N-acylation) of the 2-deoxystreptamine moiety. All 8'-ß-glycosylated aprosamine derivatives (3a-h) showed excellent antibacterial activity against carbapenem-resistant Enterobacteriaceae and 16S ribosomal RNA methyltransferase-producing multidrug-resistant Gram-negative bacteria compared to the clinical drug, arbekacin. The antibacterial activity of 5-epi (6a-d) and 5-deoxy derivatives (8a,b and 8h) of ß-glycosylated aprosamine was further enhanced. On the other hand, the derivatives (10a,b and 10h) in which the amino group at the C-1 position was acylated with (S)-4-amino-2-hydroxybutyric acid showed excellent activity (MICs 0.25-0.5 µg/mL) against resistant bacteria that produce the aminoglycoside-modifying enzyme, aminoglycoside 3-N-acetyltransferase IV, which induces high resistance against parent apramycin (MIC > 64 µg/mL). In particular, 8b and 8h showed approximately 2- to 8-fold antibacterial activity against carbapenem-resistant Enterobacteriaceae and 8- to 16-fold antibacterial activity against resistant Gram-positive bacteria, such as methicillin-resistant Staphylococcus aureus and vancomycin-resistant enterococci, compared to apramycin. Our results showed that aprosamine derivatives have immense potential in the development of therapeutic agents for multidrug-resistant bacteria.

Measles Outbreak Response Activity in Japan, and a Discussion for a Possible Strategy of Outbreak Response Using Cycle Threshold Values of Real-Time Reverse Transcription PCR for Measles Virus in Measles Elimination Settings.

Measles is a highly contagious, but vaccine-preventable disease caused by the measles virus (MeV). Although the administration of two doses of measles vaccines is the most effective strategy to prevent and eliminate measles, MeV continues to spread worldwide, even in 2022. In measles-eliminated countries, preparedness and response to measles outbreaks originating from imported cases are required to maintain elimination status. Under these circumstances, real-time reverse transcription (RT) PCR for MeV could provide a diagnostic method capable of strengthening the subnational capacity for outbreak responses. Real-time RT-PCR can detect MeV RNA from patients with measles at the initial symptomatic stage, which can enable rapid public health responses aimed at detecting their contacts and common sources of infection. Furthermore, low cycle threshold (Ct) values (i.e., high viral load) of throat swabs indicate high infectiousness in patients with measles. The high basic reproduction number of measles suggests that patients with high infectiousness can easily become super-spreaders. This opinion proposes a possible strategy of rapid and intensive responses to counter measles outbreaks caused by super-spreader candidates showing low Ct values in throat swabs. Our strategy would make it possible to effectively prevent further measles transmission, thereby leading to the early termination of measles outbreaks.

Feasibility and applicability of locomotive syndrome risk test in elderly patients who underwent total knee arthroplasty.

OBJECTIVES: The concept of locomotive syndrome (LS) and its evaluation method, the LS risk test, have been applied in an integrated manner to capture the decline in mobility resulting from musculoskeletal disorders. The purpose of this study was to evaluate the impact of total knee arthroplasty (TKA) in the elderly with knee osteoarthritis, a common disorder found in LS. METHODS: A total of 111 patients were registered prior to TKA and postoperatively followed up for 1 year. Three components of the LS risk test (the two-step test, stand-up test, and Geriatric Locomotive Function Scale-25) were assessed pre- and postoperatively. RESULTS: After surgery, all three components of the test showed significant improvements from the baseline. The ratio of Stage 3 LS patients (progressed stage of decrease in mobility) reduced from 82.3% to 33.9% postoperatively. There was no significant difference in the degree of change in the scores between the younger (60-74 years) and older (≥75 years) age groups. CONCLUSIONS: We found that TKA has a major impact in preventing the progression of LS in patients with knee osteoarthritis. The LS risk test is a feasible tool for the longitudinal evaluation of patients with musculoskeletal diseases of varying severity and with multiple symptoms.

Frequency-response relationship between exercise and locomotive syndrome across age groups: Secondary analysis of a nationwide cross-sectional study in Japan.

OBJECTIVES: This study aimed to clarify the frequency-response relationship between exercise habits and locomotive syndrome (LS) in different age groups. METHODS: We analysed 8681 independent community dwellers who participated in a nationwide cross-sectional LS study. Participants were divided into young (20-39 years), middle (40-64 years), and old (≥ 65 years) groups. The presence of LS was evaluated using the LS risk tests. Exercise habits were assessed by exercise frequency (almost every day/sometimes/occasionally/never). The association between LS and exercise was assessed using the multiple logistic regression analysis adjusted for participant characteristics for each age category. RESULTS: Overall, 3528 (41%) participants had LS. In the old group, the risk of LS was lower with an increased exercise frequency. However, this dose-response relationship was not found in the young group. The proportion of LS was lower even in participants who occasionally exercised (1-3 days/month) than those who never exercise in the young and middle age groups. CONCLUSIONS: Regular exercise reduced the risk of LS not only in middle and older adults but also in young adults. Moreover, increasing exercise frequency may provide benefits for LS prevention, especially among older adults.

Practical guidance to handle missing values in the 25-question Geriatric Locomotive Function Scale (GLFS-25): a simulation study.

OBJECTIVES: Despite the possible large number of missing values on the 25-question Geriatric Locomotive Function Scale (GLFS-25), how we should treat them is unknown. In a simulation study, we investigated how to handle missing values in the GLFS-25. DESIGN, SETTING AND PARTICIPANTS: We used three datasets with different participant characteristics: community dwellers who could walk by themselves, outpatients of orthopaedics owing to pain, and patients who required surgery for total knee replacement or lumbar spinal canal stenosis. OUTCOME MEASURES: The missing items of the datasets were artificially created, and four statistical methods, complete case analysis, multiple imputation, single imputation using individual mean, and single imputation using individual domain average, were compared in terms of bias and mean squared error. Simulation studies were conducted to compare them under varying numbers of participants with missing values (5%-40%) and under varying numbers of missing items of GLFS-25 (4-16). RESULTS: Multiple imputation had the lowest root mean squared error. Complete case analysis showed the largest bias, and the performances of the single imputation were between those methods. The relative performances were similar across the three datasets. The absolute bias of the single imputation was<0.1. The bias and mean squared error of multiple imputation and single imputation were comparable when the number of missing items was less than or equal to eight. CONCLUSIONS: Multiple imputation is preferable, although single imputation using subject average/subject domain average can be used with practically negligible bias as long as the number of missing items is up to 8 out of 25 items in each individual of the population.

Self-assessment of pediatric basic life support among nursery teachers in a rural area of Japan.

•Serious injury and deaths can occur in childcare settings.•School teachers in Japan have low confidence about managing serious injuries.•With appropriate education we can increase confidence of childcare workers.•We evaluate the confidence of nursery teachers in pediatric basic life support.

Reliability and validity of the Japanese version of Pain Disability Index.

This study evaluated the reliability and validity of a Japanese version of Pain Disability Index (PDI). Analyses were conducted on a 7-item version (PDI-J) and a 5-item (PDI-5-J version of the PDI). Using a web-based survey system, we recruited 300 individuals with chronic low back pain (lasting ≥3 months) and 300 individuals with chronic daily headache (lasting ≥15 days per month for 3 months) aged 20-64 years. Analyses revealed a one-factor with goodness-of-fit indices assessed by confirmatory factor analysis. For concurrent validity, we calculated Pearson's correlation coefficients among the PDI-J, PDI-5-J, Pain Disability Assessment Scale, Pain numerical rating scale, and revised version of Short-Form McGill Pain Questionnaire. Internal consistency was evaluated by Cronbach's α, and test-retest reliability was assessed with intraclass correlations (ICCs) in 100 of 600 participants a week after the first response. Both Japanese adaptations of the PDI demonstrated good concurrent validity and reliability (Cronbach's α was 0.89 for PDI-J in chronic low back pain or chronic daily headache, and 0.94 and 0.93 for PDI-5-J in chronic low back pain and chronic daily headache, respectively). The PDI-J and PDI-5-J showed were highly correlated (r = 0.98). ICCs were 0.67 and 0.59 for the PDI-J and 0.59 and 0.63 for the PDI-5-J in chronic low back pain and chronic daily headache, respectively. In conclusion, these two PDI versions can be potentially used for evaluating pain-related interference with daily activities among the Japanese general population.

Prevalence and municipal variation in chronic musculoskeletal pain among independent older people: data from the Japan Gerontological Evaluation Study (JAGES).

BACKGROUND: Urbanization and population aging may affect prevalence of chronic pain from various causes. This cross-sectional study aimed to investigate the prevalence of chronic musculoskeletal pain, including some subtypes, in independent Japanese older people, and whether population density and population aging rate explained prevalence and differences in pain levels between municipalities. METHODS: We analyzed data from 12,883 independent older people living in 58 municipalities who completed mailed questionnaires and did not need support for daily living. We identified three types of pain: "chronic musculoskeletal pain" lasting ≥ 3 months (overall and in each part of the body), "chronic widespread-type pain" in the spinal and peripheral area, and "chronic multisite pain" in at least three sites. The latter two were measured using new definitions. These types of pain are correlated with depressive symptoms and we therefore examined the construct validity of the definitions by comparing the Geriatric Depression Scale score. We also used analysis of covariance to compare the prevalence of these three types of pain between municipalities. Odds ratios, median odds ratios, and the municipal variance in prevalence of chronic musculoskeletal pain were estimated by Bayesian multilevel logistic regression analysis using the Markov Chain Monte Carlo method. RESULTS: The construct validity of the definitions of chronic widespread-type pain and chronic multisite pain was confirmed. The prevalence of the three types of pain (chronic musculoskeletal, widespread, and multisite pain) was 39.0%, 13.9%, and 10.3%, respectively. Chronic musculoskeletal pain showed a higher prevalence among older people and women. Individuals in underpopulated, suburban, or metropolitan areas tended to have more pain than those in urban areas, but this was not statistically significant (odds ratio [95% credible interval] 1.15 [0.86-1.51], 1.17 [0.93-1.43], 1.17 [0.94-1.46]). Population density and population aging rate did not explain the differences between municipalities. CONCLUSIONS: The prevalence of chronic musculoskeletal pain was consistent with previous global reports. Areas with overpopulation and depopulation tended to have higher pain prevalence, but population density and population aging rate did not explain municipal variance. Further research is needed to identify other factors that contribute to regional variance.