Unraveling heterogeneous susceptibility and the evolution of breast cancer using a systems biology approach

This is an Open Access article distributed under the terms of the Creative Commons Attribution License.-- et al.[Background]: An essential question in cancer is why individuals with the same disease have different clinical outcomes. Progress toward a more personalized medicine in cancer patients requires taking into account the underlying heterogeneity at different molecular levels. [Results]: Here, we present a model in which there are complex interactions at different cellular and systemic levels that account for the heterogeneity of susceptibility to and evolution of ERBB2-positive breast cancers. Our model is based on our analyses of a cohort of mice that are characterized by heterogeneous susceptibility to ERBB2-positive breast cancers. Our analysis reveals that there are similarities between ERBB2 tumors in humans and those of backcross mice at clinical, genomic, expression, and signaling levels. We also show that mice that have tumors with intrinsically high levels of active AKT and ERK are more resistant to tumor metastasis. Our findings suggest for the first time that a site-specific phosphorylation at the serine 473 residue of AKT1 modifies the capacity for tumors to disseminate. Finally, we present two predictive models that can explain the heterogeneous behavior of the disease in the mouse population when we consider simultaneously certain genetic markers, liver cell signaling and serum biomarkers that are identified before the onset of the disease. [Conclusions]: Considering simultaneously tumor pathophenotypes and several molecular levels, we show the heterogeneous behavior of ERBB2-positive breast cancer in terms of disease progression. This and similar studies should help to better understand disease variability in patient populations.JPL was partially supported by FEDER and MICINN (PLE2009-119), FIS (PI07/0057, PI10/00328, PIE14/00066), the Junta de Castilla y León (SAN673/SA26/08; SAN126/SA66/09, SA078A09, CSI034U13), the “Fundación Eugenio Rodríguez Pascual”, the Fundación Inbiomed (Instituto Oncológico Obra Social de la Caja Guipozcoa-San Sebastian, Kutxa), and the “Fundación Sandra Ibarra de Solidaridad frente al Cáncer”. AC was supported by MICINN (PLE2009-119). SCLL is funded by a JAEdoc Fellowship (CSIC)/FSE. MMSF and ABG are funded by fellowships from the Junta de Castilla y Leon. WR was supported by a Forschungsstipendium of the Deutsche Forschungsgemeinschaft (DFG) [RE 3108/1-1]. TN, BPB and DYL acknowledge support from the US Department of Energy Low-Dose SFA Program at Berkeley Lab [DE-AC02-05CH11231], the National Institutes of Health [RC1NS069177] and the California Breast Cancer Research Program [15IB-0063]. JHM was supported by the National Institutes of Health, a National Cancer Institute grant (R01 CA116481), and the Low-Dose Scientific Focus Area, Office of Biological and Environmental Research, US Department of Energy (DE-AC02-05CH11231).Peer Reviewe

A genome-wide association study follow-up suggests a possible role for PPARG in systemic sclerosis susceptibility

Introduction: A recent genome-wide association study (GWAS) comprising a French cohort of systemic sclerosis (SSc) reported several non-HLA single-nucleotide polymorphisms (SNPs) showing a nominal association in the discovery phase. We aimed to identify previously overlooked susceptibility variants by using a follow-up strategy.<p></p> Methods: Sixty-six non-HLA SNPs showing a P value <10-4 in the discovery phase of the French SSc GWAS were analyzed in the first step of this study, performing a meta-analysis that combined data from the two published SSc GWASs. A total of 2,921 SSc patients and 6,963 healthy controls were included in this first phase. Two SNPs, PPARG rs310746 and CHRNA9 rs6832151, were selected for genotyping in the replication cohort (1,068 SSc patients and 6,762 healthy controls) based on the results of the first step. Genotyping was performed by using TaqMan SNP genotyping assays. Results: We observed nominal associations for both PPARG rs310746 (PMH = 1.90 × 10-6, OR, 1.28) and CHRNA9 rs6832151 (PMH = 4.30 × 10-6, OR, 1.17) genetic variants with SSc in the first step of our study. In the replication phase, we observed a trend of association for PPARG rs310746 (P value = 0.066; OR, 1.17). The combined overall Mantel-Haenszel meta-analysis of all the cohorts included in the present study revealed that PPARG rs310746 remained associated with SSc with a nominal non-genome-wide significant P value (PMH = 5.00 × 10-7; OR, 1.25). No evidence of association was observed for CHRNA9 rs6832151 either in the replication phase or in the overall pooled analysis.<p></p> Conclusion: Our results suggest a role of PPARG gene in the development of SSc

