Identifying candidate gene effects by restricting search space in a multivariate genetic analysis of white matter microstructure

Several genetic variants are thought to influence white matter (WM) integrity, measured with diffusion tensor imaging (DTI). Voxel based methods can test genetic associations, but heavy multiple comparisons corrections are required to adjust for searching the whole brain and for all genetic variants analyzed. Thus, genetic associations are hard to detect even in large studies. Using a recently developed multi-SNP analysis, we examined the joint predictive power of a group of 18 cholesterol-related single nucleotide polymorphisms (SNPs) on WM integrity, measured by fractional anisotropy. To boost power, we limited the analysis to brain voxels that showed significant associations with total serum cholesterol levels. From this space, we identified two genes with effects that replicated in individual voxel-wise analyses of the whole brain. Multivariate analyses of genetic variants on a reduced anatomical search space may help to identify SNPs with strongest effects on the brain from a broad panel of genes.

Exploration of genetic susceptibility to pancreatic cancer: A GWAS approach

Pancreatic cancer is the 4th most common cause for cancer death in the United States, accompanied by less than 5% five-year survival rate based on current treatments, particularly because it is usually detected at a late stage. Identifying a high-risk population to launch an effective preventive strategy and intervention to control this highly lethal disease is desperately needed. The genetic etiology of pancreatic cancer has not been well profiled. We hypothesized that unidentified genetic variants by previous genome-wide association study (GWAS) for pancreatic cancer, due to stringent statistical threshold or missing interaction analysis, may be unveiled using alternative approaches. To achieve this aim, we explored genetic susceptibility to pancreatic cancer in terms of marginal associations of pathway and genes, as well as their interactions with risk factors. We conducted pathway- and gene-based analysis using GWAS data from 3141 pancreatic cancer patients and 3367 controls with European ancestry. Using the gene set ridge regression in association studies (GRASS) method, we analyzed 197 pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Using the logistic kernel machine (LKM) test, we analyzed 17906 genes defined by University of California Santa Cruz (UCSC) database. Using the likelihood ratio test (LRT) in a logistic regression model, we analyzed 177 pathways and 17906 genes for interactions with risk factors in 2028 pancreatic cancer patients and 2109 controls with European ancestry. After adjusting for multiple comparisons, six pathways were marginally associated with risk of pancreatic cancer ( P < 0.00025): Fc epsilon RI signaling, maturity onset diabetes of the young, neuroactive ligand-receptor interaction, long-term depression (Ps < 0.0002), and the olfactory transduction and vascular smooth muscle contraction pathways (P = 0.0002; Nine genes were marginally associated with pancreatic cancer risk (P < 2.62 × 10−5), including five reported genes (ABO, HNF1A, CLPTM1L, SHH and MYC), as well as four novel genes (OR13C4, OR 13C3, KCNA6 and HNF4 G); three pathways significantly interacted with risk factors on modifying the risk of pancreatic cancer (P < 2.82 × 10−4): chemokine signaling pathway with obesity ( P < 1.43 × 10−4), calcium signaling pathway (P < 2.27 × 10−4) and MAPK signaling pathway with diabetes (P < 2.77 × 10−4). However, none of the 17906 genes tested for interactions survived the multiple comparisons corrections. In summary, our current GWAS study unveiled unidentified genetic susceptibility to pancreatic cancer using alternative methods. These novel findings provide new perspectives on genetic susceptibility to and molecular mechanisms of pancreatic cancer, once confirmed, will shed promising light on the prevention and treatment of this disease. ^

Dopamine transporter genotype predicts behavioural and neural measures of response inhibition

The ability to inhibit unwanted actions is a heritable executive function that may confer risk to disorders such as attention deficit hyperactivity disorder (ADHD). Converging evidence from pharmacology and cognitive neuroscience suggests that response inhibition is instantiated within frontostriatal circuits of the brain with patterns of activity that are modulated by the catecholamines dopamine and noradrenaline. A total of 405 healthy adult participants performed the stop-signal task, a paradigmatic measure of response inhibition that yields an index of the latency of inhibition, termed the stop-signal reaction time (SSRT). Using this phenotype, we tested for genetic association, performing high-density single-nucleotide polymorphism mapping across the full range of autosomal catecholamine genes. Fifty participants also underwent functional magnetic resonance imaging to establish the impact of associated alleles on brain and behaviour. Allelic variation in polymorphisms of the dopamine transporter gene (SLC6A3: rs37020; rs460000) predicted individual differences in SSRT, after corrections for multiple comparisons. Furthermore, activity in frontal regions (anterior frontal, superior frontal and superior medial gyri) and caudate varied additively with the T-allele of rs37020. The influence of genetic variation in SLC6A3 on the development of frontostriatal inhibition networks may represent a key risk mechanism for disorders of behavioural inhibition.

Genome-wide scan of healthy human connectome discovers SPON1 gene variant influencing dementia severity

Aberrant connectivity is implicated in many neurological and psychiatric disorders, including Alzheimer's disease and schizophrenia. However, other than a few disease-associated candidate genes, we know little about the degree to which genetics play a role in the brain networks; we know even less about specific genes that influence brain connections. Twin and family-based studies can generate estimates of overall genetic influences on a trait, but genome-wide association scans (GWASs) can screen the genome for specific variants influencing the brain or risk for disease. To identify the heritability of various brain connections, we scanned healthy young adult twins with high-field, highangular resolution diffusion MRI. We adapted GWASs to screen the brain's connectivity pattern, allowing us to discover genetic variants that affect the human brain's wiring. The association of connectivity with the SPON1 variant at rs2618516 on chromosome 11 (11p15.2) reached connectome-wide, genome-wide significance after stringent statistical corrections were enforced, and it was replicated in an independent subsample. rs2618516 was shown to affect brain structure in an elderly population with varying degrees of dementia. Older people who carried the connectivity variant had significantly milder clinical dementia scores and lower risk of Alzheimer's disease. As a posthoc analysis, we conducted GWASs on several organizational and topological network measures derived from the matrices to discover variants in and around genes associated with autism (MACROD2), development (NEDD4), and mental retardation (UBE2A) significantly associated with connectivity. Connectome-wide, genome-wide screening offers substantial promise to discover genes affecting brain connectivity and risk for brain diseases.

A population-based study of adverse physical effects following prostate cancer treatments: Results from the PiCTure (Prostate cancer treatment – your experience) study

Background:

Men and clinicians need reliable population based information when making decisions about investigation and treatment of prostate cancer. In the absence of clearly preferred treatments, differences in outcomes become more important.

Aim:

To investigate rates of adverse physical effects among prostate cancer survivors 2-15 years post diagnosis by treatment, and estimate population burden.

Methods:

A cross sectional, postal survey to 6,559 survivors (all ages) diagnosed with primary, invasive prostate cancer (ICD10-C61), identified in Northern Ireland and the Republic of Ireland via cancer registries. Questions included symptoms at diagnosis, treatments received and adverse physical effects (impotence, urinary incontinence, bowel problems, breast changes, libido loss, hot flashes, fatigue) experienced ‘ever’ and ‘current’ i.e. at questionnaire completion. Physical effect levels were weighted by age, country and time since diagnosis for all prostate cancer survivors. Bonferroni corrections were applied to account for multiple comparisons.

Results:

Adjusted response rate 54%, (n=3,348). 75% reported at least one current physical effect (90% ever), with 29% reporting at least three. These varied by treatment. Current impotence was reported by 76% post-prostatectomy, 64% post-external beam radiotherapy with hormone therapy, with average for all survivors of 57%. Urinary incontinence (overall current level: 16%) was highest post-prostatectomy (current 28%, ever 70%). 42% of brachytherapy patients reported no current adverse physical effects; however 43% reported current impotence and 8% current incontinence. Current hot flashes (41%), breast changes (18%) and fatigue (28%) were reported more commonly by patients on hormone therapy.

Conclusions:

This study provides evidence that adverse physical effects following prostate cancer represent a significant public health burden; an estimated 1.6% of men over 45 is a prostate cancer survivor with a current adverse physical effect. This information should facilitate investigation and treatment decision-making and follow-up care of patients.

Patient-reported 'ever had' and 'current' long term physical symptoms following prostate cancer treatments

OBJECTIVE: To document prostate cancer patient reported 'ever experienced' and 'current' prevalence of disease specific physical symptoms stratified by primary treatment received.
PATIENTS: 3,348 prostate cancer survivors 2-15 years post diagnosis.
METHODS: Cross-sectional, postal survey of 6,559 survivors diagnosed 2-15 years ago with primary, invasive PCa (ICD10-C61) identified via national, population based cancer registries in Northern Ireland and Republic of Ireland. Questions included symptoms at diagnosis, primary treatments and physical symptoms (impotence/urinary incontinence/bowel problems/breast changes/loss of libido/hot flashes/fatigue) experienced 'ever' and at questionnaire completion ("current"). Symptom proportions were weighted by age, country and time since diagnosis. Bonferroni corrections were applied for multiple comparisons.
RESULTS: Adjusted response rate 54%; 75% reported at least one 'current' physical symptom ('ever':90%), with 29% reporting at least three. Prevalence varied by treatment; overall 57% reported current impotence; this was highest following radical prostatectomy (RP)76% followed by external beam radiotherapy with concurrent hormone therapy (HT); 64%. Urinary incontinence (overall 'current' 16%) was highest following RP ('current'28%, 'ever'70%). While 42% of brachytherapy patients reported no 'current' symptoms; 43% reported 'current' impotence and 8% 'current' incontinence. 'Current' hot flashes (41%), breast changes (18%) and fatigue (28%) were reported more often by patients on HT.
CONCLUSION: Symptoms following prostate cancer are common, often multiple, persist long-term and vary by treatment. They represent a significant health burden. An estimated 1.6% of men over 45 is a prostate cancer survivor currently experiencing an adverse physical symptom. Recognition and treatment of physical symptoms should be prioritised in patient follow-up. This information should facilitate men and clinicians when deciding about treatment as differences in survival between radical treatments is minimal.

Investigating the role of potential genetic markers that can predict risk for steroid refractory inflammatory bowel disease

Contexte - La variation interindividuelle de la réponse aux corticostéroïdes (CS) est un problème important chez les patients atteints de maladies inflammatoires d’intestin. Ce problème est bien plus accentué chez les enfants avec la prévalence de la corticodépendance extrêmement (~40 %) élevée. La maladie réfractaire au CS a des répercussions sur le développement et le bien-être physique et psychologique des patients et impose des coûts médicaux élevés, particulièrement avec la maladie active comparativement à la maladie en rémission, le coût étant 2-3 fois plus élevé en ambulatoire et 20 fois plus élevé en hôpital. Il est ainsi primordial de déterminer les marqueurs prédictifs de la réponse aux CS. Les efforts précédents de découvrir les marqueurs cliniques et démographiques ont été équivoques, ce qui souligne davantage le besoin de marqueurs moléculaires. L'action des CS se base sur des processus complexes déterminés génétiquement. Deux gènes, le ABCB1, appartenant à la famille des transporteurs transmembraneaux, et le NR3C1, encodant le récepteur glucocorticoïde, sont des éléments importants des voies métaboliques. Nous avons postulé que les variations dans ces gènes ont un rôle dans la variabilité observée de la réponse aux CS et pourraient servir en tant que les marqueurs prédictifs. Objectifs - Nous avons visé à: (1) examiner le fardeau de la maladie réfractaire aux CS chez les enfants avec la maladie de Crohn (MC) et le rôle des caractéristiques cliniques et démographiques potentiellement liés à la réponse; (2) étudier l'association entre les variantes d'ADN de gène ABCB1 et la réponse aux CS; (3) étudier les associations entre les variantes d'ADN de gène NR3C1 et la réponse aux CS. Méthodes - Afin d’atteindre ces objectifs, nous avons mené une étude de cohorte des patients recrutés dans deux cliniques pédiatriques tertiaires de gastroentérologie à l’Ottawa (CHEO) et à Montréal (HSJ). Les patients avec la MC ont été diagnostiqués avant l'âge de 18 ans selon les critères standard radiologiques, endoscopiques et histopathologiques. La corticorésistance et la corticodépendance ont été définies en adaptant les critères reconnus. L’ADN, acquise soit du sang ou de la salive, était génotypée pour des variations à travers de gènes ABCB1 et NR3C1 sélectionnées à l’aide de la méthodologie de tag-SNP. La fréquence de la corticorésistance et la corticodépendance a été estimée assumant une distribution binomiale. Les associations entre les variables cliniques/démographiques et la réponse aux CS ont été examinées en utilisant la régression logistique en ajustant pour des variables potentielles de confusion. Les associations entre variantes génétiques de ABCB1 et NR3C1 et la réponse aux CS ont été examinées en utilisant la régression logistique assumant différents modèles de la transmission. Les associations multimarqueurs ont été examinées en utilisant l'analyse de haplotypes. Les variantes nongénotypées ont été imputées en utilisant les données de HAPMAP et les associations avec SNPs imputés ont été examinées en utilisant des méthodes standard. Résultats - Parmi 645 patients avec la MC, 364 (56.2%) ont reçu CS. La majorité de patients étaient des hommes (54.9 %); présentaient la maladie de l’iléocôlon (51.7%) ou la maladie inflammatoire (84.6%) au diagnostic et étaient les Caucasiens (95.6 %). Huit pourcents de patients étaient corticorésistants et 40.9% - corticodépendants. Le plus bas âge au diagnostic (OR=1.34, 95% CI: 1.03-3.01, p=0.040), la maladie cœxistante de la région digestive supérieure (OR=1.35, 95% CI: 95% CI: 1.06-3.07, p=0.031) et l’usage simultané des immunomodulateurs (OR=0.35, 95% CI: 0.16-0.75, p=0.007) ont été associés avec la corticodépendance. Un total de 27 marqueurs génotypés à travers de ABCB1 (n=14) et NR3C1 (n=13) ont été en l'Équilibre de Hardy-Weinberg, à l’exception d’un dans le gène NR3C1 (rs258751, exclu). Dans ABCB1, l'allèle rare de rs2032583 (OR=0.56, 95% CI: 0.34-0.95, p=0.029) et génotype hétérozygote (OR=0.52, 95% CI: 0.28-0.95 p=0.035) ont été négativement associes avec la dépendance de CS. Un haplotype à 3 marqueurs, comprenant le SNP fonctionnel rs1045642 a été associé avec la dépendance de CS (p empirique=0.004). 24 SNPs imputés introniques et six haplotypes ont été significativement associés avec la dépendance de CS. Aucune de ces associations n'a cependant maintenu la signification après des corrections pour des comparaisons multiples. Dans NR3C1, trois SNPs: rs10482682 (OR=1.43, 95% CI: 0.99-2.08, p=0.047), rs6196 (OR=0.55, 95% CI: 0.31-0.95, p=0.024), et rs2963155 (OR=0.64, 95% CI: 0.42-0.98, p=0.039), ont été associés sous un modèle additif, tandis que rs4912911 (OR=0.37, 95% CI: 0.13-1.00, p=0.03) et rs2963156 (OR=0.32, 95% CI: 0.07-1.12, p=0.047) - sous un modèle récessif. Deux haplotypes incluant ces 5 SNPs (AAACA et GGGCG) ont été significativement (p=0.006 et 0.01 empiriques) associés avec la corticodépendance. 19 SNPs imputés ont été associés avec la dépendance de CS. Deux haplotypes multimarqueurs (p=0.001), incluant les SNPs génotypés et imputés, ont été associés avec la dépendance de CS. Conclusion - Nos études suggèrent que le fardeau de la corticodépendance est élevé parmi les enfants avec le CD. Les enfants plus jeunes au diagnostic et ceux avec la maladie coexistante de la région supérieure ainsi que ceux avec des variations dans les gènes ABCB1 et NR3C1 étaient plus susceptibles de devenir corticodépendants.

Évaluation de l'acuité visuelle chez la personne âgée atteinte de troubles de la cognition

Objectif : L'évaluation de l'acuité visuelle (AV) chez la personne âgée atteinte de troubles cognitifs peut être limitée par le manque de collaboration ou les difficultés de communication du patient. Très peu d'études ont examiné l'AV chez les patients atteints de déficits sévères de la cognition. L’objectif de cette étude était d’évaluer l'AV chez la personne âgée vulnérable ayant des troubles cognitifs à l'aide d'échelles variées afin de vérifier leur capacité à répondre à ces échelles. Méthodes: Trois groupes de 30 sujets chacun ont été recrutés. Le premier groupe était constitué de sujets jeunes (Moy.±ET: 24.9±3.5ans) et le second, de sujets âgés (70.0±4.5ans), ne présentant aucun trouble de la cognition ou de la communication. Le troisième groupe, composé de sujets atteints de démence faible à sévère (85.6±6.9ans), a été recruté au sein des unités de soins de longue durée de l’Institut Universitaire de Gériatrie de Montréal. Le test du Mini-Mental State Examination (MMSE) a été réalisé pour chaque sujet afin de déterminer leur niveau cognitif. L’AV de chaque participant a été mesurée à l’aide de six échelles validées (Snellen, cartes de Teller, ETDRS-lettres,-chiffres,-Patty Pics,-E directionnel) et présentées selon un ordre aléatoire. Des tests non paramétriques ont été utilisés afin de comparer les scores d’AV entre les différentes échelles, après une correction de Bonferroni-Holm pour comparaisons multiples. Résultats: Le score moyen au MMSE chez les sujets atteints de démence était de 9.8±7.5, alors qu’il était de 17.8±3.7 et 5.2±4.6 respectivement, chez les sujets atteints de démence faible à modérée (MMSE supérieur ou égal à 13; n=11) et sévère (MMSE inférieur à 13; n=19). Tous les sujets des groupes 1 et 2 ont répondu à chacune des échelles. Une grande majorité de sujets avec démence ont répondu à toutes les échelles (n=19) alors qu’un seul sujet n’a répondu à aucune échelle d’AV. Au sein du groupe 3, les échelles d’AV fournissant les scores les plus faibles ont été les cartes de Teller (20/65) et les Patty Pics (20/62), quelque soit le niveau cognitif du sujet, alors que les meilleurs scores d’AV ont été obtenus avec les échelles de Snellen (20/35) et les lettres ETDRS (20/36). Une grande proportion de sujets avec démence sévère ont répondu aux cartes de Teller (n=18) mais le score d’AV obtenu était le plus faible (20/73). Au sein des trois groupes, l’échelle de lettres-ETDRS était la seule dont les scores d’AV ne différaient pas de ceux obtenus avec l’échelle de Snellen traditionnelle. Conclusions: L’acuité visuelle peut être mesurée chez la personne âgée atteinte de troubles cognitifs ou de la communication. Nos résultats indiquent que les échelles les plus universelles, utilisant des lettres comme optotypes, peuvent être utilisées avec de bons résultats chez les personnes âgées atteintes de démence sévère. Nos résultats suggèrent de plus que la mesure d’acuité visuelle doit être tentée chez toutes les personnes, peu importe leur niveau cognitif.

Assessment of CloudSat reflectivity measurements and ice cloud properties using ground-based and airborne cloud radar observations

A quantitative assessment of Cloudsat reflectivities and basic ice cloud properties (cloud base, top, and thickness) is conducted in the present study from both airborne and ground-based observations. Airborne observations allow direct comparisons on a limited number of ocean backscatter and cloud samples, whereas the ground-based observations allow statistical comparisons on much longer time series but with some additional assumptions. Direct comparisons of the ocean backscatter and ice cloud reflectivities measured by an airborne cloud radar and Cloudsat during two field experiments indicate that, on average, Cloudsat measures ocean backscatter 0.4 dB higher and ice cloud reflectivities 1 dB higher than the airborne cloud radar. Five ground-based sites have also been used for a statistical evaluation of the Cloudsat reflectivities and basic cloud properties. From these comparisons, it is found that the weighted-mean difference ZCloudsat − ZGround ranges from −0.4 to +0.3 dB when a ±1-h time lag around the Cloudsat overpass is considered. Given the fact that the airborne and ground-based radar calibration accuracy is about 1 dB, it is concluded that the reflectivities of the spaceborne, airborne, and ground-based radars agree within the expected calibration uncertainties of the airborne and ground-based radars. This result shows that the Cloudsat radar does achieve the claimed sensitivity of around −29 dBZ. Finally, an evaluation of the tropical “convective ice” profiles measured by Cloudsat has been carried out over the tropical site in Darwin, Australia. It is shown that these profiles can be used statistically down to approximately 9-km height (or 4 km above the melting layer) without attenuation and multiple scattering corrections over Darwin. It is difficult to estimate if this result is applicable to all types of deep convective storms in the tropics. However, this first study suggests that the Cloudsat profiles in convective ice need to be corrected for attenuation by supercooled liquid water and ice aggregates/graupel particles and multiple scattering prior to their quantitative use.

Differential effects of neuromuscular blockers on twitches and tetani in the isolated rat muscle: a multiple comparison study using simultaneous confidence intervals

In the present work a comparative quantitative evaluation of the differential effects of neuromuscular blockers on twitches and tetani was performed, encompassing: atracurium, cisatracurium, mivacurium, pancuronium, rocuronium and vecuronium. The sciatic nerve-extensor digitorum longus muscle of the rat was used, in vitro. Twitches were evoked at 0.1 Hz and tetani at 50 Hz. The differential effects of the studied compounds on twitches and tetani were statistically compared using simultaneous confidence intervals for the ratios between mean IC(50) for the block of twitches and mean IC(50) for the block of tetani. The results of ratios of mean IC(50) together with their corresponding 95% simultaneous confidence intervals were: vecuronium: 2.5 (1.8-3.5); mivacurium: 3.8 (3.0-4.9); pancuronium: 3.9 (2.0-7.6); rocuronium: 6.1 (3.8-9.9); atracurium: 9.0 (6.4-12.6); cisatracurium: 13.1 (6.0-28.4). Using the criteria that neuromuscular blockers displaying disjunct confidence intervals for the ratios of mean IC(50) differ statistically with regard to differential effects on twitches and tetani, significant differences in ratios of IC(50) were detected in the following cases: vecuronium vs. rocuronium, vs. atracurium and vs. cisatracurium and mivacurium vs: cisatracurium and vs. atracurium. The results show that the magnitude of the differential effects of neuromuscular blockers on twitches and tetani, as evaluated in the present work in the form of ratios of mean IC(50), does not depend on the chemical structure (comparing steroidal and isoquinolinic compounds), but seems to depend on differential pre- and post-synaptic effects of the compounds. It is also suggested that the greater the ability of a compound to block twitches and tetani in a differential manner, the safer is the compound from the clinical anesthesiology viewpoint.

Proinflammatory and Antiinflammatory Cytokine Levels in Complicated and Noncomplicated Parapneumonic Pleural Effusions

Objectives: This study aimed to evaluate a panel of proinflammatory and antiinflammatory cytokines in noncomplicated and complicated parapneumonic pleural effusions and to correlate their levels with pleural fluid biochemical parameters. Methods: Serum and pleural effusion were collected from 60 patients with noncomplicated (n = 26) or complicated (n = 34) parapneumonic effusions and assayed for cytologic, biochemical, and proinflammatory and antiinflammatory cytokines. Student t test was used to compare serum and pleural fluid values, Spearman correlation to analyze the relationship between pleural fluid cytokines and biochemical parameters, and accuracy of pleural fluid cytokine levels to determine the optimal cutoff value for identification of complicated effusions. Corrections for multiple comparisons were applied and a P value < .05 was accepted as significant. Results: Serum and pleural fluid cytokine levels of IL-8, vascular endothelial growth factor (VEGF), IL-10, and tumor necrosis factor (TNF) soluble receptor (sR) II were similar between groups. In contrast, complicated effusions had higher levels of pleural fluid IL-1 beta, IL-1 receptor antagonist (ra), and TNF sRI. Negative correlations were found between pleural fluid glucose with IL-1 beta and TNF sRI and positive correlations between lactic dehydrogenate (LDH) with IL-1 beta, IL-8, and VEGF. Pleural fluid levels of IL-1 beta, IL-1ra, and TNF sRI were more accurate than IL-8, VEGF, IL-10, and TNF sRII in discriminating complicated effusions. Conclusions: Both proinflammatory and antiinflammatory cytokine levels in pleural fluid are elevated in complicated in comparison with noncomplicated parapneumonic pleural effusions, and they correlate with both pleural fluid glucose and LDH levels. IL-1 beta, IL-1ra, and TNF sRI had higher sensitivity and specificity than IL-8, VEGF, IL-10, and TNF sRII in discriminating complicated effusions. CHEST 2012; 141( 1):183-189

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.