De novo design of drug-binding proteins with predictable binding energy and specificity

The de novo design of small molecule-binding proteins has seen exciting recent progress; however, high-affinity binding and tunable specificity typically require laborious screening and optimization after computational design. We developed a computational procedure to design a protein that recognizes a common pharmacophore in a series of poly(ADP-ribose) polymerase-1 inhibitors. One of three designed proteins bound different inhibitors with affinities ranging from <5 nM to low micromolar. X-ray crystal structures confirmed the accuracy of the designed protein-drug interactions. Molecular dynamics simulations informed the role of water in binding. Binding free energy calculations performed directly on the designed models were in excellent agreement with the experimentally measured affinities. We conclude that de novo design of high-affinity small molecule-binding proteins with tuned interaction energies is feasible entirely from computation.

The influence of respiratory movement on preoperative CT-guided localization of lung nodules

AIM

To evaluate the motion amplitude of lung nodules in different locations during preoperative computed tomography (CT)-guided localization, and the influence of respiratory movement on CT-guided percutaneous lung puncture.

MATERIALS AND METHODS

A consecutive cohort of 398 patients (123 men and 275 women with a mean age of 53.9 ± 10.7 years) who underwent preoperative CT-guided lung nodule localization from May 2021 to Apr 2022 were included in this retrospective study. The respiratory movement-related nodule amplitude in the cranial-caudal direction during the CT scan, characteristics of patients, lesions, and procedures were statistically analyzed. Univariate and multivariate logistic regression analyses were used to evaluate the influence of these factors on CT-guided localization.

RESULTS

The nodule motion distribution showed a statistically significant correlation within the upper/middle (lingular) and lower lobes (p<0.001). Motion amplitude was an independent risk factor for CT scan times (p=0.011) and procedure duration (p=0.016), but not for the technical failure rates or the incidence of complications. Puncture depth was an independent risk factor for the CT scan times, procedure duration, technical failure rates, and complications (p<0.01). Female, prone, and supine (as opposed to lateral) positions were significant protective factors for pneumothorax, while the supine position was an independent risk factor for parenchymal hemorrhage (p=0.025).

CONCLUSION

Respiratory-induced motion amplitude of nodules was greater in the lower lobes, resulting in more CT scan times/radiation dose and longer localization duration, but showed no statistically significant influence on the technical success rates or the incidence of complications during preoperative CT-guided localization.

De novo design of drug-binding proteins with predictable binding energy and specificity

The de novo design of small-molecule-binding proteins has seen exciting recent progress; however, the ability to achieve exquisite affinity for binding small molecules while tuning specificity has not yet been demonstrated directly from computation. Here, we develop a computational procedure that results in the highest affinity binders to date with predetermined relative affinities, targeting a series of PARP1 inhibitors. Two of four designed proteins bound with affinities ranging from < 5 nM to low μM, in a predictable manner. X-ray crystal structures confirmed the accuracy of the designed protein-drug interactions. Molecular dynamics simulations informed the role of water in binding. Binding free-energy calculations performed directly on the designed models are in excellent agreement with the experimentally measured affinities, suggesting that the de novo design of small-molecule-binding proteins with tuned interaction energies is now feasible entirely from computation. We expect these methods to open many opportunities in biomedicine, including rapid sensor development, antidote design, and drug delivery vehicles.

Hyperon Polarization along the Beam Direction Relative to the Second and Third Harmonic Event Planes in Isobar Collisions at sqrt[s_{NN}]=200 GeV

The polarization of Λ and Λ[over ¯] hyperons along the beam direction has been measured relative to the second and third harmonic event planes in isobar Ru+Ru and Zr+Zr collisions at sqrt[s_{NN}]=200  GeV. This is the first experimental evidence of the hyperon polarization by the triangular flow originating from the initial density fluctuations. The amplitudes of the sine modulation for the second and third harmonic results are comparable in magnitude, increase from central to peripheral collisions, and show a mild p_{T} dependence. The azimuthal angle dependence of the polarization follows the vorticity pattern expected due to elliptic and triangular anisotropic flow, and qualitatively disagrees with most hydrodynamic model calculations based on thermal vorticity and shear induced contributions. The model results based on one of existing implementations of the shear contribution lead to a correct azimuthal angle dependence, but predict centrality and p_{T} dependence that still disagree with experimental measurements. Thus, our results provide stringent constraints on the thermal vorticity and shear-induced contributions to hyperon polarization. Comparison to previous measurements at RHIC and the LHC for the second-order harmonic results shows little dependence on the collision system size and collision energy.

Kinome Reprogramming Is a Targetable Vulnerability in ESR1 Fusion-Driven Breast Cancer

Transcriptionally active ESR1 fusions (ESR1-TAF) are a potent cause of breast cancer endocrine therapy (ET) resistance. ESR1-TAFs are not directly druggable because the C-terminal estrogen/anti-estrogen-binding domain is replaced with translocated in-frame partner gene sequences that confer constitutive transactivation. To discover alternative treatments, a mass spectrometry (MS)-based kinase inhibitor pulldown assay (KIPA) was deployed to identify druggable kinases that are upregulated by diverse ESR1-TAFs. Subsequent explorations of drug sensitivity validated RET kinase as a common therapeutic vulnerability despite remarkable ESR1-TAF C-terminal sequence and structural diversity. Organoids and xenografts from a pan-ET-resistant patient-derived xenograft model that harbors the ESR1-e6>YAP1 TAF were concordantly inhibited by the selective RET inhibitor pralsetinib to a similar extent as the CDK4/6 inhibitor palbociclib. Together, these findings provide preclinical rationale for clinical evaluation of RET inhibition for the treatment of ESR1-TAF-driven ET-resistant breast cancer.

SIGNIFICANCE

Kinome analysis of ESR1 translocated and mutated breast tumors using drug bead-based mass spectrometry followed by drug-sensitivity studies nominates RET as a therapeutic target. See related commentary by Wu and Subbiah, p. 3159.

Tumor Treatment Response Assessed During the Concurrent Chemoradiotherapy for Nasopharyngeal Patients

PURPOSE/OBJECTIVE(S)

To evaluate intratumoral treatment response distribution with using FDG-PET/CT during the chemoradiotherapy of nasopharyngeal patients (NPC).

MATERIALS/METHODS

A total of 5 of 30 patients with stage III-IVA NPC were enrolled in the institutional protocol for induction/concurrent chemoradiotherapy with radiation dose of 70 Gy in 33 fractions. For each patient, a pre-radiation treatment FDG-PET/MRI image (SUV0) and a mid-treatment image (SUVm) at the treatment dose of 31.8 Gy were obtained. Followed by deformable PET/MRI registration between SUV0 and SUVm, the tumor voxel SUV reduction ratio was obtained to construct a tumor dose response matrix (DRM). Tumor SUVavid was also constructed by limiting tumor voxel SUVm > a given value. Spatial correlations of the tumor SUV0, SUVm, SUVavid and DRM were determined.

RESULTS

The mean and coefficient variation (CV) of the SUV0, SUVm and DRM for all tumors were 5.05 (52%), 2.72 (49%) and 0.64 (63%) (Table contains the individual data), which were smaller than those on the SUVs of head-n-neck HPV+ patients reported previously due to the induction chemotherapy, but had much larger DRM mean and CV. The inter-tumoral CVs of SUV0 and DRM were 29% and 27%, which were much lower than those of the intra-tumoral CVs 43% and 57%. Meanwhile, the intra-tumoral variations on SUV0 was smaller than the one of head-neck HPV+ patients, but the DRM intra-variation was much larger. There was a weak correlation between SUV0 and SUVm with the correlation coefficient 0.13, a medium correlation of -0.55 between SUV0 and DRM, but a strong correlation, 0.72, between SUVm and DRM. However, the spatial correlation between tumor DRM and SUVavid was getting weaker as the SUVavid value increasing and equal 0.47 with SUVavid value > 3.

CONCLUSION

The spatial dose response DRM for NPC in the concurrent chemoradiotherapy was relatively high, while had relatively low baseline tumor metabolic activity SUV0. It was most likely due to the induction chemotherapy. In addition, the tumor dose response showed vary large intra-tumoral variation. The high correlations between DRM and SUVm imply that SUVavid could be used partially to guide adaptive modification of NPC treatment with carefully selected boundary value.

Tumor Treatment Response Assessed During the Chemo-Radiotherapy for Locally Advanced NSCLC

PURPOSE/OBJECTIVE(S)

To evaluate the capability of assessing intratumoral treatment response distribution with using FDG-PET/CT during the chemoradiotherapy of locally advanced NSCLC.

MATERIALS/METHODS

Twelve of total 50 patients with stage III NSCLC were enrolled in the institutional protocol for concurrent chemoradiotherapy with treatment dose of 54-60 Gy in 27-30 fractions. For each patient, a pre-treatment FDG-PET/CT image (SUV0) and a mid-treatment image (SUVm) obtained within the treatment dose of 24 ∼ 46 Gy were obtained. Followed by deformable PET/CT registration between SUV0 and SUVm, the tumor voxel SUV reduction ratio was obtained to construct a tumor dose response matrix (DRM). Tumor SUVavid was also constructed by limiting tumor voxel SUVm > a given value. Spatial correlations of the tumor SUV0, SUVm, SUVavid and DRM were determined.

