Association among polydisperse aerodynamic size of bioaerosols, biodiversity and urbanization in kindergartens.

The aerodynamic sizes of bioaerosols may significantly affect their behaviors, respiratory deposition and biodiversity. The respirable bacterial size, biodiversity, and human-associated bacteria (HAB) related bioaerosols were evaluated at three kindergartens in Taiwan. Kindergartens A, B, and C were in urban, semi-urban, and rural areas, respectively. A six-stage viable Andersen cascade impactor was used to collect bioaerosols and to determine their size distributions. The geometric mean diameter (GMD), geometric standard deviation (GSD), heat maps, and uniformity were used to evaluate the association of bacteria characteristics. A BD Phoenix-100 automated interpretation system was used to identify the airborne bacteria species. The results revealed that 1425 colonies of the sampled airborne bacteria contained 63 species in 29 genera, and overall, 63.0% were HABs. The most abundant phylum was Actinobacteria (56.6 ± 22.2%) and Firmicutes (31.6 ± 22.3%), and from the taxonomic analysis, both airborne Micrococcus and the Staphylococcus aureus are the dominant genus. All the bacteria aerodynamic particle size distributions were polydisperse distributions. The heat map and uniformity analysis had revealed most of the sampled bioaerosols distributed between 1.1-3.3 µm, and most of the polydisperse airborne Streptococcus spp. had a size in the respirable range, due to urbanization, they have potentially contributed to respiratory risk in the kindergartens. The Shannon diversity index (H) and inverse Simpson diversity index (D) of the bioaerosols in urban kindergarten were negatively correlated with GMD and GSD. The Pearson correlations revealed that the kindergarten in the rural area, with a higher temperature, a lower relative humidity, and a lower CO2 concentration than the others, tended to have the largest H and D values (P < 0.05). Multiple and stepwise regression revealed that bioaerosol aerodynamic size was statistically significantly correlated with H (P = 0.001) and D values (P = 0.002). This study sheds light on the characteristics of bioaerosols and their associations with microbiome.

CSMed® wound dressing for prophylaxis and management of radiation dermatitis in breast and head-neck cancer patients: a single hospital prospective clinical trial.

PURPOSE: CSMed® wound dressing, a dressing with various herb extracts, was tested for its therapeutic effect in radiation dermatitis of breast and head-and-neck cancer patients. METHODS: This study included 20 breast cancer patients and 10 head-and-neck cancer patients. Half of the irradiated area was covered with CSMed® and the other half was under routine treatment. The severity of radiation dermatitis was evaluated with radiation therapy oncology group (RTOG) grade throughout the treatment and the follow-up period. The RTOG grade between the dressed and undressed area were compared to illustrate the therapeutic effect of CSMed® dressing. RESULTS: The results showed that CSMed® dressed area had significant lower RTOG score at 3-7 weeks and final record during the treatment, and 1-3 weeks during follow-up than undressed area. CONCLUSIONS: This indicated that CSMed® can delay the onset, reduce the severity, and enhance healing of radiation dermatitis. CSMed® can be used for prophylaxis and management of radiation dermatitis.

Spine, hip, and femoral neck bone mineral density in relation to vegetarian type and status among Taiwanese adults.

We determined the association of vegetarian type and status with bone mineral density (BMD) Z-scores at the spine, hip, and femoral neck. Compared to non-vegetarians, current vegetarians, especially vegans, lacto-vegetarians, and lacto-ovo-vegetarians had lower Z-scores at multiple sites. Sole reliance on a vegetarian diet might be detrimental to the bone. PURPOSE: The impact of vegetarian diets on BMD is contentious. We determined the association of vegetarian type and status with the spine, hip, and femoral neck BMD Z-scores. METHODS: We analyzed data from 20,110 Taiwan Biobank volunteers. BMD was measured using dual-energy X-ray absorptiometry (DXA). The vegetarian status (non-, former, and current vegetarians) and type (non-vegetarians, ovo-vegetarians, lacto-vegetarians, lacto-ovo-vegetarians, and vegans) were determined using questionnaires. RESULTS: The participants consisted of 12,910 women and 7200 men, with a mean age of 55.5 years. Based on vegetarian status (reference: non-vegetarians), current vegetarians had significantly lower BMD Z-scores at the spine (unstandardized regression coefficient, B = - 0.195, p = 0.006), left hip (B = - 0.125, p = 0.008), and right hip (B = - 0.100, p = 0.027), respectively. Based on vegetarian status and type (reference: non-vegetarians), current vegans and non-vegans had notably lower BMD Z-scores at specific skeletal sites. For non-vegans, the BMD Z-scores were significant at the spine (B = -0.184, p = 0.010), left hip (B = - 0.124, p = 0.010), and left femoral neck (B = - 0.125, p = 0.012). For current vegans, however, the BMD Z-scores were significant only at the right hip (B = - 0.232; p = 0.028). Nonetheless, after stratifying vegetarian diet into more subgroups, current vegans exhibited a significant reduction in BMD Z-scores at the spine and right hip, with B-coefficients of - 0.326 and - 0.238, respectively. Current lacto-vegetarians also had significantly lower Z-scores (p < 0.05) at the spine (B = - 0.459), left hip (B = - 0.313), and right hip (B = - 0.214). Moreover, current lacto-ovo-vegetarians had significantly lower Z-scores at the spine (B = - 0.175) and left hip (B = - 0.115). CONCLUSION: Current vegetarians, particularly vegans, lacto-vegetarians, and lacto-ovo-vegetarians, demonstrated significantly lower BMD Z-scores at various skeletal sites compared to non-vegetarians. Sole reliance on a vegetarian diet might be detrimental to the bone.

