Biovolatilisation: a poorly studied pathway of the arsenic biogeochemical cycle

It has been known for over a hundred years that microorganisms can produce volatile arsenic (As) species, termed “arsines”. However, this topic has received relatively little attention compared to As behaviour in soils and biotransformation through the trophic level in the marine and terrestrial environment. We believe this is due to long-standing misconceptions regarding volatile As stability and transport as well as an absence, until recently, of appropriate sampling methods. First and foremost, an attempt is made to unify arsines' designations, notations and formulas, taking into account all the different terms used in the literature. Then, the stability of As volatile species is discussed and new analytical developments are explored. Further, the special cases of diffuse low-level emissions (e.g. soil and sediment biovolatilisation), and point sources with high-level emissions (geothermal environments, landfills, and natural gas) are comprehensively reviewed. In each case, future possible areas of research and unknown mechanisms are identified and their importance towards the global As biogeochemical cycle is explored. This review gathers new information regarding mechanisms, stability, transport and sampling of the very elusive arsines and shows that more research should be conducted on this important process.

Herbrand Theorems for Substructural Logics

Herbrand and Skolemization theorems are obtained for a broad family of first-order substructural logics. These logics typically lack equivalent prenex forms, a deduction theorem, and reductions of semantic consequence to satisfiability. The Herbrand and Skolemization theorems therefore take various forms, applying either to the left or right of the consequence relation, and to restricted classes of formulas.

New perspectives in long QT syndrome

Long QT Syndrome (LQTS) is a cardiac channelopathy characterized by prolonged ventricular repolarization and increased risk to sudden death secondary to ventricular dysrrhythmias. Was the first cardiac channelopathy described and is probably the best understood. After a decade of the sentinel identification of ion channel mutation in LQTS, genotype-phenotype correlations have been developed along with important improvement in risk stratification and genetic guided-treatment. Genetic screening has shown that LQTS is more frequent than expected and interestingly, ethnic specific polymorphism conferring increased susceptibility to drug induced QT prolongation and torsades de pointes have been identified. A better understanding of ventricular arrhythmias as an adverse effect of ion channel binding drugs, allow the development of more safety formulas and better control of this public health problem. Progress in understanding the molecular basis of LQTS has been remarkable; eight different genes have been identified, however still 25% of patients remain genotype-negative. This article is an overview of the main LQTS knowledge developed during the last years.

Evaluating statistical methods used to estimate the number of postsynaptic receptors.

Calcium levels in spines play a significant role in determining the sign and magnitude of synaptic plasticity. The magnitude of calcium influx into spines is highly dependent on influx through N-methyl D-aspartate (NMDA) receptors, and therefore depends on the number of postsynaptic NMDA receptors in each spine. We have calculated previously how the number of postsynaptic NMDA receptors determines the mean and variance of calcium transients in the postsynaptic density, and how this alters the shape of plasticity curves. However, the number of postsynaptic NMDA receptors in the postsynaptic density is not well known. Anatomical methods for estimating the number of NMDA receptors produce estimates that are very different than those produced by physiological techniques. The physiological techniques are based on the statistics of synaptic transmission and it is difficult to experimentally estimate their precision. In this paper we use stochastic simulations in order to test the validity of a physiological estimation technique based on failure analysis. We find that the method is likely to underestimate the number of postsynaptic NMDA receptors, explain the source of the error, and re-derive a more precise estimation technique. We also show that the original failure analysis as well as our improved formulas are not robust to small estimation errors in key parameters.

Dilution of candidates: the case of iron-related genes in restless legs syndrome