RESULTS

The mean and coefficient variation (CV) of the SUV0, SUVm and DRM for all tumors were 6.56(64%), 2.82(59%) and 0.52(70%) (Table contains the individual data), which were like those on the SUVs and the mean DRM of head-neck HPV- patients reported previously, but much larger on the DRM variation. The inter-tumoral CVs of SUV0 and DRM were 17% and 43%, which were much smaller than those of the intra-tumoral CVs 61% and 55%. Meanwhile, the intra-tumoral variations on both SUV0 and DRM were much larger than those of head-neck HPV- patients. There was a weak correlation between SUV0 and SUVm with the correlation coefficient 0.32, a medium correlation of -0.51 between SUV0 and DRM; 0.58 between SUVm and DRM. It implies that the rule of tumor dose response DRM on treatment modification decision cannot be fully replaced by either SUV0 or SUVm. The spatial correlation between tumor DRM and SUVavid was 0.23 with SUVavid value > 3, which was getting weaker when increasing SUVavid value.

CONCLUSION

Spatial dose response for NSCLC assessed using FDG-PET/CT feedback demonstrated high treatment resistant patterns, which had a large intra-tumoral variation. In addition, the medium correlations of DRM vs SUV0 and DRM vs SUVm imply that all these factors could be used to guide adaptive modification of NSCLC treatment.

Erratum: Global Polarization of Ξ and Ω Hyperons in Au+Au Collisions at sqrt[s_{NN}]=200 GeV [Phys. Rev. Lett. 126, 162301 (2021)]

This corrects the article DOI: 10.1103/PhysRevLett.126.162301.

Measurements of the Elliptic and Triangular Azimuthal Anisotropies in Central ^{3}He+Au, d+Au and p+Au Collisions at sqrt[s_{NN}]=200 GeV

The elliptic (v_{2}) and triangular (v_{3}) azimuthal anisotropy coefficients in central ^{3}He+Au, d+Au, and p+Au collisions at sqrt[s_{NN}]=200  GeV are measured as a function of transverse momentum (p_{T}) at midrapidity (|η|<0.9), via the azimuthal angular correlation between two particles both at |η|<0.9. While the v_{2}(p_{T}) values depend on the colliding systems, the v_{3}(p_{T}) values are system independent within the uncertainties, suggesting an influence on eccentricity from subnucleonic fluctuations in these small-sized systems. These results also provide stringent constraints for the hydrodynamic modeling of these systems.

Observation of Directed Flow of Hypernuclei _{Λ}^{3}H and _{Λ}^{4}H in sqrt[s_{NN}]=3 GeV Au+Au Collisions at RHIC

We report here the first observation of directed flow (v_{1}) of the hypernuclei _{Λ}^{3}H and _{Λ}^{4}H in mid-central Au+Au collisions at sqrt[s_{NN}]=3  GeV at RHIC. These data are taken as part of the beam energy scan program carried out by the STAR experiment. From 165×10^{6} events in 5%-40% centrality, about 8400 _{Λ}^{3}H and 5200 _{Λ}^{4}H candidates are reconstructed through two- and three-body decay channels. We observe that these hypernuclei exhibit significant directed flow. Comparing to that of light nuclei, it is found that the midrapidity v_{1} slopes of _{Λ}^{3}H and _{Λ}^{4}H follow baryon number scaling, implying that the coalescence is the dominant mechanism for these hypernuclei production in the 3 GeV Au+Au collisions.

Beam Energy Dependence of Triton Production and Yield Ratio (N_{t}×N_{p}/N_{d}^{2}) in Au+Au Collisions at RHIC

We report the triton (t) production in midrapidity (|y|<0.5) Au+Au collisions at sqrt[s_{NN}]=7.7-200  GeV measured by the STAR experiment from the first phase of the beam energy scan at the Relativistic Heavy Ion Collider. The nuclear compound yield ratio (N_{t}×N_{p}/N_{d}^{2}), which is predicted to be sensitive to the fluctuation of local neutron density, is observed to decrease monotonically with increasing charged-particle multiplicity (dN_{ch}/dη) and follows a scaling behavior. The dN_{ch}/dη dependence of the yield ratio is compared to calculations from coalescence and thermal models. Enhancements in the yield ratios relative to the coalescence baseline are observed in the 0%-10% most central collisions at 19.6 and 27 GeV, with a significance of 2.3σ and 3.4σ, respectively, giving a combined significance of 4.1σ. The enhancements are not observed in peripheral collisions or model calculations without critical fluctuation, and decreases with a smaller p_{T} acceptance. The physics implications of these results on the QCD phase structure and the production mechanism of light nuclei in heavy-ion collisions are discussed.

Abstract 5011: Targeting kinome reprogramming in ESR1 fusion-driven breast cancer

Background: Transcriptionally active ESR1 gene fusions (ESR1-TAF) are a potent cause of estrogen receptor alpha-positive (ERα+) breast cancer endocrine therapy (ET) resistance. These ESR1-TAF are gain-of-function mutations, exhibiting estrogen-independent cell growth, motility and ET resistance. They are not directly druggable because the ERα C-terminal ligand binding domain (LBD) encoding sequence is replaced with a translocated in-frame partner gene sequence. Herein we utilized proteomic approaches to develop novel targeted therapies against ESR1-TAF driven tumorigenesis.

Methods: ESR1 fusion cDNA constructs were expressed in ERα+ breast cancer cell lines (T47D and MCF7). Cell growth was assayed by an Alamar blue assay. A mass spectrometry (MS)-based Kinase Inhibitor Pulldown Assay (KIPA) was employed to identify druggable kinases that are commonly upregulated by diverse ESR1-TAFs. A panel of 22 ERα+ patient-derived xenograft (PDX) models were profiled using proteomics and phosphoproteomics to identify models with sensitivity to RET kinase inhibition.

Results: KIPA detected an increased abundance of a receptor tyrosine kinase, RET, in T47D cells expressing ESR1-TAFs in an estrogen-independent manner, compared to stable cell lines expressing transcriptionally inactive fusions as well as wild-type ERα protein. Interestingly, RET was also increased when constitutive activating ERα LBD point mutants, Y537S and D538G, were expressed in breast cancer cells. Inhibition of the RET kinase in vitro by repurposing pralsetinib, an FDA-approved RET inhibitor for advanced thyroid and non-small-cell lung cancers, demonstrated a significant reduction in the growth of cells expressing ESR1-TAFs and ERα LBD mutants. These data nominate RET kinase as a common therapeutic vulnerability for ESR1-TAF expressing breast cancers. Proteomic profiling of 22 biologically heterogenous ERα+ PDX tumors defined targetable pathways and predicted tumor subsets that were responsive to RET inhibition therapy. Organoids and xenografts from the pan-ET resistant WHIM18 PDX (that expresses the ESR1-YAP1 TAF) were inhibited by pralsetinib to a similar extent as the CDK4/6 inhibitor palbociclib. These data provide key preclinical rationale for the consideration of RET inhibition for the treatment of ESR1-TAF-driven ET-resistant breast cancer. Interestingly, the growth of WHIM37 PDX (that expresses ERα D538G) that had low level of RET and high level of GFRα-1, the co-receptor of RET, was also suppressed by pralsetinib. This data suggests that either RET or GFRα-1 is a predictive biomarker for RET inhibitor efficacy.

Conclusions: Kinome analysis of ESR1 translocated breast tumors using KIPA followed by drug sensitivity studies nominated RET as a new therapeutic target for ET-resistant ERα+ breast cancer.

Citation Format: Xuxu Gou, Beom-Jun Kim, Meenakshi Anurag, Jonathan T. Lei, Meggie N. Young, Matthew V. Holt, Diana Fandino, Craig T. Vollert, Purba Singh, Mohammad A. Alzubi, Anna Malovannaya, Lacey E. Dobrolecki, Michael T. Lewis, Shunqiang Li, Matthew J. Ellis, Charles E. Foulds. Targeting kinome reprogramming in ESR1 fusion-driven breast cancer. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 5011.

Measurement of Sequential ϒ Suppression in Au+Au Collisions at sqrt[s_{NN}]=200 GeV with the STAR Experiment

We report on measurements of sequential ϒ suppression in Au+Au collisions at sqrt[s_{NN}]=200  GeV with the STAR detector at the Relativistic Heavy Ion Collider (RHIC) through both the dielectron and dimuon decay channels. In the 0%-60% centrality class, the nuclear modification factors (R_{AA}), which quantify the level of yield suppression in heavy-ion collisions compared to p+p collisions, for ϒ(1S) and ϒ(2S) are 0.40±0.03(stat)±0.03(sys)±0.09(norm) and 0.26±0.08(stat)±0.02(sys)±0.06(norm), respectively, while the upper limit of the ϒ(3S) R_{AA} is 0.17 at a 95% confidence level. This provides experimental evidence that the ϒ(3S) is significantly more suppressed than the ϒ(1S) at RHIC. The level of suppression for ϒ(1S) is comparable to that observed at the much higher collision energy at the Large Hadron Collider. These results point to the creation of a medium at RHIC whose temperature is sufficiently high to strongly suppress excited ϒ states.

Beam Energy Dependence of Fifth- and Sixth-Order Net-Proton Number Fluctuations in Au+Au Collisions at RHIC

