Valutazione dell’impiego dei test per la genotipizzazione di HPV e l’espressione degli oncogeni virali nel follow-up di donne conizzate per lesioni cervicali di alto grado nello screening del cervico-carcinoma della Regione Emilia-Romagna

Obiettivi: Valutare la prevalenza dei diversi genotipi di HPV in pazienti con diagnosi di CIN2/3 nella Regione Emilia-Romagna, la persistenza genotipo-specifica di HPV e l’espressione degli oncogeni virali E6/E7 nel follow-up post-trattamento come fattori di rischio di recidiva/persistenza o progressione di malattia; verificare l’applicabilità di nuovi test diagnostici biomolecolari nello screening del cervicocarcinoma. Metodi: Sono state incluse pazienti con citologia di screening anormale, sottoposte a trattamento escissionale (T0) per diagnosi di CIN2/3 su biopsia mirata. Al T0 e durante il follow-up a 6, 12, 18 e 24 mesi, oltre al Pap test e alla colposcopia, sono state effettuate la ricerca e la genotipizzazione dell'HPV DNA di 28 genotipi. In caso di positività al DNA dei 5 genotipi 16, 18, 31, 33 e/o 45, si è proceduto alla ricerca dell'HPV mRNA di E6/E7. Risultati preliminari: Il 95.8% delle 168 pazienti selezionate è risultato HPV DNA positivo al T0. Nel 60.9% dei casi le infezioni erano singole (prevalentemente da HPV 16 e 31), nel 39.1% erano multiple. L'HPV 16 è stato il genotipo maggiormente rilevato (57%). Il 94.3% (117/124) delle pazienti positive per i 5 genotipi di HPV DNA sono risultate mRNA positive. Abbiamo avuto un drop-out di 38/168 pazienti. A 18 mesi (95% delle pazienti) la persistenza dell'HPV DNA di qualsiasi genotipo era del 46%, quella dell'HPV DNA dei 5 genotipi era del 39%, con espressione di mRNA nel 21%. Abbiamo avuto recidiva di malattia (CIN2+) nel 10.8% (14/130) a 18 mesi. Il pap test era negativo in 4/14 casi, l'HPV DNA test era positivo in tutti i casi, l'mRNA test in 11/12 casi. Conclusioni: L'HR-HPV DNA test è più sensibile della citologia, l'mRNA test è più specifico nell'individuare una recidiva. I dati definitivi saranno disponibili al termine del follow-up programmato.

Ottimizzazione del consumo dei grassi durante esercizio fisico: studio di un test per il target mirato di allenamento individuale

Il primo studio ha verificato l'affidabilità del software Polimedicus e gli effetti indotti d'allenamento arobico all’intensità del FatMax. 16 soggetti sovrappeso, di circa 40-55anni, sono stati arruolati e sottoposti a un test incrementale fino a raggiungere un RER di 0,95, e da quel momento il carico è stato aumentato di 1 km/ h ogni minuto fino a esaurimento. Successivamente, è stato verificato se i valori estrapolati dal programma erano quelli che si possono verificare durante a un test a carico costante di 1ora. I soggetti dopo 8 settimane di allenamento hanno fatto un altro test incrementale. Il dati hanno mostrato che Polimedicus non è molto affidabile, soprattutto l'HR. Nel secondo studio è stato sviluppato un nuovo programma, Inca, ed i risultati sono stati confrontati con i dati ottenuti dal primo studio con Polimedicus. I risultati finali hanno mostrato che Inca è più affidabile. Nel terzo studio, abbiamo voluto verificare l'esattezza del calcolo del FatMax con Inca e il test FATmaxwork. 25 soggetti in sovrappeso, tra 40-55 anni, sono stati arruolati e sottoposti al FATmaxwork test. Successivamente, è stato verificato se i valori estrapolati da INCA erano quelli che possono verificarsi durante un carico di prova costante di un'ora. L'analisi ha mostrato una precisione del calcolo della FatMax durante il carico di lavoro. Conclusione: E’ emersa una certa difficoltà nel determinare questo parametro, sia per la variabilità inter-individuale che intra-individuale. In futuro bisognerà migliorare INCA per ottenere protocolli di allenamento ancora più validi.

