Search

Chromothripsis is a frequent form of genome instability, whereby a presumably single catastrophic event generates extensive genomic rearrangements of one or multiple chromosome(s). However, little is known about the heterogeneity of chromothripsis across different clones from the same tumour, as well as changes in response to treatment. Here we analyse single-cell genomic and transcriptomic alterations linked with chromothripsis in human p53-deficient medulloblastoma and neural stem cells (n = 9). We reconstruct the order of somatic events, identify early alterations likely linked to chromothripsis and depict the contribution of chromothripsis to malignancy. We characterise subclonal variation of chromothripsis and its effects on extrachromosomal circular DNA, cancer drivers and putatively druggable targets. Furthermore, we highlight the causative role and the fitness consequences of specific rearrangements in neural progenitors., Competing Interests: Competing interests: The authors declare no competing interests., (© 2024. The Author(s).)

Virtually all patients with multiple myeloma become unresponsive to treatment over time. Relapsed/refractory multiple myeloma (RRMM) is accompanied by the clonal evolution of myeloma cells with heterogeneous genomic aberrations and profound changes of the bone marrow microenvironment (BME). However, the molecular mechanisms that drive drug resistance remain elusive. Here, we analyze the heterogeneous tumor cell population and its complex interaction network with the BME of 20 RRMM patients by single cell RNA-sequencing before/after treatment. Subclones with chromosome 1q-gain express a specific transcriptomic signature and frequently expand during treatment. Furthermore, RRMM cells shape an immune suppressive BME by upregulation of inflammatory cytokines and close interaction with the myeloid compartment. It is characterized by the accumulation of PD1 + γδ T-cells and tumor-associated macrophages as well as the depletion of hematopoietic progenitors. Thus, our study resolves transcriptional features of subclones in RRMM and mechanisms of microenvironmental reprogramming with implications for clinical decision-making., (© 2021. The Author(s).)

Motivation: The use of liquid biopsies for cancer patients enables the non-invasive tracking of treatment response and tumor dynamics through single or serial blood drawn tests. Next-generation sequencing assays allow for the simultaneous interrogation of extended sets of somatic single-nucleotide variants (SNVs) in circulating cell-free DNA (cfDNA), a mixture of DNA molecules originating both from normal and tumor tissue cells. However, low circulating tumor DNA (ctDNA) fractions together with sequencing background noise and potential tumor heterogeneity challenge the ability to confidently call SNVs., Results: We present a computational methodology, called Adaptive Base Error Model in Ultra-deep Sequencing data (ABEMUS), which combines platform-specific genetic knowledge and empirical signal to readily detect and quantify somatic SNVs in cfDNA. We tested the capability of our method to analyze data generated using different platforms with distinct sequencing error properties and we compared ABEMUS performances with other popular SNV callers on both synthetic and real cancer patients sequencing data. Results show that ABEMUS performs better in most of the tested conditions proving its reliability in calling low variant allele frequencies somatic SNVs in low ctDNA levels plasma samples., Availability and Implementation: ABEMUS is cross-platform and can be installed as R package. The source code is maintained on Github at http://github.com/cibiobcg/abemus, and it is also available at CRAN official R repository., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2020. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Loss of androgen receptor (AR) signaling dependence occurs in approximately 15%-20% of advanced treatment-resistant prostate cancers, and this may manifest clinically as transformation from a prostate adenocarcinoma histology to a castration-resistant neuroendocrine prostate cancer (CRPC-NE). The diagnosis of CRPC-NE currently relies on a metastatic tumor biopsy, which is invasive for patients and sometimes challenging to diagnose due to morphologic heterogeneity. By studying whole-exome sequencing and whole-genome bisulfite sequencing of cell free DNA (cfDNA) and of matched metastatic tumor biopsies from patients with metastatic prostate adenocarcinoma and CRPC-NE, we identified CRPC-NE features detectable in the circulation. Overall, there was markedly higher concordance between cfDNA and biopsy tissue genomic alterations in patients with CRPC-NE compared with castration-resistant adenocarcinoma, supporting greater intraindividual genomic consistency across metastases. Allele-specific copy number and serial sampling analyses allowed for the detection and tracking of clonal and subclonal tumor cell populations. cfDNA methylation was indicative of circulating tumor content fraction, reflective of methylation patterns observed in biopsy tissues, and was capable of detecting CRPC-NE-associated epigenetic changes (e.g., hypermethylation of ASXL3 and SPDEF; hypomethylation of INSM1 and CDH2). A targeted set combining genomic (TP53, RB1, CYLD, AR) and epigenomic (hypo- and hypermethylation of 20 differential sites) alterations applied to ctDNA was capable of identifying patients with CRPC-NE.

