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

Haplotype Structure of FSHB, the Beta-Subunit Gene for Fertility-Associated Follicle-Stimulating Hormone: Possible Influence of Balancing Selection

Follicle-stimulating hormone (FSH) is essential for human reproduction. The unique functions of this hormone are provided by the FSH receptor-binding beta-subunit encoded by the FSHB gene. Resequencing and genotyping of FSHB in three European, two Asian and one African population, as well as in the great apes (chimpanzee, gorilla, orangutan), revealed low diversity and significant excess of polymorphisms with intermediate frequency alleles. Statistical tests for FSHB showed deviations from neutrality in all populations suggesting a possible effect of balancing selection. Two core haplotypes were identified (carried by 76-96.6% of each population's sample), the sequences of which are clearly separated from each other. As fertility most directly affects an organism's fitness, the carriers of these haplotypes have apparently had more success in human history to contribute to the next generation. There is a preliminary observation suggesting that the second most frequent FSHB haplotype may be associated with rapid conception success in females. Interestingly, the same haplotype is related to an ancestral FSHB variant shared with the ancestor of the great apes. The determination of the functional consequence of the two core FSHB variants may have implications for understanding and regulating human fertility, as well as in assisting infertility treatments

Dependence of paracentric inversion rate on tract length

BACKGROUND: We develop a Bayesian method based on MCMC for estimating the relative rates of pericentric and paracentric inversions from marker data from two species. The method also allows estimation of the distribution of inversion tract lengths. RESULTS: We apply the method to data from Drosophila melanogaster and D. yakuba. We find that pericentric inversions occur at a much lower rate compared to paracentric inversions. The average paracentric inversion tract length is approx. 4.8 Mb with small inversions being more frequent than large inversions. If the two breakpoints defining a paracentric inversion tract are uniformly and independently distributed over chromosome arms there will be more short tract-length inversions than long; we find an even greater preponderance of short tract lengths than this would predict. Thus there appears to be a correlation between the positions of breakpoints which favors shorter tract lengths. CONCLUSION: The method developed in this paper provides the first statistical estimator for estimating the distribution of inversion tract lengths from marker data. Application of this method for a number of data sets may help elucidate the relationship between the length of an inversion and the chance that it will get accepted

Fluoromycobacteriophages for rapid, specific, and sensitive antibiotic susceptibility testing of Mycobacterium tuberculosis

Rapid antibiotic susceptibility testing of Mycobacterium tuberculosis is of paramount importance as multiple- and extensively- drug resistant strains of M. tuberculosis emerge and spread. We describe here a virus-based assay in which fluoromycobacteriophages are used to deliver a GFP or ZsYellow fluorescent marker gene to M. tuberculosis, which can then be monitored by fluorescent detection approaches including fluorescent microscopy and flow cytometry. Pre-clinical evaluations show that addition of either Rifampicin or Streptomycin at the time of phage addition obliterates fluorescence in susceptible cells but not in isogenic resistant bacteria enabling drug sensitivity determination in less than 24 hours. Detection requires no substrate addition, fewer than 100 cells can be identified, and resistant bacteria can be detected within mixed populations. Fluorescence withstands fixation by paraformaldehyde providing enhanced biosafety for testing MDR-TB and XDR-TB infections. © 2009 Piuri et al

Regulation of cell survival by sphingosine-1-phosphate receptor S1P1 via reciprocal ERK-dependent suppression of bim and PI-3-kinase/protein kinase C-mediated upregulation of Mcl-1

Although the ability of bioactive lipid sphingosine-1-phosphate (S1P) to positively regulate anti-apoptotic/pro-survival responses by binding to S1P1 is well known, the molecular mechanisms remain unclear. Here we demonstrate that expression of S1P1 renders CCL39 lung fibroblasts resistant to apoptosis following growth factor withdrawal. Resistance to apoptosis was associated with attenuated accumulation of pro-apoptotic BH3-only protein Bim. However, although blockade of extracellular signal-regulated kinase (ERK) activation could reverse S1P1-mediated suppression of Bim accumulation, inhibition of caspase-3 cleavage was unaffected. Instead S1P1-mediated inhibition of caspase-3 cleavage was reversed by inhibition of phosphatidylinositol-3-kinase (PI3K) and protein kinase C (PKC), which had no effect on S1P1 regulation of Bim. However, S1P1 suppression of caspase-3 was associated with increased expression of anti-apoptotic protein Mcl-1, the expression of which was also reduced by inhibition of PI3K and PKC. A role for the induction of Mcl-1 in regulating endogenous S1P receptor-dependent pro-survival responses in human umbilical vein endothelial cells was confirmed using S1P receptor agonist FTY720-phosphate (FTY720P). FTY720P induced a transient accumulation of Mcl-1 that was associated with a delayed onset of caspase-3 cleavage following growth factor withdrawal, whereas Mcl-1 knockdown was sufficient to enhance caspase-3 cleavage even in the presence of FTY720P. Consistent with a pro-survival role of S1P1 in disease, analysis of tissue microarrays from ER+ breast cancer patients revealed a significant correlation between S1P1 expression and tumour cell survival. In these tumours, S1P1 expression and cancer cell survival were correlated with increased activation of ERK, but not the PI3K/PKB pathway. In summary, pro-survival/anti-apoptotic signalling from S1P1 is intimately linked to its ability to promote the accumulation of pro-survival protein Mcl-1 and downregulation of pro-apoptotic BH3-only protein Bim via distinct signalling pathways. However, the functional importance of each pathway is dependent on the specific cellular context

Genetic Variation in Native Americans, Inferred from Latino SNP and Resequencing Data

Analyses of genetic polymorphism data have the potential to be highly informative about the demographic history of Native American populations, but due to a combination of historical and political factors, there are essentially no autosomal sequence polymorphism data from any Native American group. However, there are many resequencing studies involving Latinos, whose genomes contain segments inherited from their Native American ancestors. In this study, we introduce a new method for estimating local ancestry across the genomes of admixed individuals and show how this method, along with dense genotyping and targeted resequencing, can be used to assay genetic variation in ancestral Native American groups. We analyze roughly 6 Mb of resequencing data from 22 Mexican Americans to provide the first large-scale view of sequence level variation in Native Americans. We observe low levels of diversity and high levels of linkage disequilibrium in the Native American–derived sequences, consistent with a recent severe population bottleneck associated with the initial peopling of the Americas. Using two different computational approaches, one novel, we estimate that this bottleneck occurred roughly 12.5 Kya; when uncertainty in the estimation process is taken into account, our results are consistent with archeological estimates for the colonization of the Americas

Mitochondrial phylogeography and demographic history of the Vicuña: implications for conservation

The vicuña (Vicugna vicugna; Miller, 1924) is a conservation success story, having recovered from near extinction in the 1960s to current population levels estimated at 275 000. However, lack of information about its demographic history and genetic diversity has limited both our understanding of its recovery and the development of science-based conservation measures. To examine the evolution and recent demographic history of the vicuña across its current range and to assess its genetic variation and population structure, we sequenced mitochondrial DNA from the control region (CR) for 261 individuals from 29 populations across Peru, Chile and Argentina. Our results suggest that populations currently designated as Vicugna vicugna vicugna and Vicugna vicugna mensalis comprise separate mitochondrial lineages. The current population distribution appears to be the result of a recent demographic expansion associated with the last major glacial event of the Pleistocene in the northern (18 to 22°S) dry Andes 14–12 000 years ago and the establishment of an extremely arid belt known as the 'Dry Diagonal' to 29°S. Within the Dry Diagonal, small populations of V. v. vicugna appear to have survived showing the genetic signature of demographic isolation, whereas to the north V. v. mensalis populations underwent a rapid demographic expansion before recent anthropogenic impacts

A Dynamic Model of Interactions of Ca^(2+), Calmodulin, and Catalytic Subunits of Ca^(2+)/Calmodulin-Dependent Protein Kinase II

During the acquisition of memories, influx of Ca^(2+) into the postsynaptic spine through the pores of activated N-methyl-D-aspartate-type glutamate receptors triggers processes that change the strength of excitatory synapses. The pattern of Ca^(2+) influx during the first few seconds of activity is interpreted within the Ca^(2+)-dependent signaling network such that synaptic strength is eventually either potentiated or depressed. Many of the critical signaling enzymes that control synaptic plasticity, including Ca^(2+)/calmodulin-dependent protein kinase II (CaMKII), are regulated by calmodulin, a small protein that can bind up to 4 Ca^(2+) ions. As a first step toward clarifying how the Ca^(2+)-signaling network decides between potentiation or depression, we have created a kinetic model of the interactions of Ca^(2+), calmodulin, and CaMKII that represents our best understanding of the dynamics of these interactions under conditions that resemble those in a postsynaptic spine. We constrained parameters of the model from data in the literature, or from our own measurements, and then predicted time courses of activation and autophosphorylation of CaMKII under a variety of conditions. Simulations showed that species of calmodulin with fewer than four bound Ca^(2+) play a significant role in activation of CaMKII in the physiological regime, supporting the notion that processing ofCa^(2+) signals in a spine involves competition among target enzymes for binding to unsaturated species of CaM in an environment in which the concentration of Ca^(2+) is fluctuating rapidly. Indeed, we showed that dependence of activation on the frequency of Ca^(2+) transients arises from the kinetics of interaction of fluctuating Ca^(2+) with calmodulin/CaMKII complexes. We used parameter sensitivity analysis to identify which parameters will be most beneficial to measure more carefully to improve the accuracy of predictions. This model provides a quantitative base from which to build more complex dynamic models of postsynaptic signal transduction during learning

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