Abemaciclib plus a nonsteroidal aromatase inhibitor as initial therapy for HR+, HER2- advanced breast cancer: final overall survival results of MONARCH 3.

BACKGROUND: In MONARCH 2, the addition of abemaciclib to fulvestrant significantly improved both progression-free survival (PFS) and overall survival (OS) in patients with hormone receptor-positive (HR+), human epidermal growth factor receptor 2-negative (HER2-) advanced breast cancer (ABC) with disease progression on prior endocrine therapy. In MONARCH 3, the addition of abemaciclib to a nonsteroidal aromatase inhibitor (NSAI) as initial therapy for HR+, HER2- ABC significantly improved PFS. Here, we present the prespecified final OS results for MONARCH 3. PATIENTS AND METHODS: MONARCH 3 is a randomized, double-blind, phase III study of abemaciclib plus NSAI (anastrozole or letrozole) versus placebo plus NSAI in postmenopausal women with HR+, HER2- ABC without prior systemic therapy in the advanced setting. The primary objective was investigator-assessed PFS; OS was a gated secondary endpoint, and chemotherapy-free survival was an exploratory endpoint. RESULTS: A total of 493 women were randomized 2 : 1 to receive abemaciclib plus NSAI (n = 328) or placebo plus NSAI (n = 165). After a median follow-up of 8.1 years, there were 198 OS events (60.4%) in the abemaciclib arm and 116 (70.3%) in the placebo arm (hazard ratio, 0.804; 95% confidence interval 0.637-1.015; P = 0.0664, non-significant). Median OS was 66.8 versus 53.7 months for abemaciclib versus placebo. In the subgroup with visceral disease, there were 113 OS events (65.3%) in the abemaciclib arm and 65 (72.2%) in the placebo arm (hazard ratio, 0.758; 95% confidence interval 0.558-1.030; P = 0.0757, non-significant). Median OS was 63.7 months versus 48.8 months for abemaciclib versus placebo. The previously demonstrated PFS benefit was sustained, and chemotherapy-free survival numerically improved with the addition of abemaciclib. No new safety signals were observed. CONCLUSIONS: Abemaciclib combined with an NSAI resulted in clinically meaningful improvement in median OS (intent-to-treat population: 13.1 months; subgroup with visceral disease: 14.9 months) in patients with HR+ HER2- ABC; however, statistical significance was not reached.