Background: Genomic characterization of prostate cancer (PCa) biopsies may improve criteria for the selection of patients suitable for active surveillance (AS)., Objective: To identify somatic genomic aberrations associated with adverse outcome as AS protocol exclusion indicators., Design, Setting and Participants: Whole-exome sequencing profiles were generated for Gleason score (GS)=3+3 biopsies obtained from 54 PCa patients enrolled in two AS protocols. Patients were selected as representative of a nonindolent population, consisting of 27 patients who dropped out from AS due to upgrading (ie, finding of GS>3+3 at a follow-up biopsy) within 2 yr, and a potentially indolent population, consisting of 27 patients in AS for ≥4 yr without any evidence of reclassification., Outcome Measurements and Statistical Analysis: The genomic alteration landscape of core biopsies was analyzed using an integrated computational pipeline and correlated with patient reclassification due to upgrading., Results and Limitations: Of all the GS=3+3 biopsies of the study cohort, 34% showed clear evidence of somatic copy number aberrations along the genome. Of these, 39% came from the potentially indolent and 61% from the nonindolent population. Single-nucleotide variants demonstrated low allelic fractions and included a common F133C mutation in the SPOP gene. The minimally altered genomic landscape of the study cohort presented a distinct set of monoallelic deletions, including on 8p, 13q, 16q, and 21q, and rare amplifications of 8q, which were observed in both AS patient populations. Concerning lesions typically associated with adverse outcome, PTEN deletions and MYC amplification, though observed in a small number of cases, were detected exclusively or preferentially, respectively, in nonindolent patients. Such molecular findings were confirmed by immunohistochemistry on the same tissue blocks. The small sample size and the retrospective nature of the analysis represent the main study limitations., Conclusions: Genomic features enriched in aggressive tumors can be detected in GS=3+3 core biopsies of AS patients., Patient Summary: PTEN and MYC alterations at the time of diagnosis would deserve investigation in larger cohorts of AS patients to assess their potential as biomarkers for a more precise/earlier identification of patients at risk of reclassification., (Copyright © 2018 The Author(s). Published by Elsevier B.V. All rights reserved.)

Purpose: Metastatic castration-resistant prostate cancer (mCRPC) is a lethal but clinically heterogeneous disease, with patients having variable benefit from endocrine and cytotoxic treatments. Intrapatient genomic heterogeneity could be a contributing factor to this clinical heterogeneity. Here, we used whole-genome sequencing (WGS) to investigate genomic heterogeneity in 21 previously treated CRPC metastases from 10 patients to investigate intrapatient molecular heterogeneity (IPMH). Experimental Design: WGS was performed on topographically separate metastases from patients with advanced metastatic prostate cancer. IPMH of the RB1 gene was identified and further evaluated by FISH and IHC assays., Results: WGS identified limited IPMH for putative driver events. However, heterogeneous genomic aberrations of RB1 were detected. We confirmed the presence of these RB1 somatic copy-number aberrations, initially identified by WGS, with FISH, and identified novel structural variants involving RB1 in 6 samples from 3 of these 10 patients (30%; 3/10). WGS uncovered a novel deleterious RB1 structural lesion constituted of an intragenic tandem duplication involving multiple exons and associating with protein loss. Using RB1 IHC in a large series of mCRPC biopsies, we identified heterogeneous expression in approximately 28% of mCRPCs., Conclusions: mCRPCs have a high prevalence of RB1 genomic aberrations, with structural variants, including rearrangements, being common. Intrapatient genomic and expression heterogeneity favors RB1 aberrations as late, subclonal events that increase in prevalence due to treatment-selective pressures., (©2018 American Association for Cancer Research.)

