Search

Abstract Extracellular vesicles (EVs) are a new mechanism of cellular communication, by delivering their cargo into target cells to modulate molecular pathways. EV-mediated crosstalk contributes to tumor survival and resistance to cellular stress. However, the role of EVs in B-cell Acute Lymphoblastic Leukaemia (B-ALL) awaits to be thoroughly investigated. We recently published that ActivinA increases intracellular calcium levels and promotes actin polymerization in B-ALL cells. These biological processes guide cytoskeleton reorganization, which is a crucial event for EV secretion and internalization. Hence, we investigated the role of EVs in the context of B-ALL and the impact of ActivinA on this phenomenon. We demonstrated that leukemic cells release a higher number of EVs in response to ActivinA treatment, and they can actively uptake EVs released by other B-ALL cells. Under culture-induced stress conditions, EVs coculture promoted cell survival in B-ALL cells in a dose-dependent manner. Direct stimulation of B-ALL cells with ActivinA or with EVs isolated from ActivinA-stimulated cells was even more effective in preventing cell death. This effect can be possibly ascribed to the increase of vesiculation and modifications of EV-associated microRNAs induced by ActivinA. These data demonstrate that ActivinA boosts EV-mediated B-ALL crosstalk, improving leukemia survival in stress conditions.

Abstract Morgana is a ubiquitous HSP90 co-chaperone protein coded by the CHORDC1 gene. Morgana heterozygous mice develop with age a myeloid malignancy resembling human atypical myeloid leukemia (aCML), now renamed MDS/MPN with neutrophilia. Patients affected by this pathology exhibit low Morgana levels in the bone marrow (BM), suggesting that Morgana downregulation plays a causative role in the human malignancy. A decrease in Morgana expression levels is also evident in the BM of a subgroup of Philadelphia-positive (Ph+) chronic myeloid leukemia (CML) patients showing resistance or an incomplete response to imatinib. Despite the relevance of these data, the mechanism through which Morgana expression is downregulated in patients’ bone marrow remains unclear. In this study, we investigated the possibility that Morgana expression is regulated by miRNAs and we demonstrated that Morgana is under the control of four miRNAs (miR-15a/b and miR-26a/b) and that miR-15a may account for Morgana downregulation in CML patients.

Abstract γδ T cells provide rapid cellular immunity against pathogens. Here, we conducted matched single-cell RNA-sequencing and γδ-TCR-sequencing to delineate the molecular changes in γδ T cells during a longitudinal study following mRNA SARS-CoV-2 vaccination. While the first dose of vaccine primes Vδ2 T cells, it is the second administration that significantly boosts their immune response. Specifically, the second vaccination uncovers memory features of Vδ2 T cells, shaped by the induction of AP-1 family transcription factors and characterized by a convergent central memory signature, clonal expansion, and an enhanced effector potential. This temporally distinct effector response of Vδ2 T cells was also confirmed in vitro upon stimulation with SARS-CoV-2 spike-peptides. Indeed, the second challenge triggers a significantly higher production of IFNγ by Vδ2 T cells. Collectively, our findings suggest that mRNA SARS-CoV-2 vaccination might benefit from the establishment of long-lasting central memory Vδ2 T cells to confer protection against SARS-CoV-2 infection.

Abstract ALK and ROS1 fusions are effectively targeted by tyrosine kinase inhibitors (TKIs), however patients inevitably relapse after an initial response, often due to kinase domain mutations. We investigated circulating DNA from TKI-relapsed NSCLC patients by deep-sequencing. New EML4::ALK substitutions, L1198R, C1237Y and L1196P, were identified in the plasma of NSCLC ALK patients and characterized in a Ba/F3 cell model. Variants C1237Y and L1196P demonstrated pan-inhibitor resistance across 5 clinical and 2 investigational TKIs.

Abstract Background Copy number alterations (CNAs) are genetic changes commonly found in cancer that involve different regions of the genome and impact cancer progression by affecting gene expression and genomic stability. Computational techniques can analyze copy number data obtained from high-throughput sequencing platforms, and various tools visualize and analyze CNAs in cancer genomes, providing insights into genetic mechanisms driving cancer development and progression. However, tools for visualizing copy number data in cancer research have some limitations. In fact, they can be complex to use and require expertise in bioinformatics or computational biology. While copy number data analysis and visualization provide insights into cancer biology, interpreting results can be challenging, and there may be multiple explanations for observed patterns of copy number alterations. Results We created Control-FREEC Viewer, a tool that facilitates effective visualization and exploration of copy number data. With Control-FREEC Viewer, experimental data can be easily loaded by the user. After choosing the reference genome, copy number data are displayed in whole genome or single chromosome view. Gain or loss on a specific gene can be found and visualized on each chromosome. Analysis parameters for subsequent sessions can be stored and images can be exported in raster and vector formats. Conclusions Control-FREEC Viewer enables users to import and visualize data analyzed by the Control-FREEC tool, as well as by other tools sharing a similar tabular output, providing a comprehensive and intuitive graphical user interface for data visualization.