Machine Learning Improves Risk Stratification in Myelofibrosis: An Analysis of the Spanish Registry of Myelofibrosis

Myelofibrosis (MF) is a myeloproliferative neoplasm (MPN) with heterogeneous clinical course. Allogeneic hematopoietic cell transplantation remains the only curative therapy, but its morbidity and mortality require careful candidate selection. Therefore, accurate disease risk prognostication is critical for treatment decision-making. We obtained registry data from patients diagnosed with MF in 60 Spanish institutions (N = 1386). These were randomly divided into a training set (80%) and a test set (20%). A machine learning (ML) technique (random forest) was used to model overall survival (OS) and leukemia-free survival (LFS) in the training set, and the results were validated in the test set. We derived the AIPSS-MF (Artificial Intelligence Prognostic Scoring System for Myelofibrosis) model, which was based on 8 clinical variables at diagnosis and achieved high accuracy in predicting OS (training set c-index, 0.750; test set c-index, 0.744) and LFS (training set c-index, 0.697; test set c-index, 0.703). No improvement was obtained with the inclusion of MPN driver mutations in the model. We were unable to adequately assess the potential benefit of including adverse cytogenetics or high-risk mutations due to the lack of these data in many patients. AIPSS-MF was superior to the IPSS regardless of MF subtype and age range and outperformed the MYSEC-PM in patients with secondary MF. In conclusion, we have developed a prediction model based exclusively on clinical variables that provides individualized prognostic estimates in patients with primary and secondary MF. The use of AIPSS-MF in combination with predictive models that incorporate genetic information may improve disease risk stratification

Effect of viral storm in patients admitted to intensive care units with severe COVID-19 in Spain: a multicentre, prospective, cohort study