Androgen receptor (AR) gene aberrations are rare in prostate cancer before primary hormone treatment but emerge with castration resistance. To determine AR gene status using a minimally invasive assay that could have broad clinical utility, we developed a targeted next-generation sequencing approach amenable to plasma DNA, covering all AR coding bases and genomic regions that are highly informative in prostate cancer. We sequenced 274 plasma samples from 97 castration-resistant prostate cancer patients treated with abiraterone at two institutions. We controlled for normal DNA in patients' circulation and detected a sufficiently high tumor DNA fraction to quantify AR copy number state in 217 samples (80 patients). Detection of AR copy number gain and point mutations in plasma were inversely correlated, supported further by the enrichment of nonsynonymous versus synonymous mutations in AR copy number normal as opposed to AR gain samples. Whereas AR copy number was unchanged from before treatment to progression and no mutant AR alleles showed signal for acquired gain, we observed emergence of T878A or L702H AR amino acid changes in 13% of tumors at progression on abiraterone. Patients with AR gain or T878A or L702H before abiraterone (45%) were 4.9 and 7.8 times less likely to have a ≥50 or ≥90% decline in prostate-specific antigen (PSA), respectively, and had a significantly worse overall [hazard ratio (HR), 7.33; 95% confidence interval (CI), 3.51 to 15.34; P = 1.3 × 10(-9)) and progression-free (HR, 3.73; 95% CI, 2.17 to 6.41; P = 5.6 × 10(-7)) survival. Evaluation of plasma AR by next-generation sequencing could identify cancers with primary resistance to abiraterone., (Copyright © 2015, American Association for the Advancement of Science.)

Somatic mutations arise and accumulate both during tumor genesis and progression. However, the order in which mutations occur is an open question and the inference of the temporal ordering at the gene level could potentially impact on patient treatment. Thus, exploiting recent observations suggesting that the occurrence of mutations is a non-memoryless process, we developed a computational approach to infer timed oncogenetic directed acyclic graphs (TO-DAGs) from human tumor mutation data. Such graphs represent the path and the waiting times of alterations during tumor evolution. The probability of occurrence of each alteration in a path is the probability that the alteration occurs when all alterations prior to it have occurred. The waiting time between an alteration and the subsequent is modeled as a stochastic function of the conditional probability of the event given the occurrence of the previous one. TO-DAG performances have been evaluated both on synthetic data and on somatic non-silent mutations from prostate cancer and melanoma patients and then compared with those of current well-established approaches. TO-DAG shows high performance scores on synthetic data and recognizes mutations in gatekeeper tumor suppressor genes as trigger for several downstream mutational events in the human tumor data.

Prion diseases are fatal neurodegenerative disorders in mammals and other animal species. In humans, about 15% of these maladies are caused by pathogenic mutations (PMs) in the gene encoding for the prion protein (PrP(C)). Seven PMs are located in the naturally unfolded PrP(C) N-terminal domain, which constitutes about half of the protein. Intriguingly and in sharp contrast to other PMs clustered in the folded domain, N-terminal PMs barely affect the conversion to the pathogenic (scrapie, or PrP(Sc)) isoform of PrP(C). Here, we hypothesize that the neurotoxicity of these PMs arises from changes in structural determinants of the N-terminal domain, affecting the protein binding with its cellular partners and/or the cotranslational translocation during the PrP(C) biosynthesis. We test this idea by predicting the conformational ensemble of the wild-type (WT) and mutated mouse PrP(C) N-terminal domain, whose sequence is almost identical to that of the human one and for which the largest number of in vivo data is available. The conformational properties of the WT are consistent with those inferred experimentally. Importantly, the PMs turn out to affect in a subtle manner the intramolecular contacts in the putative N-terminal domain binding sites for Cu(2+) ions, sulphated glycosaminoglycans, and other known PrP(C) cellular partners. The PMs also alter the local structural features of the transmembrane domain and adjacent stop transfer effector, which act together to regulate the protein topology. These results corroborate the hypothesis that N-terminal PMs affect the PrP(C) binding to functional interactors and/or the translocation.

