ABCtoolbox: a versatile toolkit for approximate Bayesian computations

BACKGROUND: The estimation of demographic parameters from genetic data often requires the computation of likelihoods. However, the likelihood function is computationally intractable for many realistic evolutionary models, and the use of Bayesian inference has therefore been limited to very simple models. The situation changed recently with the advent of Approximate Bayesian Computation (ABC) algorithms allowing one to obtain parameter posterior distributions based on simulations not requiring likelihood computations. RESULTS: Here we present ABCtoolbox, a series of open source programs to perform Approximate Bayesian Computations (ABC). It implements various ABC algorithms including rejection sampling, MCMC without likelihood, a Particle-based sampler and ABC-GLM. ABCtoolbox is bundled with, but not limited to, a program that allows parameter inference in a population genetics context and the simultaneous use of different types of markers with different ploidy levels. In addition, ABCtoolbox can also interact with most simulation and summary statistics computation programs. The usability of the ABCtoolbox is demonstrated by inferring the evolutionary history of two evolutionary lineages of Microtus arvalis. Using nuclear microsatellites and mitochondrial sequence data in the same estimation procedure enabled us to infer sex-specific population sizes and migration rates and to find that males show smaller population sizes but much higher levels of migration than females. CONCLUSION: ABCtoolbox allows a user to perform all the necessary steps of a full ABC analysis, from parameter sampling from prior distributions, data simulations, computation of summary statistics, estimation of posterior distributions, model choice, validation of the estimation procedure, and visualization of the results

A high throughput genotyping approach reveals distinctive autosomal genetic signatures for European and Near Eastern wild boar

The lack of a Near Eastern genetic signature in modern European porcine breeds indicates that, although domestic pigs from the Fertile Crescent entered Europe during the Neolithic, they were completely replaced by their European counterparts in a short window of time. Whilst the absence of such genetic signature has been convincingly demonstrated at the mitochondrial level, variation at the autosomal genomes of European and Near Eastern Sus scrofa has not been compared yet. Herewith, we have explored the genetic relationships among 43 wild boar from Europe (N = 21), Near East (N = 19) and Korea (N = 3), and 40 Iberian (N = 16), Canarian (N = 4) and Mangalitza (N = 20) pigs by using a high throughput SNP genotyping platform. After data filtering, 37,167 autosomal SNPs were used to perform population genetics analyses. A multidimensional scaling plot based on genome-wide identity-by-state pairwise distances inferred with PLINK showed that Near Eastern and European wild boar populations are genetically differentiated. Maximum likelihood trees built with TreeMix supported this conclusion i.e. an early population split between Near Eastern and European Sus scrofa was observed. Moreover, analysis of the data with Structure evidenced that the sampled Iberian, Canarian and Mangalitza pigs did not carry any autosomal signature compatible with a Near Eastern ancestry, a finding that agrees well with previous mitochondrial studies

DNA repair, genome stability and cancer: a historical perspective

The multistep process of cancer progresses over many years. The prevention of mutations by DNA repair pathways led to an early appreciation of a role for repair in cancer avoidance. However, the broader role of the DNA damage response (DDR) emerged more slowly. In this Timeline article, we reflect on how our understanding of the steps leading to cancer developed, focusing on the role of the DDR. We also consider how our current knowledge can be exploited for cancer therapy

Thiacetazone, an Antitubercular Drug that Inhibits Cyclopropanation of Cell Wall Mycolic Acids in Mycobacteria