Restless legs syndrome (RLS) is a common multifactorial disease. Some genetic risk factors have been identified. RLS susceptibility also has been related to iron. We therefore asked whether known iron-related genes are candidates for association with RLS and, vice versa, whether known RLS-associated loci influence iron parameters in serum. RLS/control samples (n = 954/1814 in the discovery step, 735/736 in replication 1, and 736/735 in replication 2) were tested for association with SNPs located within 4 Mb intervals surrounding each gene from a list of 111 iron-related genes using a discovery threshold of P = 5 × 10(-4). Two population cohorts (KORA F3 and F4 with together n = 3447) were tested for association of six known RLS loci with iron, ferritin, transferrin, transferrin-saturation, and soluble transferrin receptor. Results were negative. None of the candidate SNPs at the iron-related gene loci was confirmed significantly. An intronic SNP, rs2576036, of KATNAL2 at 18q21.1 was significant in the first (P = 0.00085) but not in the second replication step (joint nominal P-value = 0.044). Especially, rs1800652 (C282Y) in the HFE gene did not associate with RLS. Moreover, SNPs at the known RLS loci did not significantly affect serum iron parameters in the KORA cohorts. In conclusion, the correlation between RLS and iron parameters in serum may be weaker than assumed. Moreover, in a general power analysis, we show that genetic effects are diluted if they are transmitted via an intermediate trait to an end-phenotype. Sample size formulas are provided for small effect sizes.

Models of DNA sequence evolution

Models of DNA sequence evolution and methods for estimating evolutionary distances are needed for studying the rate and pattern of molecular evolution and for inferring the evolutionary relationships of organisms or genes. In this dissertation, several new models and methods are developed.^ The rate variation among nucleotide sites: To obtain unbiased estimates of evolutionary distances, the rate heterogeneity among nucleotide sites of a gene should be considered. Commonly, it is assumed that the substitution rate varies among sites according to a gamma distribution (gamma model) or, more generally, an invariant+gamma model which includes some invariable sites. A maximum likelihood (ML) approach was developed for estimating the shape parameter of the gamma distribution $(\alpha)$ and/or the proportion of invariable sites $(\theta).$ Computer simulation showed that (1) under the gamma model, $\alpha$ can be well estimated from 3 or 4 sequences if the sequence length is long; and (2) the distance estimate is unbiased and robust against violations of the assumptions of the invariant+gamma model.^ However, this ML method requires a huge amount of computational time and is useful only for less than 6 sequences. Therefore, I developed a fast method for estimating $\alpha,$ which is easy to implement and requires no knowledge of tree. A computer program was developed for estimating $\alpha$ and evolutionary distances, which can handle the number of sequences as large as 30.^ Evolutionary distances under the stationary, time-reversible (SR) model: The SR model is a general model of nucleotide substitution, which assumes (i) stationary nucleotide frequencies and (ii) time-reversibility. It can be extended to SRV model which allows rate variation among sites. I developed a method for estimating the distance under the SR or SRV model, as well as the variance-covariance matrix of distances. Computer simulation showed that the SR method is better than a simpler method when the sequence length $L>1,000$ bp and is robust against deviations from time-reversibility. As expected, when the rate varies among sites, the SRV method is much better than the SR method.^ The evolutionary distances under nonstationary nucleotide frequencies: The statistical properties of the paralinear and LogDet distances under nonstationary nucleotide frequencies were studied. First, I developed formulas for correcting the estimation biases of the paralinear and LogDet distances. The performances of these formulas and the formulas for sampling variances were examined by computer simulation. Second, I developed a method for estimating the variance-covariance matrix of the paralinear distance, so that statistical tests of phylogenies can be conducted when the nucleotide frequencies are nonstationary. Third, a new method for testing the molecular clock hypothesis was developed in the nonstationary case. ^

A nonparametric method of comparing the partial areas under two correlated receiver operating characteristic curves

A non-parametric method was developed and tested to compare the partial areas under two correlated Receiver Operating Characteristic curves. Based on the theory of generalized U-statistics the mathematical formulas have been derived for computing ROC area, and the variance and covariance between the portions of two ROC curves. A practical SAS application also has been developed to facilitate the calculations. The accuracy of the non-parametric method was evaluated by comparing it to other methods. By applying our method to the data from a published ROC analysis of CT image, our results are very close to theirs. A hypothetical example was used to demonstrate the effects of two crossed ROC curves. The two ROC areas are the same. However each portion of the area between two ROC curves were found to be significantly different by the partial ROC curve analysis. For computation of ROC curves with large scales, such as a logistic regression model, we applied our method to the breast cancer study with Medicare claims data. It yielded the same ROC area computation as the SAS Logistic procedure. Our method also provides an alternative to the global summary of ROC area comparison by directly comparing the true-positive rates for two regression models and by determining the range of false-positive values where the models differ. ^