Background: The contribution of the virus to the pathogenesis of severe COVID-19 is still unclear. We aimed to evaluate associations between viral RNA load in plasma and host response, complications, and deaths in critically ill patients with COVID-19. Methods: We did a prospective cohort study across 23 hospitals in Spain. We included patients aged 18 years or older with laboratory-confirmed SARS-CoV-2 infection who were admitted to an intensive care unit between March 16, 2020, and Feb 27, 2021. RNA of the SARS-CoV-2 nucleocapsid region 1 (N1) was quantified in plasma samples collected from patients in the first 48 h following admission, using digital PCR. Patients were grouped on the basis of N1 quantity: VIR-N1-Zero ([removed]2747 N1 copies per mL). The primary outcome was all-cause death within 90 days after admission. We evaluated odds ratios (ORs) for the primary outcome between groups using a logistic regression analysis. Findings: 1068 patients met the inclusion criteria, of whom 117 had insufficient plasma samples and 115 had key information missing. 836 patients were included in the analysis, of whom 403 (48%) were in the VIR-N1-Low group, 283 (34%) were in the VIR-N1-Storm group, and 150 (18%) were in the VIR-N1-Zero group. Overall, patients in the VIR-N1-Storm group had the most severe disease: 266 (94%) of 283 patients received invasive mechanical ventilation (IMV), 116 (41%) developed acute kidney injury, 180 (65%) had secondary infections, and 148 (52%) died within 90 days. Patients in the VIR-N1-Zero group had the least severe disease: 81 (54%) of 150 received IMV, 34 (23%) developed acute kidney injury, 47 (32%) had secondary infections, and 26 (17%) died within 90 days (OR for death 0·30, 95% CI 0·16–0·55; p<0·0001, compared with the VIR-N1-Storm group). 106 (26%) of 403 patients in the VIR-N1-Low group died within 90 days (OR for death 0·39, 95% CI 0·26–0·57; p[removed]11 página

Comprensión conceptual de estadística y métodos de decisión utilizando herramientas de gamificación para la adquisición de competencias en los grados de ciencias económicas y empresariales

El objetivo general del proyecto es mejorar la motivación y la actitud de los alumnos de los grados de ciencias sociales hacia las materias cuantitativas, así como su participación activa en el proceso de enseñanza-aprendizaje, individual y en grupo, con la finalidad última de incrementar su nivel de competencias y rendimiento

A blood microRNA classifier for the prediction of ICU mortality in COVID-19 patients: a multicenter validation study

Background: The identification of critically ill COVID-19 patients at risk of fatal outcomes remains a challenge. Here, we first validated candidate microRNAs (miRNAs) as biomarkers for clinical decision-making in critically ill patients. Second, we constructed a blood miRNA classifier for the early prediction of adverse outcomes in the ICU. Methods: This was a multicenter, observational and retrospective/prospective study including 503 critically ill patients admitted to the ICU from 19 hospitals. qPCR assays were performed in plasma samples collected within the first 48 h upon admission. A 16-miRNA panel was designed based on recently published data from our group. Results: Nine miRNAs were validated as biomarkers of all-cause in-ICU mortality in the independent cohort of critically ill patients (FDR < 0.05). Cox regression analysis revealed that low expression levels of eight miRNAs were associated with a higher risk of death (HR from 1.56 to 2.61). LASSO regression for variable selection was used to construct a miRNA classifier. A 4-blood miRNA signature composed of miR-16-5p, miR-192-5p, miR-323a-3p and miR-451a predicts the risk of all-cause in-ICU mortality (HR 2.5). Kaplan‒Meier analysis confirmed these findings. The miRNA signature provides a significant increase in the prognostic capacity of conventional scores, APACHE-II (C-index 0.71, DeLong test p-value 0.055) and SOFA (C-index 0.67, DeLong test p-value 0.001), and a risk model based on clinical predictors (C-index 0.74, DeLong test-p-value 0.035). For 28-day and 90-day mortality, the classifier also improved the prognostic value of APACHE-II, SOFA and the clinical model. The association between the classifier and mortality persisted even after multivariable adjustment. The functional analysis reported biological pathways involved in SARS-CoV infection and inflammatory, fibrotic and transcriptional pathways. Conclusions: A blood miRNA classifier improves the early prediction of fatal outcomes in critically ill COVID-19 patients.11 página

Prognostic implications of comorbidity patterns in critically ill COVID-19 patients: A multicenter, observational study

Background The clinical heterogeneity of COVID-19 suggests the existence of different phenotypes with prognostic implications. We aimed to analyze comorbidity patterns in critically ill COVID-19 patients and assess their impact on in-hospital outcomes, response to treatment and sequelae. Methods Multicenter prospective/retrospective observational study in intensive care units of 55 Spanish hospitals. 5866 PCR-confirmed COVID-19 patients had comorbidities recorded at hospital admission; clinical and biological parameters, in-hospital procedures and complications throughout the stay; and, clinical complications, persistent symptoms and sequelae at 3 and 6 months. Findings Latent class analysis identified 3 phenotypes using training and test subcohorts: low-morbidity (n=3385; 58%), younger and with few comorbidities; high-morbidity (n=2074; 35%), with high comorbid burden; and renal-morbidity (n=407; 7%), with chronic kidney disease (CKD), high comorbidity burden and the worst oxygenation profile. Renal-morbidity and high-morbidity had more in-hospital complications and higher mortality risk than low-morbidity (adjusted HR (95% CI): 1.57 (1.34-1.84) and 1.16 (1.05-1.28), respectively). Corticosteroids, but not tocilizumab, were associated with lower mortality risk (HR (95% CI) 0.76 (0.63-0.93)), especially in renal-morbidity and high-morbidity. Renal-morbidity and high-morbidity showed the worst lung function throughout the follow-up, with renal-morbidity having the highest risk of infectious complications (6%), emergency visits (29%) or hospital readmissions (14%) at 6 months (p<0.01). Interpretation Comorbidity-based phenotypes were identified and associated with different expression of in-hospital complications, mortality, treatment response, and sequelae, with CKD playing a major role. This could help clinicians in day-to-day decision making including the management of post-discharge COVID-19 sequelae. Copyright (C) 2022 The Author(s). Published by Elsevier Ltd

Cross-disease Meta-analysis of Genome-wide Association Studies for Systemic Sclerosis and Rheumatoid Arthritis Reveals IRF4 as a New Common Susceptibility Locus

Objectives: Systemic sclerosis (SSc) and rheumatoid arthritis (RA) are autoimmune diseases that share clinical and immunological characteristics. To date, several shared SSc- RA loci have been identified independently. In this study, we aimed to systematically search for new common SSc-RA loci through an inter-disease meta-GWAS strategy. Methods: We performed a meta-analysis combining GWAS datasets of SSc and RA using a strategy that allowed identification of loci with both same-direction and opposingdirection allelic effects. The top single-nucleotide polymorphisms (SNPs) were followed-up in independent SSc and RA case-control cohorts. This allowed us to increase the sample size to a total of 8,830 SSc patients, 16,870 RA patients and 43,393 controls. Results: The cross-disease meta-analysis of the GWAS datasets identified several loci with nominal association signals (P-value < 5 x 10-6), which also showed evidence of association in the disease-specific GWAS scan. These loci included several genomic regions not previously reported as shared loci, besides risk factors associated with both diseases in previous studies. The follow-up of the putatively new SSc-RA loci identified IRF4 as a shared risk factor for these two diseases (Pcombined = 3.29 x 10-12). In addition, the analysis of the biological relevance of the known SSc-RA shared loci pointed to the type I interferon and the interleukin 12 signaling pathways as the main common etiopathogenic factors. Conclusions: Our study has identified a novel shared locus, IRF4, for SSc and RA and highlighted the usefulness of cross-disease GWAS meta-analysis in the identification of common risk loci

Effects of intubation timing in patients with COVID-19 throughout the four waves of the pandemic : a matched analysis

The primary aim of our study was to investigate the association between intubation timing and hospital mortality in critically ill patients with COVID-19-associated respiratory failure. We also analysed both the impact of such timing throughout the first four pandemic waves and the influence of prior non-invasive respiratory support on outcomes. This is a secondary analysis of a multicentre, observational and prospective cohort study that included all consecutive patients undergoing invasive mechanical ventilation due to COVID-19 from across 58 Spanish intensive care units (ICU) participating in the CIBERESUCICOVID project. The study period was between 29 February 2020 and 31 August 2021. Early intubation was defined as that occurring within the first 24 h of intensive care unit (ICU) admission. Propensity score (PS) matching was used to achieve balance across baseline variables between the early intubation cohort and those patients who were intubated after the first 24 h of ICU admission. Differences in outcomes between early and delayed intubation were also assessed. We performed sensitivity analyses to consider a different timepoint (48 h from ICU admission) for early and delayed intubation. Of the 2725 patients who received invasive mechanical ventilation, a total of 614 matched patients were included in the analysis (307 for each group). In the unmatched population, there were no differences in mortality between the early and delayed groups. After PS matching, patients with delayed intubation presented higher hospital mortality (27.3% versus 37.1%, p =0.01), ICU mortality (25.7% versus 36.1%, p=0.007) and 90-day mortality (30.9% versus 40.2%, p=0.02) when compared to the early intubation group. Very similar findings were observed when we used a 48-hour timepoint for early or delayed intubation. The use of early intubation decreased after the first wave of the pandemic (72%, 49%, 46% and 45% in the first, second, third and fourth wave, respectively; first versus second, third and fourth waves p<0.001). In both the main and sensitivity analyses, hospital mortality was lower in patients receiving high-flow nasal cannula (n=294) who were intubated earlier. The subgroup of patients undergoing NIV (n=214) before intubation showed higher mortality when delayed intubation was set as that occurring after 48 h from ICU admission, but not when after 24 h. In patients with COVID-19 requiring invasive mechanical ventilation, delayed intubation was associated with a higher risk of hospital mortality. The use of early intubation significantly decreased throughout the course of the pandemic. Benefits of such an approach occurred more notably in patients who had received high-flow nasal cannul

Intermediate Molecular Phenotypes to Identify Genetic Markers of Anthracycline-Induced Cardiotoxicity Risk.

Cardiotoxicity due to anthracyclines (CDA) affects cancer patients, but we cannot predict who may suffer from this complication. CDA is a complex trait with a polygenic component that is mainly unidentified. We propose that levels of intermediate molecular phenotypes (IMPs) in the myocardium associated with histopathological damage could explain CDA susceptibility, so variants of genes encoding these IMPs could identify patients susceptible to this complication. Thus, a genetically heterogeneous cohort of mice (n = 165) generated by backcrossing were treated with doxorubicin and docetaxel. We quantified heart fibrosis using an Ariol slide scanner and intramyocardial levels of IMPs using multiplex bead arrays and QPCR. We identified quantitative trait loci linked to IMPs (ipQTLs) and cdaQTLs via linkage analysis. In three cancer patient cohorts, CDA was quantified using echocardiography or Cardiac Magnetic Resonance. CDA behaves as a complex trait in the mouse cohort. IMP levels in the myocardium were associated with CDA. ipQTLs integrated into genetic models with cdaQTLs account for more CDA phenotypic variation than that explained by cda-QTLs alone. Allelic forms of genes encoding IMPs associated with CDA in mice, including AKT1, MAPK14, MAPK8, STAT3, CAS3, and TP53, are genetic determinants of CDA in patients. Two genetic risk scores for pediatric patients (n = 71) and women with breast cancer (n = 420) were generated using machine-learning Least Absolute Shrinkage and Selection Operator (LASSO) regression. Thus, IMPs associated with heart damage identify genetic markers of CDA risk, thereby allowing more personalized patient management.J.P.L.’s lab is sponsored by Grant PID2020-118527RB-I00 funded by MCIN/AEI/10.13039/ 501100011039; Grant PDC2021-121735-I00 funded by MCIN/AEI/10.13039/501100011039 and by the “European Union Next Generation EU/PRTR”, the Regional Government of Castile and León (CSI144P20). J.P.L. and P.L.S. are supported by the Carlos III Health Institute (PIE14/00066). AGN laboratory and human patients’ studies are supported by an ISCIII project grant (PI18/01242). The Human Genotyping unit is a member of CeGen, PRB3, and is supported by grant PT17/0019 of the PE I + D + i 2013–2016, funded by ISCIII and ERDF. SCLl is supported by MINECO/FEDER research grants (RTI2018-094130-B-100). CH was supported by the Department of Defense (DoD) BCRP, No. BC190820; and the National Cancer Institute (NCI) at the National Institutes of Health (NIH), No. R01CA184476. Lawrence Berkeley National Laboratory (LBNL) is a multi-program national laboratory operated by the University of California for the DOE under contract DE AC02-05CH11231. The Proteomics Unit belongs to ProteoRed, PRB3-ISCIII, supported by grant PT17/0019/0023 of the PE I + D +i, 2017–2020, funded by ISCIII and FEDER. RCC is funded by fellowships from the Spanish Regional Government of Castile and León. NGS is a recipient of an FPU fellowship (MINECO/FEDER). hiPSC-CM studies were funded in part by the “la Caixa” Banking Foundation under the project code HR18-00304 and a Severo Ochoa CNIC Intramural Project (Exp. 12-2016 IGP) to J.J.S