Non-thermal phenomena in galaxy clusters: the LOFAR revolution

Turbulence introduced into the intra-cluster medium (ICM) through cluster merger events transfers energy to non-thermal components (relativistic particles and magnetic fields) and can trigger the formation of diffuse synchrotron radio sources. Owing to their steep synchrotron spectral index, such diffuse sources can be better studied at low radio frequencies. In this respect, the LOw Frequency ARray (LOFAR) is revolutionizing our knowledge thanks to its unprecedented resolution and sensitivity below 200 MHz. In this Thesis we focus on the study of radio halos (RHs) by using LOFAR data. In the first part of this work we analyzed the largest-ever sample of galaxy clusters observed at radio frequencies. This includes 309 Planck clusters from the Second Data Release of the LOFAR Two Metre Sky Survey (LoTSS-DR2), which span previously unexplored ranges of mass and redshift. We detected 83 RHs, half of which being new discoveries. In 140 clusters we lack a detected RH; for this sub-sample we developed new techniques to derive upper limits to their radio powers. By comparing detections and upper limits, we carried out the first statistical analysis of populations of clusters observed at low frequencies and tested theoretical formation models. In the second part of this Thesis we focused on ultra-steep spectrum radio halos. These sources are almost undetected at GHz frequencies, but are thought to be common at low frequencies. We presented LOFAR observations of two interesting clusters hosting ultra-steep spectrum radio halos. With complementary radio and X-ray observations we constrained the properties and origin of these targets.

Galaxy cluster cosmology in KiDS, Planck, and Euclid surveys

In this Thesis, we present a series of works that encompass the fundamental steps of cosmological analyses based on galaxy clusters, spanning from mass calibration to deriving cosmological constraints through counts and clustering. Firstly, we focus on the 3D two-point correlation function (2PCF) of the galaxy cluster sample by Planck Collaboration XXVII (2016). The masses of these clusters are expected to be underestimated, as they are derived from a scaling relation calibrated through X-ray observations. We derived a mass bias which disagrees with simulation predictions, consistent with what derived by Planck Collaboration VI (2020). Furthermore, in this Thesis we analyse the cluster counts and 2PCF, respectively, of the photometric galaxy cluster sample developed by Maturi et al. (2019), based on the third data release of KiDS (KiDS-DR3, de Jong et al. 2017). We derived constraints on fundamental cosmological parameters which are consistent and competitive, in terms of uncertainties, with other state-of-the-art cosmological analyses. Then, we introduce a novel approach to establish galaxy colour-redshift relations for cluster weak-lensing analyses, regardless of the specific photometric bands in use. This method optimises the selection completeness of cluster background galaxies while maintaining a defined purity threshold. Based on the galaxy sample by Bisigello et al. (2020), we calibrated two colour selections, one relying on the ground-based griz bands, and the other including the griz and Euclid YJH bands. In addition, we present the preliminary work on the weak-lensing mass calibration of the clusters detected by Maturi et al. (in prep.) in the fourth data release of KiDS (KiDS-1000, Kuijken et al. 2019). This mass calibration will enable the cosmological analyses based on cluster counts and clustering, from which we expect remarkable improvements in the results compared to those derived in KiDS-DR3.

Mass transport in polymers