We report the beam energy and collision centrality dependence of fifth and sixth order cumulants (C_{5}, C_{6}) and factorial cumulants (κ_{5}, κ_{6}) of net-proton and proton number distributions, from center-of-mass energy (sqrt[s_{NN}]) 3 GeV to 200 GeV Au+Au collisions at RHIC. Cumulant ratios of net-proton (taken as proxy for net-baryon) distributions generally follow the hierarchy expected from QCD thermodynamics, except for the case of collisions at 3 GeV. The measured values of C_{6}/C_{2} for 0%-40% centrality collisions show progressively negative trend with decreasing energy, while it is positive for the lowest energy studied. These observed negative signs are consistent with QCD calculations (for baryon chemical potential, μ_{B}≤110  MeV) which contains the crossover transition range. In addition, for energies above 7.7 GeV, the measured proton κ_{n}, within uncertainties, does not support the two-component (Poisson+binomial) shape of proton number distributions that would be expected from a first-order phase transition. Taken in combination, the hyperorder proton number fluctuations suggest that the structure of QCD matter at high baryon density, μ_{B}∼750  MeV at sqrt[s_{NN}]=3  GeV is starkly different from those at vanishing μ_{B}∼24  MeV at sqrt[s_{NN}]=200  GeV and higher collision energies.

Collision-System and Beam-Energy Dependence of Anisotropic Flow Fluctuations

Elliptic flow measurements from two-, four-, and six-particle correlations are used to investigate flow fluctuations in collisions of U+U at sqrt[s_{NN}]=193  GeV, Cu+Au at sqrt[s_{NN}]=200  GeV and Au+Au spanning the range sqrt[s_{NN}]=11.5-200  GeV. The measurements show a strong dependence of the flow fluctuations on collision centrality, a modest dependence on system size, and very little if any, dependence on particle species and beam energy. The results, when compared to similar LHC measurements, viscous hydrodynamic calculations, and trento model eccentricities, indicate that initial-state-driven fluctuations predominate the flow fluctuations generated in the collisions studied.

Proteogenomic Markers of Chemotherapy Resistance and Response in Triple-Negative Breast Cancer

Microscaled proteogenomics was deployed to probe the molecular basis for differential response to neoadjuvant carboplatin and docetaxel combination chemotherapy for triple-negative breast cancer (TNBC). Proteomic analyses of pretreatment patient biopsies uniquely revealed metabolic pathways, including oxidative phosphorylation, adipogenesis, and fatty acid metabolism, that were associated with resistance. Both proteomics and transcriptomics revealed that sensitivity was marked by elevation of DNA repair, E2F targets, G2-M checkpoint, interferon-gamma signaling, and immune-checkpoint components. Proteogenomic analyses of somatic copy-number aberrations identified a resistance-associated 19q13.31-33 deletion where LIG1, POLD1, and XRCC1 are located. In orthogonal datasets, LIG1 (DNA ligase I) gene deletion and/or low mRNA expression levels were associated with lack of pathologic complete response, higher chromosomal instability index (CIN), and poor prognosis in TNBC, as well as carboplatin-selective resistance in TNBC preclinical models. Hemizygous loss of LIG1 was also associated with higher CIN and poor prognosis in other cancer types, demonstrating broader clinical implications.

SIGNIFICANCE

Proteogenomic analysis of triple-negative breast tumors revealed a complex landscape of chemotherapy response associations, including a 19q13.31-33 somatic deletion encoding genes serving lagging-strand DNA synthesis (LIG1, POLD1, and XRCC1), that correlate with lack of pathologic response, carboplatin-selective resistance, and, in pan-cancer studies, poor prognosis and CIN. This article is highlighted in the In This Issue feature, p. 2483.

Evidence for Nonlinear Gluon Effects in QCD and Their Mass Number Dependence at STAR

The STAR Collaboration reports measurements of back-to-back azimuthal correlations of di-π^{0}s produced at forward pseudorapidities (2.6<η<4.0) in p+p, p+Al, and p+Au collisions at a center-of-mass energy of 200 GeV. We observe a clear suppression of the correlated yields of back-to-back π^{0} pairs in p+Al and p+Au collisions compared to the p+p data. The observed suppression of back-to-back pairs as a function of transverse momentum suggests nonlinear gluon dynamics arising at high parton densities. The larger suppression found in p+Au relative to p+Al collisions exhibits a dependence of the saturation scale Q_{s}^{2} on the mass number A. A linear scaling of the suppression with A^{1/3} is observed with a slope of -0.09±0.01.

Abstract 1010: LIG1 deletion predicts chemotherapy resistance, chromosomal instability, and poor prognosis in triple negative breast cancer

Introduction: Cytotoxic chemotherapy for sporadic Triple Negative Breast Cancer (TNBC) remains the standard of care and the recent approval for adjuvant PD1 therapy is not biomarker guided. Pathological complete response (pCR) is often not achieved and portends poor survival. Predictive markers for individual drugs have proven elusive.

Approach: Microscaled proteogenomics (MPG) was applied to snap-frozen TNBC clinical trial core needle biopsies obtained before treatment with carboplatin and docetaxel (WashU: NCT201404107 and BCM: NCT02544987). Clinical endpoints for discovery analysis were pathological complete response (pCR) and residual cancer burden (RCB). Standard non-parametric statistical tests were employed to identify proteogenomic features associated with these endpoints.

Results: Copy number aberrations (CNA) are a recurrent feature of TNBC and a potential driver of chemotherapy sensitivity. We therefore sought CNA with concordant changes at the mRNA and protein levels that also associate with pCR status. Genes located within a recurrent interstitial deletion at chromosomal location 19q13.3 were the most significantly down-regulated at mRNA and protein level in chemotherapy resistant cases. 19q13.3 encodes multiple DNA damage response (DDR) genes; however, only LIG1, a DNA ligase required for lagging strand synthesis and DNA repair, showed concordant changes at both the mRNA and protein level. In multiple independent TNBC data sets, LIG1 deletion was associated lack of pCR and poor metastasis-free survival. Additionally in the BrighTNess TNBC trial lower LIG1 mRNA levels were associated with increased chemotherapy resistance in the carboplatin containing arms (no pCR and residual cancer burden I-III; p=0.0008 and 0.003 respectively). In PDX-derived short-term cultures and PDXs treated with docetaxel or carboplatin, a specific association of carboplatin resistance with LIG1 deletion was observed. LIG1 depleted-tumors did not harbor elevated scores for homologous recombination defect signature, suggesting LIG1 loss is an orthogonal pathway for TNBC pathogenesis The high chromosomal instability index in LIG1 deletion tumors in our TNBC study was robustly reproduced in multiple datasets (including TCGA-BRCA ; Metastatic breast cancer project). LIG1 copy number deletion was also associated with poor progression free survival, and high chromosomal instability in multiple other cancers (including TCGA-​UCEC HR=2.23, TCGA-HNSC HR=1.46, TCGA-PRAD HR=2.07, TCGA- COAD HR=1.75 and TCGA-KIRP HR=4.00).

Conclusion: Deletion of LIG1 is associated with chromosomal instability in TNBC and occurs in tumors without genomic evidence for defects in homologous recombination. Other clinical features of LIG1 deleted TNBC and how LIG1 loss may cause chromosomal instability and tumorigenesis will be discussed.

Citation Format: Meenakshi Anurag, Eric Jaehnig, Jonathan Lei, Beom-Jun Kim, Anh Minh Tran Huynh, Yongchao Dou, Tanmayi Vashist, Erik Bergstrom, Xuxu Gou, Viktoriya Korchina, Donna Marie Muzny, Kristen Otte, Harshavardhan Doddapaneni, Lacey Dobrolecki, Gloria Vittone Echeverria, Bora Lim, Mothaffar Rimawi, Karsten Krug, Ian Hageman, Henry Rodriguez, Ana I. Robles, Tara Hiltke, Kent Osborne, Michael Gillette, George Miles, Steven Carr, Michael T Lewis, Bing Zhang, Foluso Ademuyiwa, Shankha Satpathy, Matthew J. Ellis. LIG1 deletion predicts chemotherapy resistance, chromosomal instability, and poor prognosis in triple negative breast cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 1010.

Targeting kinome reprogramming in ESR1 fusion-driven metastatic breast cancer

1085

