Longitudinal Micro CT-driven biomarkers as a tool to evaluate lung fibrosis progression and antifibrotic efficacy in a refined Bleomycin mouse model

Idiopathic pulmonary fibrosis (IPF) is a chronic progressive disease with no curative pharmacological treatment. Animal models play an essential role in revealing molecular mechanisms involved in the pathogenesis of the disease. Bleomycin (BLM)-induced lung fibrosis is the most widely used and characterized model for anti-fibrotic drugs screening. However, several issues have been reported, such as the identification of an optimal BLM dose and administration scheme as well as gender-specificity. Moreover, the balance between disease resolution, an appropriate time window for therapeutic intervention and animal welfare remains critical aspects yet to be fully elucidated. In this thesis, Micro CT imaging has been used as a tool to identify the ideal BLM dose regimen to induce sustained lung fibrosis in mice as well as to assess the anti-fibrotic effect of Nintedanib (NINT) treatment upon this BLM administration regimen. In order to select the optimal BLM dose scheme, C57bl/6 male mice were treated with BLM via oropharyngeal aspiration (OA), following either double or triple BLM administration. The triple BLM administration resulted in the most promising scheme, able to balance disease resolution, appropriate time-window for therapeutic intervention and animal welfare. The fibrosis progression was longitudinally assessed by micro-CT every 7 days for 5 weeks after BLM administration and 5 animals were sacrificed at each timepoint for the BALF and histological evaluation. The antifibrotic effect of NINT was assessed following different treatment regimens in this model. Herein, we have developed an optimized mouse model of pulmonary fibrosis, enabling three weeks of the therapeutic window to screen putative anti-fibrotic drugs. micro-CT scanning, allowed us to monitor the progression of lung fibrosis and the therapeutical response longitudinally in the same subject, drastically reducing the number of animals involved in the experiment.

Content-aware approach for improving biomedical image analysis: an interdisciplinary study series

Biomedicine is a highly interdisciplinary research area at the interface of sciences, anatomy, physiology, and medicine. In the last decade, biomedical studies have been greatly enhanced by the introduction of new technologies and techniques for automated quantitative imaging, thus considerably advancing the possibility to investigate biological phenomena through image analysis. However, the effectiveness of this interdisciplinary approach is bounded by the limited knowledge that a biologist and a computer scientist, by professional training, have of each other’s fields. The possible solution to make up for both these lacks lies in training biologists to make them interdisciplinary researchers able to develop dedicated image processing and analysis tools by exploiting a content-aware approach. The aim of this Thesis is to show the effectiveness of a content-aware approach to automated quantitative imaging, by its application to different biomedical studies, with the secondary desirable purpose of motivating researchers to invest in interdisciplinarity. Such content-aware approach has been applied firstly to the phenomization of tumour cell response to stress by confocal fluorescent imaging, and secondly, to the texture analysis of trabecular bone microarchitecture in micro-CT scans. Third, this approach served the characterization of new 3-D multicellular spheroids of human stem cells, and the investigation of the role of the Nogo-A protein in tooth innervation. Finally, the content-aware approach also prompted to the development of two novel methods for local image analysis and colocalization quantification. In conclusion, the content-aware approach has proved its benefit through building new approaches that have improved the quality of image analysis, strengthening the statistical significance to allow unveiling biological phenomena. Hopefully, this Thesis will contribute to inspire researchers to striving hard for pursuing interdisciplinarity.

An experimental and theoretical approach to correct for the scattered radiation in an X-ray computer tomography system for industrial applications

The main problem connected to cone beam computed tomography (CT) systems for industrial applications employing 450 kV X-ray tubes is the high amount of scattered radiation which is added to the primary radiation (signal). This stray radiation leads to a significant degradation of the image quality. A better understanding of the scattering and methods to reduce its effects are therefore necessary to improve the image quality. Several studies have been carried out in the medical field at lower energies, whereas studies in industrial CT, especially for energies up to 450 kV, are lacking. Moreover, the studies reported in literature do not consider the scattered radiation generated by the CT system structure and the walls of the X-ray room (environmental scatter). In order to investigate the scattering on CT projections a GEANT4-based Monte Carlo (MC) model was developed. The model, which has been validated against experimental data, has enabled the calculation of the scattering including the environmental scatter, the optimization of an anti-scatter grid suitable for the CT system, and the optimization of the hardware components of the CT system. The investigation of multiple scattering in the CT projections showed that its contribution is 2.3 times the one of primary radiation for certain objects. The results of the environmental scatter showed that it is the major component of the scattering for aluminum box objects of front size 70 x 70 mm2 and that it strongly depends on the thickness of the object and therefore on the projection. For that reason, its correction is one of the key factors for achieving high quality images. The anti-scatter grid optimized by means of the developed MC model was found to reduce the scatter-toprimary ratio in the reconstructed images by 20 %. The object and environmental scatter calculated by means of the simulation were used to improve the scatter correction algorithm which could be patented by Empa. The results showed that the cupping effect in the corrected image is strongly reduced. The developed CT simulation is a powerful tool to optimize the design of the CT system and to evaluate the contribution of the scattered radiation to the image. Besides, it has offered a basis for a new scatter correction approach by which it has been possible to achieve images with the same spatial resolution as state-of-the-art well collimated fan-beam CT with a gain in the reconstruction time of a factor 10. This result has a high economic impact in non-destructive testing and evaluation, and reverse engineering.