The study of mass transport in polymeric membranes has grown in importance due to its potential application in many processes such as separation of gases and vapors, packaging, controlled drug release. The diffusion of a low molecular weight species in a polymer is often accompanied by other phenomena like swelling, reactions, stresses, that have not been investigated in all their aspects yet. Furthermore, novel materials have been developed that include inorganic fillers, reactive functional groups or ions, that make the scenery even more complicated. The present work focused on the experimental study of systems where the diffusion is accompanied by other processes; suitable models were also developed to describe the particular circumstances in order to understand the underlying concepts and be able to design the performances of the material. The effect of solvent-induced deformation in polymeric films during sorption processes was studied since the dilation, especially in constrained membranes, can cause the development of stresses and therefore early failures of the material. The bending beam technique was used to test the effects of the dilation and the stress induced in the polymer by penetrant diffusion. A model based on the laminate theory was developed that accounts for the swelling and is able to predict the stress that raise in the material. The addition of inorganic fillers affects the transport properties of polymeric films. Mixed matrix membranes based on fluorinated, high free volume matrices show attractive performances for separation purposes but there is a need for deeper investigation of the selectivity properties towards gases and vapors. A new procedure based on the NELF model was tested on the experimental data; it allows to predict solubility of every penetrant on the basis of data for one vapor. The method has proved to be useful also for the determination of the diffusion coefficient and for an estimation of the permeability in the composite materials. Oxygen scavenging systems can overcome lack of barrier properties in common polymers that forbids their application in sensitive applications as food packaging. The final goal of obtaining a membrane almost impermeable to oxygen leads to experimental times out of reach. Hence, a simple model was developed in order to describe the transport of oxygen in a membrane with also reactive groups and analyze the experimental data collected on SBS copolymers that show attractive scavenging capacity. Furthermore, a model for predicting the oxygen barrier behavior of a film formed as a blend of OSP in a common packaging material was built, considering particles capable of reactions with oxygen embedded in a non-reactive matrix. Perfluorosulphonic acid ionomers (PFSI) are capturing attention due to a high thermal and chemical resistance coupled with very peculiar transport properties, that make them appropriate to be used in fuel cells. The possible effect of different formation procedure was studied together with the swelling due to water sorption since both water uptake and dilation can dramatically affect the fuel cells performances. The water diffusion and sorption was studied with a FTIR-ATR spectrometer that can give deeper information on the bonds between water molecules and the sulphonic hydrophilic groups and, therefore, on the microstructure of the hydrated ionomer.

Identificazione di un QTL principale per resistenza a ruggine bruna sul cromosoma 7B di frumento duro

