Multiplexed Imaging Mass Cytometry Analysis in Preclinical Models of Pancreatic Cancer.

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal cancers. PDAC is characterized by a complex tumor microenvironment (TME), that plays a pivotal role in disease progression and resistance to therapy. Investigating the spatial distribution and interaction of TME cells with the tumor is the basis for understanding the mechanisms underlying disease progression and represents a current challenge in PDAC research. Imaging mass cytometry (IMC) is the major multiplex imaging technology for the spatial analysis of tumor heterogeneity. However, there is a dearth of reports of multiplexed IMC panels for different preclinical mouse models, including pancreatic cancer. We addressed this gap by utilizing two preclinical models of PDAC: the genetically engineered, bearing KRAS-TP53 mutations in pancreatic cells, and the orthotopic, and developed a 28-marker panel for single-cell IMC analysis to assess the abundance, distribution and phenotypes of cells involved in PDAC progression and their reciprocal functional interactions. Herein, we provide an unprecedented definition of the distribution of TME cells in PDAC and compare the diversity between transplanted and genetic disease models. The results obtained represent an important and customizable tool for unraveling the complexities of PDAC and deciphering the mechanisms behind therapy resistance.

Integrating AI-Powered Digital Pathology and Imaging Mass Cytometry Identifies Key Classifiers of Tumor Cells, Stroma, and Immune Cells in Non-Small Cell Lung Cancer.

Artificial intelligence (AI)-powered approaches are becoming increasingly used as histopathologic tools to extract subvisual features and improve diagnostic workflows. On the other hand, hi-plex approaches are widely adopted to analyze the immune ecosystem in tumor specimens. Here, we aimed at combining AI-aided histopathology and imaging mass cytometry (IMC) to analyze the ecosystem of non-small cell lung cancer (NSCLC). An AI-based approach was used on hematoxylin and eosin (H&E) sections from 158 NSCLC specimens to accurately identify tumor cells, both adenocarcinoma and squamous carcinoma cells, and to generate a classifier of tumor cell spatial clustering. Consecutive tissue sections were stained with metal-labeled antibodies and processed through the IMC workflow, allowing quantitative detection of 24 markers related to tumor cells, tissue architecture, CD45+ myeloid and lymphoid cells, and immune activation. IMC identified 11 macrophage clusters that mainly localized in the stroma, except for S100A8+ cells, which infiltrated tumor nests. T cells were preferentially localized in peritumor areas or in tumor nests, the latter being associated with better prognosis, and they were more abundant in highly clustered tumors. Integrated tumor and immune classifiers were validated as prognostic on whole slides. In conclusion, integration of AI-powered H&E and multiparametric IMC allows investigation of spatial patterns and reveals tissue relevant features with clinical relevance. SIGNIFICANCE: Leveraging artificial intelligence-powered H&E analysis integrated with hi-plex imaging mass cytometry provides insights into the tumor ecosystem and can translate tumor features into classifiers to predict prognosis, genotype, and therapy response.

A Multilayered Imaging and Microfluidics Approach for Evaluating the Effect of Fibrinolysis in Staphylococcus aureus Biofilm Formation.

The recognition of microbe and extracellular matrix (ECM) is a recurring theme in the humoral innate immune system. Fluid-phase molecules of innate immunity share regulatory roles in ECM. On the other hand, ECM elements have immunological functions. Innate immunity is evolutionary and functionally connected to hemostasis. Staphylococcus aureus (S. aureus) is a major cause of hospital-associated bloodstream infections and the most common cause of several life-threatening conditions such as endocarditis and sepsis through its ability to manipulate hemostasis. Biofilm-related infection and sepsis represent a medical need due to the lack of treatments and the high resistance to antibiotics. We designed a method combining imaging and microfluidics to dissect the role of elements of the ECM and hemostasis in triggering S. aureus biofilm by highlighting an essential role of fibrinogen (FG) in adhesion and formation. Furthermore, we ascertained an important role of the fluid-phase activation of fibrinolysis in inhibiting biofilm of S. aureus and facilitating an antibody-mediated response aimed at pathogen killing. The results define FG as an essential element of hemostasis in the S. aureus biofilm formation and a role of fibrinolysis in its inhibition, while promoting an antibody-mediated response. Understanding host molecular mechanisms influencing biofilm formation and degradation is instrumental for the development of new combined therapeutic approaches to prevent the risk of S. aureus biofilm-associated diseases.