New methodologies in CT perfusion and MRI analysis to develop cancer imaging biomarkers

Quantitative imaging in oncology aims at developing imaging biomarkers for diagnosis and prediction of cancer aggressiveness and therapy response before any morphological change become visible. This Thesis exploits Computed Tomography perfusion (CTp) and multiparametric Magnetic Resonance Imaging (mpMRI) for investigating diverse cancer features on different organs. I developed a voxel-based image analysis methodology in CTp and extended its use to mpMRI, for performing precise and accurate analyses at single-voxel level. This is expected to improve reproducibility of measurements and cancer mechanisms’ comprehension and clinical interpretability. CTp has not entered the clinical routine yet, although its usefulness in the monitoring of cancer angiogenesis, due to different perfusion computing methods yielding unreproducible results. Instead, machine learning applications in mpMRI, useful to detect imaging features representative of cancer heterogeneity, are mostly limited to clinical research, because of results’ variability and difficult interpretability, which make clinicians not confident in clinical applications. In hepatic CTp, I investigated whether, and under what conditions, two widely adopted perfusion methods, Maximum Slope (MS) and Deconvolution (DV), could yield reproducible parameters. To this end, I developed signal processing methods to model the first pass kinetics and remove any numerical cause hampering the reproducibility. In mpMRI, I proposed a new approach to extract local first-order features, aiming at preserving spatial reference and making their interpretation easier. In CTp, I found out the cause of MS and DV non-reproducibility: MS and DV represent two different states of the system. Transport delays invalidate MS assumptions and, by correcting MS formulation, I have obtained the voxel-based equivalence of the two methods. In mpMRI, the developed predictive models allowed (i) detecting rectal cancers responding to neoadjuvant chemoradiation showing, at pre-therapy, sparse coarse subregions with altered density, and (ii) predicting clinically significant prostate cancers stemming from the disproportion between high- and low- diffusivity gland components.

Evaluation of bone strength: microtomographic tecniques

This thesis is a part of a larger study about the characterization of mechanical and histomorphometrical properties of bone. The main objects of this study were the bone tissue properties and its resistance to mechanical loads. Moreover, the knowledge about the equipment selected to carry out the analyses, the micro-computed tomography (micro-CT), was improved. Particular attention was given to the reliability over time of the measuring instrument. In order to understand the main characteristics of bone mechanical properties a study of the skeletal, the bones of which it is composed and biological principles that drive their formation and remodelling, was necessary. This study has led to the definition of two macro-classes describing the main components responsible for the resistance to fracture of bone: quantity and quality of bone. The study of bone quantity is the current clinical standard measure for so-called bone densitometry, and research studies have amply demonstrated that the amount of tissue is correlated with its mechanical properties of elasticity and fracture. However, the models presented in the literature, including information on the mere quantity of tissue, have often been limited in describing the mechanical behaviour. Recent investigations have underlined that also the bone-structure and the tissue-mineralization play an important role in the mechanical characterization of bone tissue. For this reason in this thesis the class defined as bone quality was mainly studied, splitting it into two sub-classes of bone structure and tissue quality. A study on bone structure was designed to identify which structural parameters, among the several presented in the literature, could be integrated with the information about quantity, in order to better describe the mechanical properties of bone. In this way, it was also possible to analyse the iteration between structure and function. It has been known for long that bone tissue is capable of remodeling and changing its internal structure according to loads, but the dynamics of these changes are still being analysed. This part of the study was aimed to identify the parameters that could quantify the structural changes of bone tissue during the development of a given disease: osteoarthritis. A study on tissue quality would have to be divided into different classes, which would require a scale of analysis not suitable for the micro-CT. For this reason the study was focused only on the mineralization of the tissue, highlighting the difference between bone density and tissue density, working in a context where there is still an ongoing scientific debate.

Prediction Nomogram for 68Ga-PSMA-11 PET/CT in different clinical settings of PSA failure after radical treatment for Prostate Cancer