Leaf rust caused by Puccinia triticina is a serious disease of durum wheat (Triticum durum) worldwide. However, genetic and molecular mapping studies aimed at characterizing leaf rust resistance genes in durum wheat have been only recently undertaken. The Italian durum wheat cv. Creso shows a high level of resistance to P. triticina that has been considered durable and that appears to be due to a combination of a single dominant gene and one or more additional factors conferring partial resistance. In this study, the genetic basis of leaf rust resistance carried by Creso was investigated using 176 recombinant inbred lines (RILs) from the cross between the cv. Colosseo (C, leaf rust resistance donor) and Lloyd (L, susceptible parent). Colosseo is a cv. directly related to Creso with the leaf rust resistance phenotype inherited from Creso, and was considered as resistance donor because of its better adaptation to local (Emilia Romagna, Italy) cultivation environment. RILs have been artificially inoculated with a mixture of 16 Italian P. triticina isolates that were characterized for virulence to seedlings of 22 common wheat cv. Thatcher isolines each carrying a different leaf rust resistance gene, and for molecular genotypes at 15 simple sequence repeat (SSR) loci, in order to determine their specialization with regard to the host species. The characterization of the leaf rust isolates was conducted at the Cereal Disease Laboratory of the University of Minnesota (St. Paul, USA) (Chapter 2). A genetic linkage map was constructed using segregation data from the population of 176 RILs from the cross CL. A total of 662 loci, including 162 simple sequence repeats (SSRs) and 500 Diversity Arrays Technology markers (DArTs), were analyzed by means of the package EasyMap 0.1. The integrated SSR-DArT linkage map consisted of 554 loci (162 SSR and 392 DArT markers) grouped into 19 linkage blocks with an average marker density of 5.7 cM/marker. The final map spanned a total of 2022 cM, which correspond to a tetraploid genome (AABB) coverage of ca. 77% (Chapter 3). The RIL population was phenotyped for their resistance to leaf rust under artificial inoculation in 2006; the percentage of infected leaf area (LRS, leaf rust susceptibility) was evaluated at three stages through the disease developmental cycle and the area under disease progress curve (AUDPC) was then calculated. The response at the seedling stage (infection type, IT) was also investigated. QTL analysis was carried out by means of the Composite Interval Mapping method based on a selection of markers from the CL map. A major QTL (QLr.ubo-7B.2) for leaf rust resistance controlling both the seedling and the adult plant response, was mapped on the distal region of chromosome arm 7BL (deletion bin 7BL10-0.78-1.00), in a gene-dense region known to carry several genes/QTLs for resistance to rusts and other major cereal fungal diseases in wheat and barley. QLr.ubo-7B.2 was identified within a supporting interval of ca. 5 cM tightly associated with three SSR markers (Xbarc340.2, Xgwm146 e Xgwm344.2), and showed an R2 and an LOD peak value for the AUDPC equal to 72.9% an 44.5, respectively. Three additional minor QTLs were also detected (QLr.ubo-7B.1 on chr. 7BS; QLr.ubo-2A on chr. 2AL and QLr.ubo-3A on chr. 3AS) (Chapter 4). The presence of the major QTL (QLr.ubo-7B.2) was validated by a linkage disequilibrium (LD)-based test using field data from two different plant materials: i) a set of 62 advanced lines from multiple crosses involving Creso and his directly related resistance derivates Colosseo and Plinio, and ii) a panel of 164 elite durum wheat accessions representative of the major durum breeding program of the Mediterranean basin. Lines and accessions were phenotyped for leaf rust resistance under artificial inoculation in two different field trials carried out at Argelato (BO, Italy) in 2006 and 2007; the durum elite accessions were also evaluated in two additional field experiments in Obregon (Messico; 2007 and 2008) and in a green-house experiment (seedling resistance) at the Cereal Disease Laboratory (St. Paul, USA, 2008). The molecular characterization involved 14 SSR markers mapping on the 7BL chromosome region found to harbour the major QTL. Association analysis was then performed with a mixed-linear-model approach. Results confirmed the presence of a major QTL for leaf rust resistance, both at adult plant and at seedling stage, located between markers Xbarc340.2, Xgwm146 and Xgwm344.2, in an interval that coincides with the supporting interval (LOD-2) of QLr.ubo-7B.2 as resulted from the RIL QTL analysis. (Chapter 5). The identification and mapping of the major QTL associated to the durable leaf rust resistance carried by Creso, together with the identification of the associated SSR markers, will enhance the selection efficiency in durum wheat breeding programs (MAS, Marker Assisted Selection) and will accelerate the release of cvs. with durable resistance through marker-assisted pyramiding of the tagged resistance genes/QTLs most effective against wheat fungal pathogens.

Kinematic description of the rupture from strong motion data: strategies for a robust inversion