Background: Genomic analysis has recently identified multiple ESR1 gene translocations in estrogen receptor-alpha positive (ERα+) metastatic breast cancer (MBC) that encode chimeric proteins whereby the ESR1 ligand binding domain is replaced by C-terminal sequences from many different gene partners. Transcriptionally active ESR1 fusions promoted hormone-independent cell growth, motility and resistance to endocrine therapy. The diversity of partner genes creates a considerable diagnostic challenge and no targeted treatments exist for ESR1 translocated tumors. Thus, we have established a transcriptional signature to diagnose the presence of an active ESR1 fusion (PMID: 34711608) and developed novel targeted therapies against ESR1 fusion-driven biology. Methods: Fifteen ESR1 fusion cDNA constructs were expressed in ER+ breast cancer cell lines by lentiviral transduction. Cell growth was assayed by Alamar blue assay. A mass spectrometry (MS)-based Kinase Inhibitor Pulldown Assay (KIPA) and tandem mass tag-based proteomics were performed to identify ESR1 fusion-driven druggable kinases for subsequent pharmacological inhibition. Results: KIPA profiling demonstrated an increase of multiple receptor tyrosine kinases including RET in T47D cells expressing active ESR1 fusions. Inhibition of RET by repurposing an FDA-approved drug significantly suppressed ESR1 fusion-driven cell growth in vitro, suggesting that despite marked diversity in the 3’ partners, common kinase activities were elevated and targetable. Proteogenomic profiling, including whole exome sequencing, RNA sequencing, and MS-based proteomics and phosphoproteomics were further performed on 22 ER+ patient-derived xenograft (PDX) tumors, which demonstrated different degrees of estradiol dependence. These integrated “omic” profiles defined targetable genes/pathways and predict tumor subsets that could be responsive to kinase inhibition therapy from this biologically heterogeneous panel of PDX tumors. WHIM18, a PDX naturally harboring the ESR1-YAP1 fusion showed elevated level of RET and CDK4/6 pathways. The tumor volumes were significantly reduced by the RET inhibitor. CDK4/6 inhibitor treatment showed similar tumor reductions to RET inhibition. Interestingly, WHIM9 PDX that expressed wild-type ESR1 conferred a comparable kinome profile to WHIM18. The tumor growth was significantly suppressed by RET or CDK4/6 inhibition. Therefore, pharmacological experiments validated proteogenomics-predicted drug response in two tested ER+ PDX models. Conclusions: Proteogenomics characterization of PDX tumors can drive clinical trial hypotheses. Here, we reveal therapeutic kinase vulnerabilities in ESR1 fusion-driven tumors as exemplified by RET inhibition, which will lay the framework for future clinical trials.

Measurements of Proton High-Order Cumulants in sqrt[s_{NN}]=3 GeV Au+Au Collisions and Implications for the QCD Critical Point

We report cumulants of the proton multiplicity distribution from dedicated fixed-target Au+Au collisions at sqrt[s_{NN}]=3.0  GeV, measured by the STAR experiment in the kinematic acceptance of rapidity (y) and transverse momentum (p_{T}) within -0.5<y<0 and 0.4<p_{T}<2.0  GeV/c. In the most central 0%-5% collisions, a proton cumulant ratio is measured to be C_{4}/C_{2}=-0.85±0.09 (stat)±0.82 (syst), which is 2σ below the Poisson baseline with respect to both the statistical and systematic uncertainties. The hadronic transport UrQMD model reproduces our C_{4}/C_{2} in the measured acceptance. Compared to higher energy results and the transport model calculations, the suppression in C_{4}/C_{2} is consistent with fluctuations driven by baryon number conservation and indicates an energy regime dominated by hadronic interactions. These data imply that the QCD critical region, if created in heavy-ion collisions, could only exist at energies higher than 3 GeV.

Measurements of _{Λ}^{3}H and _{Λ}^{4}H Lifetimes and Yields in Au+Au Collisions in the High Baryon Density Region

We report precision measurements of hypernuclei _{Λ}^{3}H and _{Λ}^{4}H lifetimes obtained from Au+Au collisions at sqrt[s_{NN}]=3.0  GeV and 7.2 GeV collected by the STAR experiment at the Relativistic Heavy Ion Collider, and the first measurement of _{Λ}^{3}H and _{Λ}^{4}H midrapidity yields in Au+Au collisions at sqrt[s_{NN}]=3.0  GeV. _{Λ}^{3}H and _{Λ}^{4}H, being the two simplest bound states composed of hyperons and nucleons, are cornerstones in the field of hypernuclear physics. Their lifetimes are measured to be 221±15(stat)±19(syst)  ps for _{Λ}^{3}H and 218±6(stat)±13(syst)  ps for _{Λ}^{4}H. The p_{T}-integrated yields of _{Λ}^{3}H and _{Λ}^{4}H are presented in different centrality and rapidity intervals. It is observed that the shape of the rapidity distribution of _{Λ}^{4}H is different for 0%-10% and 10%-50% centrality collisions. Thermal model calculations, using the canonical ensemble for strangeness, describes the _{Λ}^{3}H yield well, while underestimating the _{Λ}^{4}H yield. Transport models, combining baryonic mean-field and coalescence (jam) or utilizing dynamical cluster formation via baryonic interactions (phqmd) for light nuclei and hypernuclei production, approximately describe the measured _{Λ}^{3}H and _{Λ}^{4}H yields. Our measurements provide means to precisely assess our understanding of the fundamental baryonic interactions with strange quarks, which can impact our understanding of more complicated systems involving hyperons, such as the interior of neutron stars or exotic hypernuclei.

Probing the Gluonic Structure of the Deuteron with J/ψ Photoproduction in d+Au Ultraperipheral Collisions

Understanding gluon density distributions and how they are modified in nuclei are among the most important goals in nuclear physics. In recent years, diffractive vector meson production measured in ultraperipheral collisions (UPCs) at heavy-ion colliders has provided a new tool for probing the gluon density. In this Letter, we report the first measurement of J/ψ photoproduction off the deuteron in UPCs at the center-of-mass energy sqrt[s_{NN}]=200  GeV in d+Au collisions. The differential cross section as a function of momentum transfer -t is measured. In addition, data with a neutron tagged in the deuteron-going zero-degree calorimeter is investigated for the first time, which is found to be consistent with the expectation of incoherent diffractive scattering at low momentum transfer. Theoretical predictions based on the color glass condensate saturation model and the leading twist approximation nuclear shadowing model are compared with the data quantitatively. A better agreement with the saturation model has been observed. With the current measurement, the results are found to be directly sensitive to the gluon density distribution of the deuteron and the deuteron breakup process, which provides insights into the nuclear gluonic structure.

Search for the Chiral Magnetic Effect via Charge-Dependent Azimuthal Correlations Relative to Spectator and Participant Planes in Au+Au Collisions at sqrt[s_{NN}]=200 GeV

The chiral magnetic effect (CME) refers to charge separation along a strong magnetic field due to imbalanced chirality of quarks in local parity and charge-parity violating domains in quantum chromodynamics. The experimental measurement of the charge separation is made difficult by the presence of a major background from elliptic azimuthal anisotropy. This background and the CME signal have different sensitivities to the spectator and participant planes, and could thus be determined by measurements with respect to these planes. We report such measurements in Au+Au collisions at a nucleon-nucleon center-of-mass energy of 200 GeV at the Relativistic Heavy-Ion Collider. It is found that the charge separation, with the flow background removed, is consistent with zero in peripheral (large impact parameter) collisions. Some indication of finite CME signals is seen in midcentral (intermediate impact parameter) collisions. Significant residual background effects may, however, still be present.

Measurement of the Sixth-Order Cumulant of Net-Proton Multiplicity Distributions in Au+Au Collisions at sqrt[s_{NN}]=27, 54.4, and 200 GeV at RHIC

According to first-principle lattice QCD calculations, the transition from quark-gluon plasma to hadronic matter is a smooth crossover in the region μ_{B}≤T_{c}. In this range the ratio, C_{6}/C_{2}, of net-baryon distributions are predicted to be negative. In this Letter, we report the first measurement of the midrapidity net-proton C_{6}/C_{2} from 27, 54.4, and 200 GeV Au+Au collisions at the Relativistic Heavy Ion Collider (RHIC). The dependence on collision centrality and kinematic acceptance in (p_{T}, y) are analyzed. While for 27 and 54.4 GeV collisions the C_{6}/C_{2} values are close to zero within uncertainties, it is observed that for 200 GeV collisions, the C_{6}/C_{2} ratio becomes progressively negative from peripheral to central collisions. Transport model calculations without critical dynamics predict mostly positive values except for the most central collisions within uncertainties. These observations seem to favor a smooth crossover in the high-energy nuclear collisions at top RHIC energy.

Transcriptional reprogramming differentiates active from inactive ESR1 fusions in endocrine therapy-refractory metastatic breast cancer

Genomic analysis has recently identified multiple ESR1 gene translocations in estrogen receptor-alpha positive (ERα+) metastatic breast cancer (MBC) that encode chimeric proteins whereby the ESR1 ligand binding domain (LBD) is replaced by C-terminal sequences from many different gene partners. Here we functionally screened 15 ESR1 fusions and identified 10 that promoted estradiol-independent cell growth, motility, invasion, EMT and resistance to fulvestrant. RNA sequencing identified a gene expression pattern specific to functionally active ESR1 gene fusions that was subsequently reduced to a diagnostic 24-gene signature. This signature was further examined in 20 ERα+ patient-derived xenografts (PDXs) and in 55 ERα+ MBC samples. The 24-gene signature successfully identified cases harboring ESR1 gene fusions and also accurately diagnosed the presence of activating ESR1 LBD point mutations. Therefore, the 24-gene signature represents an efficient approach to screening samples for the presence of diverse somatic ESR1 mutations and translocations that drive endocrine treatment failure in MBC.

Observation of D_{s}^{±}/D^{0} Enhancement in Au+Au Collisions at sqrt[s_{NN}]=200 GeV

