Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study.

The identification of genetically homogeneous groups of individuals is a long standing issue in population genetics. A recent Bayesian algorithm implemented in the software STRUCTURE allows the identification of such groups. However, the ability of this algorithm to detect the true number of clusters (K) in a sample of individuals when patterns of dispersal among populations are not homogeneous has not been tested. The goal of this study is to carry out such tests, using various dispersal scenarios from data generated with an individual-based model. We found that in most cases the estimated 'log probability of data' does not provide a correct estimation of the number of clusters, K. However, using an ad hoc statistic DeltaK based on the rate of change in the log probability of data between successive K values, we found that STRUCTURE accurately detects the uppermost hierarchical level of structure for the scenarios we tested. As might be expected, the results are sensitive to the type of genetic marker used (AFLP vs. microsatellite), the number of loci scored, the number of populations sampled, and the number of individuals typed in each sample.

Reliability and validity of the Alcohol Use Disorders Identification Test (AUDIT) imbedded within a general health risk screening questionnaire: results of a survey in 332 primary care patients.

BACKGROUND: Self-administered, general health risk screening questionnaires that are administered while patients wait in the doctor's office may be a reasonable and timesaving approach to address the requirements of preventive medicine in a typical 10-min medical visit. The psychometric characteristics of the Alcohol Use Disorders Identification Test (AUDIT) incorporated within a health questionnaire (H-AUDIT) have not been examined. METHODS: The reliability and validity of the self-administered AUDIT were compared between the H-AUDIT and the AUDIT used as a single scale (S-AUDIT) in 332 primary care patients. RESULTS: No major demographic or alcohol use characteristics were found between the 166 subjects who completed the H-AUDIT and the 166 individuals who completed the S-AUDIT. The test-retest reliability of the 166 subjects who completed the H-AUDIT [estimated by Spearman correlation coefficient at a 6-week interval (0.88), internal consistency (total correlation coefficients for all items ranged from 0.38 to 0.69; Cronbach alpha index 0.85), and the sensitivity and specificity of the H-AUDIT were used to identify at-risk drinkers' areas under receiver operating characteristic (0.77) and alcohol-dependent subjects' areas under receiver operating characteristic (0.89)] was similar to the same measurements obtained with the 166 individuals who completed the S-AUDIT. CONCLUSIONS: The AUDIT incorporated in a health risk screening questionnaire is a reliable and valid self-administered instrument to identify at-risk drinkers and alcohol-dependent individuals in primary care settings.

An empirical hierarchical Bayesian unification of occupational exposure assessment methods.

In occupational exposure assessment of airborne contaminants, exposure levels can either be estimated through repeated measurements of the pollutant concentration in air, expert judgment or through exposure models that use information on the conditions of exposure as input. In this report, we propose an empirical hierarchical Bayesian model to unify these approaches. Prior to any measurement, the hygienist conducts an assessment to generate prior distributions of exposure determinants. Monte-Carlo samples from these distributions feed two level-2 models: a physical, two-compartment model, and a non-parametric, neural network model trained with existing exposure data. The outputs of these two models are weighted according to the expert's assessment of their relevance to yield predictive distributions of the long-term geometric mean and geometric standard deviation of the worker's exposure profile (level-1 model). Bayesian inferences are then drawn iteratively from subsequent measurements of worker exposure. Any traditional decision strategy based on a comparison with occupational exposure limits (e.g. mean exposure, exceedance strategies) can then be applied. Data on 82 workers exposed to 18 contaminants in 14 companies were used to validate the model with cross-validation techniques. A user-friendly program running the model is available upon request.

Multi-ethnic fine-mapping of 14 central adiposity loci.

The Genetic Investigation of Anthropometric Traits (GIANT) consortium identified 14 loci in European Ancestry (EA) individuals associated with waist-to-hip ratio (WHR) adjusted for body mass index. These loci are wide and narrowing the signals remains necessary. Twelve of 14 loci identified in GIANT EA samples retained strong associations with WHR in our joint EA/individuals of African Ancestry (AA) analysis (log-Bayes factor >6.1). Trans-ethnic analyses at five loci (TBX15-WARS2, LYPLAL1, ADAMTS9, LY86 and ITPR2-SSPN) substantially narrowed the signals to smaller sets of variants, some of which are in regions that have evidence of regulatory activity. By leveraging varying linkage disequilibrium structures across different populations, single-nucleotide polymorphisms (SNPs) with strong signals and narrower credible sets from trans-ethnic meta-analysis of central obesity provide more precise localizations of potential functional variants and suggest a possible regulatory role. Meta-analysis results for WHR were obtained from 77 167 EA participants from GIANT and 23 564 AA participants from the African Ancestry Anthropometry Genetics Consortium. For fine mapping we interrogated SNPs within ± 250 kb flanking regions of 14 previously reported index SNPs from loci discovered in EA populations by performing trans-ethnic meta-analysis of results from the EA and AA meta-analyses. We applied a Bayesian approach that leverages allelic heterogeneity across populations to combine meta-analysis results and aids in fine-mapping shared variants at these locations. We annotated variants using information from the ENCODE Consortium and Roadmap Epigenomics Project to prioritize variants for possible functionality.