Interrogating colorectal cancer metastasis to liver: a search for clinically viable compounds and mechanistic insights in colorectal cancer Patient Derived Organoids.

BACKGROUND: Approximately 20-50% of patients presenting with localized colorectal cancer progress to stage IV metastatic disease (mCRC) following initial treatment and this is a major prognostic determinant. Here, we have interrogated a heterogeneous set of primary colorectal cancer (CRC), liver CRC metastases and adjacent liver tissue to identify molecular determinants of the colon to liver spreading. Screening Food and Drug Administration (FDA) approved drugs for their ability to interfere with an identified colon to liver metastasis signature may help filling an unmet therapeutic need. METHODS: RNA sequencing of primary colorectal cancer specimens vs adjacent liver tissue vs synchronous and asynchronous liver metastases. Pathways enrichment analyses. The Library of Integrated Network-based Cellular Signatures (LINCS)-based and Connectivity Map (CMAP)-mediated identification of FDA-approved compounds capable to interfere with a 22 gene signature from primary CRC and liver metastases. Testing the identified compounds on CRC-Patient Derived Organoid (PDO) cultures. Microscopy and Fluorescence Activated Cell Sorting (FACS) based analysis of the treated PDOs. RESULTS: We have found that liver metastases acquire features of the adjacent liver tissue while partially losing those of the primary tumors they derived from. We have identified a 22-gene signature differentially expressed among primary tumors and metastases and validated in public databases. A pharmacogenomic screening for FDA-approved compounds capable of interfering with this signature has been performed. We have validated some of the identified representative compounds in CRC-Patient Derived Organoid cultures (PDOs) and found that pentoxyfilline and, to a minor extent, dexketoprofen and desloratadine, can variably interfere with number, size and viability of the CRC -PDOs in a patient-specific way. We explored the pentoxifylline mechanism of action and found that pentoxifylline treatment attenuated the 5-FU elicited increase of ALDHhigh cells by attenuating the IL-6 mediated STAT3 (tyr705) phosphorylation. CONCLUSIONS: Pentoxifylline synergizes with 5-Fluorouracil (5-FU) in attenuating organoid formation. It does so by interfering with an IL-6-STAT3 axis leading to the emergence of chemoresistant ALDHhigh cell subpopulations in 5-FU treated PDOs. A larger cohort of CRC-PDOs will be required to validate and expand on the findings of this proof-of-concept study.

Quantitative analysis of disease-related metabolic dysregulation of human microbiota.

The metabolic activity of all the micro-organism composing the human microbiome interacts with the host metabolism contributing to human health and disease in a way that is not fully understood. Here, we introduce STELLA, a computational method to derive the spectrum of metabolites associated with the microbiome of an individual. STELLA integrates known information on metabolic pathways associated with each bacterial species and extracts from these the list of metabolic products of each singular reaction by means of automatic text analysis. By comparing the result obtained on a single subject with the metabolic profile data of a control set of healthy subjects, we are able to identify individual metabolic alterations. To illustrate the method, we present applications to autism spectrum disorder and multiple sclerosis.

Nuclear biophysical changes during human melanoma plasticity.

Tumor plasticity is an emerging property of tumor cells which allows them to change their phenotype in dependence on the environment. The epithelial-mesenchymal transition plays a crucial role in helping cells to acquire a more aggressive phenotype when they are in the mesenchymal state. Herein we investigated the biophysical changes occurring during phenotypic switching in human melanoma cells considering the blebbines of the nuclei, their stiffness and the involvement of polycombs with lamins. We show that the formation of cellular heterogeneity involves many crucial nuclear changes including the interaction between different types of polycombs with lamins and chromosome accessibility. All together our results shed new light on the molecular mechanisms involved in the formation of an heterogeneous cell population during phenotypic switching. In particular, our results show that phenotypic switching in melanoma involves chromatin remodeling changing the transcriptional activity of cells and consequently their phenotype.

TRPM2 Oxidation Activates Two Distinct Potassium Channels in Melanoma Cells through Intracellular Calcium Increase.

