ABSTRACT
BACKGROUND: Nursing personnel are critical for enabling access to health service in primary health care. However, the State of the World's Nursing 2020 report showed important inequalities in nurse availability between countries. METHODS: The purpose of this study/analysis was to describe the differences in nurse-to-population density in 58 countries from six regional areas and the relationship between differences in access to nurses and other indicators of health equity. RESULTS: All countries and income groups showed subnational inequalities in the distribution of nursing personnel with Gini coefficients ranging from 1 to 39. The latter indicated situation such as 13% of the population having access to 45% of nurses in a country. The average max-to-min ratio was on average of 11-fold. In our sample, the African region had the highest level of subnational inequalities with the average Gini coefficient of 19.6. The European Region had the lowest level of within-country inequalities with the average Gini coefficient being 5.6. A multivariate analysis showed a clustering of countries in three groups: (1) high Gini coefficients comprised mainly African countries; (2) moderate Gini coefficients comprised mainly South-East Asian, Central and South American countries; (3) low Gini coefficients comprised mainly Western countries, Japan, and Korea. The analysis also showed that inequality in distribution of nurses was correlated with other indices of health and inequality such as the Human Development Index, maternal mortality, and life expectancy. CONCLUSIONS: This study showed that there is a high level of geographic inequality in the distribution of nurses at subnational level. Inequalities in nursing distribution are multifactorial, to improve access to nurses, policies should be bundled, tailored to the local context and tackle the various root causes for inequalities.
Subject(s)
Nurses , Nursing Staff , Africa , Humans , Income , Life Expectancy , Socioeconomic FactorsABSTRACT
BACKGROUND: Prospective cohort studies on diet and cancer report risk associations as hazard ratios. But hazard ratios do not inform on the number of people who need to alter their dietary behaviours for preventing cancer. The objective of this study is to estimate the number of people that need to alter their diet for preventing one additional case of female breast or colorectal cancer. METHODS: Based on the largest prospective studies done in the USA and in Europe, we computed the number of subjects who need to alter their diet. RESULTS: For preventing one case of breast cancer, European women should increase their fruit consumption by 100 g/day during 33 000 person-years, and US women by 60 g/day during 10 600 person-years. For vegetables, European women should increase their consumption by 160 g/day during 26 900 person-years and US women by 100 g/day during 19 000 person-years. For preventing one case of colorectal cancer, European subjects should decrease their red meat consumption by 20 g/day during 26 100 person-years, and US subjects by 30 g/day during 8170 person-years. For processed meat, European subjects should decrease their consumption by 20 g/day during 17 400 person-years, and US subjects by 10 g/day during 7940 person-years. CONCLUSIONS: Large number of subjects would need to alter their intake of fruits, vegetables, red and processed meat during many years in order to prevent one additional breast or colorectal cancer.
Subject(s)
Colorectal Neoplasms , Vegetables , Colorectal Neoplasms/epidemiology , Colorectal Neoplasms/prevention & control , Diet , Europe , Female , Fruit , Humans , Prospective Studies , Risk FactorsABSTRACT
The ecological impacts of increasing global temperatures are evident in most ecosystems on Earth, but our understanding of how climatic variation influences natural selection and adaptive resilience across latitudes remains largely unknown. Latitudinal gradients allow testing general ecosystem-level theories relevant to climatic adaptation. We assessed differences in adaptive diversity of populations along a latitudinal region spanning highly variable temperate to subtropical climates. We generated and integrated information from environmental mapping, phenotypic variation and genome-wide data from across the geographical range of the rainbowfish Melanotaenia duboulayi, an emerging aquatic system for studies of climate change. We detected, after controlling for spatial population structure, strong interactions between genotypes and environment associated with variation in stream flow and temperature. Some of these hydroclimate-associated genes were found to interact within functional protein networks that contain genes of adaptive significance for projected future climates in rainbowfish. Hydroclimatic selection was also associated with variation in phenotypic traits, including traits known to affect fitness of rainbowfish exposed to different flow environments. Consistent with predictions from the "climatic variability hypothesis," populations exposed to extremes of important environmental variables showed stronger adaptive divergence and less variation in climate-associated genes compared to populations at the centre of the environmental gradient. Our findings suggest that populations that evolved at environmental range margins and at geographical range edges may be more vulnerable to changing climates, a finding with implications for predicting adaptive resilience and managing biodiversity under climate change.
Subject(s)
Climate Change , Ecosystem , Animals , Fishes , Genotype , Selection, GeneticABSTRACT
Risk assessment for mixtures of chemicals requires to investigate the magnitude of their potential adverse effects on living organisms. This is usually done by assessing how experimental toxicological mixture data depart from the model of Loewe additivity. Several recent scientific studies propose to perform this task using an ad hoc method known as model deviation ratio (MDR) method. Moreover, the first official European regulatory document for the study of combined exposures explicitly recommends the use of the MDR method (EFSA Scientific Committee et al. Guidance on harmonised methodologies for human health, animal health and ecological risk assessment of combined exposure to multiple chemicals. EFSA Journal, 2019). We show here that the MDR method is not rooted in statistical principles and can lead to erroneous claims. We show however that the distribution of the MDR can be evaluated by simulations and show how this allows us to devise and carry out a bona fide statistical test. The proposed method accounts for uncertainty in the estimation of ED/EC50 and does not require a minimum sample size. The computer code developped in this study is made available as an R package called MDR.
Subject(s)
Complex Mixtures/toxicity , Dose-Response Relationship, Drug , Models, Statistical , Toxicology/methods , Animals , Humans , Risk Assessment , Toxicology/statistics & numerical dataABSTRACT
Bees are exposed to a wide range of multiple chemicals "chemical mixtures" from anthropogenic (e.g. plant protection products or veterinary products) or natural origin (e.g. mycotoxins, plant toxins). Quantifying the relative impact of multiple chemicals on bee health compared with other environmental stressors (e.g. varroa, viruses, and nutrition) has been identified as a priority to support the development of holistic risk assessment methods. Here, extensive literature searches and data collection of available laboratory studies on combined toxicity data for binary mixtures of pesticides and non-chemical stressors has been performed for honey bees (Apis mellifera), wild bees (Bombus spp.) and solitary bee species (Osmia spp.). From 957 screened publications, 14 publications provided 218 binary mixture toxicity data mostly for acute mortality (lethal dose: LD50) after contact exposure (61%), with fewer studies reporting chronic oral toxicity (20%) and acute oral LC50 values (19%). From the data collection, available dose response data for 92 binary mixtures were modelled using a Toxic Unit (TU) approach and the MIXTOX modelling tool to test assumptions of combined toxicity i.e. concentration addition (CA), and interactions (i.e. synergism, antagonism). The magnitude of interactions was quantified as the Model Deviation Ratio (MDR). The CA model applied to 17% of cases while synergism and antagonism were observed for 72% (MDRâ¯>â¯1.25) and 11% (MDRâ¯<â¯0.83) respectively. Most synergistic effects (55%) were observed as interactions between sterol-biosynthesis-inhibiting (SBI) fungicides and insecticide/acaricide. The mechanisms behind such synergistic effects of binary mixtures in bees are known to involve direct cytochrome P450 (CYP) inhibition, resulting in an increase in internal dose and toxicity of the binary mixture. Moreover, bees are known to have the lowest number of CYP copies and other detoxification enzymes in the insect kingdom. In the light of these findings, occurrence of these binary mixtures in relevant crops (frequency and concentrations) would need to be investigated. Addressing this exposure dimension remains critical to characterise the likelihood and plausibility of such interactions to occur under field realistic conditions. Finally, data gaps and further work for the development of risk assessment methods to assess multiple stressors in bees including chemicals and non-chemical stressors in bees are discussed.
Subject(s)
Bees , Fungicides, Industrial/toxicity , Pesticides/toxicity , Animals , Lethal Dose 50 , Risk AssessmentABSTRACT
In a cross-sectional field study involving 51 cattle herds in Belgium, 3159 serum samples and 557 individual milk samples were collected and tested by four different commercial antibody (Ab) ELISAs on serum and two Ab ELISAs on milk. A virus neutralization test (VNT) was performed on serum samples with discording ELISA results and on all samples from non-vaccinating herds. An epidemiological survey was carried out in the same herds to collect information about herd characteristics, management practices, BVD vaccination and BVD infection status. The objective of the study was to evaluate the performances of the Ab ELISAs relatively to the VNT, to assess the possibility of using pooled samples and to give recommendations regarding serological monitoring of BVD-free herds in the context of the Belgian national BVD eradication program which started early 2015. Depending on the assays, for ELISAs on serum, the diagnostic sensitivity (DSe) was estimated to be between 93.0 and 98.7% and the diagnostic specificity (DSp) between 94.3% and 99.1%. For the two ELISAs on milk, the DSe were 91.3% and 96.7% and the DSp 94.0% and 100% respectively and the Cohen's agreement coefficients between serum and milk samples were 0.75 and 0.85. Positive serum and milk samples diluted in negative samples to mimic different pool sizes were not detected by all ELISAs at dilutions above 1:5 or 1:10, leading to the conclusion that the testing of pooled samples should be used cautiously for serological monitoring and only with ELISAs with high sensitivity. The epidemiological analysis and the seroprevalence study, based on a general estimating equation model, showed that several factors had a significant influence on overall animal seroprevalence and within-herd seroprevalence such as age class, herd size, BVD herd infection status, BVD vaccination of young and/or adult cattle and the number of stables in the farm. This study showed that the best performances obtained with commercial Ab ELISAs are observed on individual serum samples, which should therefore be the preferred matrix to monitor BVD-free herds in the context of the Belgian eradication program. By regularly testing a limited number of samples from young (6-18 months) unvaccinated cattle it is possible to confirm the BVD-free herd status or to detect a recent infection.
Subject(s)
Antibodies, Viral/blood , Bovine Virus Diarrhea-Mucosal Disease/prevention & control , Diarrhea Viruses, Bovine Viral , Disease Eradication/methods , Enzyme-Linked Immunosorbent Assay/veterinary , Milk/virology , Animals , Antibodies, Viral/immunology , Belgium/epidemiology , Bovine Virus Diarrhea-Mucosal Disease/epidemiology , Bovine Virus Diarrhea-Mucosal Disease/virology , Cattle/blood , Cattle/virology , Cross-Sectional Studies , Female , Male , Neutralization Tests/veterinaryABSTRACT
Purging of solutions to enhance sonochemical reactions is a common practice. A fundamental study combining sonoluminescence spectroscopy and sonochemical activity is adopted to study the effects of continuous Ar gas flow in the solution and of the position of the gas inlet tube on high-frequency sonolysis of aqueous solutions. It has been observed that neither sonochemical activity nor sonoluminescence intensity is controlled by the gas solubility only. Besides, the change in position of the gas inlet tube leads to opposite effects in sonoluminescence intensity and sonochemical activity: while the former increases, the latter decreases. Such an observation has never been reported despite sonochemical reactions have been carried out under different gas environments. Sonoluminescence spectroscopy indicates that more extreme conditions are reached at collapse with the gas inlet on the side, which could be explained by a more symmetrical collapse. Finally, it is shown in certain conditions that it is possible to favor the formation of some sonochemical products simply by positioning the gas inlet at different positions, which has practical significance in designing large scale sonochemical reactors for industrial applications.
ABSTRACT
This guidance on the assessment of dermal absorption has been developed to assist notifiers, users of test facilities and Member State authorities on critical aspects related to the setting of dermal absorption values to be used in risk assessments of active substances in Plant Protection Products (PPPs). It is based on the 'scientific opinion on the science behind the revision of the guidance document on dermal absorption' issued in 2011 by the EFSA Panel on Plant Protection Products and their Residues (PPR). The guidance refers to the EFSA PPR opinion in many instances. In addition, the first version of this guidance, issued in 2012 by the EFSA PPR Panel, has been revised in 2017 on the basis of new available data on human in vitro dermal absorption for PPPs and wherever clarifications were needed. Basic details of experimental design, available in the respective test guidelines and accompanying guidance for the conduct of studies, have not been addressed but recommendations specific to performing and interpreting dermal absorption studies with PPPs are given. Issues discussed include a brief description of the skin and its properties affecting dermal absorption. To facilitate use of the guidance, flow charts are included. Guidance is also provided, for example, when there are no data on dermal absorption for the product under evaluation. Elements for a tiered approach are presented including use of default values, data on closely related products, in vitro studies with human skin (regarded to provide the best estimate), data from experimental animals (rats) in vitro and in vivo, and the so called 'triple pack' approach. Various elements of study design and reporting that reduce experimental variation and aid consistent interpretation are presented. A proposal for reporting data for assessment reports is also provided. The issue of nanoparticles in PPPs is not addressed. Data from volunteer studies have not been discussed since their use is not allowed in EU for risk assessment of PPPs.
ABSTRACT
MOTIVATION: Large-scale genotype datasets can help track the dispersal patterns of epidemiological outbreaks and predict the geographic origins of individuals. Such genetically-based geographic assignments also show a range of possible applications in forensics for profiling both victims and criminals, and in wildlife management, where poaching hotspot areas can be located. They, however, require fast and accurate statistical methods to handle the growing amount of genetic information made available from genotype arrays and next-generation sequencing technologies. RESULTS: We introduce a novel statistical method for geopositioning individuals of unknown origin from genotypes. Our method is based on a geostatistical model trained with a dataset of georeferenced genotypes. Statistical inference under this model can be implemented within the theoretical framework of Integrated Nested Laplace Approximation, which represents one of the major recent breakthroughs in statistics, as it does not require Monte Carlo simulations. We compare the performance of our method and an alternative method for geospatial inference, SPA in a simulation framework. We highlight the accuracy and limits of continuous spatial assignment methods at various scales by analyzing genotype datasets from a diversity of species, including Florida Scrub-jay birds Aphelocoma coerulescens, Arabidopsis thaliana and humans, representing 41-197,146 SNPs. Our method appears to be best suited for the analysis of medium-sized datasets (a few tens of thousands of loci), such as reduced-representation sequencing data that become increasingly available in ecology. AVAILABILITY AND IMPLEMENTATION: http://www2.imm.dtu.dk/â¼gigu/Spasiba/ CONTACT: gilles.b.guillot@gmail.com SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Subject(s)
Data Interpretation, Statistical , Genotype , High-Throughput Nucleotide Sequencing , Polymorphism, Single Nucleotide , Humans , Models, Theoretical , Monte Carlo MethodABSTRACT
Recognition of evolutionary units (species, populations) requires integrating several kinds of data, such as genetic or phenotypic markers or spatial information in order to get a comprehensive view concerning the differentiation of the units. We propose a statistical model with a double original advantage: (i) it incorporates information about the spatial distribution of the samples, with the aim to increase inference power and to relate more explicitly observed patterns to geography and (ii) it allows one to analyze genetic and phenotypic data within a unified model and inference framework, thus opening the way to robust comparisons between markers and possibly combined analyses. We show from simulated data as well as real data that our method estimates parameters accurately and is an improvement over alternative approaches in many situations. The power of this method is exemplified using an intricate case of inter- and intraspecies differentiation based on an original data set of georeferenced genetic and morphometric markers obtained on Myodes voles from Sweden. A computer program is made available as an extension of the R package Geneland.
Subject(s)
Genetic Variation , Models, Genetic , Phenotype , Animals , Arvicolinae/classification , Arvicolinae/genetics , Geography , Models, Statistical , Software , Species SpecificityABSTRACT
We propose a new model to make use of georeferenced genetic data for inferring the location and shape of a hybrid zone. The model output includes the posterior distribution of a parameter that quantifies the width of the hybrid zone. The model proposed is implemented in the GUI and command-line versions of the Geneland program versions ≥ 3.3.0. Information about the program can be found on http://www2.imm.dtu.dk/gigu/Geneland/.
Subject(s)
Algorithms , Demography , Hybridization, Genetic , Models, Genetic , Software , Markov Chains , Monte Carlo MethodABSTRACT
CONTEXT: Changes in energy intake have marked and rapid effects on metabolic functions, and some of these effects may be due to changes in adipocyte gene expression that precede alterations in body weight. OBJECTIVE: The aim of the study was to identify adipocyte genes regulated by changes in caloric intake independent of alterations in body weight. RESEARCH DESIGN AND METHODS: Obese subjects given a very low-caloric diet followed by gradual reintroduction of ordinary food and healthy subjects subjected to overfeeding were investigated. Adipose tissue biopsies were taken at multiple time-points, and gene expression was measured by DNA microarray. Genes regulated in the obese subjects undergoing caloric restriction followed by refeeding were identified using two-way ANOVA corrected with Bonferroni. From these, genes regulated by caloric restriction and oppositely during the weight-stable refeeding phase were identified in the obese subjects. The genes that were also regulated, in the same direction as the refeeding phase, in the healthy subjects after overfeeding were defined as being regulated by caloric intake. Results were confirmed using real-time PCR or immunoassay. RESULTS: Using a significance level of P < 0.05 for all comparisons, 52 genes were down-regulated, and 50 were up-regulated by caloric restriction and regulated in the opposite direction by refeeding and overfeeding. Among these were genes involved in lipogenesis (ACLY, ACACA, FASN, SCD), control of protein synthesis (4EBP1, 4EBP2), ß-oxidation (CPT1B), and insulin resistance (PEDF, SPARC). CONCLUSIONS: Metabolic genes involved in lipogenesis, protein synthesis, and insulin resistance are central in the transcriptional response of adipocytes to changes in caloric intake.
Subject(s)
Adipocytes/physiology , Body Weight/genetics , Body Weight/physiology , Energy Intake/genetics , Gene Expression Regulation/physiology , Adult , Biopsy , Caloric Restriction , Energy Intake/physiology , Eye Proteins/blood , Fast Foods , Female , Humans , Immunoassay , Insulin Resistance/genetics , Lipogenesis/genetics , Male , Microarray Analysis , Middle Aged , Nerve Growth Factors/blood , Obesity/diet therapy , Obesity/genetics , Oxidation-Reduction , Reverse Transcriptase Polymerase Chain Reaction , Serpins/bloodABSTRACT
The use of dominant markers such as amplified fragment length polymorphism (AFLP) for population genetics analyses is often impeded by the lack of appropriate computer programs and rarely motivated by objective considerations. The point of the present note is twofold: (i) we describe how the computer program Geneland designed to infer population structure has been adapted to deal with dominant markers; and (ii) we use Geneland for numerical comparison of dominant and codominant markers to perform clustering. AFLP markers lead to less accurate results than bi-allelic codominant markers such as single nucleotide polymorphisms (SNP) markers but this difference becomes negligible for data sets of common size (number of individuals n≥100, number of markers L≥200). The latest Geneland version (3.2.1) handling dominant markers is freely available as an R package with a fully clickable graphical interface. Installation instructions and documentation can be found on http://www2.imm.dtu.dk/~gigu/Geneland.
ABSTRACT
The joint analysis of spatial and genetic data is rapidly becoming the norm in population genetics. More and more studies explicitly describe and quantify the spatial organization of genetic variation and try to relate it to underlying ecological processes. As it has become increasingly difficult to keep abreast with the latest methodological developments, we review the statistical toolbox available to analyse population genetic data in a spatially explicit framework. We mostly focus on statistical concepts but also discuss practical aspects of the analytical methods, highlighting not only the potential of various approaches but also methodological pitfalls.
Subject(s)
Genetic Variation , Genetics, Population , Models, Genetic , Models, Statistical , Cluster AnalysisABSTRACT
MOTIVATION: In a series of recent papers, Tess, a computer program based on the concept of hidden Markov random field, has been proposed to infer the number and locations of panmictic population units from the genotypes and spatial locations of these individuals. The method seems to be of broad appeal as it is conceptually much simpler than other competing methods and it has been reported by its authors to be fast and accurate. However, this methodology is not grounded in a formal statistical inference method and seems to rely to a large extent on arbitrary choices regarding the parameters used. The present article is an investigation of the accuracy of this method and an attempt to assess whether recent results reported on the basis of this method are genuine features of the genetic process or artefacts of the method. METHOD: I analyse simulated data consisting of populations at Hardy-Weinberg and linkage equilibrium and also data simulated under a scenario of isolation-by-distance at mutation-migration-drift equilibrium. Arabidopsis thaliana data previously analysed with this method are also reconsidered. RESULTS: Using the Tess program under the no-admixture model to analyse data consisting of several genuine HWLE populations with individuals of pure ancestries leads to highly inaccurate results; Using the Tess program under the admixture model to analyse data consisting of a continuous isolation-by-distance population leads to the inference of spurious HWLE populations whose number and features depend on the parameters used. Results previously reported about the A.thaliana using Tess seem to a large extent to be artefacts of the statistical methodology used. The findings go beyond population clustering models and can be an help to design more efficient algorithms based on graphs. AVAILABILITY: The data analysed in the present article are available from http://folk.uio.no/gillesg/Bioinformatics-HMRF.
Subject(s)
Algorithms , Genetics, Population , Arabidopsis/genetics , Linkage Disequilibrium , Markov ChainsABSTRACT
BACKGROUND: The ascertainment process of molecular markers amounts to disregard loci carrying alleles with low frequencies. This can result in strong biases in inferences under population genetics models if not properly taken into account by the inference algorithm. Attempting to model this censoring process in view of making inference of population structure (i.e.identifying clusters of individuals) brings up challenging numerical difficulties. METHOD: These difficulties are related to the presence of intractable normalizing constants in Metropolis-Hastings acceptance ratios. This can be solved via an Markov chain Monte Carlo (MCMC) algorithm known as single variable exchange algorithm (SVEA). RESULT: We show how this general solution can be implemented for a class of clustering models of broad interest in population genetics that includes the models underlying the computer programs STRUCTURE, GENELAND and GESTE. We also implement the method proposed for a simple example and show that it allows us to reduce the bias substantially. AVAILABILITY: Further details and a computer program implementing the method are available from http://folk.uio.no/gillesg/AscB/.
Subject(s)
Algorithms , Computer Simulation , Genetics, Population/methods , Alleles , Gene Frequency , Markov Chains , Monte Carlo Method , SoftwareABSTRACT
Many models for inference of population genetic parameters are based on the assumption that the data set at hand consists of groups displaying within-group Hardy-Weinberg equilibrium at individual loci and linkage equilibrium between loci. This assumption is commonly violated by the presence of within-group spatial structure arising from nonrandom mating of individuals due to isolation by distance (IBD). This paper proposes a model and simulation method implemented in a computer program to flexibly simulate data displaying such patterns. The program permits displaying of smooth spatial variations of allele frequencies due to IBD and more abrupt variations due to presence of strong barriers to gene flow. It is useful in assessing performance of various statistical inference methods and in designing spatial sampling schemes. This is shown by a simulation study aimed at assessing the extent to which IBD patterns affect accuracy of cluster inferences performed in models assuming panmixia. The program is also used to study the effects of spatial sampling scheme (e.g. sampling individuals in clumps or uniformly across the spatial domain). The accuracy of such inferences is assessed in terms of number of inferred populations, assignment of individuals to populations and location of borders between populations. The effect of spatial sampling was weak while the effect of IBD may be substantial, leading to the inference of spurious populations, especially when IBD was strong with respect to the size of the sampling domain. The model and program are new and have been embedded in the R package Geneland, for user convenience and compliance with existing data formats.
ABSTRACT
MOTIVATION: This article considers the problem of estimating population genetic subdivision from multilocus genotype data. A model is considered to make use of genotypes and possibly of spatial coordinates of sampled individuals. A particular attention is paid to the case of low genetic differentiation with the help of a previously described Bayesian clustering model where allele frequencies are assumed to be a priori correlated. Under this model, various problems of inference are considered, in particular the common and difficult, but still unaddressed, situation where the number of populations is unknown. RESULTS: A Markov chain Monte Carlo algorithm and a new post-processing scheme are proposed. It is shown that they significantly improve the accuracy of previously existing algorithms in terms of estimated number of populations and estimated population membership. This is illustrated numerically with data simulated from the prior-likelihood model used in inference and also with data simulated from a Wright-Fisher model. Improvements are also illustrated on a real dataset of eighty-eight wolverines (Gulo gulo) genotyped at 10 microsatellites loci. The interest of the solutions presented here are not specific to any clustering model and are hence relevant to many settings in populations genetics where weakly differentiated populations are assumed or sought. AVAILABILITY: The improvements implemented will be made available in version 3.0.0 of the R package Geneland. Informations on how to get and use the software are available from http://folk.uio.no/gillesg/Geneland.html. SUPPLEMENTARY INFORMATION: http://folk.uio.no/gillesg/CFM/SuppMat.pdf.
Subject(s)
Algorithms , Gene Frequency , Genetic Drift , Animals , Genetics, Population , Genotype , Markov Chains , Microsatellite Repeats , Monte Carlo Method , Mustelidae/genetics , SoftwareABSTRACT
BACKGROUND: Despite several thousands of years of close contacts, there are genetic differences between the neighbouring countries of Finland and Sweden. Within Finland, signs of an east-west duality have been observed, whereas the population structure within Sweden has been suggested to be more subtle. With a fine-scale substructure like this, inferring the cluster membership of individuals requires a large number of markers. However, some studies have suggested that this number could be reduced if the individual spatial coordinates are taken into account in the analysis. RESULTS: We genotyped 34 unlinked autosomal single nucleotide polymorphisms (SNPs), originally designed for zygosity testing, from 2044 samples from Sweden and 657 samples from Finland, and 30 short tandem repeats (STRs) from 465 Finnish samples. We saw significant population structure within Finland but not between the countries or within Sweden, and isolation by distance within Finland and between the countries. In Sweden, we found a deficit of heterozygotes that we could explain by simulation studies to be due to both a small non-random genotyping error and hidden substructure caused by immigration. Geneland, a model-based Bayesian clustering algorithm, clustered the individuals into groups that corresponded to Sweden and Eastern and Western Finland when spatial coordinates were used, whereas in the absence of spatial information, only one cluster was inferred. CONCLUSION: We show that the power to cluster individuals based on their genetic similarity is increased when including information about the spatial coordinates. We also demonstrate the importance of estimating the size and effect of genotyping error in population genetics in order to strengthen the validity of the results.
Subject(s)
Genetic Linkage , Genetic Markers/genetics , Genetics, Population , Genotype , Polymorphism, Single Nucleotide , Adult , Analysis of Variance , Cluster Analysis , Finland , Founder Effect , Genetic Variation , Homozygote , Humans , Male , Middle Aged , Models, Genetic , Principal Component Analysis , SwedenABSTRACT
UNLABELLED: We introduce a new algorithm to account for the presence of null alleles in inferences of populations clusters from individual multilocus genetic data. We show by simulations that the presence of null alleles can affect the accuracy of inferences if not properly accounted for and that our algorithm improve signficantly their accuracy. AVAILABILITY: This new algorithm is implemented in the program Geneland. It is freely available under GNU public license as an R package on the Comprehensive R Archive Network. It now includes a fully clickable graphical interface. Informations on how to get the software are available on folk.uio.no/gillesg/Geneland.html