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.

Extending private-collective innovation: a case study

The private-collective innovation model proposes incentives for individuals and firms to privately invest resources to create public goods innovations. Such innovations are characterized by non-rivalry and non-exclusivity in consumption. Examples include open source software, user-generated media products, drug formulas, and sport equipment designs. There is still limited empirical research on private-collective innovation. We present a case study to (1) provide empirical evidence of a case of private-collective innovation, showing specific benefits, and (2) to extend the private-collective innovation model by analyzing the hidden costs for the company involved. We examine the development of the Nokia Internet Tablet, which builds on both proprietary and open source software development, and that involves both Nokia developers and volunteers who are not employed by the company. Seven benefits for Nokia are identified, as are five hidden costs: difficulty to differentiate, guarding business secrets, reducing community entry barriers, giving up control, and organizational inertia. We examine the actions taken by the management to mitigate these costs throughout the development period.

A Proof-Theoretic Analysis of Theories for Stratified Inductive Definitions

In this article we study subsystems SIDᵥ of the theory ID₁ in which fixed point induction is restricted to properly stratified formulas.

VIS-NIR reflectance of water ice/regolith analogue mixtures and implications for the detectability of ice mixed within planetary regoliths

Permanently shadowed regions at the poles of the Moon and Mercury have been pointed out as candidates for hosting water ice at their surface. We have measured in the laboratory the visible and near infrared spectral range (VIS-NIR) bidirectional reflectance of intimate mixtures of water ice and the JSC-1AF lunar simulant for different ice concentrations, particle sizes, and measurement geometries. The nonlinearity between the measured reflectance and the amount of ice in the mixture can be reproduced to some extent by the mixing formulas of standard reflectance models, in particular, those of Hapke and Hiroi, which are tested here. Estimating ice concentrations from reflectance data without knowledge of the mixing coefficientsstrongly dependent on the size/shape of the grainscan result in large errors. According to our results, it is possible that considerable amounts of water ice might be intimately mixed in the regolith of the Moon and Mercury without producing noticeable photometric signatures.

Nitric Oxide Equilibrium

The explicit formulas for following the passage to chemical equilibrium are derived.