We present a non linear technique to invert strong motion records with the aim of obtaining the final slip and rupture velocity distributions on the fault plane. In this thesis, the ground motion simulation is obtained evaluating the representation integral in the frequency. The Green’s tractions are computed using the discrete wave-number integration technique that provides the full wave-field in a 1D layered propagation medium. The representation integral is computed through a finite elements technique, based on a Delaunay’s triangulation on the fault plane. The rupture velocity is defined on a coarser regular grid and rupture times are computed by integration of the eikonal equation. For the inversion, the slip distribution is parameterized by 2D overlapping Gaussian functions, which can easily relate the spectrum of the possible solutions with the minimum resolvable wavelength, related to source-station distribution and data processing. The inverse problem is solved by a two-step procedure aimed at separating the computation of the rupture velocity from the evaluation of the slip distribution, the latter being a linear problem, when the rupture velocity is fixed. The non-linear step is solved by optimization of an L2 misfit function between synthetic and real seismograms, and solution is searched by the use of the Neighbourhood Algorithm. The conjugate gradient method is used to solve the linear step instead. The developed methodology has been applied to the M7.2, Iwate Nairiku Miyagi, Japan, earthquake. The estimated magnitude seismic moment is 2.6326 dyne∙cm that corresponds to a moment magnitude MW 6.9 while the mean the rupture velocity is 2.0 km/s. A large slip patch extends from the hypocenter to the southern shallow part of the fault plane. A second relatively large slip patch is found in the northern shallow part. Finally, we gave a quantitative estimation of errors associates with the parameters.

Inferences on earthquake kinematic properties from data inversion: two different approaches

During my PhD, starting from the original formulations proposed by Bertrand et al., 2000 and Emolo & Zollo 2005, I developed inversion methods and applied then at different earthquakes. In particular large efforts have been devoted to the study of the model resolution and to the estimation of the model parameter errors. To study the source kinematic characteristics of the Christchurch earthquake we performed a joint inversion of strong-motion, GPS and InSAR data using a non-linear inversion method. Considering the complexity highlighted by superficial deformation data, we adopted a fault model consisting of two partially overlapping segments, with dimensions 15x11 and 7x7 km2, having different faulting styles. This two-fault model allows to better reconstruct the complex shape of the superficial deformation data. The total seismic moment resulting from the joint inversion is 3.0x1025 dyne.cm (Mw = 6.2) with an average rupture velocity of 2.0 km/s. Errors associated with the kinematic model have been estimated of around 20-30 %. The 2009 Aquila sequence was characterized by an intense aftershocks sequence that lasted several months. In this study we applied an inversion method that assumes as data the apparent Source Time Functions (aSTFs), to a Mw 4.0 aftershock of the Aquila sequence. The estimation of aSTFs was obtained using the deconvolution method proposed by Vallée et al., 2004. The inversion results show a heterogeneous slip distribution, characterized by two main slip patches located NW of the hypocenter, and a variable rupture velocity distribution (mean value of 2.5 km/s), showing a rupture front acceleration in between the two high slip zones. Errors of about 20% characterize the final estimated parameters.

Search for the MSSM Neutral Higgs Boson in the mu+mu- final state with the CMS experiment at LHC

In this thesis, my work in the Compact Muon Solenoid (CMS) experiment on the search for the neutral Minimal Supersymmetric Standard Model (MSSM) Higgs decaying into two muons is presented. The search is performed on the full data collected during the years 2011 and 2012 by CMS in proton-proton collisions at CERN Large Hadron Collider (LHC). The MSSM is explored within the most conservative benchmark scenario, m_h^{max}, and within its modified versions, m_h^{mod +} and m_h^{mod -}. The search is sensitive to MSSM Higgs boson production in association with a b\bar{b} quark pair and to the gluon-gluon fusion process. In the m_h^{max} scenario, the results exclude values of tanB larger than 15 in the m_A range 115-200 GeV, and values of tanB greater than 30 in the m_A range up to 300 GeV. There are no significant differences in the results obtained within the three different scenarios considered. Comparisons with other neutral MSSM Higgs searches are shown.

Development of new tools and devices for CMB and foreground data analysis and future experiments

