Statistics of selectively neutral genetic variation

Random models of evolution are instrumental in extracting rates of microscopic evolutionary mechanisms from empirical observations on genetic variation in genome sequences. In this context it is necessary to know the statistical properties of empirical observables (such as the local homozygosity for instance). Previous work relies on numerical results or assumes Gaussian approximations for the corresponding distributions. In this paper we give an analytical derivation of the statistical properties of the local homozygosity and other empirical observables assuming selective neutrality. We find that such distributions can be very non-Gaussian.Comment: 4 pages, 4 figure

Reversal Distances for Strings with Few Blocks or Small Alphabets

International audienceWe study the String Reversal Distance problem, an extension of the well-known Sorting by Reversals problem. String Reversal Distance takes two strings S and T as input, and asks for a minimum number of reversals to obtain T from S. We consider four variants: String Reversal Distance, String Prefix Reversal Distance (in which any reversal must include the first letter of the string), and the signed variants of these problems, namely Signed String Reversal Distance and Signed String Prefix Reversal Distance. We study algorithmic properties of these four problems, in connection with two parameters of the input strings: the number of blocks they contain (a block being maximal substring such that all letters in the substring are equal), and the alphabet size Σ. For instance, we show that Signed String Reversal Distance and Signed String Prefix Reversal Distance are NP-hard even if the input strings have only one letter

A model of weak selection in the infinite alleles framework

Ewens (1972) proposed a model in the infinite allele framework for populations with neutrality of all alleles at a particular locus. This paper proposes a generalisation of Ewens' result for situations where there is a form of weak selection. The models considered here are continuous time, discrete state space Markov processes.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/46945/1/285_2004_Article_BF00275643.pd

Evidence for Pervasive Adaptive Protein Evolution in Wild Mice

The relative contributions of neutral and adaptive substitutions to molecular evolution has been one of the most controversial issues in evolutionary biology for more than 40 years. The analysis of within-species nucleotide polymorphism and between-species divergence data supports a widespread role for adaptive protein evolution in certain taxa. For example, estimates of the proportion of adaptive amino acid substitutions (alpha) are 50% or more in enteric bacteria and Drosophila. In contrast, recent estimates of alpha for hominids have been at most 13%. Here, we estimate alpha for protein sequences of murid rodents based on nucleotide polymorphism data from multiple genes in a population of the house mouse subspecies Mus musculus castaneus, which inhabits the ancestral range of the Mus species complex and nucleotide divergence between M. m. castaneus and M. famulus or the rat. We estimate that 57% of amino acid substitutions in murids have been driven by positive selection. Hominids, therefore, are exceptional in having low apparent levels of adaptive protein evolution. The high frequency of adaptive amino acid substitutions in wild mice is consistent with their large effective population size, leading to effective natural selection at the molecular level. Effective natural selection also manifests itself as a paucity of effectively neutral nonsynonymous mutations in M. m. castaneus compared to humans

The sampling theory of neutral alleles and an urn model in population genetics

The behaviour of a Pólya-like urn which generates Ewens' sampling formula in population genetics is investigated. Connections are made with work of Watterson and Kingman and to the Poisson-Dirichlet distribution. The order in which novel types occur in the urn is shown to parallel the age distribution of the infinitely many alleles diffusion model and consequences of this property are explored. Finally the urn process is related to Kingman's coalescent with mutation to provide a rigorous basis for this parallel.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/46946/1/285_2004_Article_BF00276386.pd

Genomic signatures of population decline in the malaria mosquito Anopheles gambiae

Population genomic features such as nucleotide diversity and linkage disequilibrium are expected to be strongly shaped by changes in population size, and might therefore be useful for monitoring the success of a control campaign. In the Kilifi district of Kenya, there has been a marked decline in the abundance of the malaria vector Anopheles gambiae subsequent to the rollout of insecticide-treated bed nets. To investigate whether this decline left a detectable population genomic signature, simulations were performed to compare the effect of population crashes on nucleotide diversity, Tajima's D, and linkage disequilibrium (as measured by the population recombination parameter ρ). Linkage disequilibrium and ρ were estimated for An. gambiae from Kilifi, and compared them to values for Anopheles arabiensis and Anopheles merus at the same location, and for An. gambiae in a location 200 km from Kilifi. In the first simulations ρ changed more rapidly after a population crash than the other statistics, and therefore is a more sensitive indicator of recent population decline. In the empirical data, linkage disequilibrium extends 100-1000 times further, and ρ is 100-1000 times smaller, for the Kilifi population of An. gambiae than for any of the other populations. There were also significant runs of homozygosity in many of the individual An. gambiae mosquitoes from Kilifi. These results support the hypothesis that the recent decline in An. gambiae was driven by the rollout of bed nets. Measuring population genomic parameters in a small sample of individuals before, during and after vector or pest control may be a valuable method of tracking the effectiveness of interventions