We report on the first measurement of charm-strange meson D_{s}^{±} production at midrapidity in Au+Au collisions at sqrt[s_{NN}]=200  GeV from the STAR experiment. The yield ratio between strange (D_{s}^{±}) and nonstrange (D^{0}) open-charm mesons is presented and compared to model calculations. A significant enhancement, relative to a pythia simulation of p+p collisions, is observed in the D_{s}^{±}/D^{0} yield ratio in Au+Au collisions over a large range of collision centralities. Model calculations incorporating abundant strange-quark production in the quark-gluon plasma and coalescence hadronization qualitatively reproduce the data. The transverse-momentum integrated yield ratio of D_{s}^{±}/D^{0} at midrapidity is consistent with a prediction from a statistical hadronization model with the parameters constrained by the yields of light and strange hadrons measured at the same collision energy. These results suggest that the coalescence of charm quarks with strange quarks in the quark-gluon plasma plays an important role in D_{s}^{±}-meson production in heavy-ion collisions.

Abstract 742: The integration of a structure-function rule and a transcriptional signature to assign ESR1 fusion activity in metastatic breast cancer

Background: In-frame chimeric ESR1 fusion proteins can cause endocrine therapy (ET) resistance in estrogen receptor-positive (ER+) metastatic breast cancer (MBC) (PMC6171747). Here, we further investigated multiple new examples of ESR1 fusions in order to differentiate between transcriptionally active “driver” ESR1 fusions from inactive, likely “passenger” events.

Methods: ESR1 fusions were identified by RNA-sequencing (RNA-seq) or whole genome sequencing in MBC (PMC5913625, PMC6872491 and unpublished). ESR1 fusion cDNA constructs were expressed in ER+ breast cancer cell lines. Cell proliferation and cell motility were assessed. RNA-seq followed by qPCR validation was conducted to examine individual transcriptional profiles.

Results: All ESR1 fusions studied contained the first six exons of ESR1 (e6) fused in-frame to diverse partner gene sequences that replaced the ESR1 drug/ligand binding domain. Fusions with a transcription factor (TF) or coactivator (CoA) gene partner, for example ESR1-YAP1, ESR1-ARNT2-e18 and ESR1-SOX9 activated fulvestrant-resistant cell growth and hormone-independent cell motility. Other ESR1-e6 fusions, including ESR1-DAB2, ESR1-GYG1, ESR1-PCMT1 and ESR1-ARID1B did not show these phenotypes. We thus conclude that in-frame ESR1-e6 fusions arising from inter-chromosomal translocations with 3' TF/CoA partners are likely to be active. However, the active ESR1-PCDH11X fusion involves a proto-cadherin and is therefore an exception to this rule. This outlier emphasizes the need to establish additional approaches to determine ESR1 fusion activity. RNA-seq of T47D cells expressing active and inactive ESR1 fusions defined a gene signature associated with active ESR1 fusions, with activation of estrogen response and epithelial-to-mesenchymal transition (EMT) genes. This transcriptional signature was present in a patient-derived xenograft bearing the ESR1-YAP1 fusion, and thus potentially identifies the presence of an active ESR1 fusion protein. We subsequently identified two new ESR1 fusions involving recurrent TF/CoA partners, ESR1-ARNT2-e2 and ESR1-LPP. Both demonstrated tumor cell growth and motility activation, as predicted by the functional rule. The gene activation patterns were similar to the three other active fusions suggesting that despite marked diversity in the 3' partners, the transcriptional activities were similar and potentially diagnostic.

Conclusion: The integration of the structure-function rule and the stereotypic transcriptional signature distinguishes pathogenic ESR1 fusions from non-active passenger events, thus prioritizing patients bearing ESR1 translocation-driven tumors for targeted therapeutic approaches.

Citation Format: Xuxu Gou, Meenakshi Anurag, Jonathan T. Lei, Sinem Seker, Adrian V. Lee, Dan R. Robinson, Matthew J. Ellis. The integration of a structure-function rule and a transcriptional signature to assign ESR1 fusion activity in metastatic breast cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 742.

Global Polarization of Ξ and Ω Hyperons in Au+Au Collisions at sqrt[s_{NN}]=200 GeV

Global polarization of Ξ and Ω hyperons has been measured for the first time in Au+Au collisions at sqrt[s_{NN}]=200  GeV. The measurements of the Ξ^{-} and Ξ[over ¯]^{+} hyperon polarization have been performed by two independent methods, via analysis of the angular distribution of the daughter particles in the parity violating weak decay Ξ→Λ+π, as well as by measuring the polarization of the daughter Λ hyperon, polarized via polarization transfer from its parent. The polarization, obtained by combining the results from the two methods and averaged over Ξ^{-} and Ξ[over ¯]^{+}, is measured to be ⟨P_{Ξ}⟩=0.47±0.10(stat)±0.23(syst)% for the collision centrality 20%-80%. The ⟨P_{Ξ}⟩ is found to be slightly larger than the inclusive Λ polarization and in reasonable agreement with a multiphase transport model. The ⟨P_{Ξ}⟩ is found to follow the centrality dependence of the vorticity predicted in the model, increasing toward more peripheral collisions. The global polarization of Ω, ⟨P_{Ω}⟩=1.11±0.87(stat)±1.97(syst)% was obtained by measuring the polarization of daughter Λ in the decay Ω→Λ+K, assuming the polarization transfer factor C_{ΩΛ}=1.

Nonmonotonic Energy Dependence of Net-Proton Number Fluctuations

Nonmonotonic variation with collision energy (sqrt[s_{NN}]) of the moments of the net-baryon number distribution in heavy-ion collisions, related to the correlation length and the susceptibilities of the system, is suggested as a signature for the quantum chromodynamics critical point. We report the first evidence of a nonmonotonic variation in the kurtosis times variance of the net-proton number (proxy for net-baryon number) distribution as a function of sqrt[s_{NN}] with 3.1  σ significance for head-on (central) gold-on-gold (Au+Au) collisions measured solenoidal tracker at Relativistic Heavy Ion Collider. Data in noncentral Au+Au collisions and models of heavy-ion collisions without a critical point show a monotonic variation as a function of sqrt[s_{NN}].

Abstract PS17-03: Recurrent active ESR1 fusions render a diagnostic transcriptional signature in metastatic breast cancer

Background: We recently reported two ESR1 fusions (ESR1-YAP1 and ESR1-PCDH11X) that drive endocrine therapy (ET) resistance and metastasis in estrogen receptor positive (ER+) metastatic breast cancer (MBC) (PMC6171747). Here, we report the functional properties of additional ESR1 fusions in ET-resistant MBC with an emphasis on the identification of a transcriptional signature designed to diagnose the presence of an active ESR1 fusion for targeted therapies directed against ESR1 fusion-driven biology. Methods: ESR1 fusions were detected by RNA-seq in ER+ MBC samples. ESR1 fusions were reproduced as cDNA constructs and expressed in ER+ breast cancer cell lines. Hormone-independent cell growth was detected by an Alamar blue assay and activated cell motility by a scratch wound assay. The transcriptional properties of ESR1 fusions was studied by RNA-seq followed by qPCR-based validation. Signature performance was evaluated using a ROC analysis on ER+ patient derived xenografts (PDX) harboring a variety of ESR1 somatic events. Results: All ESR1 fusions studied encoded the first six exons of ESR1 fused in-frame to diverse partner genes, thus replacing the ESR1 drug/ligand binding domain (LBD). Fusions involving a known transcription factor (TF) or coactivator (CoA) gene, including ESR1-YAP1, ESR1-SOX9 and ESR1-ARNT2 drove fulvestrant-resistant cell growth and hormone-independent cell motility. Other ESR1-e6 fusions, including ESR1-DAB2, ESR1-GYG1, ESR1-PCMT1 and ESR1-ARID1B did not induce these properties. From these examples, a functional rule is emerging whereby inter-chromosomal ESR1 translocations fused in-frame to 3’ partner genes with a positive role in transcription are active. Intra-chromosomal fusions with genes with no transcriptional roles are likely inactive. The ESR1-PCDH11X fusion is an exception, suggesting the need for continued functional study of non-TF/CoA partner ESR1-e6 fusions. RNA-seq of T47D cells expressing the full panel of gene fusions demonstrated an overlapping pattern of transcriptional activation focused on estrogen response and epithelial-to-mesenchymal transition (EMT) genes driven by active fusions. This gene signature was well-preserved in a PDX naturally expressing the ESR1-YAP1 fusion. Interestingly, further study showed that a series of ET-resistant PDXs bearing a variety of ESR1 LBD point mutations induced a similar pattern to the active ESR1 fusion signature suggesting overlapping transcriptional regulatory events between ESR1 fusions and ESR1 LBD mutations. The ESR1-D538G mutation conferred the most comparable gene dysregulation to ESR1 fusions. The Y537S/N and E380Q mutations also reproduced the signature driving hormone-independent growth but with exceptions. Two PDX lines bearing either a fully heterozygous Y537S or L536P mutations were surprisingly completely estrogen-dependent. Neither of these examples exhibited the ESR1 fusion gene signature, suggesting an unknown secondary event needed to fully express the phenotype of some ESR1 mutants. The gene signature distinguished ESR1 mutations (constitutively active fusions and point mutations) from wild-type ESR1, with a 92.0% Area Under Curve. Conclusion: Here, we show that ESR1 fusions are recurrent somatic mutations that lead to drug resistance and metastasis by transcriptional reprogramming. We describe a fusion gene signature that may be useful to determine whether an ESR1 fusion or mutation is transcriptionally active and is capable of driving hormone-independent growth and endocrine therapy resistance.