Prescriptions of Chinese medicinal herbs in Switzerland: the example of suan zao ren (Ziziphi Spinosae Semen)

Introduction According to the Swiss Health Survey 2007, 1.7% of the adult population use traditional Chinese medicine (including Chinese herbal medicine, but excluding acupuncture). In contrast to conventional drugs, that contain single chemically defined substances, prescriptions of Chinese herbs are mixtures of up to 40 ingredients (parts of plants, fungi, animal substances and minerals). Originally they were taken in the form of decoctions, but nowadays granules are more popular. Medium daily dosages of granules range between 8 to 12g. In a recent work we identified the most commonly used Chinese herbs (all ingredients are referred to as herbs for reasons of simplicity) and classical formulas (mixtures). Here we present a short overview and the example of suan zao ren (Ziziphi Spinosae Semen), which is used in the treatment of insomnia and anxiety and contains saponins that have been shown to increase sleep in animal studies. Material and Methods A random sample of 1,053 prescriptions was drawn from the database of Lian Chinaherb AG, Switzerland, and analysed according to the most frequently used individual herbs and classical formulas. Cluster analysis (Jaccard similarity coefficient, complete linkage method) was applied to identify common combinations of herbs. Results The most frequently used herbs were dang gui (Angelicae Sinensis Radix), fu ling (Poria), bai shao (Paeoniae Radix Alba), and gan cao (Glycyrrhizae Radix et Rhizoma); the most frequently used classical formulas were gui pi tang (Restore the Spleen Decoction) and xiao yao san (Rambling Powder). The average number of herbs per prescription was 12.0, and the average daily dosage of granules was 8.7g. 74.3% of the prescriptions were for female, 24.8% for male patients. Suan zao ren was present in 14.2% of all prescriptions. These prescriptions contained on average 13.7 herbs, and the daily dosage of granules was 8.9g. Suan zao ren was more frequently prescribed by practitioners of non-Asian than of Asian origin but equally often for female and male patients. Cluster analysis grouped suan zao ren with yuan zhi (Polygalae Radix), bai zi ren (Platycladi Semen), sheng di huang (Rehmanniae Radix) and dan shen (Salviae Miltiorrhizae Radix et Rhizoma). Discussion Prescriptions including suan zao ren contained on average slightly more herbs than other prescriptions. This might be due to the fact that two of the three most popular classical formulas with suan zao ren are composed of 13 and 12 herbs with the possibility of adding more ingredients when necessary. Cluster analysis resulted in the clustering of suan zao ren with other herbs of the classical formula tian wang bu xin dan (Emperor of Heaven’s Special Pill to Tonify the Heart), indicating the use of suan zao ren for the treatment of insomnia and irritability. Unfortunately, the diagnoses of the patients were unavailable and thus correlations between use of suan zao ren and diseases could not be analysed.

Therapeutic drug monitoring of once daily aminoglycoside dosing: comparison of two methods and investigation of the optimal blood sampling strategy.