Analyzing Financial Efficiency of Public-Private Partnerships Hospitals in Taiwan.

This study evaluates the management capacity ability and profitable capacity of eight public-private partnership hospitals in Taiwan from 2015 to 2020. By conducting various ratio analyses of the financial statement, this study found these hospitals have achieved a balance between management efficiency and profitability, thereby confirming the viability of the PPP model for hospital management. In addition, the subject hospitals play a vital role as isolation hospitals during the COVID-19 pandemic. Beyond offering medical assistance to infected individuals, these hospitals contribute to the integrity of Taiwan's medical network, mitigating the impact of the pandemic. Overall, establishing and managing hospitals with PPP partnership is a feasible solution as it alleviates governmental financial burdens related to medical welfare and achieves profitability.

Synergistic Effect of Ginsenoside Rh2 Combines with Ionizing Radiation on CT26/luc Colon Carcinoma Cells and Tumor-Bearing Animal Model.

BACKGROUND: The local tumor control rate of colon cancer by radiotherapy is unsatisfactory due to recurrence and radioresistance. Ginsenoside Rh2 (Rh2), a panoxadiol saponin, possesses various antitumor effects. METHODS: CT26/luc murine colon carcinoma cells and a CT26/luc tumor-bearing animal model were used to investigate the therapeutic efficacy of Rh2 combined with ionizing radiation and the underlying mechanisms. RESULTS: Rh2 caused cell cycle arrest at the G1 phase in CT26/luc cells; however, when combined with ionizing radiation, the cells were arrested at the G2/M phase. Rh2 was found to suppress the activity of NF-κB induced by radiation by inhibiting the MAPK pathway, consequently affecting the expression of effector proteins. In an in vivo study, the combination treatment significantly increased tumor growth delay time and overall survival. Furthermore, the combination treatment significantly reduced NF-κB and NF-κB-related effector proteins, along with PD-1 receptor expression. Additionally, Rh2 administration led to increased levels of interleukin-12, -18, and interferon-γ in the mice's sera. Importantly, biochemical analysis revealed no toxicities associated with Rh2 alone or combined with radiation. CONCLUSIONS: The combination of Rh2 with radiation may have potential as an alternative to improve the therapeutic efficacy of colorectal cancer.

Effectiveness of multimedia courses in improving self-care among patients with breast cancer undergoing radiotherapy.

BACKGROUND: Ninety percent of patients receiving radiation therapy experience side effects. Busy schedules and intensive health education programs may lead to incomplete education content delivery and inaccurate patient self-care implementation. This study investigated whether multimedia health education improves the accuracy of patient self-care implementation compared with paper-based education. METHODS: From March 11, 2020 to February 28, 2021, 110 patients were randomly divided into experimental and control groups, each comprising 55 participants. Paper-based materials were used along with multimedia materials. Radiology self-care awareness questionnaires were administered to both groups before the first treatment and on day 10. The differences in radiology self-care awareness between the two groups was analyzed with inferential statistics, independent t tests, categorical data, and Pearson's chi-squared test. Differences between the two groups were considered significant at a p value of < 0.05. RESULTS: The treatment accuracy rate improved from 10.9 to 79.1% in the control group and from 24.8 to 98.5% in the experimental group, indicating an improvement in both groups. The difference was significant. These results indicate that the intervention could improve the effectiveness of self-care. CONCLUSIONS: Participants who used pretreatment multimedia health education exhibited a higher rate of having a correct understanding of treatment self-care than did the control group. These findings can inform the development of a patient-centered cancer treatment knowledge base for improved quality of care.

