Multiple testing in spatial epidemiology: a Bayesian approach

In this work we aim to propose a new approach for preliminary epidemiological studies on Standardized Mortality Ratios (SMR) collected in many spatial regions. A preliminary study on SMRs aims to formulate hypotheses to be investigated via individual epidemiological studies that avoid bias carried on by aggregated analyses. Starting from collecting disease counts and calculating expected disease counts by means of reference population disease rates, in each area an SMR is derived as the MLE under the Poisson assumption on each observation. Such estimators have high standard errors in small areas, i.e. where the expected count is low either because of the low population underlying the area or the rarity of the disease under study. Disease mapping models and other techniques for screening disease rates among the map aiming to detect anomalies and possible high-risk areas have been proposed in literature according to the classic and the Bayesian paradigm. Our proposal is approaching this issue by a decision-oriented method, which focus on multiple testing control, without however leaving the preliminary study perspective that an analysis on SMR indicators is asked to. We implement the control of the FDR, a quantity largely used to address multiple comparisons problems in the eld of microarray data analysis but which is not usually employed in disease mapping. Controlling the FDR means providing an estimate of the FDR for a set of rejected null hypotheses. The small areas issue arises diculties in applying traditional methods for FDR estimation, that are usually based only on the p-values knowledge (Benjamini and Hochberg, 1995; Storey, 2003). Tests evaluated by a traditional p-value provide weak power in small areas, where the expected number of disease cases is small. Moreover tests cannot be assumed as independent when spatial correlation between SMRs is expected, neither they are identical distributed when population underlying the map is heterogeneous. The Bayesian paradigm oers a way to overcome the inappropriateness of p-values based methods. Another peculiarity of the present work is to propose a hierarchical full Bayesian model for FDR estimation in testing many null hypothesis of absence of risk.We will use concepts of Bayesian models for disease mapping, referring in particular to the Besag York and Mollié model (1991) often used in practice for its exible prior assumption on the risks distribution across regions. The borrowing of strength between prior and likelihood typical of a hierarchical Bayesian model takes the advantage of evaluating a singular test (i.e. a test in a singular area) by means of all observations in the map under study, rather than just by means of the singular observation. This allows to improve the power test in small areas and addressing more appropriately the spatial correlation issue that suggests that relative risks are closer in spatially contiguous regions. The proposed model aims to estimate the FDR by means of the MCMC estimated posterior probabilities b i's of the null hypothesis (absence of risk) for each area. An estimate of the expected FDR conditional on data (\FDR) can be calculated in any set of b i's relative to areas declared at high-risk (where thenull hypothesis is rejected) by averaging the b i's themselves. The\FDR can be used to provide an easy decision rule for selecting high-risk areas, i.e. selecting as many as possible areas such that the\FDR is non-lower than a prexed value; we call them\FDR based decision (or selection) rules. The sensitivity and specicity of such rule depend on the accuracy of the FDR estimate, the over-estimation of FDR causing a loss of power and the under-estimation of FDR producing a loss of specicity. Moreover, our model has the interesting feature of still being able to provide an estimate of relative risk values as in the Besag York and Mollié model (1991). A simulation study to evaluate the model performance in FDR estimation accuracy, sensitivity and specificity of the decision rule, and goodness of estimation of relative risks, was set up. We chose a real map from which we generated several spatial scenarios whose counts of disease vary according to the spatial correlation degree, the size areas, the number of areas where the null hypothesis is true and the risk level in the latter areas. In summarizing simulation results we will always consider the FDR estimation in sets constituted by all b i's selected lower than a threshold t. We will show graphs of the\FDR and the true FDR (known by simulation) plotted against a threshold t to assess the FDR estimation. Varying the threshold we can learn which FDR values can be accurately estimated by the practitioner willing to apply the model (by the closeness between\FDR and true FDR). By plotting the calculated sensitivity and specicity (both known by simulation) vs the\FDR we can check the sensitivity and specicity of the corresponding\FDR based decision rules. For investigating the over-smoothing level of relative risk estimates we will compare box-plots of such estimates in high-risk areas (known by simulation), obtained by both our model and the classic Besag York Mollié model. All the summary tools are worked out for all simulated scenarios (in total 54 scenarios). Results show that FDR is well estimated (in the worst case we get an overestimation, hence a conservative FDR control) in small areas, low risk levels and spatially correlated risks scenarios, that are our primary aims. In such scenarios we have good estimates of the FDR for all values less or equal than 0.10. The sensitivity of\FDR based decision rules is generally low but specicity is high. In such scenario the use of\FDR = 0:05 or\FDR = 0:10 based selection rule can be suggested. In cases where the number of true alternative hypotheses (number of true high-risk areas) is small, also FDR = 0:15 values are well estimated, and \FDR = 0:15 based decision rules gains power maintaining an high specicity. On the other hand, in non-small areas and non-small risk level scenarios the FDR is under-estimated unless for very small values of it (much lower than 0.05); this resulting in a loss of specicity of a\FDR = 0:05 based decision rule. In such scenario\FDR = 0:05 or, even worse,\FDR = 0:1 based decision rules cannot be suggested because the true FDR is actually much higher. As regards the relative risk estimation, our model achieves almost the same results of the classic Besag York Molliè model. For this reason, our model is interesting for its ability to perform both the estimation of relative risk values and the FDR control, except for non-small areas and large risk level scenarios. A case of study is nally presented to show how the method can be used in epidemiology.

