RESUMEN
Cytometry enables precise single-cell phenotyping within heterogeneous populations. These cell types are traditionally annotated via manual gating, but this method lacks reproducibility and sensitivity to batch effect. Also, the most recent cytometers-spectral flow or mass cytometers-create rich and high-dimensional data whose analysis via manual gating becomes challenging and time-consuming. To tackle these limitations, we introduce Scyan https://github.com/MICS-Lab/scyan, a Single-cell Cytometry Annotation Network that automatically annotates cell types using only prior expert knowledge about the cytometry panel. For this, it uses a normalizing flow-a type of deep generative model-that maps protein expressions into a biologically relevant latent space. We demonstrate that Scyan significantly outperforms the related state-of-the-art models on multiple public datasets while being faster and interpretable. In addition, Scyan overcomes several complementary tasks, such as batch-effect correction, debarcoding and population discovery. Overall, this model accelerates and eases cell population characterization, quantification and discovery in cytometry.
Asunto(s)
Biología , Reproducibilidad de los Resultados , Citometría de Flujo/métodosRESUMEN
Single-cell RNA sequencing (scRNA-seq) technology produces an unprecedented resolution at the level of a unique cell, raising great hopes in medicine. Nevertheless, scRNA-seq data suffer from high variations due to the experimental conditions, called batch effects, preventing any aggregated downstream analysis. Adversarial Information Factorization provides a robust batch-effect correction method that does not rely on prior knowledge of the cell types nor a specific normalization strategy while being adapted to any downstream analysis task. It compares to and even outperforms state-of-the-art methods in several scenarios: low signal-to-noise ratio, batch-specific cell types with few cells, and a multi-batches dataset with imbalanced batches and batch-specific cell types. Moreover, it best preserves the relative gene expression between cell types, yielding superior differential expression analysis results. Finally, in a more complex setting of a Leukemia cohort, our method preserved most of the underlying biological information for each patient while aligning the batches, improving the clustering metrics in the aggregated dataset.
Asunto(s)
Análisis de la Célula Individual , Análisis de Expresión Génica de una Sola Célula , Humanos , Análisis de la Célula Individual/métodos , Análisis por Conglomerados , Secuenciación del Exoma , Benchmarking , Análisis de Secuencia de ARN/métodos , Perfilación de la Expresión Génica , AlgoritmosRESUMEN
The developmental history of blood cancer begins with mutation acquisition and the resulting malignant clone expansion. The two most prevalent driver mutations found in myeloproliferative neoplasms-JAK2V617F and CALRm-occur in hematopoietic stem cells, which are highly complex to observe in vivo. To circumvent this difficulty, we propose a method relying on mathematical modeling and statistical inference to determine disease initiation and dynamics. Our findings suggest that CALRm mutations tend to occur later in life than JAK2V617F. Our results confirm the higher proliferative advantage of the CALRm malignant clone compared to JAK2V617F. Furthermore, we illustrate how mathematical modeling and Bayesian inference can be used for setting up early screening strategies.
Asunto(s)
Calreticulina , Janus Quinasa 2 , Trastornos Mieloproliferativos , Teorema de Bayes , Calreticulina/genética , Humanos , Janus Quinasa 2/genética , Modelos Biológicos , Mutación , Trastornos Mieloproliferativos/genéticaRESUMEN
The availability of patient cohorts with several types of omics data opens new perspectives for exploring the disease's underlying biological processes and developing predictive models. It also comes with new challenges in computational biology in terms of integrating high-dimensional and heterogeneous data in a fashion that captures the interrelationships between multiple genes and their functions. Deep learning methods offer promising perspectives for integrating multi-omics data. In this paper, we review the existing integration strategies based on autoencoders and propose a new customizable one whose principle relies on a two-phase approach. In the first phase, we adapt the training to each data source independently before learning cross-modality interactions in the second phase. By taking into account each source's singularity, we show that this approach succeeds at taking advantage of all the sources more efficiently than other strategies. Moreover, by adapting our architecture to the computation of Shapley additive explanations, our model can provide interpretable results in a multi-source setting. Using multiple omics sources from different TCGA cohorts, we demonstrate the performance of the proposed method for cancer on test cases for several tasks, such as the classification of tumor types and breast cancer subtypes, as well as survival outcome prediction. We show through our experiments the great performances of our architecture on seven different datasets with various sizes and provide some interpretations of the results obtained. Our code is available on (https://github.com/HakimBenkirane/CustOmics).
Asunto(s)
Neoplasias de la Mama , Aprendizaje Profundo , Femenino , Humanos , Neoplasias de la Mama/genética , Biología Computacional/métodos , MultiómicaRESUMEN
Myeloproliferative neoplasms (MPNs) are hematologic malignancies that result from acquired driver mutations in hematopoietic stem cells (HSCs), causing overproduction of blood cells and an increased risk of thrombohemorrhagic events. The most common MPN driver mutation affects the JAK2 gene (JAK2V617F ). Interferon alpha (IFNα) is a promising treatment against MPNs by inducing a hematologic response and molecular remission for some patients. Mathematical models have been proposed to describe how IFNα targets mutated HSCs, indicating that a minimal dose is necessary for long-term remission. This study aims to determine a personalized treatment strategy. First, we show the capacity of an existing model to predict cell dynamics for new patients from data that can be easily obtained in clinic. Then, we study different treatment scenarios in silico for three patients, considering potential IFNα dose-toxicity relations. We assess when the treatment should be interrupted depending on the response, the patient's age, and the inferred development of the malignant clone without IFNα We find that an optimal strategy would be to treat patients with a constant dose so that treatment could be interrupted as quickly as possible. Higher doses result in earlier discontinuation but also higher toxicity. Without knowledge of the dose-toxicity relationship, trade-off strategies can be found for each patient. A compromise strategy is to treat patients with medium doses (60-120 µg/week) for 10-15 years. Altogether, this work demonstrates how a mathematical model calibrated from real data can help build a clinical decision-support tool to optimize long-term IFNα therapy for MPN patients. SIGNIFICANCE STATEMENT: Myeloproliferative neoplasms (MPNs) are chronic blood cancers. Interferon alpha (IFNα) is a promising treatment with the potential to induce a molecular response by targeting mutated hematopoietic stem cells. MPN patients are treated over several years, and there is a lack of knowledge concerning the posology strategy and the best timing for interrupting therapy. The study opens avenues for rationalizing how to treat MPN patients with IFNα over several years, promoting a more personalized approach to treatment.
Asunto(s)
Trastornos Mieloproliferativos , Neoplasias , Humanos , Trastornos Mieloproliferativos/tratamiento farmacológico , Trastornos Mieloproliferativos/genética , Trastornos Mieloproliferativos/patología , Interferón-alfa/uso terapéutico , Interferón-alfa/farmacología , Células Madre Hematopoyéticas , Inmunoterapia , Neoplasias/patología , Janus Quinasa 2/genética , MutaciónRESUMEN
OBJECTIVES: Around 30% of patients with RA have an inadequate response to MTX. We aimed to use routine clinical and biological data to build machine learning models predicting EULAR inadequate response to MTX and to identify simple predictive biomarkers. METHODS: Models were trained on RA patients fulfilling the 2010 ACR/EULAR criteria from the ESPOIR and Leiden EAC cohorts to predict the EULAR response at 9 months (± 6 months). Several models were compared on the training set using the AUROC. The best model was evaluated on an external validation cohort (tREACH). The model's predictions were explained using Shapley values to extract a biomarker of inadequate response. RESULTS: We included 493 therapeutic sequences from ESPOIR, 239 from EAC and 138 from tREACH. The model selected DAS28, Lymphocytes, Creatininemia, Leucocytes, AST, ALT, swollen joint count and corticosteroid co-treatment as predictors. The model reached an AUROC of 0.72 [95% CI (0.63, 0.80)] on the external validation set, where 70% of patients were responders to MTX. Patients predicted as inadequate responders had only 38% [95% CI (20%, 58%)] chance to respond and using the algorithm to decide to initiate MTX would decrease inadequate-response rate from 30% to 23% [95% CI: (17%, 29%)]. A biomarker was identified in patients with moderate or high activity (DAS28 > 3.2): patients with a lymphocyte count superior to 2000 cells/mm3 are significantly less likely to respond. CONCLUSION: Our study highlights the usefulness of machine learning in unveiling subgroups of inadequate responders to MTX to guide new therapeutic strategies. Further work is needed to validate this approach.
Asunto(s)
Antirreumáticos , Artritis Reumatoide , Humanos , Metotrexato/uso terapéutico , Antirreumáticos/uso terapéutico , Resultado del Tratamiento , Artritis Reumatoide/tratamiento farmacológico , Biomarcadores , Quimioterapia CombinadaRESUMEN
Classical BCR-ABL-negative myeloproliferative neoplasms (MPNs) are clonal disorders of hematopoietic stem cells (HSCs) caused mainly by recurrent mutations in genes encoding JAK2 (JAK2), calreticulin (CALR), or the thrombopoietin receptor (MPL). Interferon α (IFNα) has demonstrated some efficacy in inducing molecular remission in MPNs. To determine factors that influence molecular response rate, we evaluated the long-term molecular efficacy of IFNα in patients with MPN by monitoring the fate of cells carrying driver mutations in a prospective observational and longitudinal study of 48 patients over more than 5 years. We measured the clonal architecture of early and late hematopoietic progenitors (84 845 measurements) and the global variant allele frequency in mature cells (409 measurements) several times per year. Using mathematical modeling and hierarchical Bayesian inference, we further inferred the dynamics of IFNα-targeted mutated HSCs. Our data support the hypothesis that IFNα targets JAK2V617F HSCs by inducing their exit from quiescence and differentiation into progenitors. Our observations indicate that treatment efficacy is higher in homozygous than heterozygous JAK2V617F HSCs and increases with high IFNα dose in heterozygous JAK2V617F HSCs. We also found that the molecular responses of CALRm HSCs to IFNα were heterogeneous, varying between type 1 and type 2 CALRm, and a high dose of IFNα correlates with worse outcomes. Our work indicates that the long-term molecular efficacy of IFNα implies an HSC exhaustion mechanism and depends on both the driver mutation type and IFNα dose.
Asunto(s)
Células Madre Hematopoyéticas/efectos de los fármacos , Factores Inmunológicos/uso terapéutico , Interferón-alfa/uso terapéutico , Mutación/efectos de los fármacos , Trastornos Mieloproliferativos/tratamiento farmacológico , Calreticulina/genética , Células Madre Hematopoyéticas/metabolismo , Células Madre Hematopoyéticas/patología , Humanos , Factores Inmunológicos/farmacología , Interferón-alfa/farmacología , Janus Quinasa 2/genética , Estudios Longitudinales , Trastornos Mieloproliferativos/genética , Trastornos Mieloproliferativos/patología , Estudios Prospectivos , Receptores de Trombopoyetina/genética , Células Tumorales CultivadasRESUMEN
Standard imaging cannot reliably predict the nature of renal tumors. Among malignant renal tumors, clear cell renal cell carcinoma (ccRCC) is the most common histological subtype, in which the vascular endothelial growth factor 2 (VEGFR-2) is highly expressed in the vascular endothelium. BR55, a contrast agent for ultrasound imaging, consists of gas-core lipid microbubbles that specifically target and bind to the extracellular portion of the VEGFR-2. The specific information provided by ultrasound molecular imaging (USMI) using BR55 was compared with the vascular tumor expression of the VEGFR-2 by immunohistochemical (IHC) staining in a preclinical model of ccRCC. Patients' ccRCCs were orthotopically grafted onto Nod-Scid-Gamma (NSG) mice to generate patient-derived xenografts (PdX). Mice were divided into four groups to receive either vehicle or axitinib an amount of 2, 7.5 or 15 mg/kg twice daily. Perfusion parameters and the BR55 ultrasound contrast signal on PdX renal tumors were analyzed at D0, D1, D3, D7 and D11, and compared with IHC staining for the VEGFR-2 and CD34. Significant Pearson correlation coefficients were observed between the area under the curve (AUC) and the CD34 (0.84, p < 10-4), and between the VEGFR-2-specific signal obtained by USMI and IHC (0.72, p < 10-4). USMI with BR55 could provide instant, quantitative information on tumor VEGFR-2 expression to characterize renal masses non-invasively.
Asunto(s)
Carcinoma de Células Renales , Neoplasias Renales , Humanos , Ratones , Animales , Carcinoma de Células Renales/diagnóstico por imagen , Carcinoma de Células Renales/genética , Receptor 2 de Factores de Crecimiento Endotelial Vascular/metabolismo , Factor A de Crecimiento Endotelial Vascular , Xenoinjertos , Ultrasonografía/métodos , Imagen Molecular/métodos , Medios de Contraste , Neoplasias Renales/diagnóstico por imagenRESUMEN
Clear cell renal cell carcinoma (ccRCC) constitutes the most common renal cell carcinoma subtype and has long been recognized as an immunogenic cancer. As such, significant attention has been directed toward optimizing immune-checkpoints (IC)-based therapies. Despite proven benefits, a substantial number of patients remain unresponsive to treatment, suggesting that yet unreported, immunosuppressive mechanisms coexist within tumors and their microenvironment. Here, we comprehensively analyzed and ranked forty-four immune-checkpoints expressed in ccRCC on the basis of in-depth analysis of RNAseq data collected fromâ¯the TCGA database and advanced statistical methods designed to obtain the group of checkpoints that best discriminates tumor from healthy tissues. Immunohistochemistry and flow cytometry confirmed and enlarged the bioinformatics results. In particular, by using the recursive feature elimination method, we show that HLA-G, B7H3, PDL-1 and ILT2 are the most relevant genes that characterize ccRCC. Notably, ILT2 expression was detected for the first time on tumor cells. The levels of other ligand-receptor pairs such as CD70:CD27; 4-1BB:4-1BBL; CD40:CD40L; CD86:CTLA4; MHC-II:Lag3; CD200:CD200R; CD244:CD48 were also found highly expressed in tumors compared to adjacent non-tumor tissues. Collectively, our approach provides a comprehensible classification of forty-four IC expressed in ccRCC, some of which were never reported before to be co-expressed in ccRCC. In addition, the algorithms used allowed identifying the most relevant group that best discriminates tumor from healthy tissues. The data can potentially assist on the choice of valuable immune-therapy targets which hold potential for the development of more effective anti-tumor treatments.
Asunto(s)
Antígenos CD/inmunología , Biomarcadores de Tumor/inmunología , Carcinoma de Células Renales/inmunología , Antígenos HLA-G/inmunología , Neoplasias Renales/inmunología , Receptor Leucocitario Tipo Inmunoglobulina B1/inmunología , Glicoproteínas de Membrana/inmunología , Receptores Inmunológicos/inmunología , Adulto , Anciano , Anciano de 80 o más Años , Carcinoma de Células Renales/patología , Estudios de Casos y Controles , Femenino , Estudios de Seguimiento , Perfilación de la Expresión Génica , Humanos , Neoplasias Renales/patología , Masculino , Persona de Mediana Edad , Pronóstico , Estudios RetrospectivosRESUMEN
OBJECTIVES: The aim of our study was to investigate the association between driver oncogene alterations and metastatic patterns on imaging assessment, in a large cohort of metastatic lung adenocarcinoma patients. METHODS: From January 2010 to May 2017, 550 patients with stage IV lung adenocarcinoma with molecular analysis were studied retrospectively including 135 EGFR-mutated, 81 ALK-rearrangement, 47 BRAF-mutated, 141 KRAS-mutated, and 146 negative tumors for these 4 mutations (4N). After review of the complete imaging report by two radiologists (junior and senior) to identify metastatic sites, univariate correlation analyzes were performed. RESULTS: We found differences in metastatic tropism depending on the molecular alteration type when compared with the non-mutated 4N group: in the EGFR group, pleural metastases were more frequent (32% versus 20%; p = 0.021), and adrenal and node metastases less common (6% versus 23%; p < 0.001 and 11% versus 23%; p = 0.011). In the ALK group, there were more brain and lung metastases (respectively 42% versus 29%; p = 0.043 and 37% versus 24%; p = 0.037). In the BRAF group, pleural and pericardial metastases were more common (respectively 47% versus 20%; p < 0.001 and 11% versus 3%; p = 0.04) and bone metastases were rarer (21% versus 42%; p = 0.011). Lymphangitis was more frequent in EGFR, ALK, and BRAF groups (respectively 6%, 7%, and 15% versus 1%); p = 0.016; p = 0.009; and p < 0.001. CONCLUSION: The application of these correlations between molecular status and metastatic tropism in clinical practice may lead to earlier and more accurate identification of patients for targeted therapy. KEY POINTS: ⢠Bone and brain metastasis are the most common organs involved in lung adenocarcinoma but the relative incidence of each metastatic site depends on the molecular alteration. ⢠EGFR-mutated tumors preferentially spread to the pleura and less commonly to adrenals, ALK-rearrangement tumors usually spread to the brain and the lungs, whereas BRAF-mutated tumors are unlikely to spread to bones and have a serous (pericardial ad pleural) tropism. ⢠These correlations could help in the clinical management of patients with metastatic lung adenocarcinoma.
Asunto(s)
Neoplasias Óseas/secundario , Carcinoma de Pulmón de Células no Pequeñas/diagnóstico , ADN de Neoplasias/genética , Receptores ErbB/genética , Neoplasias Pulmonares/diagnóstico , Mutación , Estadificación de Neoplasias , Adulto , Anciano , Biomarcadores de Tumor/genética , Biomarcadores de Tumor/metabolismo , Neoplasias Óseas/diagnóstico , Neoplasias Óseas/genética , Carcinoma de Pulmón de Células no Pequeñas/genética , Carcinoma de Pulmón de Células no Pequeñas/secundario , Receptores ErbB/metabolismo , Femenino , Humanos , Neoplasias Pulmonares/genética , Neoplasias Pulmonares/metabolismo , Masculino , Persona de Mediana Edad , Estudios RetrospectivosRESUMEN
Background and Aims: Functional-structural plant models (FSPMs) describe explicitly the interactions between plants and their environment at organ to plant scale. However, the high level of description of the structure or model mechanisms makes this type of model very complex and hard to calibrate. A two-step methodology to facilitate the calibration process is proposed here. Methods: First, a global sensitivity analysis method was applied to the calibration loss function. It provided first-order and total-order sensitivity indexes that allow parameters to be ranked by importance in order to select the most influential ones. Second, the Akaike information criterion (AIC) was used to quantify the model's quality of fit after calibration with different combinations of selected parameters. The model with the lowest AIC gives the best combination of parameters to select. This methodology was validated by calibrating the model on an independent data set (same cultivar, another year) with the parameters selected in the second step. All the parameters were set to their nominal value; only the most influential ones were re-estimated. Key Results: Sensitivity analysis applied to the calibration loss function is a relevant method to underline the most significant parameters in the estimation process. For the studied winter oilseed rape model, 11 out of 26 estimated parameters were selected. Then, the model could be recalibrated for a different data set by re-estimating only three parameters selected with the model selection method. Conclusions: Fitting only a small number of parameters dramatically increases the efficiency of recalibration, increases the robustness of the model and helps identify the principal sources of variation in varying environmental conditions. This innovative method still needs to be more widely validated but already gives interesting avenues to improve the calibration of FSPMs.
Asunto(s)
Interacción Gen-Ambiente , Modelos Estadísticos , Fenómenos Fisiológicos de las Plantas , Plantas/genética , CalibraciónRESUMEN
MOTIVATION: Likelihood-free methods, like Approximate Bayesian Computation (ABC), have been extensively used in model-based statistical inference with intractable likelihood functions. When combined with Sequential Monte Carlo (SMC) algorithms they constitute a powerful approach for parameter estimation and model selection of mathematical models of complex biological systems. A crucial step in the ABC-SMC algorithms, significantly affecting their performance, is the propagation of a set of parameter vectors through a sequence of intermediate distributions using Markov kernels. RESULTS: In this article, we employ Dirichlet process mixtures (DPMs) to design optimal transition kernels and we present an ABC-SMC algorithm with DPM kernels. We illustrate the use of the proposed methodology using real data for the canonical Wnt signaling pathway. A multi-compartment model of the pathway is developed and it is compared to an existing model. The results indicate that DPMs are more efficient in the exploration of the parameter space and can significantly improve ABC-SMC performance. In comparison to alternative sampling schemes that are commonly used, the proposed approach can bring potential benefits in the estimation of complex multimodal distributions. The method is used to estimate the parameters and the initial state of two models of the Wnt pathway and it is shown that the multi-compartment model fits better the experimental data. AVAILABILITY AND IMPLEMENTATION: Python scripts for the Dirichlet Process Gaussian Mixture model and the Gibbs sampler are available at https://sites.google.com/site/kkoutroumpas/software CONTACT: konstantinos.koutroumpas@ecp.fr.
Asunto(s)
Teorema de Bayes , Modelos Estadísticos , Vía de Señalización Wnt , Funciones de Verosimilitud , Método de MontecarloRESUMEN
Water stress is a major cause of tree mortality. In response to drought, leaves wilt due to an increase in petiole flexibility. We present an analytical model coupling petiole mechanics, thermal balance, and xylem hydraulics to investigate the role of petiole flexibility in protecting a tree from water stress. Our model suggests that turgidity-dependent petiole flexibility can significantly attenuate the minimal xylem pressure and thus reduce the risk of cavitation. Moreover, we show that petiole flexibility increases water use efficiency by trees under water stress.
Asunto(s)
Hojas de la Planta/fisiología , Árboles/fisiología , Agua/fisiología , Fenómenos Biomecánicos , Dióxido de Carbono/metabolismo , Deshidratación , Modelos Biológicos , Fotosíntesis/fisiología , Estomas de Plantas/fisiología , SueloRESUMEN
We consider a plant's local leaf area index as a spatially continuous variable, subject to particular reaction-diffusion dynamics of allocation, senescence and spatial propagation. The latter notably incorporates the plant's tendency to form new leaves in bright rather than shaded locations. Applying a generalized Beer-Lambert law allows to link existing foliage to production dynamics. The approach allows for inter-individual variability and competition for light while maintaining robustness-a key weakness of comparable existing models. The analysis of the single plant case leads to a significant simplification of the system's key equation when transforming it into the well studied porous medium equation. Confronting the theoretical model to experimental data of sugar beet populations, differing in configuration density, demonstrates its accuracy.
Asunto(s)
Modelos Biológicos , Plantas/efectos de la radiación , Beta vulgaris/crecimiento & desarrollo , Beta vulgaris/efectos de la radiación , Luz , Conceptos Matemáticos , Fototropismo , Hojas de la Planta/crecimiento & desarrollo , Hojas de la Planta/efectos de la radiaciónRESUMEN
Myeloproliferative neoplasms (MPN) are blood cancers that appear after acquiring a driver mutation in a hematopoietic stem cell. These hematological malignancies result in the overproduction of mature blood cells and, if not treated, induce a risk of cardiovascular events and thrombosis. Pegylated IFN$\alpha $ is commonly used to treat MPN, but no clear guidelines exist concerning the dose prescribed to patients. We applied a model selection procedure and ran a hierarchical Bayesian inference method to decipher how dose variations impact the response to the therapy. We inferred that IFN$\alpha $ acts on mutated stem cells by inducing their differentiation into progenitor cells; the higher the dose, the higher the effect. We found that the treatment can induce long-term remission when a sufficient (patient-dependent) dose is reached. We determined this minimal dose for individuals in a cohort of patients and estimated the most suitable starting dose to give to a new patient to increase the chances of being cured.
Asunto(s)
Teorema de Bayes , Interferón-alfa , Trastornos Mieloproliferativos , Humanos , Trastornos Mieloproliferativos/tratamiento farmacológico , Interferón-alfa/administración & dosificación , Conceptos Matemáticos , Relación Dosis-Respuesta a DrogaRESUMEN
Spatial omics data allow in-depth analysis of tissue architectures, opening new opportunities for biological discovery. In particular, imaging techniques offer single-cell resolutions, providing essential insights into cellular organizations and dynamics. Yet, the complexity of such data presents analytical challenges and demands substantial computing resources. Moreover, the proliferation of diverse spatial omics technologies, such as Xenium, MERSCOPE, CosMX in spatial-transcriptomics, and MACSima and PhenoCycler in multiplex imaging, hinders the generality of existing tools. We introduce Sopa ( https://github.com/gustaveroussy/sopa ), a technology-invariant, memory-efficient pipeline with a unified visualizer for all image-based spatial omics. Built upon the universal SpatialData framework, Sopa optimizes tasks like segmentation, transcript/channel aggregation, annotation, and geometric/spatial analysis. Its output includes user-friendly web reports and visualizer files, as well as comprehensive data files for in-depth analysis. Overall, Sopa represents a significant step toward unifying spatial data analysis, enabling a more comprehensive understanding of cellular interactions and tissue organization in biological systems.
Asunto(s)
Programas Informáticos , Humanos , Procesamiento de Imagen Asistido por Computador/métodos , Análisis de la Célula Individual/métodos , Biología Computacional/métodos , Transcriptoma , AnimalesRESUMEN
Assessing donor/recipient HLA compatibility at the eplet level requires second field DNA typings but these are not always available. These can be estimated from lower-resolution data either manually or with computational tools currently relying, at best, on data containing typing ambiguities. We gathered NGS typing data from 61,393 individuals in 17 French laboratories, for loci A, B, and C (100% of typings), DRB1 and DQB1 (95.5%), DQA1 (39.6%), DRB3/4/5, DPB1, and DPA1 (10.5%). We developed HaploSFHI, a modified iterative maximum likelihood algorithm, to impute second field HLA typings from low- or intermediate-resolution ones. Compared with the reference tools HaploStats, HLA-EMMA, and HLA-Upgrade, HaploSFHI provided more accurate predictions across all loci on two French test sets and four European-independent test sets. Only HaploSFHI could impute DQA1, and solely HaploSFHI and HaploStats provided DRB3/4/5 imputations. The improved performance of HaploSFHI was due to our local and nonambiguous data. We provided explanations for the most common imputation errors and pinpointed the variability of a low number of low-resolution haplotypes. We thus provided guidance to select individuals for whom sequencing would optimize incompatibility assessment and cost-effectiveness of HLA typing, considering not only well-imputed second field typing(s) but also well-imputed eplets.
Asunto(s)
Secuenciación de Nucleótidos de Alto Rendimiento , Donantes de Tejidos , Humanos , Alelos , Haplotipos , Prueba de Histocompatibilidad , Antígenos HLA/genética , Frecuencia de los GenesRESUMEN
BACKGROUND: This retrospective study determined survival responses to immune checkpoint inhibitors (ICIs), comparing mono- (mono) and combo-immunotherapy (combo) in patients with microsatellite instability-high (MSI-H) metastatic colorectal cancer (mCRC) by analyzing quantitative imaging data and clinical factors. METHODS: One hundred fifty patients were included from two centers and divided into training (n = 105) and validation (n = 45) cohorts. Radiologists manually annotated chest-abdomen-pelvis computed tomography and calculated tumor burden. Progression-free survival (PFS) was assessed, and variables were selected through Recursive Feature Elimination. Cutoff values were determined using maximally selected rank statistics to binarize features, forming a risk score with hazard ratio-derived weights. RESULTS: In total, 2258 lesions were annotated with excellent reproducibility. Key variables in the training cohort included: total tumor volume (cutoff: 73 cm3), lesion count (cutoff: 20), age (cutoff: 60) and the presence of peritoneal carcinomatosis. Their respective weights were 1.13, 0.96, 0.91, and 0.38, resulting in a risk score cutoff of 1.36. Low-score patients showed similar overall survival and PFS regardless of treatment, while those with a high-score had significantly worse survivals with mono vs combo (P = 0.004 and P = 0.0001). In the validation set, low-score patients exhibited no significant difference in overall survival and PFS with mono or combo. However, patients with a high-score had worse PFS with mono (P = 0.046). CONCLUSIONS: A score based on total tumor volume, lesion count, the presence of peritoneal carcinomatosis, and age can guide MSI-H mCRC treatment decisions, allowing oncologists to identify suitable candidates for mono and combo ICI therapies.
Asunto(s)
Neoplasias del Colon , Neoplasias Colorrectales , Neoplasias Peritoneales , Humanos , Inhibidores de Puntos de Control Inmunológico/uso terapéutico , Pronóstico , Neoplasias Colorrectales/diagnóstico por imagen , Neoplasias Colorrectales/tratamiento farmacológico , Neoplasias Peritoneales/tratamiento farmacológico , Estudios Retrospectivos , Reproducibilidad de los Resultados , Neoplasias del Colon/tratamiento farmacológico , Inestabilidad de Microsatélites , Reparación de la Incompatibilidad de ADNRESUMEN
BACKGROUND: The DCE-US (Dynamic Contrast-Enhanced Ultrasonography) imaging protocol predicts the vascular modifications compared with Response Evaluation Criteria in Solid Tumors (RECIST) based mainly on morphological changes. A quantitative biomarker has been validated through the DCE-US multi-centric study for early monitoring of the efficiency of anti-angiogenic cancer treatments. In this context, the question of transposing the use of this biomarker to other types of ultrasound scanners, probes and settings has arisen to maintain the follow-up of patients under anti-angiogenic treatments. As a consequence, radiologists encounter standardization issues between the different generations of ultrasound scanners to perform quantitative imaging protocols. PURPOSE: The aim of this study was to develop a new calibration setup to transpose the DCE-US imaging protocol to the new generation of ultrasound scanners using both abdominal and linear probes. METHODS: This calibration method has been designed to be easily reproducible and optimized, reducing the time required and cost incurred. It is based on an original set-up that includes using a concentration splitter to measure the variation of the harmonic signal intensity, obtained from the Area Under the time-intensity Curve (AUC) as a function of various contrast-agent concentrations. The splitter provided four different concentrations simultaneously ranging from 12.5% to 100% of the initial concentration of the SonoVue contrast agent (Bracco Imaging S.p.A., Milan, Italy), therefore, measuring four AUCs in a single injection. The plot of the AUC as a function of the four contrast agent concentrations represents the intensity variation of the harmonic signal: the slope being the calibration parameter. The standardization through this method implied that both generations of ultrasound scanners had to have the same slopes to be considered as calibrated. This method was tested on two ultrasound scanners from the same manufacturer (Aplio500, Aplioi900, Canon Medical Systems, Tokyo, Japan). The Aplio500 used the settings defined by the initial multicenter DCE-US study. The Mechanical Index (MI) and the Color Gain (CG) of the Aplioi900 have been adjusted to match those of the Aplio500. The reliability of the new setup was evaluated in terms of measurement repeatability, and reproducibility with the agreement between the measurements obtained once the two ultrasound scanners were calibrated. RESULTS: The new setup provided excellent repeatability measurements with a value of 96.8%. Once the two ultrasound scanners have been calibrated for both types of probes, the reproducibility was excellent with the agreement between their respective quantitative measurement was at the lowest 95.4% and at the best 98.8%. The settings of the Aplioi900 (Canon Medical Systems) were adjusted to match those of the Aplio500 (Canon Medical Systems) and these validated settings were for the abdominal probe: MI = 0.13 and CG = 34 dB; and for the linear probe: MI = 0.10 and CG = 38 dB. CONCLUSION: This new calibration setup provided reliable measurements and enabled the rapid transfer and the use of the DCE-US imaging protocol on new ultrasound scanners, thus permitting a continuation of the therapeutic evaluation of patients through quantitative imaging.
Asunto(s)
Medios de Contraste , Humanos , Reproducibilidad de los Resultados , Calibración , Ultrasonografía/métodos , Estándares de Referencia , Estudios Multicéntricos como AsuntoRESUMEN
PURPOSE: The objective of the study is to propose the immunotherapy progression decision (iPD) score, a practical tool based on patient features that are available at the first evaluation of immunotherapy treatment, to help oncologists decide whether to continue the treatment or switch rapidly to another therapeutic line when facing a progressive disease patient at the first evaluation. EXPERIMENTAL DESIGN: This retrospective study included 107 patients with progressive disease at first evaluation according to RECIST 1.1. Clinical, radiological, and biological data at baseline and first evaluation were analyzed. An external validation set consisting of 31 patients with similar baseline characteristics was used for the validation of the score. RESULTS: Variables were analyzed in a univariate study. The iPD score was constructed using only independent variables, each considered as a worsening factor for the survival of patients. The patients were stratified in three groups: good prognosis (GP), poor prognosis (PP), and critical prognosis (CP). Each group showed significantly different survivals (GP: 11.4, PP: 4.4, CP: 2.3 months median overall survival, P < 0.001, log-rank test). Moreover, the iPD score was able to detect the pseudoprogressors better than other scores. On the validation set, CP patients had significantly worse survival than PP and GP patients (P < 0.05, log-rank test). CONCLUSIONS: The iPD score provides oncologists with a new evaluation, computable at first progression, to decide whether treatment should be continued (for the GP group), or immediately changed for the PP and CP groups. Further validation on larger cohorts is needed to prove its efficacy in clinical practice.