Human mesenchymal stem cells modulate the immune response and are good candidates for cell therapy in neuroinflammatory brain disorders affecting both adult and premature infants. Recent evidence indicates that through their secretome, mesenchymal stem cells direct microglia, brain-resident immune cells, toward pro-regenerative functions, but the mechanisms underlying microglial phenotypic transition are still under investigation. Using an in vitro coculture approach combined with transcriptomic analysis, we identified the extracellular matrix as the most relevant pathway altered by the human mesenchymal stem cell secretome in the response of microglia to inflammatory cytokines. We confirmed extracellular matrix remodeling in microglia exposed to the mesenchymal stem cell secretome via immunofluorescence analysis of the matrix component fibronectin and the extracellular crosslinking enzyme transglutaminase-2. Furthermore, an analysis of hallmark microglial functions revealed that changes in the extracellular matrix enhance ruffle formation by microglia and cell motility. These findings point to extracellular matrix changes, associated plasma membrane remodeling, and enhanced microglial migration as novel mechanisms by which mesenchymal stem cells contribute to the pro-regenerative microglial transition.

Abstract Recurring sequences of genomic alterations occurring across patients can highlight repeated evolutionary processes with significant implications for predicting cancer progression. Leveraging the ever-increasing availability of cancer omics data, here we unveil cancer’s evolutionary signatures tied to distinct disease outcomes, representing “favored trajectories” of acquisition of driver mutations detected in patients with similar prognosis. We present a framework named ASCETIC (Agony-baSed Cancer EvoluTion InferenCe) to extract such signatures from sequencing experiments generated by different technologies such as bulk and single-cell sequencing data. We apply ASCETIC to (i) single-cell data from 146 myeloid malignancy patients and bulk sequencing from 366 acute myeloid leukemia patients, (ii) multi-region sequencing from 100 early-stage lung cancer patients, (iii) exome/genome data from 10,000+ Pan-Cancer Atlas samples, and (iv) targeted sequencing from 25,000+ MSK-MET metastatic patients, revealing subtype-specific single-nucleotide variant signatures associated with distinct prognostic clusters. Validations on several datasets underscore the robustness and generalizability of the extracted signatures.

Abstract Within the chromatin, distal elements interact with promoters to regulate specific transcriptional programs. Histone acetylation, interfering with the net charges of the nucleosomes, is a key player in this regulation. Here, we report that the oncoprotein SET is a critical determinant for the levels of histone acetylation within enhancers. We disclose that a condition in which SET is accumulated, the severe Schinzel-Giedion Syndrome (SGS), is characterized by a failure in the usage of the distal regulatory regions typically employed during fate commitment. This is accompanied by the usage of alternative enhancers leading to a massive rewiring of the distal control of the gene transcription. This represents a (mal)adaptive mechanism that, on one side, allows to achieve a certain degree of differentiation, while on the other affects the fine and corrected maturation of the cells. Thus, we propose the differential in cis-regulation as a contributing factor to the pathological basis of SGS and possibly other the SET-related disorders in humans.

Abstract Background Longitudinal single-cell sequencing experiments of patient-derived models are increasingly employed to investigate cancer evolution. In this context, robust computational methods are needed to properly exploit the mutational profiles of single cells generated via variant calling, in order to reconstruct the evolutionary history of a tumor and characterize the impact of therapeutic strategies, such as the administration of drugs. To this end, we have recently developed the LACE framework for the Longitudinal Analysis of Cancer Evolution. Results The LACE 2.0 release aimed at inferring longitudinal clonal trees enhances the original framework with new key functionalities: an improved data management for preprocessing of standard variant calling data, a reworked inference engine, and direct connection to public databases. Conclusions All of this is accessible through a new and interactive Shiny R graphical interface offering the possibility to apply filters helpful in discriminating relevant or potential driver mutations, set up inferential parameters, and visualize the results. The software is available at: github.com/BIMIB-DISCo/LACE.