Body temperature regulation and outcome after cardiac arrest and therapeutic hypothermia.

OBJECTIVE: Therapeutic temperature modulation is recommended after cardiac arrest (CA). However, body temperature (BT) regulation has not been extensively studied in this setting. We investigated BT variation in CA patients treated with therapeutic hypothermia (TH) and analyzed its impact on outcome. METHODS: A prospective cohort of comatose CA patients treated with TH (32-34°C, 24h) at the medical/surgical intensive care unit of the Lausanne University Hospital was studied. Spontaneous BT was recorded on hospital admission. The following variables were measured during and after TH: time to target temperature (TTT=time from hospital admission to induced BT target <34°C), cooling rate (spontaneous BT-induced BT target/TTT) and time of passive rewarming to normothermia. Associations of spontaneous and induced BT with in-hospital mortality were examined. RESULTS: A total of 177 patients (median age 61 years; median time to ROSC 25 min) were studied. Non-survivors (N=90, 51%) had lower spontaneous admission BT than survivors (median 34.5 [interquartile range 33.7-35.9]°C vs. 35.1 [34.4-35.8]°C, p=0.04). Accordingly, time to target temperature was shorter among non-survivors (200 [25-363]min vs. 270 [158-375]min, p=0.03); however, when adjusting for admission BT, cooling rates were comparable between the two outcome groups (0.4 [0.2-0.5]°C/h vs. 0.3 [0.2-0.4]°C/h, p=0.65). Longer duration of passive rewarming (600 [464-744]min vs. 479 [360-600]min, p<0.001) was associated with mortality. CONCLUSIONS: Lower spontaneous admission BT and longer time of passive rewarming were associated with in-hospital mortality after CA and TH. Impaired thermoregulation may be an important physiologic determinant of post-resuscitation disease and CA prognosis. When assessing the benefit of early cooling on outcome, future trials should adjust for patient admission temperature and use the cooling rate rather than the time to target temperature.

Associer et évaluer

Le forensicien joue fréquemment le rôle d'expert judiciaire. Dans cette responsabilité, il doit par exemple évaluer le lien qu'est susceptible d'entretenir une trace avec une source putative (une personne, un objet) donnée. L'expert examine aussi parfois la relation entre deux traces (proviennent-elles de la même source ?), combine les indices, ou interprète les activités possibles en regard des traces recueillies. Un mouvement dominant et intégrateur préconise d'effectuer cette évaluation au moyen de probabilités qui expriment dans un cadre formel les incertitudes inhérentes à ce genre de questions. Basé sur le théorème découvert par le révérend Bayes sur la probabilité des causes, le rapport de vraisemblance équilibre le point de vue de l'accusation (par exemple, probabilité de la trace, sachant que x en est à l'origine) et de la défense (la probabilité de la trace si la source n'est pas x). L'exploitation de ce modèle fait appel à l'expérience de l'expert, mais aussi à de grandes quantités de données représentatives qui doivent encore être rassemblées. Une tendance consiste aussi à exploiter pour cela les données stockées par l'utilisation de systèmes automatisés (empreintes digitales, projectiles, etc.). De nouvelles possibilités d'extraire des caractéristiques encore inexplorées à partir des traces sont susceptibles d'augmenter considérablement les moyens pour les experts judiciaires d'apporter des réponses probabilistes dans des circonstances toujours plus variées aux autorités judiciaires qui les mandatent.

The perception of touch and the ventral somatosensory pathway.