Evidence of a recent decline in UK emissions of hydrofluorocarbons determined by the InTEM inverse model and atmospheric measurements

National greenhouse gas inventories (GHGIs) are submitted annually to the United Nations Framework Convention on Climate Change (UNFCCC). They are estimated in compliance with Intergovernmental Panel on Climate Change (IPCC) methodological guidance using activity data, emission factors and facility-level measurements. For some sources, the outputs from these calculations are very uncertain. Inverse modelling techniques that use high-quality, long-term measurements of atmospheric gases have been developed to provide independent verification of national GHGIs. This is considered good practice by the IPCC as it helps national inventory compilers to verify reported emissions and to reduce emission uncertainty. Emission estimates from the InTEM (Inversion Technique for Emission Modelling) model are presented for the UK for the hydrofluorocarbons (HFCs) reported to the UNFCCC (HFC-125, HFC-134a, HFC-143a, HFC-152a, HFC-23, HFC-32, HFC-227ea, HFC-245fa, HFC-43-10mee and HFC-365mfc). These HFCs have high global warming potentials (GWPs), and the global background mole fractions of all but two are increasing, thus highlighting their relevance to the climate and a need for increasing the accuracy of emission estimation for regulatory purposes. This study presents evidence that the long-term annual increase in growth of HFC-134a has stopped and is now decreasing. For HFC-32 there is an early indication, its rapid global growth period has ended, and there is evidence that the annual increase in global growth for HFC-125 has slowed from 2018. The inverse modelling results indicate that the UK implementation of European Union regulation of HFC emissions has been successful in initiating a decline in UK emissions from 2018. Comparison of the total InTEM UK HFC emissions in 2020 with the average from 2009–2012 shows a drop of 35 %, indicating progress toward the target of a 79 % decrease in sales by 2030. The total InTEM HFC emission estimates (2008–2018) are on average 73 (62–83) % of, or 4.3 (2.7–5.9) Tg CO2-eq yr−1 lower than, the total HFC emission estimates from the UK GHGI. There are also significant discrepancies between the two estimates for the individual HFCs.</p

Multichromosomal median and halving problems under different genomic distances

<p>Abstract</p> <p>Background</p> <p>Genome median and genome halving are combinatorial optimization problems that aim at reconstructing ancestral genomes as well as the evolutionary events leading from the ancestor to extant species. Exploring complexity issues is a first step towards devising efficient algorithms. The complexity of the median problem for unichromosomal genomes (permutations) has been settled for both the breakpoint distance and the reversal distance. Although the multichromosomal case has often been assumed to be a simple generalization of the unichromosomal case, it is also a relaxation so that complexity in this context does not follow from existing results, and is open for all distances.</p> <p>Results</p> <p>We settle here the complexity of several genome median and halving problems, including a surprising polynomial result for the breakpoint median and guided halving problems in genomes with circular and linear chromosomes, showing that the multichromosomal problem is actually easier than the unichromosomal problem. Still other variants of these problems are NP-complete, including the DCJ double distance problem, previously mentioned as an open question. We list the remaining open problems.</p> <p>Conclusion</p> <p>This theoretical study clears up a wide swathe of the algorithmical study of genome rearrangements with multiple multichromosomal genomes.</p

Measuring the effectiveness of in-hospital and on-base Prevent Alcohol and Risk-related Trauma in Youth (P.A.R.T.Y.) programs on reducing alcohol related harms in naval trainees: P.A.R.T.Y. Defence study protocol

Abstract Background Reducing alcohol related harms in Australian Defence Force (ADF) trainees has been identified as a priority, but there are few evidence-based prevention programs available for the military setting. The study aims to test whether the P.A.R.T.Y. program delivered in-hospital or on-base, can reduce harmful alcohol consumption among ADF trainees. Methods/design The study is a 3-arm randomized controlled trial, involving 953 Royal Australian Navy trainees from a single base. Trainees, aged 18 to 30 years, will be randomly assigned to the study arms: i. in-hospital P.A.R.T.Y.; ii. On-base P.A.R.T.Y.; and iii. Control group. All groups will receive the routine ADF annual alcohol awareness training. The primary outcome is the proportion of participants reporting an Alcohol Use Disorders Identification Test (AUDIT) score of 8 or above at 12 months’ post-intervention. The secondary outcome is the number of alcohol related incidents reported to the Royal Australian Navy (RAN) in the 12 months’ post-intervention. Discussion This is the first trial of the use of the P.A.R.T.Y. program in the military. If the proposed intervention proves efficacious, it may be a useful program in the early education of RAN trainees. Trial registration Australian New Zealand Clinical Trials Registry (ANZCTR): ACTRN12614001332617 , date of registration: 18/12/2014 ‘retrospectively registered’