Tumor microenvironments are often characterized by an increase in oxidative stress levels. We studied the response to oxidative stimulation in human primary (IGR39) or metastatic (IGR37) cell lines obtained from the same patient, performing patch-clamp recordings, intracellular calcium ([Ca2+]i) imaging, and RT-qPCR gene expression analysis. In IGR39 cells, chloramine-T (Chl-T) activated large K+ currents (KROS) that were partially sensitive to tetraethylammonium (TEA). A large fraction of KROS was inhibited by paxilline-a specific inhibitor of large-conductance Ca2+-activated BK channels. The TEA-insensitive component was inhibited by senicapoc-a specific inhibitor of the Ca2+-activated KCa3.1 channel. Both BK and KCa3.1 activation were mediated by an increase in [Ca2+]i induced by Chl-T. Both KROS and [Ca2+]i increase were inhibited by ACA and clotrimazole-two different inhibitors of the calcium-permeable TRPM2 channel. Surprisingly, IGR37 cells did not exhibit current increase upon the application of Chl-T. Expression analysis confirmed that the genes encoding BK, KCa3.1, and TRPM2 are much more expressed in IGR39 than in IGR37. The potassium currents and [Ca2+]i increase observed in response to the oxidizing agent strongly suggest that these three molecular entities play a major role in the progression of melanoma. Pharmacological targeting of either of these ion channels could be a new strategy to reduce the metastatic potential of melanoma cells, and could complement classical radio- or chemotherapeutic treatments.

Synergistic effects of contaminants in Lombardy waters.

Quantifying synergistic environmental effects in water contamination is still an open issue. Here, we have analyzed geolocalized data of pollutants recorded in 2018 in surface and groundwater of Lombardy, one of the areas with the highest agricultural production rates, not only in Italy, but also in Europe. Both herbicides and insecticides are present at concentration levels above the legal limit, mainly in surface waters. Geolocalized analysis allows us to identify interesting areas particularly affected by a combination of multiple pesticides. We thus investigated possible synergistic effects of these compounds on the environment, using the alga C. reinhardtii as a biosensor. Our results show that exposure for 7 days to four compounds, that we found present together at high concentration in surface waters, was able to induce a stress in the algae, as indicated by the presence of palmelloids. Our work results in a pipeline that could easily be exported to monitor other territories in Italy and abroad.

Role of body temperature variations in bat immune response to viral infections.

The ability of bats to coexist with viruses without being harmed is an interesting issue that is still under investigation. Here we use a mathematical model to show that the pattern of body temperature variations observed in bats between day and night is responsible for their ability to keep viruses in check. From the dynamical systems point of view, our model displays an intriguing quasi-periodic behaviour that might be relevant in making the system robust by avoiding viral escape due to perturbations in the body temperature cycle.

A Model-Driven Quantitative Analysis of Retrotransposon Distributions in the Human Genome.

Retrotransposons, DNA sequences capable of creating copies of themselves, compose about half of the human genome and played a central role in the evolution of mammals. Their current position in the host genome is the result of the retrotranscription process and of the following host genome evolution. We apply a model from statistical physics to show that the genomic distribution of the two most populated classes of retrotransposons in human deviates from random placement, and that this deviation increases with time. The time dependence suggests a major role of the host genome dynamics in shaping the current retrotransposon distributions. Focusing on a neutral scenario, we show that a simple model based on random placement followed by genome expansion and sequence duplications can reproduce the empirical retrotransposon distributions, even though more complex and possibly selective mechanisms can have contributed. Besides the inherent interest in understanding the origin of current retrotransposon distributions, this work sets a general analytical framework to analyze quantitatively the effects of genome evolutionary dynamics on the distribution of genomic elements.

MicroRNA-222 Regulates Melanoma Plasticity.

Melanoma is one of the most aggressive and highly resistant tumors. Cell plasticity in melanoma is one of the main culprits behind its metastatic capabilities. The detailed molecular mechanisms controlling melanoma plasticity are still not completely understood. Here we combine mathematical models of phenotypic switching with experiments on IgR39 human melanoma cells to identify possible key targets to impair phenotypic switching. Our mathematical model shows that a cancer stem cell subpopulation within the tumor prevents phenotypic switching of the other cancer cells. Experiments reveal that hsa-mir-222 is a key factor enabling this process. Our results shed new light on melanoma plasticity, providing a potential target and guidance for therapeutic studies.

Comparative analysis of metabolic and transcriptomic features of Nothobranchius furzeri.