IntroductionHigher frequencies of mucosal-associated invariant T (MAIT) cells were associated with an increased adaptive response to mRNA BNT162b2 SARS-CoV-2 vaccine, however, the mechanistic insights into this relationship are unknown. In the present study, we hypothesized that the TNF response of MAIT cells supports B cell activation following SARS-CoV-2 immunization.MethodsTo investigate the effects of repeated SARS-CoV-2 vaccinations on the peripheral blood mononuclear cells (PBMCs), we performed a longitudinal single cell (sc)RNA-seq and scTCR-seq analysis of SARS-CoV-2 vaccinated healthy adults with two doses of the Pfizer-BioNTech BNT162b2 mRNA vaccine. Collection of PBMCs was performed 1 day before, 3 and 17 days after prime vaccination, and 3 days and 3 months following vaccine boost. Based on scRNA/TCR-seq data related to regulatory signals induced by the vaccine, we used computational approaches for the functional pathway enrichment analysis (Reactome), dynamics of the effector cell-polarization (RNA Velocity and CellRank), and cell-cell communication (NicheNet).ResultsWe identified MAIT cells as an important source of TNF across circulating lymphocytes in response to repeated SARS-CoV-2 BNT162b2 vaccination. The TNFhigh signature of MAIT cells was induced by the second administration of the vaccine. Notably, the increased TNF expression was associated with MAIT cell proliferation and efficient anti-SARS-CoV-2 antibody production. Finally, by decoding the ligand-receptor interactions and incorporating intracellular signaling, we predicted TNFhigh MAIT cell interplay with different B cell subsets. In specific, predicted TNF-mediated activation was selectively directed to conventional switched memory B cells, which are deputed to high-affinity long-term memory.DiscussionOverall, our results indicate that SARS-CoV-2 BNT162b2 vaccination influences MAIT cell frequencies and their transcriptional effector profile with the potential to promote B cell activation. This research also provides a blueprint for the promising use of MAIT cells as cellular adjuvants in mRNA-based vaccines.

Gastrointestinal graft-versus-host disease (GvHD) is a major cause of mortality and morbidity following allogeneic bone marrow transplantation (allo-BMT). Chemerin is a chemotactic protein that recruits leukocytes to inflamed tissues by interacting with ChemR23/CMKLR1, a chemotactic receptor expressed by leukocytes, including macrophages. During acute GvHD, chemerin plasma levels were strongly increased in allo-BM-transplanted mice. The role of the chemerin/CMKLR1 axis in GvHD was investigated using Cmklr1-KO mice. WT mice transplanted with an allogeneic graft from Cmklr1-KO donors (t-KO) had worse survival and more severe GvHD. Histological analysis demonstrated that the gastrointestinal tract was the organ mostly affected by GvHD in t-KO mice. The severe colitis of t-KO mice was characterized by massive neutrophil infiltration and tissue damage associated with bacterial translocation and exacerbated inflammation. Similarly, Cmklr1-KO recipient mice showed increased intestinal pathology in both allogeneic transplant and dextran sulfate sodium–induced colitis. Notably, the adoptive transfer of WT monocytes into t-KO mice mitigated GvHD manifestations by decreasing gut inflammation and T cell activation. In patients, higher chemerin serum levels were predictive of GvHD development. Overall, these results suggest that CMKLR1/chemerin may be a protective pathway for the control of intestinal inflammation and tissue damage in GvHD.

Cyclin-dependent kinase (CDK) 4/6 inhibitors have significantly improved progression-free survival in hormone-receptor-positive (HR+), human-epidermal-growth-factor-receptor-type-2-negative (HER2−) metastatic luminal breast cancer (mLBC). Several studies have shown that in patients with endocrine-sensitive or endocrine-resistant LBC, the addition of CDK4/6 inhibitors to endocrine therapy significantly prolongs progression-free survival. However, the percentage of patients who are unresponsive or refractory to these therapies is as high as 40%, and no reliable and reproducible biomarkers have been validated to select a priori responders or refractory patients. The selection of mutant clones in the target oncoprotein is the main cause of resistance. Other mechanisms such as oncogene amplification/overexpression or mutations in other pathways have been described in several models. In this study, we focused on palbociclib, a selective CDK4/6 inhibitor. We generated a human MCF-7 luminal breast cancer cell line that was able to survive and proliferate at different concentrations of palbociclib and also showed cross-resistance to abemaciclib. The resistant cell line was characterized via RNA sequencing and was found to strongly activate the epithelial-to-mesenchymal transition. Among the top deregulated genes, we found a dramatic downregulation of the CDK4 inhibitor CDKN2B and an upregulation of the TWIST1 transcription factor. TWIST1 was further validated as a target for the reversal of palbociclib resistance. This study provides new relevant information about the mechanisms of resistance to CDK4/6 inhibitors and suggests potential new markers for patients’ follow-up care during treatment.

Summary: We outline the features of the R package SparseSignatures and its application to determine the signatures contributing to mutation profiles of tumor samples. We describe installation details and illustrate a step-by-step approach to (1) prepare the data for signature analysis, (2) determine the optimal parameters, and (3) employ them to determine the signatures and related exposure levels in the point mutation dataset.For complete details on the use and execution of this protocol, please refer to Lal et al. (2021). : Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

The pathological consequences of type 2 diabetes mellitus (T2DM) also involve the central nervous system; indeed, T2DM patients suffer from learning and memory disabilities with a higher risk of developing dementia. Although several factors have been proposed as possible contributors, how neuroactive steroids and the gut microbiome impact brain pathophysiology in T2DM remain unexplored. On this basis, in male Zucker diabetic fatty (ZDF) rats, we studied whether T2DM alters memory abilities using the novel object recognition test, neuroactive steroid levels by liquid chromatography–tandem mass spectrometry, hippocampal parameters using molecular assessments, and gut microbiome composition using 16S next-generation sequencing. Results obtained reveal that T2DM worsens memory abilities and that these are correlated with increased levels of corticosterone in plasma and with a decrease in allopregnanolone in the hippocampus, where neuroinflammation, oxidative stress, and mitochondrial dysfunction were reported. Interestingly, our analysis highlighted a small group of taxa strictly related to both memory impairment and neuroactive steroid levels. Overall, the data underline an interesting role for allopregnanolone and microbiota that may represent candidates for the development of therapeutic strategies.

