Development of machine learning methods for multi-modal biomarkers detection and integration

In medicine, innovation depends on a better knowledge of the human body mechanism, which represents a complex system of multi-scale constituents. Unraveling the complexity underneath diseases proves to be challenging. A deep understanding of the inner workings comes with dealing with many heterogeneous information. Exploring the molecular status and the organization of genes, proteins, metabolites provides insights on what is driving a disease, from aggressiveness to curability. Molecular constituents, however, are only the building blocks of the human body and cannot currently tell the whole story of diseases. This is why nowadays attention is growing towards the contemporary exploitation of multi-scale information. Holistic methods are then drawing interest to address the problem of integrating heterogeneous data. The heterogeneity may derive from the diversity across data types and from the diversity within diseases. Here, four studies conducted data integration using customly designed workflows that implement novel methods and views to tackle the heterogeneous characterization of diseases. The first study devoted to determine shared gene regulatory signatures for onco-hematology and it showed partial co-regulation across blood-related diseases. The second study focused on Acute Myeloid Leukemia and refined the unsupervised integration of genomic alterations, which turned out to better resemble clinical practice. In the third study, network integration for artherosclerosis demonstrated, as a proof of concept, the impact of network intelligibility when it comes to model heterogeneous data, which showed to accelerate the identification of new potential pharmaceutical targets. Lastly, the fourth study introduced a new method to integrate multiple data types in a unique latent heterogeneous-representation that facilitated the selection of important data types to predict the tumour stage of invasive ductal carcinoma. The results of these four studies laid the groundwork to ease the detection of new biomarkers ultimately beneficial to medical practice and to the ever-growing field of Personalized Medicine.

The role of interleuchin 6 (IL-6) in the promotion of an aggressive and stem cell-like phenotype in human breast cancer cells and in stem/progenitor cells expanded in vitro as mammospheres

High serum levels of Interleukin-6 (IL-6) correlate with poor outcome in breast cancer patients. However no data are available on the relationship between IL-6 and stem/progenitor cells which may fuel the genesis of breast cancer in vivo. Herein, we address this issue in mammospheres (MS), multi-cellular structures enriched in stem/progenitor cells of the mammary gland, and also in MCF-7 breast cancer cells. We show that MS from node invasive breast carcinoma tissues express IL-6 mRNA at higher levels than MS from matched non-neoplastic mammary glands. We find that IL-6 mRNA is detectable only in basal-like breast carcinoma tissues, an aggressive variant showing stem cell features. Our results reveal that IL-6 triggers a Notch-3-dependent up-regulation of the Notch ligand Jagged-1, whose interaction with Notch-3 promotes the growth of MS and MCF-7 derived spheroids. Moreover, IL-6 induces a Notch-3-dependent up-regulation of the carbonic anhydrase IX gene, which promotes a hypoxia-resistant/invasive phenotype in MCF-7 cells and MS. Finally, an autocrine IL-6 loop relies upon Notch-3 activity to sustain the aggressive features of MCF-7-derived hypoxia-selected cells. In conclusion, our data support the hypothesis that IL-6 induces malignant features in Notch-3 expressing, stem/progenitor cells from human ductal breast carcinoma and normal mammary gland.

Molecular analysis of special type breast carcinomas