Some species have a longer lifespan than others, but usually lifespan is correlated with typical body weight. Here, we study the lifetime evolution of the metabolic behaviour of Nothobranchius furzeri, a killifish with an extremely short lifespan with respect to other fishes, even when taking into account rescaling by body weight. Comparison of the gene expression patterns of N. furzeri with those of zebrafish Danio rerio and mouse (Mus musculus) shows that a broad set of metabolic genes and pathways are affected in N. furzeri during ageing in a way that is consistent with a global deregulation of chromatin. Computational analysis of the glycolysis pathway for the three species highlights a rapid increase in the metabolic activity during the lifetime of N. furzeri with respect to the other species. Our results highlight that the unusually short lifespan of N. furzeri is associated with peculiar patterns in the metabolic activities and in chromatin dynamics.

Unjamming of active rotators.

Active particle assemblies can exhibit a wide range of interesting dynamical phases depending on internal parameters such as density, adhesion strength or self-propulsion. Active self-rotations are rarely studied in this context, although they can be relevant for active matter systems, as we illustrate by analyzing the motion of Chlamydomonas reinhardtii algae under different experimental conditions. Inspired by this example, we simulate the dynamics of a system of interacting active disks endowed with active torques and self-propulsive forces. At low packing fractions, adhesion causes the formation of small rotating clusters, resembling those observed when algae are stressed. At higher densities, the model shows a jamming to unjamming transition promoted by active torques and hindered by adhesion. We also study the interplay between self-propulsion and self-rotation and derive a phase diagram. Our results yield a comprehensive picture of the dynamics of active rotators, providing useful guidance to interpret experimental results in cellular systems where rotations might play a role.

Chromatin and Cytoskeletal Tethering Determine Nuclear Morphology in Progerin-Expressing Cells.

The nuclear morphology of eukaryotic cells is determined by the interplay between the lamina forming the nuclear skeleton, the chromatin inside the nucleus, and the coupling with the cytoskeleton. Nuclear alterations are often associated with pathological conditions as in Hutchinson-Gilford progeria syndrome, in which a mutation in the lamin A gene yields an altered form of the protein, named progerin, and an aberrant nuclear shape. Here, we introduce an inducible cellular model of Hutchinson-Gilford progeria syndrome in HeLa cells in which increased progerin expression leads to alterations in the coupling of the lamin shell with cytoskeletal or chromatin tethers as well as with polycomb group proteins. Furthermore, our experiments show that progerin expression leads to enhanced nuclear shape fluctuations in response to cytoskeletal activity. To interpret the experimental results, we introduce a computational model of the cell nucleus that explicitly includes chromatin fibers, the nuclear shell, and coupling with the cytoskeleton. The model allows us to investigate how the geometrical organization of the chromatin-lamin tether affects nuclear morphology and shape fluctuations. In sum, our findings highlight the crucial role played by lamin-chromatin and lamin-cytoskeletal alterations in determining nuclear shape morphology and in affecting cellular functions and gene regulation.

Cancer stem cells, plasticity, and drug resistance.

Melanoma is a highly aggressive tumor and almost always fatal when metastatic. Herein, we discuss recent findings on the mechanisms of resistance of human cutaneous melanoma. To achieve a precision medicine approach, the heterogeneity and plasticity of tumor cells are two crucial aspects to be investigated in depth. In fact, to understand the mechanisms that cells use to acquire a resistant phenotype after chemotherapy or how resistant cells inside the tumor are selected, it is the most important issue for a successful therapy. Since new therapeutic strategies are trying to go in this direction, we discuss here the state of the art of the research and the clinical impact of these strategies. We also discuss and suggest future research directions to develop approaches able to define the best concentration and time of exposure of the drug or the cocktails of drugs for each specific patient based on his/her biological features.

Cross-Talk Between circRNAs and mRNAs Modulates MiRNA-mediated Circuits and Affects Melanoma Plasticity.

CircularRNAs (circRNAs) are non-coding RNAs which compete for microRNA (miRNA) binding, influencing the abundance and stability of other RNA species. Herein we have investigated the effect of circRNAs on the mir200-ZEB1 feedback loop in relationship with the aggressiveness of human melanoma cells. We first compared the level of expression of key factors in the mir200-ZEB1 feedback loop in primary human melanoma cells compared with their matching metastatic one and found a correlation between the aggressiveness of the cells and the level of expression of ZEB1 and SNAI1. We also analyzed factors in the mir200-ZEB1 feedback loop, including circZEB1, during the phenotypic switching of human melanoma cells. Our results showed a correlation between the level of ZEB1 and SNAI1 and the fraction of cancer stem cells in the population. The level of circZEB1 was, however, consistently high during the entire phenotypic transformation. To understand this result we propose a mathematical model of the regulatory circuit. According to the model, the experimental observations can be explained by the presence of a back-splicing factor limiting circRNA production.