Schinzel-Giedion syndrome (SGS) is a fatal developmental syndrome characterized by severe intellectual and physical deficits due, at least in part, to early neurodegeneration. Here the authors introduce a human SGS model that displays disease-relevant phenotypes to demonstrate that neuronal death in SGS originates from developmental alterations mainly in safeguarding cell identity and homeostasis.

Background More than 50% of all patients with colorectal cancer (CRC) develop liver metastases (CLM), a clinical condition characterized by poor prognosis and lack of reliable prognostic markers. Vδ1 cells are a subset of tissue-resident gamma delta (γδ) T lymphocytes endowed with a broad array of antitumor functions and showing a natural high tropism for the liver. However, little is known about their impact in the clinical outcomes of CLM.Methods We isolated human γδ T cells from peripheral blood (PB) and peritumoral (PT) tissue of 93 patients undergone surgical procedures to remove CLM. The phenotype of freshly purified γδ T cells was assessed by multiparametric flow cytometry, the transcriptional profiles by single cell RNA-sequencing, the functional annotations by Gene Ontology enrichment analyses and the clonotype by γδ T cell receptor (TCR)-sequencing.Results The microenvironment of CLM is characterized by a heterogeneous immune infiltrate comprising different subsets of γδ tumor-infiltrating lymphocytes (TILs) able to egress the liver and re-circulate in PB. Vδ1 T cells represent the largest population of γδ TILs within the PT compartment of CLM that is greatly enriched in Vδ1 T effector (TEF) cells expressing constitutive high levels of CD69. These Vδ1 CD69+ TILs express a distinct phenotype and transcriptional signature, show high antitumor potential and correlate with better patient clinical outcomes in terms of lower numbers of liver metastatic lesions and longer overall survival (OS). Moreover, intrahepatic CD69+ Vδ1 TILs can egress CLM tissue to re-circulate in PB, where they retain a phenotype, transcriptional signature and TCR clonal repertoires resembling their liver origin. Importantly, even the increased frequencies of the CD69+ terminally differentiated (TEMRA) Vδ1 cells in PB of patients with CLM significantly correlate with longer OS. The positive prognostic score of high frequencies of CD69+ TEMRA Vδ1 cells in PB is independent from the neoadjuvant chemotherapy and immunotherapy regimens administered to patients with CLM prior surgery.Conclusions The enrichment of tissue-resident CD69+ Vδ1 TEMRA cells re-circulating at high frequencies in PB of patients with CLM limits tumor progression and represents a new important clinical tool to either predict the natural history of CLM or develop alternative therapeutic protocols of cellular therapies.

We identified Pycard and BC017158 genes as putative effectors of the Quantitative Trait locus (QTL) that we mapped at distal chromosome 7 named Irm1 for Inflammatory response modulator 1, controlling acute inflammatory response (AIR) and the production of IL-1β, dependent on the activation of the NLRP3 inflammasome. We obtained the mapping through genome-wide linkage analysis of Single Nucleotide Polymorphisms (SNPs) in a cross between High (AIRmax) and Low (AIRmin) responder mouse lines that we produced by several generations of bidirectional selection for Acute Inflammatory Response. A highly significant linkage signal (LOD score peak of 72) for ex vivo IL-1β production limited a 4 Mbp interval to chromosome 7. Sequencing of the locus region revealed 14 SNPs between “High” and “Low” responders that narrowed the locus to a 420 Kb interval. Variants were detected in non-coding regions of Itgam, Rgs10 and BC017158 genes and at the first exon of Pycard gene, resulting in an E19K substitution in the protein ASC (apoptosis associated speck-like protein containing a CARD) an adaptor molecule in the inflammasome complex. Silencing of BC017158 inhibited IL1-β production by stimulated macrophages and the E19K ASC mutation carried by AIRmin mice impaired the ex vivo IL-1β response and the formation of ASC specks in stimulated cells. IL-1β and ASC specks play major roles in inflammatory reactions and in inflammation-related diseases. Our results delineate a novel genetic factor and a molecular mechanism affecting the acute inflammatory response.