Association between ESR1 rs2234693 single nucleotide polymorphism and uterine fibroids in Taiwanese premenopausal and postmenopausal women.

BACKGROUND: Uterine fibroids (UFs) are uterine smooth muscle neoplasms that affect women, especially during the reproductive stage. Both genetic and lifestyle factors affect the onset of the disease. We examined the association between the estrogen receptor 1 (ESR1) rs2234693 variant (whose genotypes are TT, TC, and CC) and UFs in Taiwanese premenopausal and postmenopausal women. METHODS: We linked individual-level data of 3588 participants from the Taiwan Biobank to the National Health Insurance Research Database at the Health and Welfare Data Science Center. The association of the ESR1 rs2234693 variant and other variables with UFs was determined by multiple logistic regression, and the results were presented as odds ratios and 95% confidence intervals (CIs). RESULTS: The 3588 participants comprised 622 cases and 2966 controls. In all the participants, the ESR1 rs2234693 TC and CC genotypes compared to the reference genotype (TT) were associated with a lower risk of UFs. However, the results were significant only for the CC genotype (OR; 95% CI = 0.70; 0.52-0.93). Noteworthy, the association of TC and CC with UFs was dose-dependent (p-trend = 0.012). Based on menopausal status, both TC and CC were significantly and dose-dependently associated with a lower risk of UFs in premenopausal women (OR; 95% CI = 0.76; 0.59-0.98 for TC and 0.64; 0.43-0.95 for CC: p-trend = 0.010). CONCLUSION: The TC and CC genotypes of the ESR1 rs2234693 variant may reduce susceptibility to UFs, especially in premenopausal women.

Interval changes of histogram-derived diffusion indices predict treatment response to induction chemotherapy in head and neck cancer: a feasibility study.

Background: This retrospective study investigated whether the interval change of apparent diffusion coefficient (∆ADC) [baseline and after the first cycle of induction chemotherapy (ICT)] can be used as a valid predictive imaging biomarker of the treatment response to ICT in head and neck cancer (HNC). Methods: A total of 19 consecutive patients with HNC who underwent diffusion-weighted magnetic resonance imaging (DWI) at baseline and after the first cycle of ICT were included. Whole-tumor ADC histogram parameters (mean, median, kurtosis, skewness, entropy, minimal, maximum, 25th percentile, and 75th percentile) were obtained. The correlations of ∆ADC histogram parameters, volume, T-stage, N-stage, and age with the treatment response were examined using the Mann-Whitney U test. The predictive value of histogram parameters was examined using receiver operating characteristic (ROC) curves. Results: Responders showed significantly higher values of ∆ADC25 (0.19±0.23) and ∆ADCmin (1.78±2.98) than non-responders (-0.09±0.15 and -0.73±0.36; P=0.035 and 0.009, respectively). When ∆ADC25 and ∆ADCmin were used for predicting the treatment response, the area under the ROC curve was 0.850/0.933 with a sensitivity of 73.3%/80.0% and specificity of 100%/100% (P=0.036 and 0.009, respectively). Conclusions: ∆ADC25 and ∆ADCmin derived from whole-tumor histogram analysis are valuable imaging biomarkers for the early prediction of the ICT response in HNC.

High vs. low radiotherapy dose in locally advanced esophageal squamous cell carcinoma patients treated with neoadjuvant concurrent chemoradiotherapy: an endemic area population-based study.

BACKGROUND: The optimal radiotherapy dose for locally advanced esophageal squamous cell carcinoma in endemic areas treated with neoadjuvant concurrent chemoradiotherapy is unclear. METHODS: Eligible patients diagnosed between 2010 and 2019 were identified via the Taiwan Cancer Registry. We used propensity score (PS) weighting to balance observable potential confounders. The hazard ratio (HR) of death was compared between high dose (50-50.4 Gy) and low dose (40-41.4 Gy) radiotherapy. We also evaluated other outcomes and performed supplementary analyses via an alternative approach. RESULTS: Our study population consisted of 644 patients. The PS weight-adjusted HR of death was 0.92 (95% confidence interval: 0.7-1.19, p = 0.51). There were no statistically significant differences for other outcomes or supplementary analyses. CONCLUSIONS: In this population-based study from an endemic area, we found no significant difference in overall survival between high vs. low radiotherapy doses.