In humans, touching the skin is known to activate, among others, the contralateral primary somatosensory cortex on the postcentral gyrus together with the bilateral parietal operculum (i.e. the anatomical site of the secondary somatosensory cortex). But which brain regions beyond the postcentral gyrus specifically contribute to the perception of touch remains speculative. In this study we collected structural magnetic resonance imaging scans and neurological examination reports of patients with brain injuries or stroke in the left or right hemisphere, but not in the postcentral gyrus as the entry site of cortical somatosensory processing. Using voxel-based lesion-symptom mapping, we compared patients with impaired touch perception (i.e. hypoaesthesia) to patients without such touch impairments. Patients with hypoaesthesia as compared to control patients differed in one single brain cluster comprising the contralateral parietal operculum together with the anterior and posterior insular cortex, the putamen, as well as subcortical white matter connections reaching ventrally towards prefrontal structures. This finding confirms previous speculations on the 'ventral pathway of somatosensory perception' and causally links these brain structures to the perception of touch.

Role of the chemokine receptor CX3CR1 in the mobilization of phagocytic retinal microglial cells.

We recently showed that subretinal CX3CR1-dependent microglial cell (MC) accumulation may lead to age-related macular degeneration. The fate of MC after engulfing retinal debris is poorly understood. Severe photoreceptor degeneration was observed 40days after exposure to bright light in CX3CR1-deficient but not control mice, and more MCs accumulated in the subretinal space of the former than the latter. To study the fate of subretinal MCs in CX3CR1 competent animals, we used a dystrophic rat model in which abundant subretinal MC accumulation is observed secondary to retinal degeneration. In dystrophic rats, MCs containing rhodopsin or rod outer segment (ROS) debris were found outside the outer retina at sites suggesting choroidal and ciliary egress. In conclusion, our data indicate that MC accumulation at injury sites is independent of CX3CR1 and precedes photoreceptor degeneration. The ectopic presence of rhodopsin-positive MCs suggests that CX3CR1 participates in MC egress from the outer retina.

Primary percutaneous coronary intervention for unprotected left main disease in patients with acute ST-segment elevation myocardial infarction the AMIS (Acute Myocardial Infarction in Switzerland) plus registry experience.

OBJECTIVES: This study sought to assess outcomes in patients with ST-segment elevation myocardial infarction undergoing primary percutaneous coronary intervention (PCI) for unprotected left main (LM) disease. BACKGROUND: Limited data are available on outcomes in patients with ST-segment elevation myocardial infarction undergoing LM PCI. METHODS: Of 9,075 patients with ST-segment elevation myocardial infarction enrolled in the AMIS (Acute Myocardial Infarction in Switzerland) Plus registry between 2005 and June 30, 2010, 6,666 underwent primary PCI. Of them, 348 (5.2%; mean age: 63.5 ± 12.6 years) underwent LM PCI, either isolated (n = 208) or concomitant to PCI for other vessel segments (n = 140). They were compared with 6,318 patients (94.8%; mean age: 61.9 ± 12.5 years) undergoing PCI of non-LM vessel segments only. RESULTS: The LM patients had higher rates of cardiogenic shock (12.2% vs. 3.5%; p < 0.001), cardiac arrest (10.6% vs. 6.3%; p < 0.01), in-hospital mortality (10.9% vs. 3.8%; p < 0.001), and major adverse cardiac and cerebrovascular events (12.4% vs. 5.0%; p < 0.001) than non-LM PCI. Rates of mortality and major adverse cardiac and cerebrovascular events were highest for concurrent LM and non-LM PCI (17.9% and 18.6%, respectively), intermediate for isolated LM PCI (6.3% and 8.3%, respectively), and lowest for non-LM PCI (3.8% and 5.0%, respectively). Rates of mortality and major adverse cardiac and cerebrovascular events for LM PCI were higher than for non-LM multivessel PCI (10.9% vs. 4.9%, p < 0.001, and 12.4% vs. 6.4%, p < 0.001, respectively). LM disease independently predicted in-hospital death (odds ratio: 2.36; 95% confidence interval: 1.34 to 4.17; p = 0.003). CONCLUSIONS: Emergent LM PCI in the context of acute myocardial infarction, even including 12% cardiogenic shock, appears to have a remarkably high (89%) in-hospital survival. Concurrent LM and non-LM PCI has worse outcomes than isolated LM PCI.

Spatio-temporal trends of mortality in small areas of Southern Spain