Summary: A key task of genomic surveillance of infectious viral diseases lies in the early detection of dangerous variants. Unexpected help to this end is provided by the analysis of deep sequencing data of viral samples, which are typically discarded after creating consensus sequences. Such analysis allows one to detect intra-host low-frequency mutations, which are a footprint of mutational processes underlying the origination of new variants. Their timely identification may improve public-health decision-making with respect to traditional approaches exploiting consensus sequences. We present the analysis of 220,788 high-quality deep sequencing SARS-CoV-2 samples, showing that many spike and nucleocapsid mutations of interest associated to the most circulating variants, including Beta, Delta, and Omicron, might have been intercepted several months in advance. Furthermore, we show that a refined genomic surveillance system leveraging deep sequencing data might allow one to pinpoint emerging mutation patterns, providing an automated data-driven support to virologists and epidemiologists.

Abstract Background Oxidative stress is a hallmark of many cancers. The increment in reactive oxygen species (ROS), resulting from an increased mitochondrial respiration, is the major cause of oxidative stress. Cell fate is known to be intricately linked to the amount of ROS produced. The direct generation of ROS is also one of the mechanisms exploited by common anticancer therapies, such as chemotherapy. Methods We assessed the role of NFKBIA with various approaches, including in silico analyses, RNA-silencing and xenotransplantation. Western blot analyses, immunohistochemistry and RT-qPCR were used to detect the expression of specific proteins and genes. Immunoprecipitation and pull-down experiments were used to evaluate protein-protein interactions. Results Here, by using an in silico approach, following the identification of NFKBIA (the gene encoding IκBα) amplification in various cancers, we described an inverse correlation between IκBα, oxidative metabolism, and ROS production in lung cancer. Furthermore, we showed that novel IκBα targeting compounds combined with cisplatin treatment promote an increase in ROS beyond the tolerated threshold, thus causing death by oxytosis. Conclusions NFKBIA amplification and IκBα overexpression identify a unique cancer subtype associated with specific expression profile and metabolic signatures. Through p65-NFKB regulation, IκBα overexpression favors metabolic rewiring of cancer cells and distinct susceptibility to cisplatin. Lastly, we have developed a novel approach to disrupt IκBα/p65 interaction, restoring p65-mediated apoptotic responses to cisplatin due to mitochondria deregulation and ROS-production.

ETNK1 mutations are recurrent in leukemia but how they contribute to oncogenesis is still unclear. Here, the authors show that ETNK1 mutations increase mitochondrial activity, ROS and H2AX levels and that these effects can be rescued upon phosphoethanolamine supplementation.

Abstract Leu387Trp mutation, aroused in an imatinib‐non‐responsive CML patient, was selected by imatinib treatment along with other unknown factors responsible for resistance, and then it was overcome by bosutinib. These results will be useful for treating patients with this rare mutation and will advise against automatically considering a new mutation as the cause of TKI resistance.

Summary: We describe the procedures to perform the following: (1) the de novo discovery of mutational signatures from raw sequencing data of viral samples and (2) the association of existing viral mutational signatures to the samples of a given dataset. The goal is to identify and characterize the nucleotide substitution patterns related to the mutational processes that underlie the origination of variants in viral genomes. The VirMutSig protocol is available at this link: https://github.com/BIMIB-DISCo/VirMutSig.For complete information on the theoretical aspects of this protocol, please refer to Graudenzi et al. (2021).

We present a large-scale analysis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) substitutions, considering 1,585,456 high-quality raw sequencing samples, aimed at investigating the existence and quantifying the effect of mutational processes causing mutations in SARS-CoV-2 genomes when interacting with the human host. As a result, we confirmed the presence of three well-differentiated mutational processes likely ruled by reactive oxygen species (ROS), apolipoprotein B editing complex (APOBEC), and adenosine deaminase acting on RNA (ADAR). We then evaluated the activity of these mutational processes in different continental groups, showing that some samples from Africa present a significantly higher number of substitutions, most likely due to higher APOBEC activity. We finally analyzed the activity of mutational processes across different SARS-CoV-2 variants, and we found a significantly lower number of mutations attributable to APOBEC activity in samples assigned to the Omicron variant.

Atypical chronic myeloid leukemia is a rare disease whose pathogenesis has long been debated. It currently belongs to the group of myelodysplastic/myeloproliferative disorders. In this review, an overview on the current knowledge about diagnosis, prognosis, and genetics is presented, with a major focus on the recent molecular findings. We describe here the molecular pathogenesis of the disease, focusing on the mechanisms of action of the main mutations as well as on gene expression profiling. We also present the treatment options focusing on emerging targeted therapies.

Summary: Human Vδ2 cells are innate-like γδ T effectors performing potent immune surveillance against tumors. The constitutive expression of NKG2A identifies a subset of Vδ2 T cells licensed with an intrinsic hyper-responsiveness against cancer. Indeed, the transcriptomic profiles of NKG2A+ and NKG2A− cells characterize two distinct “intralineages” of Vδ2 T lymphocytes that appear early during development, keep their phenotypes, and show self-renewal capabilities in adult life. The hyper-responsiveness of NKG2A+ Vδ2 T cells is counterbalanced by the inhibitory signaling delivered by human leukocyte antigen E (HLA-E) expressed on malignant cells as a tumor-escape mechanism. However, either masking or knocking out NKG2A restores the capacity of Vδ2 T cells to exert the highest effector functions even against HLA-E+ tumors. This is highly relevant in the clinic, as the different degrees of engagement of the NKG2A-HLA-E checkpoint in hepatocellular carcinoma, glioblastoma, and non-small cell lung cancer directly impact patients’ overall survival. These findings open avenues for developing combined cellular and immunologic anticancer therapies.