PURPOSE Therapeutic drug monitoring of patients receiving once daily aminoglycoside therapy can be performed using pharmacokinetic (PK) formulas or Bayesian calculations. While these methods produced comparable results, their performance has never been checked against full PK profiles. We performed a PK study in order to compare both methods and to determine the best time-points to estimate AUC0-24 and peak concentrations (C max). METHODS We obtained full PK profiles in 14 patients receiving a once daily aminoglycoside therapy. PK parameters were calculated with PKSolver using non-compartmental methods. The calculated PK parameters were then compared with parameters estimated using an algorithm based on two serum concentrations (two-point method) or the software TCIWorks (Bayesian method). RESULTS For tobramycin and gentamicin, AUC0-24 and C max could be reliably estimated using a first serum concentration obtained at 1 h and a second one between 8 and 10 h after start of the infusion. The two-point and the Bayesian method produced similar results. For amikacin, AUC0-24 could reliably be estimated by both methods. C max was underestimated by 10-20% by the two-point method and by up to 30% with a large variation by the Bayesian method. CONCLUSIONS The ideal time-points for therapeutic drug monitoring of once daily administered aminoglycosides are 1 h after start of a 30-min infusion for the first time-point and 8-10 h after start of the infusion for the second time-point. Duration of the infusion and accurate registration of the time-points of blood drawing are essential for obtaining precise predictions.

Fast parallel kriging-based stepwise uncertainty reduction with application to the identification of an excursion set

Stepwise uncertainty reduction (SUR) strategies aim at constructing a sequence of points for evaluating a function  f in such a way that the residual uncertainty about a quantity of interest progressively decreases to zero. Using such strategies in the framework of Gaussian process modeling has been shown to be efficient for estimating the volume of excursion of f above a fixed threshold. However, SUR strategies remain cumbersome to use in practice because of their high computational complexity, and the fact that they deliver a single point at each iteration. In this article we introduce several multipoint sampling criteria, allowing the selection of batches of points at which f can be evaluated in parallel. Such criteria are of particular interest when f is costly to evaluate and several CPUs are simultaneously available. We also manage to drastically reduce the computational cost of these strategies through the use of closed form formulas. We illustrate their performances in various numerical experiments, including a nuclear safety test case. Basic notions about kriging, auxiliary problems, complexity calculations, R code, and data are available online as supplementary materials.

Explicit mathematics and operational set theory: some ontological comparisons

We discuss several ontological properties of explicit mathematics and operational set theory: global choice, decidable classes, totality and extensionality of operations, function spaces, class and set formation via formulas that contain the definedness predicate and applications.

The Utilization of Classical Spin Monte Carlo Methods to Simulate the Magnetic Behavior of Extended Three-Dimensional Cubic Networks Incorporating M(II) Ions with an S = 5/2 Ground State Spin