Designing Mobile Epidemic Prevention Medical Stations for the COVID-19 Pandemic and International Medical Aid.

The demand for mobile epidemic prevention medical stations originated from the rapid spread of the COVID-19 pandemic. In order to reduce the infection risk of medical practitioners and provide flexible medical facilities in response to the variable needs of the pandemic, this research aimed to design mobile medical stations for COVID-19 epidemic prevention, the emergence of which began in February 2020. The mobile medical stations include a negative pressure isolation ward, a positive pressure swabbing station, a fever clinic and a laboratory. In Taiwan, many medical institutions used the mobile swabbing station design of this study to practice COVID-19 screening pre-tests. Internationally, this study assisted Palau in setting up medical stations to provide anti-epidemic goods and materials. The design of this study not only provides a highly flexible and safe medical environment but the benefits of screening can also be used as resources for medical research, forming an economic circulation for operation sustainability. In addition, the design of this study can also be used during the non-epidemic period as a healthcare station for rural areas or as a long-term community medical station.

Association of Peripheral Arterial Occlusive Disease and Deep Venous Thrombosis with Risk of Consequent Sepsis Event: A Retrospective Population-Based Cohort Study.

Peripheral artery occlusive disease (PAOD) and deep vein thrombosis (DVT) can cause a variety of acute and chronic vascular complications and put patients at risk of subsequent sepsis. This study aimed to determine whether DVT compared with PAOD patients would increase the risk of sepsis. This study recruited 43,535 patients newly diagnosed as having PAOD and 6932 patients who were newly diagnosed as having DVT from a population of 2 million patients from the Longitudinal Health Insurance Database. Propensity score matching (PSM) between the PAOD and DVT groups was performed for age, sex, comorbidities, and prior antibiotic administration. A total of 4383 patients with PAOD and 4383 patients with DVT were analyzed for risk of sepsis. The incidence density of sepsis per 1000 person years for patients with PAOD was 25.75 (95% CI = 23.90 to 27.74) and 35.61 (95% CI = 33.29 to 38.09) for patients with DVT. After age, sex, associated comorbidities, and antibiotic administration were adjusted for, the risk of sepsis for the DVT group was 1.46-fold (95% CI = 1.32-1.62) higher than that for the PAOD group. In conclusion, patients with DVT were associated with a higher risk of subsequent sepsis than patients with PAOD. Aging was another risk factor.

Effects of aerobic and resistance exercise on glycosylated hemoglobin (HbA1c) concentrations in non-diabetic Taiwanese individuals based on the waist-hip ratio.

BACKGROUND: Glycosylated hemoglobin (HbA1c) reflects the average blood sugar over the past eight to twelve weeks. Several demographic and lifestyle factors are known to affect HbA1c levels. We evaluated the association of HbA1c with aerobic and resistance exercise in non-diabetic Taiwanese adults based on the waist-hip ratio (WHR). METHODS: We conducted this study based on TWB data collected from 90,958 individuals between 2008 and 2019. We estimated the Beta (ß) coefficient and 95% confidence intervals (CI) for HbA1c using multivariate regression models. RESULTS: Based on the multivariate analysis, lower HbA1c levels were associated with both resistance exercise (ß-coefficient = -0.027, 95% CI -0.037 to -0.017) and aerobic exercise (ß-coefficient = 0.018, 95% CI, -0.023 to -0.013). Higher HbA1c levels were associated with abnormal WHR compared to normal WHR (ß-coefficient = 0.091, 95% CI, 0.086 to 0.096). We detected an interaction between exercise and WHR (p for interaction = 0.0181). To determine the magnitude of the interaction, we performed additional analyses (with the reference group being 'abnormal WHR with no exercise') and observed substantial decreases in HbA1c regardless of the WHR and exercise category. However, the largest reduction occurred in the 'normal WHR and resistance exercise' group (ß = -0.121, 95% CI, -0.132 to -0.109). CONCLUSIONS: We found that normal resistance exercise, coupled with a normal WHR was significantly associated with lower HbA1c levels among non-diabetic individuals in Taiwan.

Impact of Genetic Variants of Long Noncoding RNA Metastasis-Associated Lung Adenocarcinoma Transcript 1 on Uterine Cervical Cancer.

