Polydispersity index from linear viscoelastic data: Unimodal and bimodal linear polymer melts

This article describes a method for determining the polydispersity index Ip2=Mz/Mw of the molecular weight distribution (MWD) of linear polymeric materials from linear viscoelastic data. The method uses the Mellin transform of the relaxation modulus of a simple molecular rheological model. One of the main features of this technique is that it enables interesting MWD information to be obtained directly from dynamic shear experiments. It is not necessary to achieve the relaxation spectrum, so the ill-posed problem is avoided. Furthermore, a determinate shape of the continuous MWD does not have to be assumed in order to obtain the polydispersity index. The technique has been developed to deal with entangled linear polymers, whatever the form of the MWD is. The rheological information required to obtain the polydispersity index is the storage G′(ω) and loss G″(ω) moduli, extending from the terminal zone to the plateau region. The method provides a good agreement between the proposed theoretical approach and the experimental polydispersity indices of several linear polymers for a wide range of average molecular weights and polydispersity indices. It is also applicable to binary blends.

Determinants of methicillin-susceptible Staphylococcus aureus native bone and joint infection treatment failure: a retrospective cohort study.

BACKGROUND: Although methicillin-susceptible Staphylococcus aureus (MSSA) native bone and joint infection (BJI) constitutes the more frequent clinical entity of BJI, prognostic studies mostly focused on methicillin-resistant S. aureus prosthetic joint infection. We aimed to assess the determinants of native MSSA BJI outcomes. METHODS: Retrospective cohort study (2001-2011) of patients admitted in a reference hospital centre for native MSSA BJI. Treatment failure determinants were assessed using Kaplan-Meier curves and binary logistic regression. RESULTS: Sixty-six patients (42 males [63.6%]; median age 61.2 years; interquartile range [IQR] 45.9-71.9) presented an acute (n = 38; 57.6%) or chronic (n = 28; 42.4%) native MSSA arthritis (n = 15; 22.7%), osteomyelitis (n = 19; 28.8%) or spondylodiscitis (n = 32; 48.5%), considered as "difficult-to-treat" in 61 cases (92.4%). All received a prolonged (27.1 weeks; IQR, 16.9-36.1) combined antimicrobial therapy, after surgical management in 37 cases (56.1%). Sixteen treatment failures (24.2%) were observed during a median follow-up period of 63.3 weeks (IQR, 44.7-103.1), including 13 persisting infections, 1 relapse after treatment disruption, and 2 super-infections. Independent determinants of treatment failure were the existence of a sinus tract (odds ratio [OR], 5.300; 95% confidence interval [CI], 1.166-24.103) and a prolonged delay to infectious disease specialist referral (OR, 1.134; 95% CI 1.013-1.271). CONCLUSIONS: The important treatment failure rate pinpointed the difficulty of cure encountered in complicated native MSSA BJI. An early infectious disease specialist referral is essential, especially in debilitated patients or in presence of sinus tract.

Validation of real-time methylation-specific PCR to determine O6-methylguanine-DNA methyltransferase gene promoter methylation in glioma.