The numerical simulations of the magnetic properties of extended three-dimensional networks containing M(II) ions with an S = 5/2 ground-state spin have been carried out within the framework of the isotropic Heisenberg model. Analytical expressions fitting the numerical simulations for the primitive cubic, diamond, together with (10−3) cubic networks have all been derived. With these empirical formulas in hands, we can now extract the interaction between the magnetic ions from the experimental data for these networks. In the case of the primitive cubic network, these expressions are directly compared with those from the high-temperature expansions of the partition function. A fit of the experimental data for three complexes, namely [(N(CH3)4][Mn(N3)] 1, [Mn(CN4)]n 2, and [FeII(bipy)3][MnII2(ox)3] 3, has been carried out. The best fits were those obtained using the following parameters, J = −3.5 cm-1, g = 2.01 (1); J = −8.3 cm-1, g = 1.95 (2); and J = −2.0 cm-1, g = 1.95 (3).

Factors influencing the accuracy of fetal weight estimation with a focus on preterm birth at the limit of viability: a systematic literature review

BACKGROUND Fetal weight estimation (FWE) is an important factor for clinical management decisions, especially in imminent preterm birth at the limit of viability between 23(0/7) and 26(0/7) weeks of gestation. It is crucial to detect and eliminate factors that have a negative impact on the accuracy of FWE. DATA SOURCES In this systematic literature review, we investigated 14 factors that may influence the accuracy of FWE, in particular in preterm neonates born at the limit of viability. RESULTS We found that gestational age, maternal body mass index, amniotic fluid index and ruptured membranes, presentation of the fetus, location of the placenta and the presence of multiple fetuses do not seem to have an impact on FWE accuracy. The influence of the examiner's grade of experience and that of fetal gender were discussed controversially. Fetal weight, time interval between estimation and delivery and the use of different formulas seem to have an evident effect on FWE accuracy. No results were obtained on the impact of active labor. DISCUSSION This review reveals that only few studies investigated factors possibly influencing the accuracy of FWE in preterm neonates at the limit of viability. Further research in this specific age group on potential confounding factors is needed.

Body Height Estimation from Femur Measurements in Postmortem Computed Tomography

After attending this presentation, attendees will: (1) understand how body height from computed tomography data can be estimated; and, (2) gain knowledge about the accuracy of estimated body height and limitations. The presentation will impact the forensic science community by providing knowledge and competence which will enable attendees to develop formulas for single bones to reconstruct body height using postmortem Computer Tomography (p-CT) data. The estimation of Body Height (BH) is an important component of the identification of corpses and skeletal remains. Stature can be estimated with relative accuracy via the measurement of long bones, such as the femora. Compared to time-consuming maceration procedures, p-CT allows fast and simple measurements of bones. This study undertook four objectives concerning the accuracy of BH estimation via p-CT: (1) accuracy between measurements on native bone and p-CT imaged bone (F1 according to Martin 1914); (2) intra-observer p-CT measurement precision; (3) accuracy between formula-based estimation of the BH and conventional body length measurement during autopsy; and, (4) accuracy of different estimation formulas available.1 In the first step, the accuracy of measurements in the CT compared to those obtained using an osteometric board was evaluated on the basis of eight defleshed femora. Then the femora of 83 female and 144 male corpses of a Swiss population for which p-CTs had been performed, were measured at the Institute of Forensic Medicine in Bern. After two months, 20 individuals were measured again in order to assess the intraobserver error. The mean age of the men was 53±17 years and that of the women was 61±20 years. Additionally, the body length of the corpses was measured conventionally. The mean body length was 176.6±7.2cm for men and 163.6±7.8cm for women. The images that were obtained using a six-slice CT were reconstructed with a slice thickness of 1.25mm. Analysis and measurements of CT images were performed on a multipurpose workstation. As a forensic standard procedure, stature was estimated by means of the regression equations by Penning & Riepert developed on a Southern German population and for comparison, also those referenced by Trotter & Gleser “American White.”2,3 All statistical tests were performed with a statistical software. No significant differences were found between the CT and osteometric board measurements. The double p-CT measurement of 20 individuals resulted in an absolute intra-observer difference of 0.4±0.3mm. For both sexes, the correlation between the body length and the estimated BH using the F1 measurements was highly significant. The correlation coefficient was slightly higher for women. The differences in accuracy of the different formulas were small. While the errors of BH estimation were generally ±4.5–5.0cm, the consideration of age led to an increase in accuracy of a few millimetres to about 1cm. BH estimations according to Penning & Riepert and Trotter & Gleser were slightly more accurate when age-at-death was taken into account.2,3 That way, stature estimations in the group of individuals older than 60 years were improved by about 2.4cm and 3.1cm.2,3 The error of estimation is therefore about a third of the common ±4.7cm error range. Femur measurements in p-CT allow very accurate BH estimations. Estimations according to Penning led to good results that (barely) come closer to the true value than the frequently used formulas by Trotter & Gleser “American White.”2,3 Therefore, the formulas by Penning & Riepert are also validated for this substantial recent Swiss population.

Skolemization for Substructural Logics

The usual Skolemization procedure, which removes strong quantifiers by introducing new function symbols, is in general unsound for first-order substructural logics defined based on classes of complete residuated lattices. However, it is shown here (following similar ideas of Baaz and Iemhoff for first-order intermediate logics in [1]) that first-order substructural logics with a semantics satisfying certain witnessing conditions admit a “parallel” Skolemization procedure where a strong quantifier is removed by introducing a finite disjunction or conjunction (as appropriate) of formulas with multiple new function symbols. These logics typically lack equivalent prenex forms. Also, semantic consequence does not in general reduce to satisfiability. The Skolemization theorems presented here therefore take various forms, applying to the left or right of the consequence relation, and to all formulas or only prenex formulas.