Citation Format: Xuxu Gou, Meenakshi Anurag, Jonathan T Lei, Purba Singh, Sinem Seker, Adrian V Lee, Dan R Robinson, Matthew J Ellis. Recurrent active ESR1 fusions render a diagnostic transcriptional signature in metastatic breast cancer [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PS17-03.

Abstract 5677: Recurrent transcriptionally active ESR1 fusions render therapeutic vulnerabilities to kinase inhibition in advanced breast cancer

Background: We recently reported two ESR1 fusions (ESR1-YAP1 and ESR1-PCDH11X) that drive endocrine therapy (ET) resistance and metastasis in estrogen receptor positive (ER+) metastatic breast cancer (MBC) (PMC6171747). Here we report additional ESR1 fusions with diverse C-terminal partner genes - ESR1-DAB2, ESR1-GYG1, ESR1-SOX9, ESR1-ARNT2, ESR1-PCMT1 and ESR1-ARID1B. Their functional characteristics and effects on kinase biology will be described.

Methods: ESR1 fusions were identified by RNA-seq in MBC using ChimeraScan to detect fusion junction reads. ESR1 fusion cDNA constructs were expressed in ER+ breast cancer cell lines. An alamarBlue assay assessed cell proliferation. RNA-seq followed by mRNA-qPCR assessed the transcriptional properties. A scratch wound assay assessed cell motility. A Kinase Inhibitor Pulldown (KIP) mass spectrometry-based assay was conducted to examine ESR1 fusion-driven druggable kinases.

Results: All fusions retained the first 6 exons (e6) of ESR1, fused in-frame to C-terminal sequences of diverse partner genes. In addition to ESR1-YAP1 and ESR1-PCDH11X fusions, ESR1-SOX9 and ESR1-ARNT2 drove fulvestrant-resistant growth. RNA-seq revealed fusion-specific transcriptional signatures indicating enrichment of estrogen responsive and epithelial-to-mesenchymal transition (EMT) pathways that were confirmed by mRNA-qPCR. Transcriptionally active ESR1 fusions also promoted hormone-independent cell motility. KIP profiling demonstrated an increase in protein abundance of multiple receptor tyrosine kinases including RET and insulin like growth factor 1 receptor (IGF1R) in T47D cells expressing active ESR1 fusion proteins, both of which have been previously implicated in driving ET and MBC. Both proteins were also elevated in a patient-derived xenograft naturally harboring the ESR1-YAP1 fusion. Combinatorial inhibition of RET and IGF1R significantly suppressed ESR1 fusion-driven cell growth in vitro. Neither transcriptional activation nor kinase upregulation was observed in other ESR1-e6 fusions. Lastly, we report recurrent examples of specific known active ESR1 fusions: ESR1-PCDH11X (3 examples) (TEMPUS unpublished data) and ESR1-ARNT2 (2 examples) (PMC6872491).

Conclusion: A subset of ESR1-e6 fusions are active, drive ET resistance and metastasis/EMT in experimental models. When ESR1 is fused in-frame with another transcriptional regulator, activity is predictable. However, recurrent partners that have non-transcriptional roles (PCDH11X) suggests cDNA expression-based functional screens should be continued. A common pattern of kinase activation indicates that ESR1 fusion specific therapeutic strategies could be devised.

Citation Format: Xuxu Gou, Jonathan Lei, Beom-Jun Kim, Meenakshi Anurag, Sinem Seker, Saif Rehman, Adrian V. Lee, Kevin White, Michael Caldwell, Jonathan Ball, Dan R. Robinson, Matthew J. Ellis. Recurrent transcriptionally active ESR1 fusions render therapeutic vulnerabilities to kinase inhibition in advanced breast cancer [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 5677.

Long Noncoding RNA FENDRR Exhibits Antifibrotic Activity in Pulmonary Fibrosis

Abnormal activation of lung fibroblasts contributes to the initiation and progression of idiopathic pulmonary fibrosis (IPF). The objective of the present study was to investigate the role of fetal-lethal noncoding developmental regulatory RNA (FENDRR) in the activation of lung fibroblasts. Dysregulated long noncoding RNAs in IPF lungs were identified by next-generation sequencing analysis from the two online datasets. FENDRR expression in lung tissues from patients with IPF and mice with bleomycin-induced pulmonary fibrosis was determined by quantitative real-time PCR. IRP1 (iron-responsive element-binding protein 1), a protein partner of FENDRR, was identified by RNA pulldown-coupled mass spectrometric analysis and confirmed by RNA immunoprecipitation. The interaction region between FENDRR and IRP1 was determined by cross-linking immunoprecipitation. The in vivo role of FENDRR in pulmonary fibrosis was studied using adenovirus-mediated gene transfer in mice. The expression of FENDRR was downregulated in fibrotic human and mouse lungs as well as in primary lung fibroblasts isolated from bleomycin-treated mice. TGF-β1 (transforming growth factor-β1)-SMAD3 signaling inhibited FENDRR expression in lung fibroblasts. FENDRR was preferentially localized in the cytoplasm of adult lung fibroblasts and bound IRP1, suggesting its role in iron metabolism. FENDRR reduced pulmonary fibrosis by inhibiting fibroblast activation by reducing iron concentration and acting as a competing endogenous RNA of the profibrotic microRNA-214. Adenovirus-mediated FENDRR gene transfer in the mouse lung attenuated bleomycin-induced lung fibrosis and improved lung function. Our data suggest that FENDRR is an antifibrotic long noncoding RNA and a potential therapeutic target for pulmonary fibrosis.

Abstract P6-04-04: Functional characterization of ESR1 fusions in breast cancer

Background: We recently reported two in-frame estrogen receptor alpha gene (ESR1) fusion transcripts, ESR1-YAP1 and ESR1-PCDH11X, in endocrine therapy (ET)-refractory, metastatic estrogen receptor positive (ER+) breast cancer patients. Both fusions were causal in driving ET resistance and metastasis in experimental models. CDK4/6 inhibition suppressed tumor growth driven by ESR1 fusions, but did not induce tumor regression (Lei et al., Cell Rep., 2018). Additional ESR1 fusions with diverse partner genes have been recently reported (ESR1-DAB2, ESR1-GYG1, and ESR1-SOX9) (Hartmaier et al., Ann Oncol., 2018), or communicated to us (ESR1-ARNT2, ESR1-PCMT1 and ESR1-ARID1B) (Dr. Dan R. Robinson, University of Michigan). All were identified in ET-refractory tumors. We therefore undertook detailed functional characterization of these new ESR1 fusions to examine their role as disease drivers. We also used mass spectrometry-based kinomics to identify therapeutic targets induced by individual ESR1 fusions.

Methods: Unreported ESR1 fusions were identified by RNA sequencing followed by informatics designed to identify fusion gene reads. Stable ER+ breast cancer cell lines (T47D and MCF7) expressing in-frame ESR1 fusion proteins were generated by lentiviral transduction. ET-resistant cell proliferation was assessed in low estradiol medium and/or in the presence of anti-estrogen. Expression of estrogen responsive and epithelial-to-mesenchymal transition (EMT) genes were measured by mRNA-qPCR. Cell motility driven by these fusions was examined using scratch wound assay. Kinase expression patterns were examined using a Kinase Inhibitor Pulldown (KIP) mass spectrometry-based assay.

Results: All the fusions we studied contain the first 6 exons (ESR1-e6) of ESR1 that provides the first 4 coding exons. In-frame C-terminal sequences were provided by different partner genes. In addition to previously described ESR1-YAP1 and ESR1-PCDH11X fusions, ESR1-SOX9 and ESR1-ARNT2 conferred estrogen-independent and fulvestrant-resistant growth in stable ER+ breast cancer cells (T47D and MCF7). In T47D cells, these ESR1 fusions upregulated not only classical ER targets, including GREB1, PGR, and TFF1, but also EMT-associated genes, SNAI1 (Snail) and VCAN in an estradiol-independent, fulvestrant-resistant manner. In addition, transcriptionally active ESR1 fusions promoted cell motility. KIP profiling of T47D cells expressing transcriptionally active ESR1 fusions and a patient-derived xenograft naturally harboring the ESR1-YAP1 fusion revealed an increase in abundance of multiple receptor tyrosine kinases including RET and insulin like growth factor 1 receptor (IGF1R). Responses to targeting these kinases have been observed in vitro. Other ESR1-e6 fusions produced stable proteins (ESR1-DAB2, ESR1-GYG1, ESR1-PCMT1, ESR1-ARID1B) but had much less clear transcriptional phenotypes and fewer upregulated kinases identified by KIP.

Conclusion: A classification of ESR1 fusions in ET-refractory breast cancer is emerging from this study. About half of ESR1-e6 fusions studied are transcriptionally active, presumably because the 3’ partner in these examples is encoded by transcription factor gene (TF) (SOX9 or ARNT2) or co-activator gene (CoA) (YAP1). These can be classified as ESR1-e6 TF/CoA fusions. The phenotype of the remaining fusions remains under investigation because of the functionally diverse nature of the 3’ partners and the low activity in functional assays. While some of these may be passengers, the example of PCDH11X suggests a non-TF/CoA partner ESR1-e6 fusion can still be active and further studies on additional fusions will be necessary. Downstream targeting of fusion-dependent kinases may be effective.

Citation Format: Xuxu Gou, Jonathan T Lei, Sinem Seker, Beom-Jun Kim, Saif Rehman, Adrian V Lee, Dan R Robinson, Matthew J Ellis. Functional characterization of ESR1 fusions in breast cancer [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr P6-04-04.

Abstract 3479: Functional significance of ESR1 fusions with diverse gene partners in endocrine therapy resistant breast cancer

Emerging evidence suggests that ESR1 fusion transcripts that produce functional fusion proteins in estrogen receptor positive (ER+) breast cancer play a role in acquired endocrine therapy resistance. We recently reported two in-frame ESR1 fusion transcripts, ESR1-YAP1 and ESR1-PCDH11X, identified in patients with endocrine therapy resistant breast cancer. Both ESR1 fusions generated stable in-frame, fusion proteins that were transcriptionally active, driving endocrine therapy resistant proliferation and also promoting an epithelial-to-mesenchymal transition (EMT) transcriptional program leading to metastasis in experimental models. Here, we extend our studies by examining properties of additional ESR1 fusions with diverse partner genes, ESR1-DAB2, ESR1-GYG1, ESR1-SOX9, ESR1-ARNT2, ESR1-PCMT1, and ESR1-ARID1B, all identified in endocrine-refractory breast tumors, to better understand the role of ESR1fusion gene formation in endocrine therapy resistance. Structurally, ESR1-ARNT2, ESR1-PCMT1, and ESR1-ARID1B fusions retain the first 6 exons of ESR1, therefore lacking ESR1 exons encoding the ligand-binding domain (LBD) that endocrine therapies recognize, but are instead fused in-frame to C-terminal sequences from partner genes, which follows the same fusion pattern of previously identified ESR1-YAP1, ESR1-PCDH11X, ESR1-DAB2, ESR1-GYG1, and ESR1-SOX9 fusions. ESR1-ARNT2, ESR1-PCMT1, and ESR1-ARID1B constructs, along with ESR1-DAB2, ESR1-GYG1, and ESR1-SOX9 constructs were stably expressed in an ER+ breast cancer cell line, T47D, and all produced stable fusion proteins. ESR1-ARNT2 and ESR1-SOX9 promoted hormone-independent and fulvestrant-resistant cell proliferation that was sensitive to cyclin-dependent kinase (CDK) 4/6 inhibitors, palbociclib and abemaciclib. In addition, ESR1-ARNT2 induced hormone-independent activation of estrogen responsive genes (TFF1, GREB1, and PGR), and EMT genes (SNAI1 and VCAN) along with upregulation of Snail protein. These results indicate that ESR1-ARNT2 and ESR1-SOX9 have similar functional and pharmacological properties to our previously described ESR1-YAP1 and ESR1-PCDH11X fusions. These data further support a role for ESR1 fusions in driving not only endocrine therapy resistance but also activating the metastatic process. Ongoing studies are examining the ability of ESR1 fusions to induce EMT phenotypes and to determine structure function relationships of ESR1 fusions. Although the formation of ESR1 fusions likely confers resistance to all endocrine therapies that target the LBD, ESR1 fusion driven growth remained sensitive to CDK4/6 inhibitor treatment therefore providing rationale for testing ESR1 fusion detection in guiding treatment with CDK4/6 monotherapy.

Citation Format: Jonathan T. Lei, Xuxu Gou, Sinem Seker, Vaishnaivi Devorakonda, Kimberly R. Holloway, Adrian V. Lee, Dan R. Robinson, Matthew J. Ellis. Functional significance of ESR1 fusions with diverse gene partners in endocrine therapy resistant breast cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 3479.

Repeated Exposure to Streptococcus pneumoniae Exacerbates Chronic Obstructive Pulmonary Disease

Streptococcus pneumoniae is commonly found in patients with chronic obstructive pulmonary disease (COPD) and is linked to acute exacerbation of COPD. However, current clinical therapy neglects asymptomatic insidious S. pneumoniae colonization. We studied the roles of repeated exposure to S. pneumoniae in COPD progression using a mouse model. C57BL/6J mice were intranasally inoculated with S. pneumoniae ST262 every 4 weeks with or without cigarette smoke (CS) exposure up to 20 weeks to maintain persistent S. pneumoniae presence in the lower airways. Streptococcus pneumoniae enhanced CS-induced inflammatory cell infiltration at 12 to 20 weeks of exposure. Streptococcus pneumoniae also increased CS-induced release of inflammatory cytokines, including IL-1β, tumor necrosis factor-α, IL-12 (p70), and IL-5 at 20 weeks of exposure. Moreover, a combination of CS and S. pneumoniae caused alveolar epithelial injury, a decline in lung function, and an increased expression of platelet-activating factor receptor and bacterial load. Our results suggest that repeated exposure to S. pneumoniae in lower airways exacerbates CS-induced COPD.

ESR1 alterations and metastasis in estrogen receptor positive breast cancer

Endocrine therapy is essential for the treatment of patients with estrogen receptor positive (ER+) breast cancer, however, resistance and the development of metastatic disease is common. Understanding how ER+ breast cancer metastasizes is critical since the major cause of death in breast cancer is metastasis to distant organs. Results from many studies suggest dysregulation of the estrogen receptor alpha gene (ESR1 ) contributes to therapeutic resistance and metastatic biology. This review covers both pre-clinical and clinical evidence on the spectrum of ESR1 alterations including amplification, point mutations, and genomic rearrangement events driving treatment resistance and metastatic potential of ER+ breast cancer. Importantly, we describe how these ESR1 alterations may provide therapeutic opportunities to improve outcomes in patients with lethal, metastatic breast cancer.

Tea Polyphenol inhibits autophagy to sensitize Epirubicin-induced apoptosis in human bladder cancer cells

Resistance to anticancer agents such as Epirubicin (EPI) becomes a great challenge for treating bladder cancer. However, the mechanism by which chemoresistance arised is still elusive. In the present study, we showed evidence that EPI induced cytoprotective autophagy in bladder cancer cell lines T24 and BIU87. In addition, EPI robustly activated JNK-mediated phosphorylation of Bcl-2 and disruption of Bcl-2/Beclin-1 complex. Furthermore, the green tea derivative tea polyphenol (TP) inhibited EPI-induced autophagy and promoted apoptosis induced by EPI in bladder cancer cells. These results revealed a pathway for EPI-induced autophagy that involved in JNK/Bcl-2/Beclin-1 in bladder cancer cells, and that TP synergistically promoted EPI-induced apoptosis at least partly through autophagy inhibition. Thus, TP could be utilized in combination with EPI to improve EPI-based bladder cancer therapy.

Abstract P5-04-01: Functional and therapeutic significance of ESR1 fusions in metastatic ER+ breast cancer

Background. Next-generation sequencing methods have identified several ESR1 fusion genes in treatment refractory ER+ breast cancer, however detailed functional studies in experimental models are lacking and how they might be targeted remains poorly understood. We recently reported two transcriptionally active, in-frame ESR1 fusions, ESR1-YAP1 and ESR1-PCDH11X, identified in a small cohort of metastatic ER+ cases, that induce not only pan-endocrine therapy resistance but also metastatic disease progression (Lei et al., Cell Reports, in press). Limited characterization of ESR1-DAB2 and ESR1-GYG1, also identified in metastatic ER+ disease from a recent study, suggests these two ESR1 fusions also drive estrogen-independent gene activation (Hartmaier et al., Annals of Oncology, 2018). Here, we functionally characterize ESR1-DAB2 and ESR1-GYG1 along with additional ESR1 fusions discovered in metastatic ER+ breast tumors to further support a causal role for in-frame ESR1 fusions in driving endocrine therapy resistance and promoting metastasis-associated biology, and explore therapeutic vulnerabilities induced by ESR1 fusion gene formation.

Methods. RNA-seq identified ESR1 fusions from treatment refractory, ER+ metastatic breast tumors. In-frame ESR1 fusions constructs were generated and stably expressed in ER+ breast cancer cell lines: T47D, MCF7, and ZR75-1. Estrogen-independent and fulvestrant-resistant growth was monitored in hormone-deprived stable cell lines. mRNA-qPCR was performed to examine expression of estrogen responsive and epithelial-to-mesenchymal transition (EMT) genes. In vitro sensitivity to CDK4/6 inhibition was tested with palbociclib and abemaciclib.

Results. In addition to previously described ESR1-YAP1, ESR1-PCDH11X, ESR1-DAB2, and ESR1-GYG1, that follow a pattern retaining the first 6 exons of ESR1 (ESR1-e6) fused in-frame to C-terminal sequences provided by the partner gene, additional in-frame ESR1-e6 fusions, ESR1-PCMT1, ESR1-ARNT2, and ESR1-ARID1B, all identified in metastatic ER+ samples, were found to follow the same fusion pattern. ESR1-DAB2 and ESR1-GYG1 produced stable ESR1 fusion proteins in ER+ breast cancer cell lines. In T47D, these two fusions drove estrogen-independent and fulvestrant-resistant growth. In addition, T47D and ZR75-1 models revealed that ESR1-DAB2 drove estrogen-independent expression of estrogen responsive genes and also EMT genes, including SNAI1, suggesting this fusion, like ESR1-YAP1 and ESR1-PCDH11X, could also drive metastasis. Treatment with CDK4/6 inhibitors suppressed growth induced by ESR1-DAB2 and ESR1-GYG1.

Conclusion. The majority of in-frame ESR1 exon 6 fusions found in metastatic ER+ breast are transcriptionally active, drive endocrine therapy resistant proliferation, and induce an EMT-like transcriptional program. The ability to block ESR1 fusion induced growth with a CDK4/6 inhibitor is clinically significant as ESR1 fusion gene formation renders ER insensitive to all endocrine therapies that target the ligand binding domain. Furthermore, clinical diagnosis of an active ESR1 fusion could potentially stratify patients for CDK4/6 inhibitor treatment. This presentation is the most complete description of the role for ESR1 fusions in endocrine therapy resistance and metastasis described to date.

Citation Format: Lei JT, Gou X, Seker S, Haricharan S, Lee AV, Robinson DR, Ellis MJ. Functional and therapeutic significance of ESR1 fusions in metastatic ER+ breast cancer [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr P5-04-01.

[Primary mucinous adenocarcinoma of the larynx: a case report]

Mucinous adenocarcinoma is most commonly found in the digestive tract,and the prognosis is poor.We present here a case of primary mucinous adenocarcinoma of the larynx, which is an extremely rare entity and very few have been reported in the literature.A 72-year-old male complained of intermittent hoarseness for over 2 months of duration.Fiberoptic laryngoscope showed negative result for tumor.While CT scan showed laryngeal space occupying lesion,and thyroid cartilage damage was observed.The patient underwent total laryngectomy.Histopathological examination and immunohistochemistry(IHC) analysis supported the diagnosis of mucinous adenocarcinoma.

ESR1 fusions drive endocrine therapy resistance and metastasis in breast cancer

Estrogen receptor alpha gene (ESR1) fusion transcripts have been identified in breast cancer but their role in breast cancer is not completely understood. Here, we report a causal role for ESR1 fusions in driving both endocrine therapy resistance and metastasis, and describe a therapeutic strategy to target ESR1 fusion-induced growth.

EZH2 enhances the differentiation of fibroblasts into myofibroblasts in idiopathic pulmonary fibrosis

The accumulation of fibroblasts/myofibroblasts in fibrotic foci is one of the characteristics of idiopathic pulmonary fibrosis (IPF). Enhancer of zeste homolog 2 (EZH2) is the catalytic component of a multiprotein complex, polycomb repressive complex 2, which is involved in the trimethylation of histone H3 at lysine 27. In this study, we investigated the role and mechanisms of EZH2 in the differentiation of fibroblasts into myofibroblasts. We found that EZH2 was upregulated in the lungs of patients with IPF and in mice with bleomycin-induced lung fibrosis. The upregulation of EZH2 occurred in myofibroblasts. The inhibition of EZH2 by its inhibitor 3-deazaneplanocin A (DZNep) or an shRNA reduced the TGF-β1-induced differentiation of human lung fibroblasts into myofibroblasts, as demonstrated by the expression of the myofibroblast markers α-smooth muscle actin and fibronectin, and contractility. DZNep inhibited Smad2/3 nuclear translocation without affecting Smad2/3 phosphorylation. DZNep treatment attenuated bleomycin-induced pulmonary fibrosis in mice. We conclude that EZH2 induces the differentiation of fibroblasts to myofibroblasts by enhancing Smad2/3 nuclear translocation.

MicroRNA-26b Modulates the NF-κB Pathway in Alveolar Macrophages by Regulating PTEN

NF-κB is one of the best-characterized transcription factors, providing the link between early membrane-proximal signaling events and changes in many inflammatory genes. MicroRNAs are small noncoding RNAs that regulate gene expression at the posttranscriptional level. In this study, we evaluated the role of miR-26b in the LPS-induced inflammatory response in bovine alveolar macrophages (bAMs). LPS stimulation of bAMs upregulated miR-26b at 1 h and downregulated it at 6 and 36 h. Overexpression of miR-26b in bAMs enhanced the LPS-induced mRNA expression of proinflammatory cytokines and chemokines, including TNF-α, IL-1β, IL-8, and IL-10, but it directly inhibited that of IL-6. A similar trend was observed for the release of these cytokines and chemokines from bAMs. miR-26b directly bound the 3'-untranslated region of PTEN, leading to the reduction of PTEN protein in bAMs. miR-26b also enhanced the LPS-induced NF-κB signaling pathway, as revealed by increased NF-κB transcriptional activity and phosphorylation of p65, IκBα, IκB kinase, and Akt. Moreover, PTEN silencing increased the LPS-induced mRNA expression of TNF-α, IL-1β, IL-6, IL-8, and IL-10 and upregulated the NF-κB pathway. Taken together, we conclude that miR-26b participates in the inflammatory response of LPS-stimulated bAMs by modulating the NF-κB pathway through targeting PTEN.

Expression of somatostatin receptor subtype 2 in growth hormone-secreting pituitary adenoma and the regulation of miR-185

INTRODUCTION

Long-acting somatostatin analogs (SSAs) are most widely used to treat growth hormone (GH)-secreting pituitary adenoma. However, approximately 30 % of treated patients show resistance to SSAs, which may be associated with the reduction of somatostatin receptor subtype 2 (SSTR2) mRNA and protein expression.

MATERIALS AND METHODS

The present study used immunohistochemistry to detect the expression of SSTR2 and SSTR5 in twenty human GH-secreting adenoma samples treated with SSAs and seven normal pituitary samples.

RESULTS

The staining intensities of SSTR2 and SSTR5 were stronger in most adenoma samples than in normal pituitary. The expression of SSTR2 tended to be lower in the SSA non-responder group than in responders. A search of the Bioinformatics data bank and the miRCURY™ LNA array confirmed miR-185 as the putative mircoRNA (miRNA) regulating the expression of SSTR2. An in vitro study using Dual Luciferase reporter assay demonstrated that miR-185 likely targets the 3'-UTR of SSTR2 mRNA in the rat pituitary adenoma GH3 cell line. MiR-185 also downregulated or upregulated the expression of SSTR2 mRNA and SSTR2 protein, following transfection with miR-185 mimics or inhibitors, respectively.

CONCLUSION

MiR-185 enhanced the cell proliferation and inhibited the apoptosis of GH3 cells.

Regulation of myofibroblast differentiation by miR-424 during epithelial-to-mesenchymal transition

Idiopathic pulmonary fibrosis (IPF) is one of the most common and severe interstitial lung diseases. Epithelial-to-mesenchymal transition (EMT) is a process whereby epithelial cells undergo transition to a mesenchymal phenotype. This process has been shown to contribute to IPF. MicroRNAs (miRNAs) are small non-coding RNAs of 18-24 nucleotides in length which regulate gene expression. Several studies have implicated miRNAs in EMT; however, specific miRNAs that regulate EMT in IPF have not yet been identified. In this study, we identified 6 up-regulated and 3 down-regulated miRNAs in a human lung epithelial cell EMT model using miRNA microarray and real-time PCR. Overexpression of one of these up-regulated miRNAs, miR-424, increased the expression of α-smooth muscle actin, an indicator of myofibroblast differentiation, but had no effects on the epithelial or mesenchymal cell markers. miR-424 enhanced the activity of the TGF-β signaling pathway, as demonstrated by a luciferase reporter assay. Further experiments showed that miR-424 decreased the protein expression of Smurf2, a negative regulator of TGF-β signaling, indicating that miR-424 exerts a forward regulatory loop in the TGF-β signaling pathway. Our results suggest that miR-424 regulates the myofibroblast differentiation during EMT by potentiating the TGF-β signaling pathway, likely through Smurf2.

Prognostic value of 1p deletion for multiple myeloma: a meta-analysis

INTRODUCTION

Previous studies have indicated that a deletion on the short arm of chromosome 1 negatively predicts survival in patients with multiple myeloma (MM). Due to the small sample size in each study, we performed this meta-analysis to comprehensively investigate the association between the 1p deletion and survival in patients with MM.

METHODS

A literature search was conducted in both foreign and Chinese databases, including SinoMed, CNKI, Wanfang, PubMed, EMBASE, and Scopus. Hazard ratios (HR) with 95% confidence intervals (CI) for overall survival (OS) and progression-free survival (PFS) in 11 eligible articles were extracted or calculated to analyze the pooled HR, which was estimated by fixed-effect or random-effect models based on the heterogeneity between included articles. A subgroup analysis and a meta-regression were conducted, and Galbraith plots were generated to examine any possible heterogeneity.

RESULTS

The HRs for OS were available in nine articles, whereas five articles discussed HRs for PFS. The HR with 95%CI was 1.989 (95%CI 1.522-2.600, P = 0.017, I(2)  = 57.1%) when comparing the OS of patients with 1p deletion with the OS of those without this deletion. For PFS, 1p deletion still predicted a poor prognosis (HR 2.11, 95%CI 1.54-2.88, P = 0.292, I(2)  = 19.3%). Moreover, the subgroup analysis suggested that either the deleted gene on 1p or techniques for detecting chromosome abnormalities contributed to the heterogeneity, which was partially consistent with the results derived from a meta-regression analysis and the Galbraith plot method.

CONCLUSION

Our meta-analysis provides globally quantifiable confirmation of the adverse prognostic role of 1p deletion in OS and PFS for patients with MM.