DOCK8 deficiency is a combined immunodeficiency due to biallelic variants in dedicator of cytokinesis 8 (DOCK8) gene. The disease has a wide clinical spectrum encompassing recurrent infections (candidiasis, viral and bacterial infections), virally driven malignancies and immune dysregulatory features, including autoimmune (cytopenia and vasculitis) as well as allergic disorders (eczema, asthma, and food allergy). Hypomorphic function and somatic reversion of DOCK8 has been reported to result in incomplete phenotype without IgE overproduction. Here we describe a case of DOCK8 deficiency in a 8-year-old Caucasian girl. The patient’s disease was initially classified as autoimmune thrombocytopenia, which then evolved toward a combined immunodeficiency phenotype with recurrent infections, persistent EBV infection and lymphoproliferation. Two novel variants (one deletion and one premature stop codon) were characterized, resulting in markedly reduced, but not absent, DOCK8 expression. Somatic reversion of the DOCK8 deletion was identified in T cells. Hypomorphic function and somatic reversion were associated with restricted T cell repertoire, decreased STAT5 phosphorylation and impaired immune synapse functioning in T cells. Although the patient presented with incomplete phenotype (absence of markedly increase IgE and eosinophil count), sclerosing cholangitis was incidentally detected, thus indicating that hypomorphic function and somatic reversion of DOCK8 may delay disease progression but do not necessarily prevent from severe complications.

Recently we showed that Ethanolamine Kinase 1 (ETNK1) mutations cause a decreased synthesis of phosphoethanolamine, and that phosphoethanolamine is able to modulate mitochondrial activity through competition with succinate for complex II. The decreased phosphoethanolamine concentration leads to increased mitochondria activity and reactive oxygen species production, which causes the accumulation of new mutations.

Summary: To dissect the mechanisms underlying the inflation of variants in the Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2) genome, we present a large-scale analysis of intra-host genomic diversity, which reveals that most samples exhibit heterogeneous genomic architectures, due to the interplay between host-related mutational processes and transmission dynamics. The decomposition of minor variants profiles unveils three non-overlapping mutational signatures related to nucleotide substitutions and likely ruled by APOlipoprotein B Editing Complex (APOBEC), Reactive Oxygen Species (ROS), and Adenosine Deaminase Acting on RNA (ADAR), highlighting heterogeneous host responses to SARS-CoV-2 infections. A corrected-for-signatures dN/dS analysis demonstrates that such mutational processes are affected by purifying selection, with important exceptions. In fact, several mutations appear to transit toward clonality, defining new clonal genotypes that increase the overall genomic diversity. Furthermore, the phylogenomic analysis shows the presence of homoplasies and supports the hypothesis of transmission of minor variants. This study paves the way for the integrated analysis of intra-host genomic diversity and clinical outcomes of SARS-CoV-2 infections.

Abstract. Atypical chronic myeloid leukemia (aCML) is a BCR-ABL1-negative clonal disorder, which belongs to the myelodysplastic/myeloproliferative group. This disease is characterized by recurrent somatic mutations in SETBP1, ASXL1 and ETNK1 genes, as well as high genetic heterogeneity, thus posing a great therapeutic challenge. To provide a comprehensive genomic characterization of aCML we applied a high-throughput sequencing strategy to 43 aCML samples, including both whole-exome and RNA-sequencing data. Our dataset identifies ASXL1, SETBP1, and ETNK1 as the most frequently mutated genes with a total of 43.2%, 29.7 and 16.2%, respectively. We characterized the clonal architecture of 7 aCML patients by means of colony assays and targeted resequencing. The results indicate that ETNK1 variants occur early in the clonal evolution history of aCML, while SETBP1 mutations often represent a late event. The presence of actionable mutations conferred both ex vivo and in vivo sensitivity to specific inhibitors with evidence of strong in vitro synergism in case of multiple targeting. In one patient, a clinical response was obtained. Stratification based on RNA-sequencing identified two different populations in terms of overall survival, and differential gene expression analysis identified 38 significantly overexpressed genes in the worse outcome group. Three genes correctly classified patients for overall survival.