Molecular mechanisms of heterogeneous oligomerization of huntingtin proteins.

There is still no successful strategy to treat Huntington's disease, an inherited autosomal disorder associated with the aggregation of mutated forms of the huntingtin protein containing polyglutamine tracts with more than 36 repeats. Recent experimental evidence is challenging the conventional view of the disease by revealing transcellular transfer of mutated huntingtin proteins which are able to seed oligomers involving wild type forms of the protein. Here we decipher the molecular mechanism of this unconventional heterogeneous oligomerization by performing discrete molecular dynamics simulations. We identify the most probable oligomer conformations and the molecular regions that can be targeted to destabilize them. Our computational findings are complemented experimentally by fluorescence-lifetime imaging microscopy/fluorescence resonance energy transfer (FLIM-FRET) of cells co-transfected with huntingtin proteins containing short and large polyglutamine tracts. Our work clarifies the structural features responsible for heterogeneous huntingtin aggregation with possible implications to contrast the prion-like spreading of Huntington's disease.

Sulforaphane Cannot Protect Human Fibroblasts From Repeated, Short and Sublethal Treatments with Hydrogen Peroxide.

A delicate balance of reactive oxygen species (ROS) exists inside the cell: when the mechanisms that control the level of ROS fail, the cell is in an oxidative stress state, a condition that can accelerate aging processes. To contrast the pro-aging effect of ROS, the supplementation of antioxidants has been recently proposed. Sulforaphane (SFN) is an isothiocyanate isolated from Brassica plants that has been shown to modulate many critical factors inside the cells helping to counteract aging processes. In the present work, we exposed human dermal fibroblast to short, sublethal and repeated treatments with hydrogen peroxide for eight days, without or in combination with low concentration of SFN. Hydrogen peroxide treatments did not affect the oxidative status of the cells, without any significant change of the intracellular ROS levels or the number of mitochondria or thiols in total proteins. However, our regime promoted cell cycle progression and cell viability, increased the anti-apoptotic factor survivin and increased DNA damage, measured as number of foci positive for γ -H2AX. On the other hand, the treatment with SFN alone seemed to exert a protective effect, increasing the level of p53, which can block the expansion of possible DNA damaged cells. However, continued exposure to SFN at this concentration could not protect the cells from stress induced by hydrogen peroxide.

Impact of the cross-talk between circular and messenger RNAs on cell regulation.

Circular RNAs (circRNA) are non-coding RNAs characterized by a closed-loop structure providing increased stability and enhancing their ability to compete with messenger RNAs (mRNA) for microRNA (miRNA) binding. An important open question is if circRNA plays a physiological role in the cells regulating critical cellular functions. We address this question with a combination of theoretical modeling and experimental analysis. Theoretically, we investigated two possible scenarios, one in which circRNAs act as sponges for miRNAs but there is no other relation between the two RNAs and the other one where circRNAs are co-generated with their corresponding mRNAs. We thus compared the results with data reporting the level of circRNAs (ZEB1, CANCX, ABCC1) modulating specific miRNAs in 27 cell lines. Due to the high stability of circRNAs, we show that a miRNA mediated cross-talk between circRNA and mRNA appears for a broad range of physiological parameters. This confirms the relevance of circRNAs in cell regulation, suggesting that they could be used as biomarkers.

Regeneration in distantly related species: common strategies and pathways.

While almost all animals are able to at least partially replace some lost parts, regeneration abilities vary considerably across species. Here we study gene expression patterns in distantly related species to investigate conserved regeneration strategies. To this end, we collect from the literature transcriptomic data obtained during the regeneration of three species (Hydra magnipapillata, Schmidtea mediterranea, and Apostichopus japonicus), and compare them with gene expression during regeneration in vertebrates and mammals. This allows us to identify a common set of differentially expressed genes and relevant shared pathways that are conserved across species during the early stage of the regeneration process. We also find a set of differentially expressed genes that in mammals are associated to the presence of macrophages and to the epithelial-mesenchymal transition. This suggests that features of the sophisticated wound healing strategy of mammals are already observable in earlier emerging metazoans.