Background. Mycolic acids are a complex mixture of branched, long-chain fatty acids, representing key components of the highly hydrophobic mycobacterial cell wall. Pathogenic mycobacteria carry mycolic acid sub-types that contain cyclopropane rings. Double bonds at specific sites on mycolic acid precursors are modified by the action of cyclopropane mycolic acid synthases (CMASs). The latter belong to a family of S-adenosyl-methionine-dependent methyl transferases, of which several have been well studied in Mycobacterium tuberculosis, namely, MmaA1 through A4, PcaA and CmaA2. Cyclopropanated mycolic acids are key factors participating in cell envelope permeability, host immunomodulation and persistence of M. tuberculosis. While several antitubercular agents inhibit mycolic acid synthesis, to date, the CMASs have not been shown to be drug targets. Methodology/Principle Findings. We have employed various complementary approaches to show that the antitubercular drug, thiacetazone (TAC), and its chemical analogues, inhibit mycolic acid cyclopropanation. Dramatic changes in the content and ratio of mycolic acids in the vaccine strainMycobacterium bovis BCG, as well as in the related pathogenic speciesMycobacterium marinum were observed after treatment with the drugs. Combination of thin layer chromatography, mass spectrometry and Nuclear Magnetic Resonance (NMR) analyses of mycolic acids purified fromdrug-treated mycobacteria showed a significant loss of cyclopropanation in both the a- and oxygenated mycolate sub-types. Additionally, High-Resolution Magic Angle Spinning (HR-MAS) NMR analyses on whole cells was used to detect cell wall-associated mycolates and to quantify the cyclopropanation status of the cell envelope. Further, overexpression of cmaA2, mmaA2 or pcaA in mycobacteria partially reversed the effects of TAC and its analogue on mycolic acid cyclopropanation, suggesting that the drugs act directly on CMASs. Conclusions/Significance. This is a first report on them echanism of action of TAC, demonstrating the CMASs as its cellular targets in mycobacteria. The implications of this study may be important for the design of alternative strategies for tuberculosis treatment

DNA Repair in Human Pluripotent Stem Cells Is Distinct from That in Non-Pluripotent Human Cells

The potential for human disease treatment using human pluripotent stem cells, including embryonic stem cells and induced pluripotent stem cells (iPSCs), also carries the risk of added genomic instability. Genomic instability is most often linked to DNA repair deficiencies, which indicates that screening/characterization of possible repair deficiencies in pluripotent human stem cells should be a necessary step prior to their clinical and research use. In this study, a comparison of DNA repair pathways in pluripotent cells, as compared to those in non-pluripotent cells, demonstrated that DNA repair capacities of pluripotent cell lines were more heterogeneous than those of differentiated lines examined and were generally greater. Although pluripotent cells had high DNA repair capacities for nucleotide excision repair, we show that ultraviolet radiation at low fluxes induced an apoptotic response in these cells, while differentiated cells lacked response to this stimulus, and note that pluripotent cells had a similar apoptotic response to alkylating agent damage. This sensitivity of pluripotent cells to damage is notable since viable pluripotent cells exhibit less ultraviolet light-induced DNA damage than do differentiated cells that receive the same flux. In addition, the importance of screening pluripotent cells for DNA repair defects was highlighted by an iPSC line that demonstrated a normal spectral karyotype, but showed both microsatellite instability and reduced DNA repair capacities in three out of four DNA repair pathways examined. Together, these results demonstrate a need to evaluate DNA repair capacities in pluripotent cell lines, in order to characterize their genomic stability, prior to their pre-clinical and clinical use

Mycobacterium tuberculosis Lineage Influences Innate Immune Response and Virulence and Is Associated with Distinct Cell Envelope Lipid Profiles