Background: Most mortality atlases show static maps from count data aggregated over time. This procedure has several methodological problems and serious limitations for decision making in Public Health. The evaluation of health outcomes, including mortality, should be approached from a dynamic time perspective that is specific for each gender and age group. At the moment, researches in Spain do not provide a dynamic image of the population’s mortality status from a spatio-temporal point of view. The aim of this paper is to describe the spatial distribution of mortality from all causes in small areas of Andalusia (Southern Spain) and evolution over time from 1981 to 2006. Methods: A small-area ecological study was devised using the municipality as the unit for analysis. Two spatiotemporal hierarchical Bayesian models were estimated for each age group and gender. One of these was used to estimate the specific mortality rate, together with its time trends, and the other to estimate the specific rate ratio for each municipality compared with Spain as a whole. Results: More than 97% of the municipalities showed a diminishing or flat mortality trend in all gender and age groups. In 2006, over 95% of municipalities showed male and female mortality specific rates similar or significantly lower than Spanish rates for all age groups below 65. Systematically, municipalities in Western Andalusia showed significant male and female mortality excess from 1981 to 2006 only in age groups over 65. Conclusions: The study shows a dynamic geographical distribution of mortality, with a different pattern for each year, gender and age group. This information will contribute towards a reflection on the past, present and future of mortality in Andalusia.

Area deprivation and mortality in the provincial capital cities of Andalusia and Catalonia (Spain)

Objective: To study the linkage between material deprivation and mortality from all causes, for men and women separately, in the capital cities of the provinces in Andalusia and Catalonia (Spain). Methods: A small-area ecological study was devised using the census section as the unit for analysis. 188 983 Deaths occurring in the capital cities of the Andalusian provinces and 109 478 deaths recorded in the Catalan capital cities were examined. Principal components factorial analysis was used to devise a material deprivation index comprising the percentage of manual labourers, unemployment and illiteracy. A hierarchical Bayesian model was used to study the relationship between mortality and area deprivation. Main results: In most cities, results show an increased male mortality risk in the most deprived areas in relation to the least depressed. In Andalusia, the relative risks between the highest and lowest deprivation decile ranged from 1.24 (Malaga) to 1.40 (Granada), with 95% credibility intervals showing a significant excess risk. In Catalonia, relative risks ranged between 1.08 (Girona) and 1.50 (Tarragona). No evidence was found for an excess of female mortality in most deprived areas in either of the autonomous communities. Conclusions: Within cities, gender-related differences were revealed when deprivation was correlated geographically with mortality rates. These differences were found from an ecological perspective. Further research is needed in order to validate these results from an individual approach. The idea to be analysed is to identify those factors that explain these differences at an individual level.

Cancer mortality inequalities in urban areas: a Bayesian small area analysis in Spanish cities

Background Intra-urban inequalities in mortality have been infrequently analysed in European contexts. The aim of the present study was to analyse patterns of cancer mortality and their relationship with socioeconomic deprivation in small areas in 11 Spanish cities. Methods It is a cross-sectional ecological design using mortality data (years 1996-2003). Units of analysis were the census tracts. A deprivation index was calculated for each census tract. In order to control the variability in estimating the risk of dying we used Bayesian models. We present the RR of the census tract with the highest deprivation vs. the census tract with the lowest deprivation. Results In the case of men, socioeconomic inequalities are observed in total cancer mortality in all cities, except in Castellon, Cordoba and Vigo, while Barcelona (RR = 1.53 95%CI 1.42-1.67), Madrid (RR = 1.57 95%CI 1.49-1.65) and Seville (RR = 1.53 95%CI 1.36-1.74) present the greatest inequalities. In general Barcelona and Madrid, present inequalities for most types of cancer. Among women for total cancer mortality, inequalities have only been found in Barcelona and Zaragoza. The excess number of cancer deaths due to socioeconomic deprivation was 16,413 for men and 1,142 for women. Conclusion This study has analysed inequalities in cancer mortality in small areas of cities in Spain, not only relating this mortality with socioeconomic deprivation, but also calculating the excess mortality which may be attributed to such deprivation. This knowledge is particularly useful to determine which geographical areas in each city need intersectorial policies in order to promote a healthy environment.

Atlas interactivo de mortalidad de Andalucía (AIMA)