What do rendering options tell us about the translating mind? Testing the choice network analysis hypothesis

Frame. Assessing the difficulty of source texts and parts thereof is important in CTIS, whether for research comparability, for didactic purposes or setting price differences in the market. In order to empirically measure it, Campbell & Hale (1999) and Campbell (2000) developed the Choice Network Analysis (CNA) framework. Basically, the CNA’s main hypothesis is that the more translation options (a group of) translators have to render a given source text stretch, the higher the difficulty of that text stretch will be. We will call this the CNA hypothesis. In a nutshell, this research project puts the CNA hypothesis to the test and studies whether it does actually measure difficulty. Data collection. Two groups of participants (n=29) of different profiles and from two universities in different countries had three translation tasks keylogged with Inputlog, and filled pre- and post-translation questionnaires. Participants translated from English (L2) into their L1s (Spanish or Italian), and worked—first in class and then at home—using their own computers, on texts ca. 800–1000 words long. Each text was translated in approximately equal halves in two 1-hour sessions, in three consecutive weeks. Only the parts translated at home were considered in the study. Results. A very different picture emerged from data than that which the CNA hypothesis might predict: there was no prevalence of disfluent task segments when there were many translation options, nor was a prevalence of fluent task segments associated to fewer translation options. Indeed, there was no correlation between the number of translation options (many and few) and behavioral fluency. Additionally, there was no correlation between pauses and both behavioral fluency and typing speed. The discussed theoretical flaws and the empirical evidence lead to the conclusion that the CNA framework does not and cannot measure text and translation difficulty.

Mechanisms of cell senescence: p21 and DNA damage response in aging and longevity

Theory of aging postulates that aging is a remodeling process where the body of survivors progressively adapts to internal and external damaging agents they are exposed to during several decades. Thus , stress response and adaptation mechanisms play a fundamental role in the aging process where the capability of adaptating effects, certainly, also is related the lifespan of each individual. A key gene linking aging to stress response is indeed p21, an induction of cyclin-dependent kinase inhibitor which triggers cell growth arrest associated with senescence and damage response and notably is involved in the up-regulation of multiple genes that have been associated with senescence or implicated in age-related . This PhD thesis project that has been performed in collaboration with the Roninson Lab at Ordway Research Institute in Albany, NY had two main aims: -the testing the hypothesis that p21 polymorphisms are involved in longevity -Evaluating age-associated differences in gene expression and transcriptional response to p21 and DNA damage In the first project, trough PCR-sequencing and Sequenom strategies, we we found out that there are about 30 polymorphic variants in the p21 gene. In addition, we found an haplotpype located in -5kb region of the p21 promoter whose frequency is ~ 2 fold higher in centenarians than in the general population (Large-scale analysis of haplotype frequencies is currently in progress). Functional studies I carried out on the promoter highilighted that the ―centenarian‖ haplotype doesn’t affect the basal p21 promoter activity or its response to p53. However, there are many other possible physiological conditions in which the centenarian allele of the p21 promoter may potentially show a different response (IL6, IFN,progesterone, vitamin E, Vitamin D etc). In the second part, project #2, trough Microarrays we seeked to evaluate the differences in gene expression between centenarians, elderly, young in dermal fibroblast cultures and their response to p21 and DNA damage. Microarray analysis of gene expression in dermal fibroblast cultures of individuals of different ages yielded a tentative "centenarian signature". A subset of genes that were up- or downregulated in centenarians showed the same response to ectopic expression of p21, yielding a putative "p21-centenarian" signature. Trough RQ-PCR (as well Microarrays studies whose analysis is in progress) we tested the DNA damage response of the p21-centenarian signature genes showing a correlation stress/aging in additional sets of young and old samples treated with p21-inducing drug doxorubicin thus finding for a subset of of them , a response to stress age-related.