The six major genetic lineages of Mycobacterium tuberculosis are strongly associated with specific geographical regions, but their relevance to bacterial virulence and the clinical consequences of infection are unclear. Previously, we found that in Vietnam, East Asian/Beijing and Indo-Oceanic strains were significantly more likely to cause disseminated tuberculosis with meningitis than those from the Euro-American lineage. To investigate this observation we characterised 7 East Asian/Beijing, 5 Indo-Oceanic and 6 Euro-American Vietnamese strains in bone-marrow-derived macrophages, dendritic cells and mice. East Asian/Beijing and Indo-Oceanic strains induced significantly more TNF-α and IL-1β from macrophages than the Euro-American strains, and East Asian/Beijing strains were detectable earlier in the blood of infected mice and grew faster in the lungs. We hypothesised that these differences were induced by lineage-specific variation in cell envelope lipids. Whole lipid extracts from East Asian/Beijing and Indo-Oceanic strains induced higher concentrations of TNF-α from macrophages than Euro-American lipids. The lipid extracts were fractionated and compared by thin layer chromatography to reveal a distinct pattern of lineage-associated profiles. A phthiotriol dimycocerosate was exclusively produced by East Asian/Beijing strains, but not the phenolic glycolipid previously associated with the hyper-virulent phenotype of some isolates of this lineage. All Indo-Oceanic strains produced a unique unidentified lipid, shown to be a phenolphthiocerol dimycocerosate dependent upon an intact pks15/1 for its production. This was described by Goren as the ‘attenuation indictor lipid’ more than 40 years ago, due to its association with less virulent strains from southern India. Mutation of pks15/1 in a representative Indo-Oceanic strain prevented phenolphthiocerol dimycocerosate synthesis, but did not alter macrophage cytokine induction. Our findings suggest that the early interactions between M. tuberculosis and host are determined by the lineage of the infecting strain; but we were unable to show these differences are driven by lineage-specific cell-surface expressed lipids

Evidence of Latitudinal Migration in Tri-colored Bats, Perimyotis subflavus

Background: Annual movements of tri-colored bats (Perimyotis subflavus) are poorly understood. While this species has been considered a regional migrant, some evidence suggests that it may undertake annual latitudinal migrations, similar to other long distance North American migratory bat species. Methodology/Principal Findings: We investigated migration in P. subflavus by conducting stable hydrogen isotope analyses of 184 museum specimen fur samples and comparing these results (dDfur) to published interpolated dD values of collection site growing season precipitation (dDprecip). Results suggest that the male molt period occurred between June 23 and October 16 and 33 % of males collected during the presumed non-molt period were south of their location of fur growth. For the same time period, 16 % of females were south of their location of fur growth and in general, had not travelled as far as migratory males. There were strong correlations between dDfur from the presumed molt period and both growing season dD precip (males – r 2 = 0.86; p,0.01; females – r 2 = 0.75; p,0.01), and latitude of collection (males – r 2 = 0.85; p,0.01; females – r 2 = 0.73; p,0.01). Most migrants were collected at the northern (.40uN; males and females) and southern (,35uN; males only) extents of the species ’ range. Conclusions/Significance: These results indicate a different pattern of migration for this species than previously documented, suggesting that some P. subflavus engage in annual latitudinal migrations and that migratory tendency varie

Ecological approaches in veterinary epidemiology: mapping the risk of bat-borne rabies using vegetation indices and night-time light satellite imagery

Rabies remains a disease of significant public health concern. In the Americas, bats are an important source of rabies for pets, livestock, and humans. For effective rabies control and prevention, identifying potential areas for disease occurrence is critical to guide future research, inform public health policies, and design interventions. To anticipate zoonotic infectious diseases distribution at coarse scale, veterinary epidemiology needs to advance via exploring current geographic ecology tools and data using a biological approach. We analyzed bat-borne rabies reports in Chile from 2002 to 2012 to establish associations between rabies occurrence and environmental factors to generate an ecological niche model (ENM). The main rabies reservoir in Chile is the bat species Tadarida brasiliensis; we mapped 726 occurrences of rabies virus variant AgV4 in this bat species and integrated them with contemporary Normalized Difference Vegetation Index (NDVI) data from the Moderate Resolution Imaging Spectroradiometer (MODIS). The correct prediction of areas with rabies in bats and the reliable anticipation of human rabies in our study illustrate the usefulness of ENM for mapping rabies and other zoonotic pathogens. Additionally, we highlight critical issues with selection of environmental variables, methods for model validation, and consideration of sampling bias. Indeed, models with weak or incorrect validation approaches should be interpreted with caution. In conclusion, ecological niche modeling applications for mapping disease risk at coarse geographic scales have a promising future, especially with refinement and enrichment of models with additional information, such as night-time light data, which increased substantially the model’s ability to anticipate human rabies