Objective The objective of this study was to develop a clinical nomogram to predict gallium-68 prostate-specific membrane antigen positron emission tomography/computed tomography (68Ga-PSMA-11-PET/CT) positivity in different clinical settings of PSA failure. Materials and methods Seven hundred three (n = 703) prostate cancer (PCa) patients with confirmed PSA failure after radical therapy were enrolled. Patients were stratified according to different clinical settings (first-time biochemical recurrence [BCR]: group 1; BCR after salvage therapy: group 2; biochemical persistence after radical prostatectomy [BCP]: group 3; advanced stage PCa before second-line systemic therapies: group 4). First, we assessed 68Ga-PSMA-11-PET/CT positivity rate. Second, multivariable logistic regression analyses were used to determine predictors of positive scan. Third, regression-based coefficients were used to develop a nomogram predicting positive 68Ga-PSMA-11-PET/CT result and 200 bootstrap resamples were used for internal validation. Fourth, receiver operating characteristic (ROC) analysis was used to identify the most informative nomogram’s derived cut-off. Decision curve analysis (DCA) was implemented to quantify nomogram’s clinical benefit. Results 68Ga-PSMA-11-PET/CT overall positivity rate was 51.2%, while it was 40.3% in group 1, 54% in group 2, 60.5% in group 3, and 86.9% in group 4 (p < 0.001). At multivariable analyses, ISUP grade, PSA, PSA doubling time, and clinical setting were independent predictors of a positive scan (all p ≤ 0.04). A nomogram based on covariates included in the multivariate model demonstrated a bootstrap-corrected accuracy of 82%. The nomogram-derived best cut-off value was 40%. In DCA, the nomogram revealed clinical net benefit of > 10%. Conclusions This novel nomogram proved its good accuracy in predicting a positive scan, with values ≥ 40% providing the most informative cut-off in counselling patients to 68Ga-PSMA-11-PET/CT. This tool might be important as a guide to clinicians in the best use of PSMA-based PET imaging.

Inside standard and hyperspectral low-dose CT variational imaging: parameter identification and material decomposition

The main contribution of this thesis is the proposal of novel strategies for the selection of parameters arising in variational models employed for the solution of inverse problems with data corrupted by Poisson noise. In light of the importance of using a significantly small dose of X-rays in Computed Tomography (CT), and its need of using advanced techniques to reconstruct the objects due to the high level of noise in the data, we will focus on parameter selection principles especially for low photon-counts, i.e. low dose Computed Tomography. For completeness, since such strategies can be adopted for various scenarios where the noise in the data typically follows a Poisson distribution, we will show their performance for other applications such as photography, astronomical and microscopy imaging. More specifically, in the first part of the thesis we will focus on low dose CT data corrupted only by Poisson noise by extending automatic selection strategies designed for Gaussian noise and improving the few existing ones for Poisson. The new approaches will show to outperform the state-of-the-art competitors especially in the low-counting regime. Moreover, we will propose to extend the best performing strategy to the hard task of multi-parameter selection showing promising results. Finally, in the last part of the thesis, we will introduce the problem of material decomposition for hyperspectral CT, which data encodes information of how different materials in the target attenuate X-rays in different ways according to the specific energy. We will conduct a preliminary comparative study to obtain accurate material decomposition starting from few noisy projection data.

Efficacia delle griglie in titanio con osso particolato nella ricostruzione dei difetti alveolari tridimensionali dei mascellari

Questa tesi valuta l’efficacia della tecnica delle griglie in titanio con osso particolato nella ricostruzione dei difetti alveolari tridimensionali ai fini della riabilitazione dentale implanto-protesica. Il primo studio ha considerato la metodica in termini di complicanze post-operatorie e di risultati implanto-protesici. Sono stati considerati 24 pazienti con difetti tridimensionali trattati con l’applicazione di 34 griglie di titanio e osso particolato e riabilitati protesicamente dopo circa 8-9 mesi. 4 su 34 griglie sono state rimosse prima dell’inserimento implantare (11.76% di fallimento totale); 20 su 34 griglie si sono esposte per deiscenza dei tessuti molli (58.82% di complicanze): 4 (11.77%) prima e 16 (47.05%) dopo le prime 4-6 settimane dall’intervento; in nessun caso il piano di trattamento implanto-protesico ha subito variazioni. Dopo un follow-up medio di 20 (3-48) mesi dal carico protesico, nessuno degli 88 impianti ha perso la propria osteo-integrazione (100% di sopravvivenza implantare), con un valore complessivo di successo implantare di 82.9%. Il secondo studio ha calcolato in termini volumetrici la ricostruzione ossea ottenuta con griglie e la sua corre-lazione con l’estensione dell’esposizione e la tempistica del suo verificarsi. Sono stati valutati 12 pazienti con 15 difetti alveolari. Per ciascun sito sono state studiate le immagini TC con un software dedicato per misurare i volumi in tre dimensioni: il volume di osso non formatosi rispetto a quanto pianificato, lacking bone volume (LBV), è stato calcolato sottraendo il volume di osso ricostruito, reconstructed bone volume (RBV) in fase di ri-entro chirurgico dal volume di osso pianificato pre-operativamente, planned bone volume (PBV). LBV è risultato direttamente proporzionale all’area di esposizione della griglia, con un valore del 16.3% di LBV per ogni cm2 di griglia esposta. Si sono evidenziate, inoltre, correlazioni positive tra LBV , la tempistica precoce di esposizione e il valore di PBV.