Tumori primitivi multipli del polmone. Profilo clinico e biologico in pazienti affetti da neoplasie polmonari multiple. Analisi di elementi clinici e marcatori biologici come possibili fattori di differenziazione dei tumori polmonari multipli

Background: MPLC represents a diagnostic challenge. Topic of the discussion is how to distinguish these patients as a metastatic or a multifocal disease. While in case of the different histology there are less doubt on the opposite in case of same histology is mandatory to investigate on other clinical features to rule out this question. Matherials and Methods: A retrospective review identified all patients treated surgically for a presumed diagnosis of SPLC. Pre-operative staging was obtained with Total CT scan and fluoro-deoxy positron emission tomography and mediastinoscopy. Patients with nodes interest or extra-thoracic location were excluded from this study. Epidermal growth factor receptor (EGFR) expression with complete immunohistochemical analisis was evaluated. Survival was estimated using Kaplan-Meyer method, and clinical features were estimated using a long-rank test or Cox proportional hazards model for categorical and continuous variable, respectively. Results: According to American College Chest Physician, 18 patients underwent to surgical resection for a diagnosis of MPLC. Of these, 8 patients had 3 or more nodules while 10 patients had less than 3 nodules. Pathologic examination demonstrated that 13/18(70%) of patients with multiple histological types was Adenocarcinoma, 2/18(10%) Squamous carcinoma, 2/18(10%) large cell carcinoma and 1/18(5%) Adenosquamosu carcinoma. Expression of EGFR has been evaluated in all nodules: in 7 patients of 18 (38%) the percentage of expression of each nodule resulted different. Conclusions: MPLC represent a multifocal disease where interactions of clinical informations with biological studies reinforce the diagnosis. EGFR could contribute to differentiate the nodules. However, further researches are necessary to validate this hypothesis.

Modulation of hypoxia-inducible factors as a therapeutic target for multiple myeloma

Hypoxia-inducible factor-1 alpha (HIF-1α) plays a critical role in survival and is associated with poor prognosis in solid tumors. The role of HIF-1α in multiple myeloma is not completely known. In the present study, we explored the effect of EZN2968, an locked nucleic acid antisense oligonucleotide against HIF-1α, as a molecular target in MM. A panel of MM cell lines and primary samples from MM patients were cultured in vitro in the presence of EZN2968 . Under normoxia culture condition, HIF-1α mRNA and protein expression was detectable in all MM cell lines and in CD138+ cells from newly diagnosed MM patients samples. Significant up-regulation of HIF-1α protein expression was observed after incubation with IL6 or IGF-I, confirming that HIF-1α can be further induced by biological stimuli. EZN2968 efficiently induces a selective and stable down-modulation of HIF-1α and decreased the secretion of VEGF released by MM cell. Treatment with EZN2968 gave rise to a progressive accumulation of cells in the S and subG0 phase. The analysis of p21, a cyclin-dependent kinase inhibitors controlling cell cycle check point, shows upregulation of protein levels. These results suggest that HIF-1α inhibition is sufficient for cell cycle arrest in normoxia, and for inducing an apoptotic pathways.. In the presence of bone marrow microenvironment, HIF-1α inhibition blocks MAPK kinase pathway and secretion of pro-surviaval cytokines ( IL6,VEGF,IL8) In this study we provide evidence that HIF-1α, even in the absence of hypoxia signal, is expressed in MM plasma cells and further inducible by bone marrow milieu stimuli; moreover its inhibition is sufficient to induce a permanent cell cycle arrest. Our data support the hypothesis that HIF-1α inhibition may suppress tumor growth by preventing proliferation of plasma cells through p21 activation and blocking pro-survival stimuli from bone marrow microenvironment.

The impact of polymorphisms in P-gp, DNA repair and folic acid metabolism genes in newly diagnosed multiple myeloma patients treated with thalidomide plus dexamethasone, with or without bortezomib