Observation of the Anomalous Shape of X(1840) in J/ψ→γ3(π

Using a sample of (10087±44)×10^{6} J/ψ events, which is about 45 times larger than that was previously analyzed, a further investigation on the J/ψ→γ3(π^{+}π^{-}) decay is performed. A significant distortion at 1.84 GeV/c^{2} in the line shape of the 3(π^{+}π^{-}) invariant mass spectrum is observed for the first time, which could be resolved by two overlapping resonant structures, X(1840) and X(1880). The new state X(1880) is observed with a statistical significance larger than 10σ. The mass and width of X(1880) are determined to be 1882.1±1.7±0.7 MeV/c^{2} and 30.7±5.5±2.4 MeV, respectively, which indicates the existence of a pp[over ¯] bound state.

Investigation of the ΔI=1/2 Rule and Test of CP Symmetry through the Measurement of Decay Asymmetry Parameters in Ξ

Using (10087±44)×10^{6} J/ψ events collected with the BESIII detector, numerous Ξ^{-} and Λ decay asymmetry parameters are simultaneously determined from the process J/ψ→Ξ^{-}Ξ[over ¯]^{+}→Λ(pπ^{-})π^{-}Λ[over ¯](n[over ¯]π^{0})π^{+} and its charge-conjugate channel. The precisions of α_{Λ0} for Λ→nπ^{0} and α[over ¯]_{Λ0} for Λ[over ¯]→n[over ¯]π^{0} compared to world averages are improved by factors of 4 and 1.7, respectively. The ratio of decay asymmetry parameters of Λ→nπ^{0} to that of Λ→pπ^{-}, ⟨α_{Λ0}⟩/⟨α_{Λ-}⟩, is determined to be 0.873±0.012_{-0.010}^{+0.011}, where the first and the second uncertainties are statistical and systematic, respectively. The ratio is smaller than unity more than 5σ, which signifies the existence of the ΔI=3/2 transition in Λ for the first time. Besides, we test for CP symmetry in Ξ^{-}→Λπ^{-} and in Λ→nπ^{0} with the best precision to date.

Observation of Significant Flavor-SU(3) Breaking in the Kaon Wave Function at 12

We present cross sections for the reaction e^{+}e^{-}→K_{S}^{0}K_{L}^{0} at center-of-mass energies ranging from 3.51 to 4.95 GeV using data samples collected in the BESIII experiment, corresponding to a total integrated luminosity of 26.5 fb^{-1}. The ratio of neutral-to-charged kaon form factors at large momentum transfers (12

Measurements of All-Particle Energy Spectrum and Mean Logarithmic Mass of Cosmic Rays from 0.3 to 30 PeV with LHAASO-KM2A.

We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at 3.67±0.05±0.15 PeV. Below the knee, the spectral index is found to be -2.7413±0.0004±0.0050, while above the knee, it is -3.128±0.005±0.027, with the sharpness of the transition measured with a statistical error of 2%. The mean logarithmic mass of cosmic rays is almost heavier than helium in the whole measured energy range. It decreases from 1.7 at 0.3 PeV to 1.3 at 3 PeV, representing a 24% decline following a power law with an index of -0.1200±0.0003±0.0341. This is equivalent to an increase in abundance of light components. Above the knee, the mean logarithmic mass exhibits a power law trend towards heavier components, which is reversal to the behavior observed in the all-particle energy spectrum. Additionally, the knee position and the change in power-law index are approximately the same. These findings suggest that the knee observed in the all-particle spectrum corresponds to the knee of the light component, rather than the medium-heavy components.

Observation of D

We perform for the first time an amplitude analysis of the decay D^{+}→K_{S}^{0}π^{+}η and report the observation of the decay D^{+}→K_{S}^{0}a_{0}(980)^{+} using 2.93 fb^{-1} of e^{+}e^{-} collision data taken at a center-of-mass energy of 3.773 GeV with the BESIII detector. As the only W-annihilation-free decay among D to a_{0}(980) pseudoscalar, D^{+}→K_{S}^{0}a_{0}(980)^{+} is the ideal decay in extracting the contributions of the W-emission amplitudes involving a_{0}(980) and to study the final-state interactions. The absolute branching fraction of D^{+}→K_{S}^{0}π^{+}η is measured to be (1.27±0.04_{stat}±0.03_{syst})%. The branching fractions of intermediate processes D^{+}→K_{S}^{0}a_{0}(980)^{+} with a_{0}(980)^{+}→π^{+}η and D^{+}→π^{+}K[over ¯]_{0}^{*}(1430)^{0} with K[over ¯]_{0}^{*}(1430)^{0}→K_{S}^{0}η are measured to be (1.33±0.05_{stat}±0.04_{syst})% and (0.14±0.03_{stat}±0.01_{syst})%, respectively.

Observation of Structures in the Processes e

We present measurements of the Born cross sections for the processes e^{+}e^{-}→ωχ_{c1} and ωχ_{c2} at center-of-mass energies sqrt[s] from 4.308 to 4.951 GeV. The measurements are performed with data samples corresponding to an integrated luminosity of 11.0 fb^{-1} collected with the BESIII detector operating at the Beijing Electron Positron Collider storage ring. Assuming the e^{+}e^{-}→ωχ_{c2} signals come from a single resonance, the mass and width are determined to be M=(4413.6±9.0±0.8) MeV/c^{2} and Γ=(110.5±15.0±2.9) MeV, respectively, which is consistent with the parameters of the well-established resonance ψ(4415). In addition, we also use one single resonance to describe the e^{+}e^{-}→ωχ_{c1} line shape and determine the mass and width to be M=(4544.2±18.7±1.7) MeV/c^{2} and Γ=(116.1±33.5±1.7) MeV, respectively. The structure of this line shape, observed for the first time, requires further understanding.

First Study of Antihyperon-Nucleon Scattering Λ[over ¯]p→Λ[over ¯]p and Measurement of Λp→Λp Cross Section.

Using (10.087±0.044)×10^{9} J/ψ events collected with the BESIII detector at the BEPCII storage ring, the processes Λp→Λp and Λ[over ¯]p→Λ[over ¯]p are studied, where the Λ/Λ[over ¯] baryons are produced in the process J/ψ→ΛΛ[over ¯] and the protons are the hydrogen nuclei in the cooling oil of the beam pipe. Clear signals are observed for the two reactions. The cross sections in -0.9≤cosθ_{Λ/Λ[over ¯]}≤0.9 are measured to be σ(Λp→Λp)=(12.2±1.6_{stat}±1.1_{syst}) and σ(Λ[over ¯]p→Λ[over ¯]p)=(17.5±2.1_{stat}±1.6_{syst}) mb at the Λ/Λ[over ¯] momentum of 1.074 GeV/c within a range of ±0.017 GeV/c, where the θ_{Λ/Λ[over ¯]} are the scattering angles of the Λ/Λ[over ¯] in the Λp/Λ[over ¯]p rest frames. Furthermore, the differential cross sections of the two reactions are also measured, where there is a slight tendency of forward scattering for Λp→Λp, and a strong forward peak for Λ[over ¯]p→Λ[over ¯]p. We present an approach to extract the total elastic cross sections by extrapolation. The study of Λ[over ¯]p→Λ[over ¯]p represents the first study of antihyperon-nucleon scattering, and these new measurements will serve as important inputs for the theoretical understanding of the (anti)hyperon-nucleon interaction.

Determination of Spin-Parity Quantum Numbers of X(2370) as 0

Based on (10087±44)×10^{6} J/ψ events collected with the BESIII detector, a partial wave analysis of the decay J/ψ→γK_{S}^{0}K_{S}^{0}η^{'} is performed. The mass and width of the X(2370) are measured to be 2395±11(stat)_{-94}^{+26}(syst) MeV/c^{2} and 188_{-17}^{+18}(stat)_{-33}^{+124}(syst) MeV, respectively. The corresponding product branching fraction is B[J/ψ→γX(2370)]×B[X(2370)→f_{0}(980)η^{'}]×B[f_{0}(980)→K_{S}^{0}K_{S}^{0}]=(1.31±0.22(stat)_{-0.84}^{+2.85}(syst))×10^{-5}. The statistical significance of the X(2370) is greater than 11.7σ and the spin parity is determined to be 0^{-+} for the first time. The measured mass and spin parity of the X(2370) are consistent with the predictions of the lightest pseudoscalar glueball.

Measurements of Normalized Differential Cross Sections of Inclusive η Production in e

Using data samples collected with the BESIII detector operating at the BEPCII storage ring, the cross section of the inclusive process e^{+}e^{-}→η+X, normalized by the total cross section of e^{+}e^{-}→hadrons, is measured at eight center-of-mass energy points from 2.0000 to 3.6710 GeV. These are the first measurements with momentum dependence in this energy region. Our measurement shows a significant discrepancy compared to the existing fragmentation functions. To address this discrepancy, a new QCD analysis is performed at the next-to-next-to-leading order with hadron mass corrections and higher twist effects, which can explain both the established high-energy data and our measurements reasonably well.

Observation of D_{s}

By analyzing 7.33 fb^{-1} of e^{+}e^{-} annihilation data collected at center-of-mass energies between 4.128 and 4.226 GeV with the BESIII detector, we report the observation of the semileptonic decay D_{s}^{+}→η^{'}µ^{+}ν_{µ}, with a statistical significance larger than 10σ, and the measurements of the D_{s}^{+}→ηµ^{+}ν_{µ} and D_{s}^{+}→η^{'}µ^{+}ν_{µ} decay dynamics for the first time. The branching fractions of D_{s}^{+}→ηµ^{+}ν_{µ} and D_{s}^{+}→η^{'}µ^{+}ν_{µ} are determined to be (2.235±0.051_{stat}±0.052_{syst})% and (0.801±0.055_{stat}±0.028_{syst})%, respectively, with precision improved by factors of 6.0 and 6.6 compared to the previous best measurements. Combined with the results for the decays D_{s}^{+}→ηe^{+}ν_{e} and D_{s}^{+}→η^{'}e^{+}ν_{e}, the ratios of the decay widths are examined both inclusively and in several ℓ^{+}ν_{ℓ} four-momentum transfer ranges. No evidence for lepton flavor universality violation is found within the current statistics. The products of the hadronic form factors f_{+,0}^{η^{(')}}(0) and the c→s Cabibbo-Kobayashi-Maskawa matrix element |V_{cs}| are determined. The results based on the two-parameter series expansion are f_{+,0}^{η}(0)|V_{cs}|=0.452±0.010_{stat}±0.007_{syst} and f_{+,0}^{η^{'}}(0)|V_{cs}|=0.504±0.037_{stat}±0.012_{syst}, which help to constrain present models on f_{+,0}^{η^{(')}}(0). The forward-backward asymmetries are determined to be ⟨A_{FB}^{η}⟩=-0.059±0.031_{stat}±0.005_{syst} and ⟨A_{FB}^{η^{'}}⟩=-0.064±0.079_{stat}±0.006_{syst} for the first time, which are consistent with the theoretical calculation.

Determination of the Σ

Based on data samples collected with the BESIII detector at the BEPCII collider, the process e^{+}e^{-}→Σ^{+}Σ[over ¯]^{-} is studied at center-of-mass energies sqrt[s]=2.3960, 2.6454, and 2.9000 GeV. Using a fully differential angular description of the final state particles, both the relative magnitude and phase information of the Σ^{+} electromagnetic form factors in the timelike region are extracted. The relative phase between the electric and magnetic form factors is determined to be sinΔΦ=-0.67±0.29(stat)±0.18(syst) at sqrt[s]=2.3960 GeV, ΔΦ=55°±19°(stat)±14°(syst) at sqrt[s]=2.6454 GeV, and 78°±22°(stat)±9°(syst) at sqrt[s]=2.9000 GeV. For the first time, the phase of the hyperon electromagnetic form factors is explored in a wide range of four-momentum transfer. The evolution of the phase along with four-momentum transfer is an important input for understanding its asymptotic behavior and the dynamics of baryons.

Coupled-Channel Analysis of the χ_{c1}(3872) Line Shape with BESIII Data.

We perform a study of the χ_{c1}(3872) line shape using the data samples of e^{+}e^{-}→γχ_{c1}(3872), χ_{c1}(3872)→D^{0}D[over ¯]^{0}π^{0}, and π^{+}π^{-}J/ψ collected with the BESIII detector. The effects of the coupled channels and the off-shell D^{*0} are included in the parametrization of the line shape. The line shape mass parameter is obtained to be M_{X}=(3871.63±0.13_{-0.05}^{+0.06}) MeV. Two poles are found on the first and second Riemann sheets corresponding to the D^{*0}D[over ¯]^{0} branch cut. The pole location on the first sheet is much closer to the D^{*0}D[over ¯]^{0} threshold than the other, and is determined to be 7.04±0.15_{-0.08}^{+0.07} MeV above the D^{0}D[over ¯]^{0}π^{0} threshold with an imaginary part -0.19±0.08_{-0.19}^{+0.14} MeV.

Study of the f_{0}(980) and f_{0}(500) Scalar Mesons through the Decay D_{s}

Using e^{+}e^{-} collision data corresponding to an integrated luminosity of 7.33 fb^{-1} recorded by the BESIII detector at center-of-mass energies between 4.128 and 4.226 GeV, we present an analysis of the decay D_{s}^{+}→π^{+}π^{-}e^{+}ν_{e}, where the D_{s}^{+} is produced via the process e^{+}e^{-}→D_{s}^{*±}D_{s}^{∓}. We observe the f_{0}(980) in the π^{+}π^{-} system and the branching fraction of the decay D_{s}^{+}→f_{0}(980)e^{+}ν_{e} with f_{0}(980)→π^{+}π^{-} measured to be (1.72±0.13_{stat}±0.10_{syst})×10^{-3}, where the uncertainties are statistical and systematic, respectively. The dynamics of the D_{s}^{+}→f_{0}(980)e^{+}ν_{e} decay are studied with the simple pole parametrization of the hadronic form factor and the Flatté formula describing the f_{0}(980) in the differential decay rate, and the product of the form factor f_{+}^{f_{0}}(0) and the c→s Cabibbo-Kobayashi-Maskawa matrix element |V_{cs}| is determined for the first time to be f_{+}^{f_{0}}(0)|V_{cs}|=0.504±0.017_{stat}±0.035_{syst}. Furthermore, the decay D_{s}^{+}→f_{0}(500)e^{+}ν_{e} is searched for the first time but no signal is found. The upper limit on the branching fraction of D_{s}^{+}→f_{0}(500)e^{+}ν_{e}, f_{0}(500)→π^{+}π^{-} decay is set to be 3.3×10^{-4} at 90% confidence level.

First Observation of a Three-Resonance Structure in e

We report the measurement of the inclusive cross sections for e^{+}e^{-}→nOCH (where nOCH denotes non-open charm hadrons) with improved precision at center-of-mass (c.m.) energies from 3.645 to 3.871 GeV. We observe three resonances: R(3760), R(3780), and R(3810) with significances of 8.1σ, 13.7σ, and 8.8σ, respectively. The R(3810) state is observed for the first time, while the R(3760) and R(3780) states are observed for the first time in the nOCH cross sections. Two sets of resonance parameters describe the energy-dependent line shape of the cross sections well. In set I [set II], the R(3810) state has mass (3805.7±1.1±2.7) [(3805.7±1.1±2.7)] MeV/c^{2}, total width (11.6±2.9±1.9) [(11.5±2.8±1.9)] MeV, and an electronic width multiplied by the nOCH decay branching fraction of (10.9±3.8±2.5) [(11.0±3.4±2.5)] eV. In addition, we measure the branching fractions B[R(3760)→nOCH]=(25.2±16.1±30.4)%[(6.4±4.8±7.7)%] and B[R(3780)→nOCH]=(12.3±6.6±8.3)%[(10.4±4.8±7.0)%] for the first time. The R(3760) state can be interpreted as an open-charm (OC) molecular state, but containing a simple four-quark state component. The R(3810) state can be interpreted as a hadrocharmonium state.

First Measurement of the Decay Asymmetry in the Pure W-Boson-Exchange Decay Λ_{c}

Based on 4.4 fb^{-1} of e^{+}e^{-} annihilation data collected at the center-of-mass energies between 4.60 and 4.70 GeV with the BESIII detector at the BEPCII collider, the pure W-boson-exchange decay Λ_{c}^{+}→Ξ^{0}K^{+} is studied with a full angular analysis. The corresponding decay asymmetry is measured for the first time to be α_{Ξ^{0}K^{+}}=0.01±0.16(stat)±0.03(syst). This result reflects the noninterference effect between the S- and P-wave amplitudes. The phase shift between S- and P-wave amplitudes has two solutions, which are δ_{p}-δ_{s}=-1.55±0.25(stat)±0.05(syst) rad or 1.59±0.25(stat)±0.05(syst) rad.

Search for a Sub-eV Sterile Neutrino Using Daya Bay's Full Dataset.

This Letter presents results of a search for the mixing of a sub-eV sterile neutrino with three active neutrinos based on the full data sample of the Daya Bay Reactor Neutrino Experiment, collected during 3158 days of detector operation, which contains 5.55×10^{6} reactor ν[over ¯]_{e} candidates identified as inverse beta-decay interactions followed by neutron capture on gadolinium. The analysis benefits from a doubling of the statistics of our previous result and from improvements of several important systematic uncertainties. No significant oscillation due to mixing of a sub-eV sterile neutrino with active neutrinos was found. Exclusion limits are set by both Feldman-Cousins and CLs methods. Light sterile neutrino mixing with sin^{2}2θ_{14}≳0.01 can be excluded at 95% confidence level in the region of 0.01 eV^{2}≲|Δm_{41}^{2}|≲0.1 eV^{2}. This result represents the world-leading constraints in the region of 2×10^{-4} eV^{2}≲|Δm_{41}^{2}|≲0.2 eV^{2}.

Quantification of left atrial strain in patients with idiopathic inflammatory myopathy using cardiovascular magnetic resonance feature tracking.

BACKGROUND: Left atrial (LA) dysfunction is involved in idiopathic inflammatory myopathy (IIM). Multiparametric cardiovascular magnetic resonance (CMR) strain imaging is a feasible and reproducible tool for examining global and regional LA functions, as well as left ventricular (LV) function in IIM patients. AIM: The aim of this study was to evaluate the feasibility and reproducibility of LA strain occurrence and strain rate for LA function assessment using CMR in IIM cases. MATERIALS AND METHODS: A total of 36 IIM and 42 healthy control cases were included. Baseline ventricular function was comparatively assessed in both groups. LA strain occurrence and strain rate were examined by cine cardiac magnetic resonance imaging [MRI] utilizing an in-house semiautomated technique. LA global function indexes were quantitated, including reservoir, conduit, and booster-pump functions. RESULTS: A total of 78 participants were enrolled in this study. There was no significant difference in left/right ventricular routine functions between IIM patients and control individuals (p>0.05); the same results (p>0.05) was also observed between patients with high hs-cTnI and normal. However, LV mass index had significant difference (p1=0.003, p2<0.01). Compared with IIM patients and control individuals, only total strain (εs) (p4=0.046) and passive strain (εe) (p4=0.002) showed significant difference, and in cases with high hs-cTnI and normal hs-cTnI, there are differences for εs (p3=0.012) and εe (p4=0.047). The strongest association was found between εe and LV ejection fraction (LVEF) (r=0.581, p<0.01). CONCLUSION: IIM cases have altered LA reservoir and conduit functions, and LA strain could reflect LA function.

A novel strategy for predicting the efficacy of temozolomide treatment for metastatic pheochromocytomas/paragangliomas.

BACKGROUND: There are few studies on the efficacy of temozolomide (TMZ) in the treatment of Metastatic pheochromocytoma / paraganglioma (MPP) patients. And it remains unclear which MPP patients may benefit from TMZ treatment. METHODS: This was a prospective study. MPP patients were enrolled. Patients were treated with TMZ until disease progression or intolerable toxicities. The primary endpoints were disease control rate (DCR) and objective response rate (ORR). Secondary endpoints included biochemical response rate progression-free survival (PFS) and safety. We compared the difference between effective and ineffective groups, to explore which patients are more suitable for TMZ treatment. RESULTS: 62 patients with MPP were enrolled and tumor response were evaluated in 54 patients. The DCR was 83% (35/42), and the ORR was 24% (10/41) among the progressive patients. PFS was 25.2 ± 3.1 months. The most common adverse event was nausea (41/55). We found that 92.9% (13/14) of patients with MGMT methylation greater than 7% respond to treatment. For the patients with MGMT methylation less than 7%, Ki-67 index could be used to guide the use of TMZ in these patients. Among the patients with Ki-67 index less than 5%, 66% (8/12) patients showed respond to treatment, and only 33% (4/12) patients with Ki-67 index more than 5% showed respond to TMZ. CONCLUSIONS: This study indicated that TMZ is a potential choice for the treatment of MPP with the high ability on disease control and well tolerability. We recommended to MGMT methylation analysis test and Ki-67 index to guide TMZ application.

Novel alternative tools for metastatic pheochromocytomas/paragangliomas prediction.

OBJECTIVE: The existing prediction models for metastasis in pheochromocytomas/paragangliomas (PPGLs) showed high heterogeneity in different centers. Therefore, this study aimed to establish new prediction models integrating multiple variables based on different algorithms. DESIGN AND METHODS: Data of patients with PPGLs undergoing surgical resection at the Peking Union Medical College Hospital from 2007 to 2022 were collected retrospectively. Patients were randomly divided into the training and testing sets in a ratio of 7:3. Subsequently, decision trees, random forest, and logistic models were constructed for metastasis prediction with the training set and Cox models for metastasis-free survival (MFS) prediction with the total population. Additionally, Ki-67 index and tumor size were transformed into categorical variables for adjusting models. The testing set was used to assess the discrimination and calibration of models and the optimal models were visualized as nomograms. Clinical characteristics and MFS were compared between patients with and without risk factors. RESULTS: A total of 198 patients with 59 cases of metastasis were included and classified into the training set (n = 138) and testing set (n = 60). Among all models, the logistic regression model showed the best discrimination for metastasis prediction with an AUC of 0.891 (95% CI, 0.793-0.990), integrating SDHB germline mutations [OR: 96.72 (95% CI, 16.61-940.79)], S-100 (-) [OR: 11.22 (95% CI, 3.04-58.51)], ATRX (-) [OR: 8.42 (95% CI, 2.73-29.24)] and Ki-67 ≥ 3% [OR: 7.98 (95% CI, 2.27-32.24)] evaluated through immunohistochemistry (IHC), and tumor size ≥ 5 cm [OR: 4.59 (95% CI, 1.34-19.13)]. The multivariate Cox model including the above risk factors also showed a high C-index of 0.860 (95% CI, 0.810-0.911) in predicting MFS after surgery. Furthermore, patients with the above risk factors showed a significantly poorer MFS (P ≤ 0.001). CONCLUSIONS: Models established in this study provided alternative and reliable tools for clinicians to predict PPGLs patients' metastasis and MFS. More importantly, this study revealed for the first time that IHC of ATRX could act as an independent predictor of metastasis in PPGLs.