One of the key challenges in current cancer research is the development of computational strategies to support clinicians in the identification of successful personalized treatments. Control theory might be an effective approach to this end, as proven by the long-established application to therapy design and testing. In this respect, we here introduce the Control Theory for Therapy Design (CT4TD) framework, which employs optimal control theory on patient-specific pharmacokinetics (PK) and pharmacodynamics (PD) models, to deliver optimized therapeutic strategies. The definition of personalized PK/PD models allows to explicitly consider the physiological heterogeneity of individuals and to adapt the therapy accordingly, as opposed to standard clinical practices. CT4TD can be used in two distinct scenarios. At the time of the diagnosis, CT4TD allows to set optimized personalized administration strategies, aimed at reaching selected target drug concentrations, while minimizing the costs in terms of toxicity and adverse effects. Moreover, if longitudinal data on patients under treatment are available, our approach allows to adjust the ongoing therapy, by relying on simplified models of cancer population dynamics, with the goal of minimizing or controlling the tumor burden. CT4TD is highly scalable, as it employs the efficient dCRAB/RedCRAB optimization algorithm, and the results are robust, as proven by extensive tests on synthetic data. Furthermore, the theoretical framework is general, and it might be applied to any therapy for which a PK/PD model can be estimated, and for any kind of administration and cost. As a proof of principle, we present the application of CT4TD to Imatinib administration in Chronic Myeloid leukemia, in which we adopt a simplified model of cancer population dynamics. In particular, we show that the optimized therapeutic strategies are diversified among patients, and display improvements with respect to the current standard regime.

SETBP1 variants occur as somatic mutations in several malignancies and as de novo germline mutations in developmental disorders. Here the authors provide evidence that SETBP1 binds to gDNA in AT-rich promoter regions to promote target gene upregulation, indicating SETBP1 functions directly to regulate transcription.

BRAF is the most frequently mutated gene in melanoma. Constitutive activation of mutant BRAFV600E leads to aberrant Ras-independent MAPK signaling and cell transformation. Inhibition of mutant BRAF is a current frontline therapy for such cases, with improved survival compared with chemotherapy. Unfortunately, reactivation of MAPK signaling by several mechanisms has been shown to cause drug resistance and disease recurrence. In this work, we describe the co-occurrence of an in-frame deletion within an amplified BRAFV600E locus and a missense point mutation of the transcriptional repressor BCORL1 in vemurafenib-resistant A375 melanoma cells. Functional data confirmed that truncated p47BRAFV600E and mutant BCORL1Q1076H both contribute to resistance. Interestingly, either endogenous BCORL1 silencing or ectopic BCORL1Q1076H expression mimicked the effects of a CRISPR/Cas9-edited BCORL1Q1076H locus, suggesting a complex mixture of loss- and gain-of-function effects caused by the mutation. Transcriptomic data confirmed this hypothesis. Finally, we show that the pan-RAF inhibitor sorafenib is not affected by expression of BRAF deletion variant and effectively synergizes with vemurafenib to block resistant cells, suggesting a possible intervention for this class of mutants.

Despite the advent of tyrosine kinase inhibitors, a proportion of chronic myeloid leukemia patients in chronic phase fail to respond to imatinib or to second-generation inhibitors and progress to blast crisis. Until now, improvements in the understanding of the molecular mechanisms responsible for chronic myeloid leukemia transformation from chronic phase to the aggressive blast crisis remain limited. Here we present a large parallel sequencing analysis of 10 blast crisis samples and of the corresponding autologous chronic phase controls that reveals, for the first time, recurrent mutations affecting the ubiquitin-conjugating enzyme E2A gene (UBE2A, formerly RAD6A). Additional analyses on a cohort of 24 blast crisis, 41 chronic phase as well as 40 acute myeloid leukemia and 38 atypical chronic myeloid leukemia patients at onset confirmed that UBE2A mutations are specifically acquired during chronic myeloid leukemia progression, with a frequency of 16.7% in advanced phases. In vitro studies show that the mutations here described cause a decrease in UBE2A activity, leading to an impairment of myeloid differentiation in chronic myeloid leukemia cells.

The presence of thousands of repetitive sequences makes the centromere a fragile region subject to breakage. In this study we collected 31 cases of rearrangements of chromosome 18, of which 16 involved an acrocentric chromosome, during genetic screening done in three centers. We noticed a significant enrichment of reciprocal translocations between the centromere of chromosome 18 and the centromeric or pericentromeric regions of the acrocentrics. We describe five cases with translocation between chromosome 18 and an acrocentric chromosome, and one case involving the common telomere regions of chromosomes 18p and 22p. In addition, we bring evidence to support the hypothesis that chromosome 18 preferentially recombines with acrocentrics: (i) the presence on 18p11.21 of segmental duplications highly homologous to acrocentrics, that can justify a NAHR mechanism; (ii) the observation by 2D-FISH of the behavior of the centromeric regions of 18 respect to the centromeric regions of acrocentrics in the nuclei of normal subjects; (iii) the contact analysis among these regions on published Hi-C data from the human lymphoblastoid cell line (GM12878).