We investigated the interaction between hypoxia and hypercapnia on ventilation and on cerebro-cardio-vascular control. A group of 12 healthy subjects performed rebreathing tests to determine the ventilatory response to hypoxia, at different levels of carbon dioxide (CO(2)), and to normoxic hypercapnia. Oxygen saturation (SaO(2)), end-tidal CO(2) (et-CO(2)), minute ventilation, blood pressure, R-R interval and mid-cerebral artery flow velocity (MCFV) were continuously recorded. The hypoxic ventilatory response significantly increased under hypercapnia and decreased under hypocapnia (slopes L/min/% Sa O(2): -0.33 +/- 0.05, -0.74 +/- 0.02 and -1.59 +/- 0.3, p < 0.0001, in hypocapnia, normocapnia and hypercapnia, respectively). At similar degrees of ventilation, MCFV increased more markedly during normocapnic hypoxia than normoxic hypercapnia; the slopes linking MCFV to hypoxia remained unchanged at increasing levels of et-CO(2), whereas the regression lines were shifted upward. The R-R interval decreased more markedly during normocapnic hypoxia than normoxic hypercapnia and the arterial baroreflex sensitivity was decreased only by hypoxia. Cardiovascular responses to hypoxia were not affected by different levels of et-CO(2). We conclude that concomitant hypoxia and hypercapnia, while increasing ventilation synergistically, exert an additive effect on cerebral blood flow. Increased sympathetic activity (and reduced baroreflex sensitivity) is one of the mechanisms by which hypoxia stimulates cardiac sympathetic activity.

Sympathetic hyperactivity and parasympathetic withdrawal may cause and sustain hypertension. This autonomic imbalance is in turn related to a reduced or reset arterial baroreflex sensitivity and chemoreflex-induced hyperventilation. Slow breathing at 6 breaths/min increases baroreflex sensitivity and reduces sympathetic activity and chemoreflex activation, suggesting a potentially beneficial effect in hypertension. We tested whether slow breathing was capable of modifying blood pressure in hypertensive and control subjects and improving baroreflex sensitivity. Continuous noninvasive blood pressure, RR interval, respiration, and end-tidal CO2 (CO2-et) were monitored in 20 subjects with essential hypertension (56.4+/-1.9 years) and in 26 controls (52.3+/-1.4 years) in sitting position during spontaneous breathing and controlled breathing at slower (6/min) and faster (15/min) breathing rate. Baroreflex sensitivity was measured by autoregressive spectral analysis and "alpha angle" method. Slow breathing decreased systolic and diastolic pressures in hypertensive subjects (from 149.7+/-3.7 to 141.1+/-4 mm Hg, P<0.05; and from 82.7+/-3 to 77.8+/-3.7 mm Hg, P<0.01, respectively). Controlled breathing (15/min) decreased systolic (to 142.8+/-3.9 mm Hg; P<0.05) but not diastolic blood pressure and decreased RR interval (P<0.05) without altering the baroreflex. Similar findings were seen in controls for RR interval. Slow breathing increased baroreflex sensitivity in hypertensives (from 5.8+/-0.7 to 10.3+/-2.0 ms/mm Hg; P<0.01) and controls (from 10.9+/-1.0 to 16.0+/-1.5 ms/mm Hg; P<0.001) without inducing hyperventilation. During spontaneous breathing, hypertensive subjects showed lower CO2 and faster breathing rate, suggesting hyperventilation and reduced baroreflex sensitivity (P<0.001 versus controls). Slow breathing reduces blood pressure and enhances baroreflex sensitivity in hypertensive patients. These effects appear potentially beneficial in the management of hypertension.

Polycythemia is one of the key factors involved in the chronic mountain sickness syndrome, a condition frequent in Andean natives but whose causes still remain unclear. In theory, polycythemia may be secondary to abnormalities in ventilation, occurring during day or night (e.g. due to sleep abnormalities) stimulating excessive erythropoietin (Epo) production, or else it may result from either autogenous production, or from co-factors like cobalt. To assess the importance of these points, we studied subjects with or without polycythemia, born and living in Cerro de Pasco (Peru, 4330m asl, CP) and evaluated the relationship between Epo and respiratory variables both in CP and sea level. We also assessed the relationship between sleep abnormalities and the circadian rhythm of Epo. Polycythemic subjects showed higher Epo in all conditions, lower SaO2 and hypoxic ventilatory response, higher physiological dead space and higher CO2, suggesting ventilatory inefficiency. Epo levels could be highly modified by the level of oxygenation, and were related to similar directional changes in SaO2. Cobalt levels were normal in all subjects and correlated poorly with hematologic variables. The diurnal variations in Epo were grossly abnormal in polycythemic subjects, with complete loss of the circadian rhythm. These abnormalities correlated with the levels of hypoxemia during the night, but not with sleep abnormalities, which were only minor even in polycythemic subjects. The increased Epo production is mainly related to a greater ventilatory inefficiency, and not to altered sensitivity to hypoxia, cobalt or sleep abnormalities. Improving oxygenation can represent a possible therapeutic option for this syndrome.