The project was developed into three parts: the analysis of p63 isoform in breast tumours; the study of intra-tumour eterogeneicity in metaplastic breast carcinoma; the analysis of oncocytic breast carcinoma. p63 is a sequence-specific DNA-binding factor, homologue of the tumour suppressor and transcription factor p53. The human p63 gene is composed of 15 exons and transcription can occur from two distinct promoters: the transactivating isoforms (TAp63) are generated by a promoter upstream of exon 1, while the alternative promoter located in intron 3 leads to the expression of N-terminal truncated isoforms (ΔNp63). It has been demonstrated that anti-p63 antibodies decorate the majority of squamous cell carcinomas of different organs; moreover tumours with myoepithelial differentiation of the breast show nuclear p63 expression. Two new isoforms have been described with the same sequence as TAp63 and ΔNp63 but lacking exon 4: d4TAp63 and ΔNp73L, respectively. Purpose of the study was to investigate the molecular expression of N-terminal p63 isoforms in benign and malignant breast tissues. In the present study 40 specimens from normal breast, benign lesions, DIN/DCIS, and invasive carcinomas were analyzed by immunohistochemistry and RT-PCR (Reverse Transcriptase-PCR) in order to disclose the patterns of p63 expression. We have observed that the full-length isoforms can be detected in non neoplastic and neoplastic lesions, while the short isoforms are only present in the neoplastic cells of invasive carcinomas. Metaplastic carcinomas of the breast are a heterogeneous group of neoplasms which exhibit varied patterns of metaplasia and differentiation. The existence of such non-modal populations harbouring distinct genetic aberrations may explain the phenotypic diversity observed within a given tumour. Intra-tumour morphological heterogeneity is not uncommon in breast cancer and it can often be appreciated in metaplastic breast carcinomas. Aim of this study was to determine the existence of intra-tumour genetic heterogeneity in metaplastic breast cancers and whether areas with distinct morphological features in a given tumour might be underpinned by distinct patterns of genetic aberrations. 47 cases of metaplastic breast carcinomas were retrieved. Out of the 47 cases, 9 had areas that were of sufficient dimensions to be independently microdissected. Our results indicate that at least some breast cancers are composed of multiple non-modal populations of clonally related cells and provide direct evidence that at least some types of metaplastic breast cancers are composed of multiple non-modal clones harbouring distinct genetic aberrations. Oncocytic tumours represent a distinctive set of lesions with typical granular cytoplasmatic eosinophilia of the neoplastic cells. Only rare example of breast oncocytic carcinomas have been reported in literature and the incidence is probably underestimated. In this study we have analysed 33 cases of oncocytic invasive breast carcinoma of the breast, selected according to morphological and immunohistochemical criteria. These tumours were morphologically classified and studied by immunohistochemistry and aCGH. We have concluded that oncocytic breast carcinoma is a morphologic entity with distinctive ultrastructural and histological features; immunohistochemically is characterized by a luminal profile, it has a frequency of 19.8%, has not distinctive clinical features and, at molecular level, shows a specific constellation of genetic aberration.

Stem cell pathways in the basal-like breast carcinoma: the role of beta-catenin and HIF-1alpha

Basal-like tumor is an aggressive breast carcinoma subtype that displays an expression signature similar to that of the basal/myoepithelial cells of the breast tissue. Basal-like carcinoma are characterized by over-expression of the Epidermal Growth Factor receptor (EGFR), high frequency of p53 mutations, cytoplasmic/nuclear localization of beta-catenin, overexpression of the Hypoxia inducible factor (HIF)-1alpha target Carbonic Anhydrase isoenzime 9 (CA9) and a gene expression pattern similar to that of normal and cancer stem cells, including the over-expression of the mammary stem cell markers CD44. In this study we investigated the role of p53, EGFR, beta-catenin and HIF-1alpha in the regulation of stem cell features and genes associated with the basal-like gene expression profile. The findings reported in this investigation indicate that p53 inactivation in ductal breast carcinoma cells leads to increased EGFR mRNA and protein levels. In our experimental model, EGFR overexpression induces beta-catenin cytoplasmatic stabilization and transcriptional activity and, by that, leads to increased aggressive features including mammosphere (MS) forming and growth capacity, invasive potential and overexpression of the mammary stem cell gene CD44. Moreover we found that EGFR/beta-catenin axis promotes hypoxia survival in breast carcinoma cells via increased CA9 expression. Indeed beta-catenin positively regulates CA9 expression upon hypoxia exposure. Interestingly we found that beta-catenin inhibits HIF-1alpha transcriptional activity. Looking for the mechanism, we found that CA9 expression is promoted by HIF-1alpha and cytoplasmatic beta-catenin further increased it post-transcriptionally, via direct mRNA binding and stabilization. These data reveal a functional beta-catenin/HIF-1alpha interplay among hallmarks of basal-like tumors and unveil a new functional role for cytoplasmic beta-catenin in the phenotype of such tumors. Therefore it can be proposed that the interplay here described among EGFR/beta-catenin and HIF-1alpha may play a role in breast cancer stem cell survival and function.

Trasduzione del segnale inositide-dipendente nel carcinoma della mammella