HFE-related hereditary hemochromatosis (HH) is characterized by marked phenotypic heterogeneity. Homozygosity for p.C282Y is a low penetrance genotype suggesting that the HFE-HH is a multifactorial disease resulting from a complex interaction involving a major gene defect, genetic background and environmental factors. We performed a targeted NGS-based gene panel to identify new candidate modifiers by using an extreme phenotype sampling study based on serum ferritin and iron removed/age ratio. We found an increased prevalence of the HIF1A p.Phe582Ser and p.Ala588Thr variants in patients with a severe iron and clinical phenotype. Accordingly, Huh-7 cells transfected with both variants showed significantly lower HAMP promoter activity by luciferase assay. The qRT-PCR assays showed a downregulation of hepcidin and an upregulation of the HIF1A target genes (VEGF, HMOX, FUR, TMPRSS6) in cells transfected with the HIF1A-P582S vector. We identified mutations in other genes (e.g., Serpina1) that might have some relevance in single cases in aggravating or mitigating disease manifestation. In conclusion, the present study identified HIF1A as a possible modifier of the HFE-HH phenotype cooperating with the genetic defect in downregulating hepcidin synthesis. In addition, this study highlights that an NGS-based approach could broaden our knowledge and help in characterizing the genetic complexity of HFE-HH patients with a severe phenotype expression.

BIM is a proapoptotic member of the Bcl-2 family. Here, we investigated the epigenetic status of the BIM locus in NPM/ALK+ anaplastic large cell lymphoma (ALCL) cell lines and in lymph node biopsies from NPM/ALK+ ALCL patients. We show that BIM is epigenetically silenced in cell lines and lymph node specimens and that treatment with the deacetylase inhibitor trichostatin A restores the histone acetylation, strongly upregulates BIM expression, and induces cell death. BIM silencing occurs through recruitment of MeCP2 and the SIN3a/histone deacetylase 1/2 (HDAC1/2) corepressor complex. This event requires BIM CpG methylation/demethylation with 5-azacytidine that leads to detachment of the MeCP2 corepressor complex and reacetylation of the histone tails. Treatment with the ALK inhibitor PF2341066 or with an inducible shRNA targeting NPM/ALK does not restore BIM locus reacetylation; however, enforced expression of NPM/ALK in an NPM/ALK-negative cell line significantly increases the methylation at the BIM locus. This study demonstrates that BIM is epigenetically silenced in NPM/ALK-positive cells through recruitment of the SIN3a/HDAC1/2 corepressor complex and that NPM/ALK is dispensable to maintain BIM epigenetic silencing but is able to act as an inducer of BIM methylation.

Copy number alterations (CNA) are common events occurring in leukaemias and solid tumors. Comparative Genome Hybridization (CGH) is actually the gold standard technique to analyze CNAs; however, CGH analysis requires dedicated instruments and is able to perform only low resolution Loss of Heterozygosity (LOH) analyses. Here we present CEQer (Comparative Exome Quantification analyzer), a new graphical, event-driven tool for CNA/allelic-imbalance (AI) coupled analysis of exome sequencing data. By using case-control matched exome data, CEQer performs a comparative digital exonic quantification to generate CNA data and couples this information with exome-wide LOH and allelic imbalance detection. This data is used to build mixed statistical/heuristic models allowing the identification of CNA/AI events. To test our tool, we initially used in silico generated data, then we performed whole-exome sequencing from 20 leukemic specimens and corresponding matched controls and we analyzed the results using CEQer. Taken globally, these analyses showed that the combined use of comparative digital exon quantification and LOH/AI allows generating very accurate CNA data. Therefore, we propose CEQer as an efficient, robust and user-friendly graphical tool for the identification of CNA/AI in the context of whole-exome sequencing data.

The ERG gene belongs to the ETS family of transcription factors and has been found to be involved in atypical chromosomal rearrangements in several cancers. To gain insight into the oncogenic activity of ERG, we compared the gene expression profile of NIH-3T3 cells stably expressing the coding regions of the three main ERG oncogenic fusions: TMPRSS2/ERG (tERG), EWS/ERG and FUS/ERG. We found that all three ERG fusions significantly up-regulate PIM1 expression in the NIH-3T3 cell line. PIM1 is a serine/threonine kinase frequently over-expressed in cancers of haematological and epithelial origin. We show here that tERG expression induces PIM1 in the non-malignant prostate cell line RWPE-1, strengthening the relation between tERG and PIM1 up-regulation in the initial stages of prostate carcinogenesis. Silencing of tERG reversed PIM1 induction. A significant association between ERG and PIM1 expression in clinical prostate carcinoma specimens was found, suggesting that such a mechanism may be relevant in vivo. Chromatin Immunoprecipitation experiments showed that tERG directly binds to PIM1 promoter in the RWPE-1 prostate cell line, suggesting that tERG could be a direct regulator of PIM1 expression. The up-regulation of PIM1 induced by tERG over-expression significantly modified Cyclin B1 levels and increased the percentage of aneuploid cells in the RWPE-1 cell line after taxane-based treatment. Here we provide the first evidence for an ERG-mediated PIM1 up-regulation in prostate cells in vitro and in vivo, suggesting a direct effect of ERG transcriptional activity in the alteration of genetic stability.