Genetic variants of long noncoding RNA metastasis-associated lung adenocarcinoma transcript 1 (lncRNA MALAT1) have been reported to be associated with several cancers. Until now, no study reveals the associations between lncRNA MALAT1 polymorphisms and cervical cancer (CC). The objectives of this study were to explore the correlations among MALAT1 polymorphisms and occurrence and clinicopathological parameters of CC, as well as patient 5 years survival in Taiwanese women. The study recruited 116 patients with cervical invasive cancer and 89 patients with cervical precancerous lesions, as well as 268 non-cancer control women. LncRNA MALAT1 polymorphisms rs3200401, rs619586 and rs1194338 were selected and their genotypic frequencies were defined by real-time polymerase chain reaction. Our results revealed that there are no relationships between lncRNA MALAT1 genetic variants and occurrence of CC. The independent factor among lncRNA MALAT1 genetic variants and clinicopathological parameters were positive pelvic lymph node metastasis (p=0.001, HR: 10.94, 95% CI: 2.65-45.23). In conclusions, lncRNA MALAT1 genetic variants are not related to occurrence and clinicopathological characteristics of CC and patient 5 years survival in Taiwanese women. Pelvic lymph node metastasis could independently predict the patient 5 years survival among various MALAT1 polymorphisms and clinicopathological factors in CC.

Risk of uterine leiomyoma based on BET1L rs2280543 single nucleotide polymorphism and vegetarian diet.

BACKGROUND: Bet1 Golgi vesicular membrane trafficking protein-like (BET1L) rs2280543 single nucleotide polymorphism (SNP) and diet have been independently associated with uterine leiomyoma (UL). However, whether the SNP and diet could jointly influence the risk of UL is yet to be assessed. Therefore, we investigated the independent and interactive effects of vegetarian diet and BET1L rs2280543 on uterine fibroids in Taiwanese women. METHODS: We linked participants' electronic data in the Taiwan Biobank (TWB) database to their medical records in the National Health Insurance Research Database (NHIRD). The TWB had genotypic, lifestyle, and biochemical data between 2008 and 2015 and the NHIRD had data on disease diagnoses between 1998 and 2015. In this study, we included 1997 premenopausal women with complete data. RESULTS: Compared  to participants with the BET1L rs2280543 CC genotype (wildtype), those with CT/CC genotype had an odds ratio (OR) of 0.69 and a 95% confidence interval (CI) of 0.51-0.93. Vegetarian diet and UL were not significantly associated: OR = 1.09 and 95% CI = 0.77-1.55. However, the test for interaction between rs2280543 and vegetarian diet was significant (p = 0.046). Compared to individuals with the CC genotype, the risk of UL was lower among vegetarians with the CT/TT genotype: OR (95% CI) = 0.15 (0.05-0.47). CONCLUSION: The BET1L rs2280543 CT/TT genotype was associated with a lower risk of UL especially among vegetarians.

Evaluating environmental radiations of the tomotherapy facility by optimizing full factorial design of the TLD technique.

BACKGROUND: In the last 40 years, the number of deaths due to cancer has been the highest in TaiwanOBJECTIVE: To optimize the readout system of the thermoluminescent dosimeter (TLD)-100H, the radiation rates among the Tomotherapy (TOMO) facility of the Department of Radiology Oncology of Chung Shan Medical University Hospital (CSMUH) were calculated with a 32 full factorial design (FFD). METHODS: A ten-month survey of the facility was employed using the sensitive and accurate TLD method. The TLD system was optimized for maximum temperature, heat rate, and preheat temperature of Harshaw 3500 reader. Eight analyzed groups with different factors were tested. RESULTS: The TOMO facility had significantly different radiation rates. The farther away from the gantry head, environmental radiation rates. The half value layer (HVL) was also determined. These results were compared with published. No significant contributions of environmental gamma radiations were detected except in the treatment room. CONCLUSIONS: Those were far below the occupational doses recommended by ICRP 60.

Interaction between a diabetes-related methylation site (TXNIP cg19693031) and variant (GLUT1 rs841853) on fasting blood glucose levels among non-diabetics.