The principle aim of this study was to investigate biological predictors of response and resistance to multiple myeloma treatment. Two hypothesis had been proposed as responsible of responsiveness: SNPs in DNA repair and Folate pathway, and P-gp dependent efflux. As a first objective, panel of SNPs in DNA repair and Folate pathway genes, were analyzed. It was a retrospective study in a group of 454, previously untreated, MM patients enrolled in a randomized phase III open-label study. Results show that some SNPs in Folate pathway are correlated with response to MM treatment. MTR genotype was associated with favorable response in the overall population of MM patients. However, this relation, disappear after adjustment for treatment response. When poor responder includes very good partial response, partial response and stable/progressive disease MTFHR rs1801131 genotype was associated with poor response to therapy. This relation - unlike in MTR – was still significant after adjustment for treatment response. Identification of this genetic variant in MM patients could be used as an independent prognostic factor for therapeutic outcome in the clinical practice. In the second objective, basic disposition characteristics of bortezomib was investigated. We demonstrated that bortezomib is a P-gp substrate in a bi-directional transport study. We obtain apparent permeability rate values that together with solubility values can have a crucial implication in better understanding of bortezomib pharmacokinetics with respect to the importance of membrane transporters. Subsequently, in view of the importance of P-gp for bortezomib responsiveness a panel of SNPs in ABCB1 gene - coding for P-gp - were analyzed. In particular we analyzed five SNPs, none of them however correlated with treatment responsiveness. However, we found a significant association between ABCB1 variants and cytogenetic abnormalities. In particular, deletion of chromosome 17 and t(4;14) translocation were present in patients harboring rs60023214 and rs2038502 variants respectively.

A threshold hypothesis of institutional change: collective action in the Italian Alps during the 13th - 19th centuries.

This dissertation is about collective action issues in common property resources. Its focus is the “threshold hypothesis,” which posits the existence of a threshold in group size that drives the process of institutional change. This hypothesis is tested using a six-century dataset concerning the management of the commons by hundreds of communities in the Italian Alps. The analysis seeks to determine the group size threshold and the institutional changes that occur when groups cross this threshold. There are five main findings. First, the number of individuals in villages remained stable for six centuries, despite the population in the region tripling in the same period. Second, the longitudinal analysis of face-to-face assemblies and community size led to the empirical identification of a threshold size that triggered the transition from informal to more formal regimes to manage common property resources. Third, when groups increased in size, gradual organizational changes took place: large groups split into independent subgroups or structured interactions into multiple layers while maintaining a single formal organization. Fourth, resource heterogeneity seemed to have had no significant impact on various institutional characteristics. Fifth, social heterogeneity showed statistically significant impacts, especially on institutional complexity, consensus, and the relative importance of governance rules versus resource management rules. Overall, the empirical evidence from this research supports the “threshold hypothesis.” These findings shed light on the rationale of institutional change in common property regimes, and clarify the mechanisms of collective action in traditional societies. Further research may generalize these conclusions to other domains of collective action and to present-day applications.

Response times in computerized adaptive testing: a method for cheating detection

In the field of educational and psychological measurement, the shift from paper-based to computerized tests has become a prominent trend in recent years. Computerized tests allow for more complex and personalized test administration procedures, like Computerized Adaptive Testing (CAT). CAT, following the Item Response Theory (IRT) models, dynamically generates tests based on test-taker responses, driven by complex statistical algorithms. Even if CAT structures are complex, they are flexible and convenient, but concerns about test security should be addressed. Frequent item administration can lead to item exposure and cheating, necessitating preventive and diagnostic measures. In this thesis a method called "CHeater identification using Interim Person fit Statistic" (CHIPS) is developed, designed to identify and limit cheaters in real-time during test administration. CHIPS utilizes response times (RTs) to calculate an Interim Person fit Statistic (IPS), allowing for on-the-fly intervention using a more secret item bank. Also, a slight modification is proposed to overcome situations with constant speed, called Modified-CHIPS (M-CHIPS). A simulation study assesses CHIPS, highlighting its effectiveness in identifying and controlling cheaters. However, it reveals limitations when cheaters possess all correct answers. The M-CHIPS overcame this limitation. Furthermore, the method has shown not to be influenced by the cheaters’ ability distribution or the level of correlation between ability and speed of test-takers. Finally, the method has demonstrated flexibility for the choice of significance level and the transition from fixed-length tests to variable-length ones. The thesis discusses potential applications, including the suitability of the method for multiple-choice tests, assumptions about RT distribution and level of item pre-knowledge. Also limitations are discussed to explore future developments such as different RT distributions, unusual honest respondent behaviors, and field testing in real-world scenarios. In summary, CHIPS and M-CHIPS offer real-time cheating detection in CAT, enhancing test security and ability estimation while not penalizing test respondents.