ST8 and ST72 Methicillin-Resistant S. aureus Bacteremia in Korea: A Comparative Analysis of Clinical and Microbiological Characteristics
Yun Woo Lee, Ji-Hun Kim, So Yun Lim, Euijin Chang, Seongman Bae, Jiwon Jung, Min Jae Kim, Yong Pil Chong, Sung-Han Kim, Sang-Ho Choi, Sang-Oh Lee, Yang Soo Kim

TL;DR
This study compares ST8 and ST72 MRSA bloodstream infections in Korea, finding similar outcomes but differences in antibiotic resistance and genetic profiles.
Contribution
The study provides the first comparative analysis of clinical and microbiological features of ST8 and ST72 MRSA bacteremia in Korea.
Findings
ST8 isolates showed higher resistance to ciprofloxacin and erythromycin compared to ST72 isolates.
Vancomycin E-test MICs were higher for ST8 isolates, though broth microdilution results were similar.
ST8 and ST72 MRSA bacteremia had comparable mortality and clinical outcomes despite differing resistance patterns.
Abstract
Methicillin-resistant Staphylococcus aureus (MRSA) remains a major cause of bloodstream infection worldwide. In Korea, sequence type (ST) 72 has predominated, whereas ST8, including the USA300 lineage, has recently emerged. Comparative data on these genotypes in MRSA bacteremia (MRSAB) are limited. We conducted a retrospective cohort study of adult patients with MRSAB admitted to a 2700-bed tertiary care hospital in Republic of Korea between July 2008 and December 2020. Clinical features and outcomes of patients with ST8 MRSA were compared with those of patients with ST72 MRSA. Among 1975 cases of S. aureus bacteremia, 998 (50.5%) were due to MRSA, including 327 (32.7%) ST72 and 23 (2.3%) ST8 isolates. Demographics and comorbidities were similar, though pneumonia appeared more frequent in ST8 cases without statistical significance. ST8 isolates exhibited greater resistance to…
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsAntimicrobial Resistance in Staphylococcus · Bacterial Identification and Susceptibility Testing · Mycobacterium research and diagnosis
1. Introduction
Staphylococcus aureus remains a major human pathogen responsible for a broad spectrum of community- and healthcare-associated infections worldwide [1]. Since the introduction of antibiotics, the organism has acquired resistance to multiple agents, and methicillin-resistant S. aureus (MRSA) has emerged as a major healthcare problem [2,3]. MRSA initially emerged as a predominant cause of hospital- and healthcare-associated infections [1,4]. In the late 1990s, however, community-associated MRSA (CA-MRSA) clones became widespread, with sequence type (ST) 8-Staphylococcal Cassette Chromosome (SCC)mec IVa (USA300) rapidly establishing itself as the predominant lineage in the United States and subsequently spreading globally [5,6,7]. USA300 strains are typically Panton–Valentine leucocidin (PVL)-positive and have been strongly associated with skin and soft tissue infections, but they are also capable of causing severe invasive infections such as necrotizing pneumonia, osteomyelitis, septic arthritis, and bacteremia [7].
In contrast, the epidemiology of MRSA in Republic of Korea has been distinct. PVL-negative ST72-SCCmec IVc has been the predominant CA-MRSA clone in both community and healthcare settings [8,9,10]. Recently, however, ST8 lineages including USA300 have begun to emerge in Korea. Sporadic cases of USA300 infection were first reported in 2008 in a patient with a travel history to Hawaii, followed by additional cases without overseas exposure, implying in-country acquisition [11,12,13,14]. A multicenter Korean study identified PVL-positive ST8-MRSA (USA300) accounting for 0.4–12% of MRSA bacteremia across three centers, a range higher than expected [15]. Longitudinal surveillance has further demonstrated a steady increase in ST8 isolates, raising concerns about their potential expansion in both community and hospital settings [8]. Similar trends have been documented in neighboring countries such as Japan, where a distinct ST8 clone (ST8 CA-MRSA/J) with close genetic similarity to USA300 but unique virulence gene content has emerged and caused severe invasive infections [16]. In Taiwan as well, recent studies have reported a rising prevalence of USA300-like clones, underscoring the need for ongoing regional and global surveillance of these lineages [17].
Despite these observations, comparative studies of ST8 and ST72 in invasive infections remain scarce in Korea. In this study, we evaluated the clinical and microbiological characteristics and outcomes of MRSA bacteremia caused by ST8 and compared them with those caused by ST72, the long-dominant community genotype in Korea.
2. Materials and Methods
2.1. Study Design and Population
This retrospective cohort study was conducted using data from a previously established institutional cohort of adult patients with S. aureus bacteremia (SAB) admitted to a 2700-bed tertiary care hospital (Asan Medical Center, Seoul, Republic of Korea) between July 2008 and December 2020. All consecutive patients with a first episode of SAB were eligible. The exclusion criteria were as follows: (i) age under 18 years, (ii) discharge before culture results were available, which precluded outcome assessment, (iii) polymicrobial bacteremia, (iv) a previous episode of SAB within the preceding 90 days, (v) receipt of ≥3 days of treatment at another hospital, (vi) clinically insignificant bacteremia, defined as a single positive blood culture without signs or symptoms of infection that resolved without anti-staphylococcal therapy, and (vii) identification in outpatient settings. Clinical data, demographics, main focus of infection, and test results of metastatic infection, were collected within one week of the initial positive culture. All SAB cases were then reviewed by infectious disease specialists. Outcomes, including duration of bacteremia, recurrence, and death, were assessed through medical record review up to 90 days after the first positive culture. This study specifically analyzed and compared the clinical and microbiological characteristics and outcomes MRSA bacteremia caused by ST8 and ST72 isolates within this institutional cohort.
2.2. Laboratory and Microbiological Data
All S. aureus isolates were identified using the standard methods recommended by the Clinical and Laboratory Standards Institute (CLSI) and confirmed by morphological and biochemical testing, including the slide coagulase assay, or by matrix-assisted laser desorption ionization–time of flight (MALDI-TOF) mass spectrometry [18]. Antimicrobial susceptibility testing was interpreted according to standard criteria using the MicroScan system (Beckman Coulter, Brea, CA, USA) in accordance with the CLSI guidelines [18]. Methicillin resistance was determined based on oxacillin minimum inhibitory concentration (MIC) values and the detection of the mecA gene. Vancomycin MICs were assessed by both the broth microdilution (BMD) method and E test. BMD was performed following CLSI recommendations, whereas E test was performed according to the manufacturer’s protocol using vancomycin E test strips (bioMérieux, Marcy-l’Étoile, France). Multilocus sequence typing (MLST), assessment of δ-hemolysin activity for evaluation of accessory gene regulator (agr) function, and staphylococcal protein A (spa) genotyping were conducted as previously described [19,20,21]. These molecular analyses were performed on all MRSA isolates collected between July 2008 and December 2020.
2.3. Data Collection and Definitions
Electronic medical records were reviewed to collect data on age, sex, mode of acquisition, medical history, site of infection, severity at admission, microbiological results, and clinical outcomes for all enrolled patients. The mode of acquisition was categorized as community-acquired, community-onset healthcare-associated, or nosocomial, in accordance with the Friedman criteria [22]. The cause (main focus) of bacteremia was determined for each patient based on clinical, microbiological, and radiological findings, and classified as catheter-related, skin and soft tissue, bone and joint, pneumonia, endocarditis, urinary tract infection, arteriovenous graft infection, or primary bacteremia. Metastatic infection was defined as a newly developed sterile-site infection that was not clinically apparent at the time of the initial blood culture and was not identified during the initial diagnostic workup. The primary focus of infection was determined by infectious disease specialists based on microbiological evidence, radiological findings, and clinical course. The duration of bacteremia was defined as the interval between the first and last positive SAB culture, regardless of antimicrobial therapy. Recurrence was defined as the reappearance of signs and symptoms of infection more than 7 days after clinical improvement with confirmed negative blood cultures.
2.4. Statistical Analysis
Data manipulation and statistical analyses were conducted using R, version 4.0.4 (R Foundation for Statistical Computing, Vienna, Austria). Continuous variables were compared using the Wilcoxon rank-sum test, and categorical variables were analyzed using Fisher’s exact test or Pearson’s chi-squared test with simulated p-values, as appropriate. Univariate and multivariate logistic regression models were applied to identify independent risk factors for mortality. Multivariate models were adjusted for the following covariates: age ≥ 65, gender, hypertension, diabetes mellitus, liver cirrhosis, end-stage renal disease, chronic lung disease, solid cancer, hematologic malignancy, solid organ transplantation, ischemic heart disease, Charlson Comorbidity Index (CCI) > 4, focus of infection, presence of metastatic infection, and status of infection focus removal. To minimize small-sample bias due to the limited number of ST8 cases and zero 30-day deaths, Firth penalized and exact logistic regression analyses were additionally performed, including the prespecified covariates (ST type, age > 60, CCI > 4, and infection severity). Post hoc power for 30-day mortality difference was calculated (α = 0.05) with effect size expressed as Cohen’s h. Temporal trends in the annual proportions of ST72 and ST8 MRSA among SAB cases were analyzed using ANOVA, Kruskal–Wallis tests, and logistic regression with calendar year as a continuous variable. Yearly trends in the proportion of isolates with vancomycin MIC > 1 mg/L (by BMD and E-test) were also evaluated using binomial logistic regression, reporting odds ratios (ORs) per year with 95% confidence intervals. Because only 8 of 17 isolates in 2020 had available E-test results, MIC trend analyses were restricted to data up to 2019. To assess the potential clinical impact of vancomycin MIC, outcomes were further compared between isolates with MIC ≤1 mg/L and >1 mg/L. All reported p-values were two-sided, and a p-value < 0.05 was considered to indicate statistical significance.
3. Results
3.1. Baseline and Clinical Characteristics
A total of 1975 cases of SAB were identified during the study period, including 998 cases of MRSA (50.5%). ST72 and ST8 accounted for 327 (32.7%) and 23 (2.3%) of the MRSA isolates, respectively. The baseline characteristics of 350 patients with ST72 MRSA and ST8 MRSA are summarized in Table 1. Patients in both groups were of similar age (median, 64 years) and sex distribution, and the mode of acquisition did not differ significantly, with the majority of infections being healthcare-associated or hospital-acquired. The prevalence of comorbidities such as hypertension, diabetes mellitus, and solid cancer was generally comparable between the two groups. CCI scores and predisposing conditions, including recent surgery, prior antibiotic treatment, and use of immunosuppressive agents, were also similar across groups. With respect to the focus of infection, pneumonia was more common in the ST8 group (30% vs. 9.2%), although the difference did not reach statistical significance. Other infection sites, including catheter-related infections, bone and joint infection, skin and soft tissue infection, and endocarditis, were distributed similarly. The occurrence of metastatic infection (18% vs. 22%) and infection severity (sepsis or septic shock) did not differ significantly between ST72 and ST8 MRSA bacteremia.
3.2. Microbiological Characteristics
The microbiological characteristics and antibiotic susceptibilities of ST72 and ST8 MRSA isolates are summarized in Table 2. Overall, resistance to most antibiotics was comparable between the two groups; however, ST8 isolates demonstrated significantly higher resistance to ciprofloxacin (78% vs. 8.9%, p < 0.001) and erythromycin (61% vs. 26%, p = 0.003) compared with ST72 isolates. Resistance to clindamycin (13–20%) and gentamicin (13–30%) was observed at moderate frequencies, while resistance to fusidic acid, rifampin, and tetracycline remained low (<5%), with no significant differences between ST72 and ST8 isolates. All MRSA isolates were susceptible to vancomycin according to CLSI criteria; no vancomycin-resistant strains were identified. Vancomycin MICs determined by BMD were similar between the two groups, with the majority of isolates showing an MIC of 1 mg/L. However, E-test results revealed a significantly different distribution (p < 0.001), with ST8 isolates more frequently exhibiting MIC values ≥1 mg/L. Four ST72 isolates (1.2%) showed E-test MICs of 3 mg/L, which may reflect inter-method variability, as all values remained well below the CLSI resistance breakpoint. The prevalence of agr dysfunction did not differ significantly between ST72 and ST8 (12% vs. 17%, p = 0.501). spa typing revealed distinct clonal distributions: ST72 isolates were dominated by t324 (48%), t148 (12%), and diverse other types, whereas ST8 isolates were largely associated with spa type t008 (65%).
Over the 12-year study period, the annual proportion of ST8 MRSA among SAB cases gradually increased from 1.3% in 2008 to 3.3% in 2020, whereas ST72 remained relatively stable (Figure S1, Table S3). When analyzed by 4-year intervals (2008–2012, 2013–2016, 2017–2020), the ST8 proportion showed a significant upward trend (ANOVA p = 0.02; Kruskal–Wallis p = 0.04), while no significant change was observed for ST72. Logistic regression demonstrated a significant year-on-year increase for ST8 (OR 1.30 per year, p < 0.001) but not for ST72 (OR 1.03, p = 0.06) (Figure S2, Table S4). Logistic regression demonstrated a decreasing proportion of isolates with MIC > 1 mg/L for both methods (BMD: OR = 0.82 per year, 95% CI 0.73–0.91, p < 0.001; E-test: OR = 0.86 per year, 95% CI 0.79–0.93, p < 0.001; Figure S3).
Clinical outcomes did not differ significantly according to vancomycin MIC category (≤1 vs. >1 mg/L) (Table S5).
3.3. Clinical Course and Outcome
Management and clinical outcomes of patients with MRSA bacteremia are summarized in Table 3. Most patients in both groups received vancomycin with similar rates of therapeutic drug monitoring and comparable trough levels. The use of alternative agents, time to appropriate therapy, and duration of intravenous antibiotics did not differ significantly. The frequencies of eradicable foci, focus removal, and persistent bacteremia were also similar. Clinical outcomes, including 30- and 90-day mortality, showed no significant differences between ST72 and ST8 cases. The 30-day mortality rate was 13% in the ST72 group and 0% in the ST8 group. However, due to complete separation arising from the absence of deaths in the ST8 group, ORs could not be estimated. Fisher’s exact test revealed no statistically significant difference between the two groups (p = 0.093).
Univariable and multivariable logistic regression analyses were performed to identify risk factors associated with 30-day mortality (Table 4). In the univariable analysis, liver cirrhosis (OR, 3.26; 95% CI, 1.52–6.73; p = 0.002), solid cancer (OR, 3.03; 95% CI, 1.58–5.93; p < 0.001), and a CCI > 4 (OR, 4.70; 95% CI, 2.44–9.21; p < 0.001) were significantly associated with increased mortality. The presence of a removable focus was protective, with a significantly lower risk of 30-day mortality when the focus was removed (OR, 0.43; 95% CI, 0.20–0.86; p = 0.021).
In the multivariable analysis, liver cirrhosis (aOR, 2.82; 95% CI, 1.20–6.44; p = 0.015) and CCI > 4 (aOR, 4.67; 95% CI, 2.00–11.5; p < 0.001) remained independent predictors of 30-day mortality. Removal of the infectious focus showed a trend toward reduced mortality but did not reach statistical significance (aOR, 0.52; 95% CI, 0.23–1.11; p = 0.10). Other factors, including age, sex, MRSA genotype (ST72 vs. ST8), agr dysfunction, and vancomycin MIC, were not significantly associated with 30-day mortality.
In Firth and exact logistic regression analyses, the OR for 30-day mortality in ST8 versus ST72 was 0.17 (95% CI, 0.00–1.34; p = 0.11), consistent across both models (Table S1). In addition, a post hoc power analysis demonstrated that the study had sufficient statistical power (0.958; Cohen’s h = 0.795, α = 0.05) to detect the observed absolute difference in 30-day mortality between the groups, suggesting that the non-significant results were less likely due to limited sample size (Table S2).
In both univariate and multivariable analyses, ST8 MRSA was not significantly associated with 90-day recurrence or mortality compared with ST72 MRSA (Table 5).
4. Discussion
In this study, we compared the clinical and microbiological characteristics and outcomes of MRSA bacteremia caused by ST72 and ST8 isolates in Korea.
Baseline characteristics, including age, sex, comorbidities, and infection sites, were largely comparable between the two groups, which may partly explain the absence of significant differences in outcomes. Notably, pneumonia was more frequently observed among ST8 cases, echoing prior reports linking USA300-related clones with pulmonary infections [13,23]; however, in our cohort this trend did not reach statistical significance.
Microbiological analyses revealed distinct differences between the two clonal groups. ST8 isolates exhibited significantly higher resistance rates to ciprofloxacin and erythromycin and showed a predominance of spa type t008, consistent with the USA300 lineage described in other regions [5]. In contrast, ST72 isolates demonstrated more diverse spa types and lower resistance profiles. Although vancomycin MICs determined by BMD were similar between the groups, E-test results showed higher MIC distributions in ST8 isolates (mostly 1 mg/L). No vancomycin-resistant MRSA isolates (MIC ≥ 16 mg/L) were identified, and all values remained within the susceptible range by CLSI criteria [18]. This subtle upward shift underscores the need for continued surveillance of potential MIC creep [24,25]. However, this difference was not associated with adverse clinical outcomes in our study.
Despite these microbiological differences, management strategies and short-term outcomes did not differ significantly. Most patients in both groups received vancomycin therapy with therapeutic drug monitoring, and the use of alternative agents, duration of intravenous therapy, and frequency of source control procedures were comparable. Likewise, persistent bacteremia and 30- and 90-day mortality rates showed no significant differences, suggesting that clonal background exerted limited influence on treatment response or outcomes in this cohort. Nevertheless, ST8/USA300 typically harbors SCCmec IVa, arginine catabolic mobile element (ACME), and PVL genes, features that may confer enhanced virulence compared to other MRSA strains [26,27,28], whereas ST72-SCCmecIVc isolates are generally PVL-negative and their pathogenic potential remains less well defined [10,29]. Further genome-wide studies and functional analyses will be necessary to clarify the clinical implications of these clonal differences.
In this cohort, most MRSA bacteremia cases were healthcare-associated (40%) or hospital-acquired (49%), indicating that infections predominantly occurred in patients with prior healthcare exposure, such as those with vascular catheters, recent surgery, or indwelling devices. This pattern reflects the tertiary care hospital setting and aligns with known risk factors for MRSA bacteremia. Epidemiologically, ST72-SCCmecIVc has been the dominant community-associated MRSA clone in Korea since the mid-2000s but has increasingly spread into healthcare environments [10,29,30]. In contrast, ST8 MRSA, including USA300 and related variants, has emerged more recently. Although USA300 originally arose in the community in North America, recent Korean studies have identified its establishment within hospital settings, even among patients without travel history, suggesting local adaptation and circulation of hospital-associated ST8 lineages [15]. Continued genomic and epidemiological surveillance is warranted to differentiate imported USA300-like strains from locally evolved ST8 variants and to monitor their clinical and public health impact.
When risk factors for 30-day mortality were examined, host-related factors rather than bacterial genotype emerged as key determinants. Higher Charlson Comorbidity Index (CCI > 4) and underlying liver cirrhosis were independently associated with mortality, consistent with previous studies demonstrating the impact of comorbidity burden and cirrhosis on adverse outcomes in SAB [31,32]. Although the protective effect of source control did not remain statistically significant in the multivariable analysis, its clinical relevance is well supported in the literature [33,34]. Together, these findings suggest that in Korea, the prognosis of MRSA bacteremia is primarily driven by host comorbidities and adequacy of management, rather than clonal lineage.
Finally, analysis of 90-day outcomes demonstrated that ST8 was not significantly associated with increased recurrence or mortality compared with ST72. This further supports the notion that host factors, rather than clonal background, are the predominant determinants of prognosis in this setting. However, given the relatively small number of ST8 cases, these findings should be interpreted with caution, and larger multicenter studies are needed to validate them. Considering the established virulence and adverse clinical impact of USA300/ST8 in other regions [26,27,28], continued molecular epidemiologic surveillance remains essential to monitor potential shifts in clinical relevance.
This study expands upon our institution’s long-term SAB cohort, providing novel comparative insights into the clinical and microbiological distinctions between ST72 and ST8 MRSA clones—an aspect not specifically addressed in earlier analyses from the same cohort [8,10,29,32]. Through this work, we contribute further evidence to the understanding of MRSA clonal diversity and its clinical implications in Korea.
Our study had some limitations. First, it was conducted in a single center with a relatively small number of ST8 cases, which limits the generalizability of our findings and reduces the statistical power to detect subtle differences between clonal groups. In particular, the absence of 30-day mortality events in the ST8 group constrained statistical estimation, highlighting the challenges of interpreting rare outcomes in small cohorts. However, supplementary analyses using Firth and exact logistic regression, along with post hoc power estimation, demonstrated that our findings were robust and not substantially affected by small-sample bias. Second, the retrospective design introduces potential information bias and confounding, particularly regarding the assessment of clinical management such as adequacy of source control. In addition, the exclusion of patients diagnosed outside the hospital or transferred after more than three days may have introduced selection bias; however, these criteria were applied to ensure consistency in microbiological testing, antimicrobial stewardship, and clinical management within a single institutional framework. Third, detailed genotypic information such as PVL, ACME, and SCCmec subtypes was not available due to the retrospective nature of the study and limited sample availability. As a result, molecular confirmation of ST8 isolates as USA300 or locally evolved ST8 variants was not possible. The absence of these molecular data, together with the lack of virulence gene profiling or whole-genome sequencing, precluded a more comprehensive evaluation of clonal virulence characteristics and their potential impact on clinical outcomes. Further genomic characterization will be important to clarify the genetic relatedness of emerging ST8 strains in Korea and their potential epidemiologic significance. Therefore, our findings should be interpreted with caution, and larger, prospective, multicenter studies incorporating advanced molecular analyses are warranted to validate and extend these results.
5. Conclusions
In conclusion, clinical outcomes of MRSA bacteremia did not differ significantly between ST72 and ST8 isolates in Korea, despite distinct microbiological profiles. Mortality was mainly driven by host comorbidities such as liver cirrhosis and high CCI rather than clonal background. Although the small number of ST8 cases limits definitive interpretation, continued molecular surveillance is warranted to monitor the potential clinical impact of emerging ST8 clones.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Lowy F.D. Staphylococcus aureus infections N. Engl. J. Med.199833952053210.1056/NEJM 1998082033908069709046 · doi ↗ · pubmed ↗
- 2Chambers H.F. Deleo F.R. Waves of resistance: Staphylococcus aureus in the antibiotic era Nat. Rev. Microbiol.2009762964110.1038/nrmicro 220019680247 PMC 2871281 · doi ↗ · pubmed ↗
- 3Lakhundi S. Zhang K. Methicillin-Resistant Staphylococcus aureus: Molecular Characterization, Evolution, and Epidemiology Clin. Microbiol. Rev.201831 e 00020-1810.1128/CMR.00020-1830209034 PMC 6148192 · doi ↗ · pubmed ↗
- 4Grundmann H. Aires-de-Sousa M. Boyce J. Tiemersma E. Emergence and resurgence of meticillin-resistant Staphylococcus aureus as a public-health threat Lancet 200636887488510.1016/S 0140-6736(06)68853-316950365 · doi ↗ · pubmed ↗
- 5Diekema D.J. Richter S.S. Heilmann K.P. Dohrn C.L. Riahi F. Tendolkar S. Mc Danel J.S. Doern G.V. Continued emergence of USA 300 methicillin-resistant Staphylococcus aureus in the United States: Results from a nationwide surveillance study Infect. Control Hosp. Epidemiol.20143528529210.1086/67528324521595 · doi ↗ · pubmed ↗
- 6Diep B.A. Otto M. The role of virulence determinants in community-associated MRSA pathogenesis Trends Microbiol.20081636136910.1016/j.tim.2008.05.00218585915 PMC 2778837 · doi ↗ · pubmed ↗
- 7Nimmo G.R. USA 300 abroad: Global spread of a virulent strain of community-associated methicillin-resistant Staphylococcus aureus Clin. Microbiol. Infect. Off. Publ. Eur. Soc. Clin. Microbiol. Infect. Dis.20121872573410.1111/j.1469-0691.2012.03822.x 22448902 · doi ↗ · pubmed ↗
- 8Choi S.H. Lee J. Jung J. Kim E.S. Kim M.J. Chong Y.P. Kim S.H. Lee S.O. Choi S.H. Woo J.H. A Longitudinal Study of Adult Patients with Staphylococcus aureus Bacteremia over 11 Years in Korea J. Korean Med. Sci.202136 e 10410.3346/jkms.2021.36.e 10433904260 PMC 8076844 · doi ↗ · pubmed ↗