Epigenetic silencing of the DNA repair protein O(6)-methylguanine-DNA methyltransferase (MGMT) by promoter methylation predicts successful alkylating agent therapy, such as with temozolomide, in glioblastoma patients. Stratified therapy assignment of patients in prospective clinical trials according to tumor MGMT status requires a standardized diagnostic test, suitable for high-throughput analysis of small amounts of formalin-fixed, paraffin-embedded tumor tissue. A direct, real-time methylation-specific PCR (MSP) assay was developed to determine methylation status of the MGMT gene promoter. Assay specificity was obtained by selective amplification of methylated DNA sequences of sodium bisulfite-modified DNA. The copy number of the methylated MGMT promoter, normalized to the beta-actin gene, provides a quantitative test result. We analyzed 134 clinical glioma samples, comparing the new test with the previously validated nested gel-based MSP assay, which yields a binary readout. A cut-off value for the MGMT methylation status was suggested by fitting a bimodal normal mixture model to the real-time results, supporting the hypothesis that there are two distinct populations within the test samples. Comparison of the tests showed high concordance of the results (82/91 [90%]; Cohen's kappa = 0.80; 95% confidence interval, 0.82-0.95). The direct, real-time MSP assay was highly reproducible (Pearson correlation 0.996) and showed valid test results for 93% (125/134) of samples compared with 75% (94/125) for the nested, gel-based MSP assay. This high-throughput test provides an important pharmacogenomic tool for individualized management of alkylating agent chemotherapy.

Polymorphism of the grasshopper Rhammatocerus schistocercoides populations revealed by RAPD

The objective of this work was to study the genetic variability of the grasshopper Rhammatocerus schistocercoides (Orthoptera: Acrididae) using RAPD analysis among individuals from three populations, one from Colombia and two from Brazil (Goiás and Mato Grosso States). Ninety scorable binary markers were obtained by fingerprinting with 11 oligonucleotide primers. Most of the polymorphism was attributed to 42 markers with variable frequency among the different populations. Although the existence of significant difference among populations (P<0.0001), most of the genetic variability was found within populations (87.7% of total variation). Pairwise distances between Colombian and Brazilian populations were 0.12 (P<0.0001) and 0.18 (P<0.0001) for Goiás and Mato Grosso, respectively. The pairwise distance between Goiás and Mato Grosso populations was 0.06 (P<0.0001). These data indicated that the phenotypic differences among populations are associated mainly with the geographical distances between the Brazilian and Colombian populations.

SCRIB and PUF60 Are Primary Drivers of the Multisystemic Phenotypes of the 8q24.3 Copy-Number Variant.

Copy-number variants (CNVs) represent a significant interpretative challenge, given that each CNV typically affects the dosage of multiple genes. Here we report on five individuals with coloboma, microcephaly, developmental delay, short stature, and craniofacial, cardiac, and renal defects who harbor overlapping microdeletions on 8q24.3. Fine mapping localized a commonly deleted 78 kb region that contains three genes: SCRIB, NRBP2, and PUF60. In vivo dissection of the CNV showed discrete contributions of the planar cell polarity effector SCRIB and the splicing factor PUF60 to the syndromic phenotype, and the combinatorial suppression of both genes exacerbated some, but not all, phenotypic components. Consistent with these findings, we identified an individual with microcephaly, short stature, intellectual disability, and heart defects with a de novo c.505C>T variant leading to a p.His169Tyr change in PUF60. Functional testing of this allele in vivo and in vitro showed that the mutation perturbs the relative dosage of two PUF60 isoforms and, subsequently, the splicing efficiency of downstream PUF60 targets. These data inform the functions of two genes not associated previously with human genetic disease and demonstrate how CNVs can exhibit complex genetic architecture, with the phenotype being the amalgam of both discrete dosage dysfunction of single transcripts and also of binary genetic interactions.

Decision curve analysis revisited: overall net benefit, relationships to ROC curve analysis, and application to case-control studies.

ABSTRACT: BACKGROUND: Decision curve analysis has been introduced as a method to evaluate prediction models in terms of their clinical consequences if used for a binary classification of subjects into a group who should and into a group who should not be treated. The key concept for this type of evaluation is the "net benefit", a concept borrowed from utility theory. METHODS: We recall the foundations of decision curve analysis and discuss some new aspects. First, we stress the formal distinction between the net benefit for the treated and for the untreated and define the concept of the "overall net benefit". Next, we revisit the important distinction between the concept of accuracy, as typically assessed using the Youden index and a receiver operating characteristic (ROC) analysis, and the concept of utility of a prediction model, as assessed using decision curve analysis. Finally, we provide an explicit implementation of decision curve analysis to be applied in the context of case-control studies. RESULTS: We show that the overall net benefit, which combines the net benefit for the treated and the untreated, is a natural alternative to the benefit achieved by a model, being invariant with respect to the coding of the outcome, and conveying a more comprehensive picture of the situation. Further, within the framework of decision curve analysis, we illustrate the important difference between the accuracy and the utility of a model, demonstrating how poor an accurate model may be in terms of its net benefit. Eventually, we expose that the application of decision curve analysis to case-control studies, where an accurate estimate of the true prevalence of a disease cannot be obtained from the data, is achieved with a few modifications to the original calculation procedure. CONCLUSIONS: We present several interrelated extensions to decision curve analysis that will both facilitate its interpretation and broaden its potential area of application.

Computational problems in perfect phylogeny haplotyping: Typing without calling the allele

A haplotype is an m-long binary vector. The XOR-genotype of two haplotypes is the m-vector of their coordinate-wise XOR. We study the following problem: Given a set of XOR-genotypes, reconstruct their haplotypes so that the set of resulting haplotypes can be mapped onto a perfect phylogeny (PP) tree. The question is motivated by studying population evolution in human genetics and is a variant of the PP haplotyping problem that has received intensive attention recently. Unlike the latter problem, in which the input is '' full '' genotypes, here, we assume less informative input and so may be more economical to obtain experimentally. Building on ideas of Gusfield, we show how to solve the problem in polynomial time by a reduction to the graph realization problem. The actual haplotypes are not uniquely determined by the tree they map onto and the tree itself may or may not be unique. We show that tree uniqueness implies uniquely determined haplotypes, up to inherent degrees of freedom, and give a sufficient condition for the uniqueness. To actually determine the haplotypes given the tree, additional information is necessary. We show that two or three full genotypes suffice to reconstruct all the haplotypes and present a linear algorithm for identifying those genotypes.

FAM-MDR: a flexible family-based multifactor dimensionality reduction technique to detect epistasis using related individuals

We propose a novel multifactor dimensionality reduction method for epistasis detection in small or extended pedigrees, FAM-MDR. It combines features of the Genome-wide Rapid Association using Mixed Model And Regression approach (GRAMMAR) with Model-Based MDR (MB-MDR). We focus on continuous traits, although the method is general and can be used for outcomes of any type, including binary and censored traits. When comparing FAM-MDR with Pedigree-based Generalized MDR (PGMDR), which is a generalization of Multifactor Dimensionality Reduction (MDR) to continuous traits and related individuals, FAM-MDR was found to outperform PGMDR in terms of power, in most of the considered simulated scenarios. Additional simulations revealed that PGMDR does not appropriately deal with multiple testing and consequently gives rise to overly optimistic results. FAM-MDR adequately deals with multiple testing in epistasis screens and is in contrast rather conservative, by construction. Furthermore, simulations show that correcting for lower order (main) effects is of utmost importance when claiming epistasis. As Type 2 Diabetes Mellitus (T2DM) is a complex phenotype likely influenced by gene-gene interactions, we applied FAM-MDR to examine data on glucose area-under-the-curve (GAUC), an endophenotype of T2DM for which multiple independent genetic associations have been observed, in the Amish Family Diabetes Study (AFDS). This application reveals that FAM-MDR makes more efficient use of the available data than PGMDR and can deal with multi-generational pedigrees more easily. In conclusion, we have validated FAM-MDR and compared it to PGMDR, the current state-of-the-art MDR method for family data, using both simulations and a practical dataset. FAM-MDR is found to outperform PGMDR in that it handles the multiple testing issue more correctly, has increased power, and efficiently uses all available information.

An efficient algorithm to perform multiple testing in epistasis screening

Background: Research in epistasis or gene-gene interaction detection for human complex traits has grown over the last few years. It has been marked by promising methodological developments, improved translation efforts of statistical epistasis to biological epistasis and attempts to integrate different omics information sources into the epistasis screening to enhance power. The quest for gene-gene interactions poses severe multiple-testing problems. In this context, the maxT algorithm is one technique to control the false-positive rate. However, the memory needed by this algorithm rises linearly with the amount of hypothesis tests. Gene-gene interaction studies will require a memory proportional to the squared number of SNPs. A genome-wide epistasis search would therefore require terabytes of memory. Hence, cache problems are likely to occur, increasing the computation time. In this work we present a new version of maxT, requiring an amount of memory independent from the number of genetic effects to be investigated. This algorithm was implemented in C++ in our epistasis screening software MBMDR-3.0.3. We evaluate the new implementation in terms of memory efficiency and speed using simulated data. The software is illustrated on real-life data for Crohn’s disease. Results: In the case of a binary (affected/unaffected) trait, the parallel workflow of MBMDR-3.0.3 analyzes all gene-gene interactions with a dataset of 100,000 SNPs typed on 1000 individuals within 4 days and 9 hours, using 999 permutations of the trait to assess statistical significance, on a cluster composed of 10 blades, containing each four Quad-Core AMD Opteron(tm) Processor 2352 2.1 GHz. In the case of a continuous trait, a similar run takes 9 days. Our program found 14 SNP-SNP interactions with a multiple-testing corrected p-value of less than 0.05 on real-life Crohn’s disease (CD) data. Conclusions: Our software is the first implementation of the MB-MDR methodology able to solve large-scale SNP-SNP interactions problems within a few days, without using much memory, while adequately controlling the type I error rates. A new implementation to reach genome-wide epistasis screening is under construction. In the context of Crohn’s disease, MBMDR-3.0.3 could identify epistasis involving regions that are well known in the field and could be explained from a biological point of view. This demonstrates the power of our software to find relevant phenotype-genotype higher-order associations.

At-risk drinking and drug use among patients seeking care in an emergency department.

OBJECTIVE: This study reports the frequency of alcohol use and associated tobacco and drug use among emergency department (ED) patients, in order to increase physician awareness and treatment of women and men seeking care in ED settings. METHOD: All adults seen in the ED at the University Hospital in Lausanne, Switzerland, between 11 AM and 11 PM were screened by direct interview for at-risk drinking, tobacco use, drug use, and depression during an 18-month period. RESULTS: A total of 8,599 patients (4,006 women and 4,593 men) participated in the screening procedure and provided full data on the variables in our analysis. The mean age was 51.9 years for women and 45.0 years for men; 57.5% (n = 2,304) of women and 58.5% (n = 2,688) of men were being treated for trauma. Based on guidelines of the National Institute on Alcohol Abuse and Alcoholism, 13.1% (n = 523) of the women were at-risk drinkers, 57.3% (n = 2,301) were low-risk drinkers, and 29.6% (n = 1,182) were abstinent. Among men, 32.8% (n = 1,507) met criteria for at-risk drinking, 51.8% (n = 2,380) met criteria for low-risk drinking, and 15.4% (n = 706) were abstinent. Younger individuals (ages 18-30) had significantly higher rates of episodic heavy drinking episodes, whereas at-risk older patients were more likely to drink on a daily basis. A binary model found that women and men who drank at at-risk levels are more likely to use tobacco (odds ratio [OR] = 2.48, 95% confidence interval [CI]: 2.0-3.08) and illicit drugs (OR = 5.91, CI: 3.32- 10.54) compared with abstinent and low-risk drinkers. CONCLUSIONS: This study supports systematic alcohol screening of women and men seen in EDs and suggests that patterns of alcohol and drug use vary by age and gender.

GUS gene expression driven by a citrus promoter in transgenic tobacco and 'Valencia' sweet orange

The objective of this work was the transformation of tobacco and 'Valencia' sweet orange with the GUS gene driven by the citrus phenylalanine ammonia-lyase (PAL) gene promoter (CsPP). Transformation was accomplished by co-cultivation of tobacco and 'Valência' sweet orange explants with Agrobacterium tumefaciens containing the binary vector CsPP-GUS/2201. After plant transformation and regeneration, histochemical analyses using GUS staining revealed that CsPP promoter preferentially, but not exclusively, conferred gene expression in xylem tissues of tobacco. Weaker GUS staining was also detected throughout the petiole region in tobacco and citrus CsPP transgenic plants.

Multislice CT for assessing in-stent dimensions after left main coronary artery stenting: a comparison with three dimensional intravascular ultrasound

Objective: To evaluate the agreement between multislice CT (MSCT) and intravascular ultrasound (IVUS) to assess the in-stent lumen diameters and lumen areas of left main coronary artery (LMCA) stents. Design: Prospective, observational single centre study. Setting: A single tertiary referral centre. Patients: Consecutive patients with LMCA stenting excluding patients with atrial fibrillation and chronic renal failure. Interventions: MSCT and IVUS imaging at 912 months follow-up were performed for all patients. Main outcome measures: Agreement between MSCT and IVUS minimum luminal area (MLA) and minimum luminal diameter (MLD). A receiver operating characteristic (ROC) curve was plotted to find the MSCT cut-off point to diagnose binary restenosis equivalent to 6 mm2 by IVUS. Results: 52 patients were analysed. PassingBablok regression analysis obtained a β coefficient of 0.786 (0.586 to 1.071) for MLA and 1.250 (0.936 to 1.667) for MLD, ruling out proportional bias. The α coefficient was −3.588 (−8.686 to −0.178) for MLA and −1.713 (−3.583 to −0.257) for MLD, indicating an underestimation trend of MSCT. The ROC curve identified an MLA ≤4.7 mm2 as the best threshold to assess in-stent restenosis by MSCT. Conclusions: Agreement between MSCT and IVUS to assess in-stent MLA and MLD for LMCA stenting is good. An MLA of 4.7 mm2 by MSCT is the best threshold to assess binary restenosis. MSCT imaging can be considered in selected patients to assess LMCA in-stent restenosis

Statistical mechanics of multi-edge networks

Statistical properties of binary complex networks are well understood and recently many attempts have been made to extend this knowledge to weighted ones. There are, however, subtle yet important considerations to be made regarding the nature of the weights used in this generalization. Weights can be either continuous or discrete magnitudes, and in the latter case, they can additionally have undistinguishable or distinguishable nature. This fact has not been addressed in the literature insofar and has deep implications on the network statistics. In this work we face this problem introducing multiedge networks as graphs where multiple (distinguishable) connections between nodes are considered. We develop a statistical mechanics framework where it is possible to get information about the most relevant observables given a large spectrum of linear and nonlinear constraints including those depending both on the number of multiedges per link and their binary projection. The latter case is particularly interesting as we show that binary projections can be understood from multiedge processes. The implications of these results are important as many real-agent-based problems mapped onto graphs require this treatment for a proper characterization of their collective behavior.

Do stacked species distribution models reflect altitudinal diversity patterns?

The objective of this study was to evaluate the performance of stacked species distribution models in predicting the alpha and gamma species diversity patterns of two important plant clades along elevation in the Andes. We modelled the distribution of the species in the Anthurium genus (53 species) and the Bromeliaceae family (89 species) using six modelling techniques. We combined all of the predictions for the same species in ensemble models based on two different criteria: the average of the rescaled predictions by all techniques and the average of the best techniques. The rescaled predictions were then reclassified into binary predictions (presence/absence). By stacking either the original predictions or binary predictions for both ensemble procedures, we obtained four different species richness models per taxa. The gamma and alpha diversity per elevation band (500 m) was also computed. To evaluate the prediction abilities for the four predictions of species richness and gamma diversity, the models were compared with the real data along an elevation gradient that was independently compiled by specialists. Finally, we also tested whether our richness models performed better than a null model of altitudinal changes of diversity based on the literature. Stacking of the ensemble prediction of the individual species models generated richness models that proved to be well correlated with the observed alpha diversity richness patterns along elevation and with the gamma diversity derived from the literature. Overall, these models tend to overpredict species richness. The use of the ensemble predictions from the species models built with different techniques seems very promising for modelling of species assemblages. Stacking of the binary models reduced the over-prediction, although more research is needed. The randomisation test proved to be a promising method for testing the performance of the stacked models, but other implementations may still be developed.

Cluster glasses of ultrasoft particles

We present molecular dynamics (MD) simulations results for dense fluids of ultrasoft, fully penetrable particles. These are a binary mixture and a polydisperse system of particles interacting via the generalized exponential model, which is known to yield cluster crystal phases for the corresponding monodisperse systems. Because of the dispersity in the particle size, the systems investigated in this work do not crystallize and form disordered cluster phases. The clusteringtransition appears as a smooth crossover to a regime in which particles are mostly located in clusters, isolated particles being infrequent. The analysis of the internal cluster structure reveals microsegregation of the big and small particles, with a strong homo-coordination in the binary mixture. Upon further lowering the temperature below the clusteringtransition, the motion of the clusters" centers-of-mass slows down dramatically, giving way to a cluster glass transition. In the cluster glass, the diffusivities remain finite and display an activated temperature dependence, indicating that relaxation in the cluster glass occurs via particle hopping in a nearly arrested matrix of clusters. Finally we discuss the influence of the microscopic dynamics on the transport properties by comparing the MD results with Monte Carlo simulations.