The discovery of the Cosmic Microwave Background (CMB) radiation in 1965 is one of the fundamental milestones supporting the Big Bang theory. The CMB is one of the most important source of information in cosmology. The excellent accuracy of the recent CMB data of WMAP and Planck satellites confirmed the validity of the standard cosmological model and set a new challenge for the data analysis processes and their interpretation. In this thesis we deal with several aspects and useful tools of the data analysis. We focus on their optimization in order to have a complete exploitation of the Planck data and contribute to the final published results. The issues investigated are: the change of coordinates of CMB maps using the HEALPix package, the problem of the aliasing effect in the generation of low resolution maps, the comparison of the Angular Power Spectrum (APS) extraction performances of the optimal QML method, implemented in the code called BolPol, and the pseudo-Cl method, implemented in Cromaster. The QML method has been then applied to the Planck data at large angular scales to extract the CMB APS. The same method has been applied also to analyze the TT parity and the Low Variance anomalies in the Planck maps, showing a consistent deviation from the standard cosmological model, the possible origins for this results have been discussed. The Cromaster code instead has been applied to the 408 MHz and 1.42 GHz surveys focusing on the analysis of the APS of selected regions of the synchrotron emission. The new generation of CMB experiments will be dedicated to polarization measurements, for which are necessary high accuracy devices for separating the polarizations. Here a new technology, called Photonic Crystals, is exploited to develop a new polarization splitter device and its performances are compared to the devices used nowadays.

Advanced modelling of time domain electromagnetic data with updated hydrogeological interpretations

Several countries have acquired, over the past decades, large amounts of area covering Airborne Electromagnetic data. Contribution of airborne geophysics has dramatically increased for both groundwater resource mapping and management proving how those systems are appropriate for large-scale and efficient groundwater surveying. We start with processing and inversion of two AEM dataset from two different systems collected over the Spiritwood Valley Aquifer area, Manitoba, Canada respectively, the AeroTEM III (commissioned by the Geological Survey of Canada in 2010) and the “Full waveform VTEM” dataset, collected and tested over the same survey area, during the fall 2011. We demonstrate that in the presence of multiple datasets, either AEM and ground data, due processing, inversion, post-processing, data integration and data calibration is the proper approach capable of providing reliable and consistent resistivity models. Our approach can be of interest to many end users, ranging from Geological Surveys, Universities to Private Companies, which are often proprietary of large geophysical databases to be interpreted for geological and\or hydrogeological purposes. In this study we deeply investigate the role of integration of several complimentary types of geophysical data collected over the same survey area. We show that data integration can improve inversions, reduce ambiguity and deliver high resolution results. We further attempt to use the final, most reliable output resistivity models as a solid basis for building a knowledge-driven 3D geological voxel-based model. A voxel approach allows a quantitative understanding of the hydrogeological setting of the area, and it can be further used to estimate the aquifers volumes (i.e. potential amount of groundwater resources) as well as hydrogeological flow model prediction. In addition, we investigated the impact of an AEM dataset towards hydrogeological mapping and 3D hydrogeological modeling, comparing it to having only a ground based TEM dataset and\or to having only boreholes data.

Differential expression analysis for sequence count data via mixtures of negative binomials

The recent advent of Next-generation sequencing technologies has revolutionized the way of analyzing the genome. This innovation allows to get deeper information at a lower cost and in less time, and provides data that are discrete measurements. One of the most important applications with these data is the differential analysis, that is investigating if one gene exhibit a different expression level in correspondence of two (or more) biological conditions (such as disease states, treatments received and so on). As for the statistical analysis, the final aim will be statistical testing and for modeling these data the Negative Binomial distribution is considered the most adequate one especially because it allows for "over dispersion". However, the estimation of the dispersion parameter is a very delicate issue because few information are usually available for estimating it. Many strategies have been proposed, but they often result in procedures based on plug-in estimates, and in this thesis we show that this discrepancy between the estimation and the testing framework can lead to uncontrolled first-type errors. We propose a mixture model that allows each gene to share information with other genes that exhibit similar variability. Afterwards, three consistent statistical tests are developed for differential expression analysis. We show that the proposed method improves the sensitivity of detecting differentially expressed genes with respect to the common procedures, since it is the best one in reaching the nominal value for the first-type error, while keeping elevate power. The method is finally illustrated on prostate cancer RNA-seq data.