BACKGROUND: Type 2 diabetes mellitus (T2DM) is caused by a combination of environmental, genetic, and epigenetic factors including, fasting blood glucose (FBG), genetic variant rs841853, and cg19693031 methylation. We evaluated the interaction between rs841853 and cg19693031 on the FBG levels of non-diabetic Taiwanese adults. METHODS: We used Taiwan Biobank (TWB) data collected between 2008 and 2016. The TWB data source contains information on basic demographics, personal lifestyles, medical history, methylation, and genotype. The study participants included 1300 people with DNA methylation data. The association of cg19693031 methylation (stratified into quartiles) with rs841853 and FBG was determined using multiple linear regression analysis. The beta-coefficients (ß) and p-values were estimated. RESULTS: The mean ± standard deviation (SD) of FBG in rs841853-CC individuals (92.07 ± 7.78) did not differ significantly from that in the CA + AA individuals (91.62 ± 7.14). However, the cg19693031 methylation levels were significantly different in the two groups (0.7716 ± 0.05 in CC individuals and 0.7631 ± 0.05 in CA + AA individuals (p = 0.002). The cg19693031 methylation levels according to quartiles were ß < 0.738592 (< Q1), 0.738592 ≤ 0.769992 (Q1-Q2), 0.769992 ≤ 0.800918 (Q2-Q3), and ß ≥ 0.800918 (≥ Q3). FBG increased with decreasing cg19693031 methylation levels in a dose-response manner (ptrend = 0.005). The ß-coefficient was - 0.0236 (p = 0.965) for Q2-Q3, 1.0317 (p = 0.058) for Q1-Q2, and 1.3336 (p = 0.019 for < Q1 compared to the reference quartile (≥ Q3). The genetic variant rs841853 was not significantly associated with FBG. However, its interaction with cg19693031 methylation was significant (p-value = 0.036). Based on stratification by rs841853 genotypes, only the CC group retained the inverse and dose-response association between FBG and cg19693031 methylation. The ß (p-value) was 0.8082 (0.255) for Q2-Q3, 1.6930 (0.022) for Q1-Q2, and 2.2190 (0.004) for < Q1 compared to the reference quartile (≥ Q3). The ptrend was 0.002. CONCLUSION: Summarily, methylation at cg19693031 was inversely associated with fasting blood glucose in a dose-dependent manner. The inverse association was more prominent in rs841853-CC individuals, suggesting that rs841853 could modulate the association between cg19693031 methylation and FBG. Our results suggest that genetic variants may be involved in epigenetic mechanisms associated with FBG, a hallmark of diabetes. Therefore, integrating genetic and epigenetic data may provide more insight into the early-onset of diabetes.

Investigating surgical smoke in otolaryngology operating rooms.

Surgical smoke is a common chemical hazard produced from the use of electrocautery, laser, or ultrasonic scalpels during surgery. It has been proved harmful to medical personnel. Thus, it is important to monitor surgical smoke concentrations in the operating room. In the past decade, many researches regarding surgical smoke were discussed in different professional healthcare fields, but few showed the correlation between surgical smoke and otolaryngology surgery. In this study, the concentrations of particulate matter and formaldehyde were measured during thirty cases of several types of otolaryngology surgery in a regional research hospital in Taiwan. The concentrations of 0.3 µm and 0.5 µm particulate matter raised rapidly in the main knife range at the beginning of the electrocautery knife used, and then decreased by half after 5-10 min of use. The concentrations of formaldehyde were ranged from 1 to 2 ppm during the surgery, which is higher than the permissible exposure limit. While many medical staffs are working in the operating room and are exposed to the smoke hazard, effective strategies for collecting and eliminating the smoke should be taken in all medical facilities.

The Effects of Frequent Coffee Drinking on Female-Dominated Healthcare Workers Experiencing Musculoskeletal Pain and a Lack of Sleep.

Previous research has demonstrated that chronic diseases can occur due to musculoskeletal (MS) pain and poor sleep. It is also worth noting that the caffeine in coffee can reduce overall sleep duration, efficiency, and quality. Thus, the present study examines the effects of frequent coffee drinking (two cups per day) on individuals experiencing MS pain and a lack of sleep during the COVID-19 period. This observational and cross-sectional study recruited 1615 individuals who completed the self-reported (Nordic musculoskeletal) questionnaire. Long-term, frequent coffee drinking and a sleep duration of less than 6 h per day were significantly associated with neck and shoulder pain among healthy individuals. The mediation model demonstrated that the shorter sleep duration and drinking multiple cups of coffee per day had a two-way relationship that worsened such pain over the long term. Specifically, individuals who experienced such pain frequently drank multiple cups of coffee per day, which, in turn, shortened their sleep durations. In summary, long-term coffee drinking creates a vicious cycle between MS pain and sleep duration. Therefore, the amount of coffee should be fewer than two cups per day for individuals who sleep less than 6 h per day or suffer from MS pain, especially neck and shoulder pain.