Breast carcinoma, one of the most frequent malignancies in women, is a complex disease in which a number of different factors combine to drive pathogenesis. The biopathological characterization of these tumors is essential to determine their aggressiveness and to find the most appropriate therapy. As in others neoplasms, the deregulation of signal transduction pathways is frequently responsible for conferring selective biological advantages to the tumor. Phosphoinositides play an essential role in diverse cellular functions, their metabolism is highly active, and is tightly controlled. Among the enzymes implicated in this pathway, phospholipase C beta 1 (PLCβ1) is one of the key regulators, both at the cytoplasmic and the nuclear level. The PLCβ1 gene maps onto the short arm of chromosome 20, a region that has been shown to be altered in several solid tumors, including breast cancer. In the present study a FISH approach was used to investigate the genetic alterations of the PLCβ1 gene in various classes of breast cancer which differ in their invasiveness and proliferation status, according to their mitotic index. The overall aim was to find out whether this enzyme could be a suitable prognostic marker for this neoplasm. Our results show that 83% of cases had aneusomies at the 20p12 level, and the most frequent alteration is a gain in this specific locus. Indeed, we found that this amplification is not related to the invasion status since there were no differences in amplified tumor frequencies between in situ and invasive breast cancer. On the contrary, the gain of PLCβ1 was significantly related to the mitotic index (p = 0.001). To verify if the change in genetic dosage influences the expression of PLCβ1 we performed Real Time PCR and Immunohystochemical analysis. Our results confirmed that amplified tumors have higher levels of PLCβ1 mRNA, which is the sum of the two splicing isoforms 1a and 1b. On the other hand, even if protein levels were higher in the majority of cases compared to the nontumoral specimens, there were no significant associations between gain and overexpression. Finally, the significant association between the amplification of PLCβ1 and others important clinicopathological parameters, such as grading and hormonal receptors status, confirmed a correlation of this enzyme with the aggressiveness of breast cancer. This suggests that PLCβ1 has the potential to be a prognostic marker in these tumors. However, further work needs to be carried out to validate these preliminary findings.

Changes in tumor stiffness for early prediction of tumor response to sorafenib: a proof-of-concept study with elastosonography in an animal model of Hepatocellular Carcinoma (HCC)

Background and aims: Sorafenib is the reference therapy for advanced Hepatocellular Carcinoma (HCC). No method exists to predict in the very early period subsequent individual response. Starting from the clinical experience in humans that subcutaneous metastases may rapidly change consistency under sorafenib and that elastosonography a new ultrasound based technique allows assessment of tissue stiffness, we investigated the role of elastonography in the very early prediction of tumor response to sorafenib in a HCC animal model. Methods: HCC (Huh7 cells) subcutaneous xenografting in mice was utilized. Mice were randomized to vehicle or treatment with sorafenib when tumor size was 5-10 mm. Elastosonography (Mylab 70XVG, Esaote, Genova, Italy) of the whole tumor mass on a sagittal plane with a 10 MHz linear transducer was performed at different time points from treatment start (day 0, +2, +4, +7 and +14) until mice were sacrified (day +14), with the operator blind to treatment. In order to overcome variability in absolute elasticity measurement when assessing changes over time, values were expressed in arbitrary units as relative stiffness of the tumor tissue in comparison to the stiffness of a standard reference stand-off pad lying on the skin over the tumor. Results: Sor-treated mice showed a smaller tumor size increase at day +14 in comparison to vehicle-treated (tumor volume increase +192.76% vs +747.56%, p=0.06). Among Sor-treated tumors, 6 mice showed a better response to treatment than the other 4 (increase in volume +177% vs +553%, p=0.011). At day +2, median tumor elasticity increased in Sor-treated group (+6.69%, range –30.17-+58.51%), while decreased in the vehicle group (-3.19%, range –53.32-+37.94%) leading to a significant difference in absolute values (p=0.034). From this time point onward, elasticity decreased in both groups, with similar speed over time, not being statistically different anymore. In Sor-treated mice all 6 best responders at day 14 showed an increase in elasticity at day +2 (ranging from +3.30% to +58.51%) in comparison to baseline, whereas 3 of the 4 poorer responders showed a decrease. Interestingly, these 3 tumours showed elasticity values higher than responder tumours at day 0. Conclusions: Elastosonography appears a promising non-invasive new technique for the early prediction of HCC tumor response to sorafenib. Indeed, we proved that responder tumours are characterized by an early increase in elasticity. The possibility to distinguish a priori between responders and non responders based on the higher elasticity of the latter needs to be validated in ad-hoc experiments as well as a confirmation of our results in humans is warranted.

Profili di espressione genetica ed epigenetica ad impatto prognostico e predittivo nel carcinoma squamoso e nelle lesioni potenzialmente maligne del cavo orale