Until now, mortality atlases have been static. Most of them describe the geographical distribution of mortality using count data aggregated over time and standardized mortality rates. However, this methodology has several limitations. Count data aggregated over time produce a bias in the estimation of death rates. Moreover, this practice difficult the study of temporal changes in geographical distribution of mortality. On the other hand, using standardized mortality hamper to check differences in mortality among groups. The Interactive Mortality Atlas in Andalusia (AIMA) is an alternative to conventional static atlases. It is a dynamic Geographical Information System that allows visualizing in web-site more than 12.000 maps and 338.00 graphics related to the spatio-temporal distribution of the main death causes in Andalusia by age and sex groups from 1981. The objective of this paper is to describe the methods used for AIMA development, to show technical specifications and to present their interactivity. The system is available from the link products in www.demap.es. AIMA is the first interactive GIS that have been developed in Spain with these characteristics. Spatio-temporal Hierarchical Bayesian Models were used for statistical data analysis. The results were integrated into web-site using a PHP environment and a dynamic cartography in Flash. Thematic maps in AIMA demonstrate that the geographical distribution of mortality is dynamic, with differences among year, age and sex groups. The information nowadays provided by AIMA and the future updating will contribute to reflect on the past, the present and the future of population health in Andalusia.

Population pharmacokinetics and effects of efavirenz in patients with human immunodeficiency virus infection.

OBJECTIVE: The reverse transcriptase inhibitor efavirenz is currently used at a fixed dose of 600 mg/d. However, dosage individualization based on plasma concentration monitoring might be indicated. This study aimed to assess the efavirenz pharmacokinetic profile and interpatient versus intrapatient variability in patients who are positive for human immunodeficiency virus, to explore the relationship between drug exposure, efficacy, and central nervous system toxicity and to build up a Bayesian approach for dosage adaptation. METHODS: The population pharmacokinetic analysis was performed by use of NONMEM based on plasma samples from a cohort of unselected patients receiving efavirenz. With the use of a 1-compartment model with first-order absorption, the influence of demographic and clinical characteristics on oral clearance and oral volume of distribution was examined. The average drug exposure during 1 dosing interval was estimated for each patient and correlated with markers of efficacy and toxicity. The population kinetic parameters and the variabilities were integrated into a Bayesian equation for dosage adaptation based on a single plasma sample. RESULTS: Data from 235 patients with a total of 719 efavirenz concentrations were collected. Oral clearance was 9.4 L/h, oral volume of distribution was 252 L, and the absorption rate constant was 0.3 h(-1). Neither the demographic covariates evaluated nor the comedications showed a clinically significant influence on efavirenz pharmacokinetics. A large interpatient variability was found to affect efavirenz relative bioavailability (coefficient of variation, 54.6%), whereas the intrapatient variability was small (coefficient of variation, 26%). An inverse correlation between average drug exposure and viral load and a trend with central nervous system toxicity were detected. This enabled the derivation of a dosing adaptation strategy suitable to bring the average concentration into a therapeutic target from 1000 to 4000 microg/L to optimize viral load suppression and to minimize central nervous system toxicity. CONCLUSIONS: The high interpatient and low intrapatient variability values, as well as the potential relationship with markers of efficacy and toxicity, support the therapeutic drug monitoring of efavirenz. However, further evaluation is needed before individualization of an efavirenz dosage regimen based on routine drug level monitoring should be recommended for optimal patient management.

Municipal distribution of bladder cancer mortality in Spain: Possible role of mining and industry

BACKGROUND Spain shows the highest bladder cancer incidence rates in men among European countries. The most important risk factors are tobacco smoking and occupational exposure to a range of different chemical substances, such as aromatic amines. METHODS This paper describes the municipal distribution of bladder cancer mortality and attempts to "adjust" this spatial pattern for the prevalence of smokers, using the autoregressive spatial model proposed by Besag, York and Molliè, with relative risk of lung cancer mortality as a surrogate. RESULTS It has been possible to compile and ascertain the posterior distribution of relative risk for bladder cancer adjusted for lung cancer mortality, on the basis of a single Bayesian spatial model covering all of Spain's 8077 towns. Maps were plotted depicting smoothed relative risk (RR) estimates, and the distribution of the posterior probability of RR>1 by sex. Towns that registered the highest relative risks for both sexes were mostly located in the Provinces of Cadiz, Seville, Huelva, Barcelona and Almería. The highest-risk area in Barcelona Province corresponded to very specific municipal areas in the Bages district, e.g., Suría, Sallent, Balsareny, Manresa and Cardona. CONCLUSION Mining/industrial pollution and the risk entailed in certain occupational exposures could in part be dictating the pattern of municipal bladder cancer mortality in Spain. Population exposure to arsenic is a matter that calls for attention. It would be of great interest if the relationship between the chemical quality of drinking water and the frequency of bladder cancer could be studied.