Il carcinoma squamoso orale (CSO) è spesso preceduto da lesioni definite potenzialmente maligne tra cui la leucoplachia e il lichen ma una diagnosi precoce avviene ancora oggi in meno della metà dei casi. Inoltre spesso un paziente trattato per CSO svilupperà secondi tumori. Scopo del lavoro di ricerca è stato: 1) Studiare, mediante metodica di next generation sequencing, lo stato di metilazione di un gruppo di geni a partire da prelievi brushing del cavo orale al fine di identificare CSO o lesioni ad alto rischio di trasformazione maligna. 2) Valurare la relazione esistente tra sovraespressione di p16INK4A e presenza di HPV in 35 pazienti affetti da lichen 3) Valutare la presenza di marker istopatologici predittivi di comparsa di seconde manifestazioni tumorali 4) valutare la relazione clonale tra tumore primitivo e metastasi linfonodale in 8 pazienti mediante 2 metodiche di clonalità differenti: l’analisi di mtDNA e delle mutazioni del gene TP53. I risultati hanno mostrato: 1) i geni ZAP70 e GP1BB hanno presentato un alterato stato di metilazione rispettivamente nel 100% e nel 90,9% di CSO e lesioni ad alto rischio, mentre non sono risultati metilati nei controlli sani; ipotizzando un ruolo come potenziali marcatori per la diagnosi precoce nel CSO. 2)Una sovraespressione di p16INK4A è risultata in 26/35 pazienti affetti da lichen ma HPV-DNA è stato identificato in soli 4 campioni. Nessuna relazione sembra essere tra sovraespressione di p16INK4A e virus HPV. 3)L’invasione perineurale è risultato un marker predittivo della comparsa di recidiva locale e metastasi linfonodale, mentre lo stato dei margini chirurgici si è rilevato un fattore predittivo per la comparsa di secondi tumori primitivi 4) Un totale accordo nei risultati c’è stato tra analisi di mtDNA e analisi di TP53 e le due metodiche hanno identificato la presenza di 4 metastasi linfonodali non clonalmente correlate al tumore primitivo.

Indicatori patologici di rilevanza diagnostica nel carcinoma squamocellulare orale felino

Il carcinoma squamocellulare è il tumore maligno orale più frequente nel gatto e si caratterizza per diagnosi spesso tardiva e prognosi infausta. Il progetto riguarda la ricerca di marker di rilevanza dia-gnostica nel carcinoma squamocellulare orale felino (FOSCC), al fine di sviluppare un test di scree-ning non invasivo. È stata condotta un’analisi retrospettiva delle disregolazioni del gene oncosoppres-sore TP53 in campioni istologici di FOSCC e di una popolazione di controllo (lesioni infiammatorie croniche orali e mucose orali normali feline). Tramite next-generation sequencing (NGS) sono state rilevate mutazioni di TP53 nel 69% dei FOSCC, ed anche l’espressione immunoistochimica della pro-teina p53 era presente nel 69% dei tumori, con una concordanza discreta (77%) fra le due alterazioni. Nella popolazione di controllo erano presenti disregolazioni di p53 solo in due lesioni infiammatorie (3%). Successivamente è stata effettuata un’analisi prospettica con NGS della metilazione del DNA di 17 geni, noti per essere disregolati nel carcinoma squamocellulare orale umano o felino, insieme all’analisi mutazionale di TP53, in campioni istologici di FOSCC e in un gruppo di controllo. Le stesse indagini molecolari sono state svolte in parallelo su campioni di cellule prelevate mediante brushing orale. Utilizzando 6 dei geni indagati differenzialmente metilati nei FOSCC (FLI1, MiR124-1, KIF1A, MAGEC2, ZAP70, MiR363) e lo stato mutazionale diTP53, è stato impostato un algoritmo diagnostico per differenziare i FOSCC dalla mucosa orale non neoplastica. Applicato ai brushing, l’algoritmo è risultato positivo (indicativo di carcinoma) in 24/35 (69%) gatti con FOSCC, contro 2/60 (3%) controlli (sensibilità: 69%; specifici-tà: 97%). La quota di FOSCC identificati era significativamente maggiore nei gatti sottoposti a prelievo in anestesia generale rispetto ai gatti svegli. Questi risultati sono incoraggianti per il riconoscimento precoce del FOSCC tramite brushing orale. Saranno necessari ulteriori studi su casistiche più ampie per validare questa metodica e migliorarne la sensibilità.