Rgtsp: a generalized top scoring pairs package for class prediction.

SUMMARY: A top scoring pair (TSP) classifier consists of a pair of variables whose relative ordering can be used for accurately predicting the class label of a sample. This classification rule has the advantage of being easily interpretable and more robust against technical variations in data, as those due to different microarray platforms. Here we describe a parallel implementation of this classifier which significantly reduces the training time, and a number of extensions, including a multi-class approach, which has the potential of improving the classification performance. AVAILABILITY AND IMPLEMENTATION: Full C++ source code and R package Rgtsp are freely available from http://lausanne.isb-sib.ch/~vpopovic/research/. The implementation relies on existing OpenMP libraries.

Comparative performance of supertree algorithms in large data sets using the soapberry family (Sapindaceae) as a case study.

For the last 2 decades, supertree reconstruction has been an active field of research and has seen the development of a large number of major algorithms. Because of the growing popularity of the supertree methods, it has become necessary to evaluate the performance of these algorithms to determine which are the best options (especially with regard to the supermatrix approach that is widely used). In this study, seven of the most commonly used supertree methods are investigated by using a large empirical data set (in terms of number of taxa and molecular markers) from the worldwide flowering plant family Sapindaceae. Supertree methods were evaluated using several criteria: similarity of the supertrees with the input trees, similarity between the supertrees and the total evidence tree, level of resolution of the supertree and computational time required by the algorithm. Additional analyses were also conducted on a reduced data set to test if the performance levels were affected by the heuristic searches rather than the algorithms themselves. Based on our results, two main groups of supertree methods were identified: on one hand, the matrix representation with parsimony (MRP), MinFlip, and MinCut methods performed well according to our criteria, whereas the average consensus, split fit, and most similar supertree methods showed a poorer performance or at least did not behave the same way as the total evidence tree. Results for the super distance matrix, that is, the most recent approach tested here, were promising with at least one derived method performing as well as MRP, MinFlip, and MinCut. The output of each method was only slightly improved when applied to the reduced data set, suggesting a correct behavior of the heuristic searches and a relatively low sensitivity of the algorithms to data set sizes and missing data. Results also showed that the MRP analyses could reach a high level of quality even when using a simple heuristic search strategy, with the exception of MRP with Purvis coding scheme and reversible parsimony. The future of supertrees lies in the implementation of a standardized heuristic search for all methods and the increase in computing power to handle large data sets. The latter would prove to be particularly useful for promising approaches such as the maximum quartet fit method that yet requires substantial computing power.

Midregional pro-atrial natriuretic peptide and procalcitonin improve survival prediction in VAP.

Ventilator-associated pneumonia (VAP) affects mortality, morbidity and cost of critical care. Reliable risk estimation might improve end-of-life decisions, resource allocation and outcome. Several scoring systems for survival prediction have been established and optimised over the last decades. Recently, new biomarkers have gained interest in the prognostic field. We assessed whether midregional pro-atrial natriuretic peptide (MR-proANP) and procalcitonin (PCT) improve the predictive value of the Simplified Acute Physiologic Score (SAPS) II and Sequential Related Organ Failure Assessment (SOFA) in VAP. Specified end-points of a prospective multinational trial including 101 patients with VAP were analysed. Death <28 days after VAP onset was the primary end-point. MR-proANP and PCT were elevated at the onset of VAP in nonsurvivors compared with survivors (p = 0.003 and p = 0.017, respectively) and their slope of decline differed significantly (p = 0.018 and p = 0.039, respectively). Patients with the highest MR-proANP quartile at VAP onset were at increased risk for death (log rank p = 0.013). In a logistic regression model, MR-proANP was identified as the best predictor of survival. Adding MR-proANP and PCT to SAPS II and SOFA improved their predictive properties (area under the curve 0.895 and 0.880). We conclude that the combination of two biomarkers, MR-proANP and PCT, improve survival prediction of clinical severity scores in VAP.

Prediction of resting energy expenditure from fat-free mass and fat mass.

On the basis of literature values, the relationship between fat-free mass (FFM), fat mass (FM), and resting energy expenditure [REE (kJ/24 h)] was determined for 213 adults (86 males, 127 females). The objectives were to develop a mathematical model to predict REE based on body composition and to evaluate the contribution of FFM and FM to REE. The following regression equations were derived: 1) REE = 1265 + (93.3 x FFM) (r2 = 0.727, P < 0.001); 2) REE = 1114 + (90.4 x FFM) + (13.2 x FM) (R2 = 0.743, P < 0.001); and 3) REE = (108 x FFM) + (16.9 x FM) (R2 = 0.986, P < 0.001). FM explained only a small part of the variation remaining after FFM was accounted for. The models that include both FFM and FM are useful in examination of the changes in REE that occur with a change in both the FFM and FM. To account for more of the variability in REE, FFM will have to be divided into organ mass and skeletal muscle mass in future analyses.

Ruling out coronary heart disease in primary care: external validation of a clinical prediction rule.

BACKGROUND: The Marburg Heart Score (MHS) aims to assist GPs in safely ruling out coronary heart disease (CHD) in patients presenting with chest pain, and to guide management decisions. AIM: To investigate the diagnostic accuracy of the MHS in an independent sample and to evaluate the generalisability to new patients. DESIGN AND SETTING: Cross-sectional diagnostic study with delayed-type reference standard in general practice in Hesse, Germany. METHOD: Fifty-six German GPs recruited 844 males and females aged ≥ 35 years, presenting between July 2009 and February 2010 with chest pain. Baseline data included the items of the MHS. Data on the subsequent course of chest pain, investigations, hospitalisations, and medication were collected over 6 months and were reviewed by an independent expert panel. CHD was the reference condition. Measures of diagnostic accuracy included the area under the receiver operating characteristic curve (AUC), sensitivity, specificity, likelihood ratios, and predictive values. RESULTS: The AUC was 0.84 (95% confidence interval [CI] = 0.80 to 0.88). For a cut-off value of 3, the MHS showed a sensitivity of 89.1% (95% CI = 81.1% to 94.0%), a specificity of 63.5% (95% CI = 60.0% to 66.9%), a positive predictive value of 23.3% (95% CI = 19.2% to 28.0%), and a negative predictive value of 97.9% (95% CI = 96.2% to 98.9%). CONCLUSION: Considering the diagnostic accuracy of the MHS, its generalisability, and ease of application, its use in clinical practice is recommended.

MGMT methylation analysis of glioblastoma on the Infinium methylation BeadChip identifies two distinct CpG regions associated with gene silencing and outcome, yielding a prediction model for comparisons across datasets, tumor grades, and CIMP-status.

The methylation status of the O(6)-methylguanine-DNA methyltransferase (MGMT) gene is an important predictive biomarker for benefit from alkylating agent therapy in glioblastoma. Recent studies in anaplastic glioma suggest a prognostic value for MGMT methylation. Investigation of pathogenetic and epigenetic features of this intriguingly distinct behavior requires accurate MGMT classification to assess high throughput molecular databases. Promoter methylation-mediated gene silencing is strongly dependent on the location of the methylated CpGs, complicating classification. Using the HumanMethylation450 (HM-450K) BeadChip interrogating 176 CpGs annotated for the MGMT gene, with 14 located in the promoter, two distinct regions in the CpG island of the promoter were identified with high importance for gene silencing and outcome prediction. A logistic regression model (MGMT-STP27) comprising probes cg1243587 and cg12981137 provided good classification properties and prognostic value (kappa = 0.85; log-rank p < 0.001) using a training-set of 63 glioblastomas from homogenously treated patients, for whom MGMT methylation was previously shown to be predictive for outcome based on classification by methylation-specific PCR. MGMT-STP27 was successfully validated in an independent cohort of chemo-radiotherapy-treated glioblastoma patients (n = 50; kappa = 0.88; outcome, log-rank p < 0.001). Lower prevalence of MGMT methylation among CpG island methylator phenotype (CIMP) positive tumors was found in glioblastomas from The Cancer Genome Atlas than in low grade and anaplastic glioma cohorts, while in CIMP-negative gliomas MGMT was classified as methylated in approximately 50 % regardless of tumor grade. The proposed MGMT-STP27 prediction model allows mining of datasets derived on the HM-450K or HM-27K BeadChip to explore effects of distinct epigenetic context of MGMT methylation suspected to modulate treatment resistance in different tumor types.

Advanced geostatistical and machine-learning models for spatial data analysis of radioactively contaminated regions

Radioactive soil-contamination mapping and risk assessment is a vital issue for decision makers. Traditional approaches for mapping the spatial concentration of radionuclides employ various regression-based models, which usually provide a single-value prediction realization accompanied (in some cases) by estimation error. Such approaches do not provide the capability for rigorous uncertainty quantification or probabilistic mapping. Machine learning is a recent and fast-developing approach based on learning patterns and information from data. Artificial neural networks for prediction mapping have been especially powerful in combination with spatial statistics. A data-driven approach provides the opportunity to integrate additional relevant information about spatial phenomena into a prediction model for more accurate spatial estimates and associated uncertainty. Machine-learning algorithms can also be used for a wider spectrum of problems than before: classification, probability density estimation, and so forth. Stochastic simulations are used to model spatial variability and uncertainty. Unlike regression models, they provide multiple realizations of a particular spatial pattern that allow uncertainty and risk quantification. This paper reviews the most recent methods of spatial data analysis, prediction, and risk mapping, based on machine learning and stochastic simulations in comparison with more traditional regression models. The radioactive fallout from the Chernobyl Nuclear Power Plant accident is used to illustrate the application of the models for prediction and classification problems. This fallout is a unique case study that provides the challenging task of analyzing huge amounts of data ('hard' direct measurements, as well as supplementary information and expert estimates) and solving particular decision-oriented problems.

Identification of human IKK-2 inhibitors of natural origin (Part II): in Silico prediction of IKK-2 inhibitors in natural extracts with known anti-inflammatory activity.

Human inhibitor NF-κB kinase 2 (hIKK-2) is the primary component responsible for activating NF-κB in response to various inflammatory stimuli. Thus, synthetic ATP-competitive inhibitors for hIKK-2 have been developed as anti-inflammatory compounds. We recently reported a virtual screening protocol (doi:10.1371/journal.pone.0016903) that is able to identify hIKK-2 inhibitors that are not structurally related to any known molecule that inhibits hIKK-2 and that have never been reported to have anti-inflammatory activity. In this study, a stricter version of this protocol was applied to an in-house database of 29,779 natural products annotated with their natural source. The search identified 274 molecules (isolated from 453 different natural extracts) predicted to inhibit hIKK-2. An exhaustive bibliographic search revealed that anti-inflammatory activity has been previously described for: (a) 36 out of these 453 extracts; and (b) 17 out of 30 virtual screening hits present in these 36 extracts. Only one of the remaining 13 hit molecules in these extracts shows chemical similarity with known synthetic hIKK-2 inhibitors. Therefore, it is plausible that a significant portion of the remaining 12 hit molecules are lead-hopping candidates for the development of new hIKK-2 inhibitors.

Official Positions for FRAX® clinical regarding rheumatoid arthritis from Joint Official Positions Development Conference of the International Society for Clinical Densitometry and International Osteoporosis Foundation on FRAX®.

Rheumatoid arthritis is the only secondary cause of osteoporosis that is considered independent of bone density in the FRAX(®) algorithm. Although input for rheumatoid arthritis in FRAX(®) is a dichotomous variable, intuitively, one would expect that more severe or active disease would be associated with a greater risk for fracture. We reviewed the literature to determine if specific disease parameters or medication use could be used to better characterize fracture risk in individuals with rheumatoid arthritis. Although many studies document a correlation between various parameters of disease activity or severity and decreased bone density, fewer have associated these variables with fracture risk. We reviewed these studies in detail and concluded that disability measures such as HAQ (Health Assessment Questionnaire) and functional class do correlate with clinical fractures but not morphometric vertebral fractures. One large study found a strong correlation with duration of disease and fracture risk but additional studies are needed to confirm this. There was little evidence to correlate other measures of disease such as DAS (disease activity score), VAS (visual analogue scale), acute phase reactants, use of non-glucocorticoid medications and increased fracture risk. We concluded that FRAX(®) calculations may underestimate fracture probability in patients with impaired functional status from rheumatoid arthritis but that this could not be quantified at this time. At this time, other disease measures cannot be used for fracture prediction. However only a few, mostly small studies addressed other disease parameters and further research is needed. Additional questions for future research are suggested.

Treatment of ruptured abdominal aortic aneurysm, a permanent challenge or a waste of resources? Prediction of outcome using a multi-organ-dysfunction score.

OBJECTIVES: in a retrospective study, attempts have been made to identify individual organ-dysfunction risk profiles influencing the outcome after surgery for ruptured abdominal aortic aneurysms. METHODS: out of 235 patients undergoing graft replacement for abdominal aortic aneurysms, 57 (53 men, four women, mean age 72 years [s.d. 8.8]) were treated for ruptured aneurysms in a 3-year period. Forty-eight preoperative, 13 intraoperative and 34 postoperative variables were evaluated statistically. A simple multi-organ dysfunction (MOD) score was adopted. RESULTS: the perioperative mortality was 32%. Three patients died intraoperatively, four within 48 h and 11 died later. A significant influence for pre-existing risk factors was identified only for cardiovascular diseases. Multiple linear-regression analysis indicated that a haemoglobin <90 g/l, systolic blood pressure <80 mmHg and ECG signs of ischaemia at admission were highly significant risk factors. The cause of death for patients, who died more than 48 h postoperatively, was mainly MOD. All patients with a MOD score >/=4 died (n=7). These patients required 27% of the intensive-care unit (ICU) days of all patients and 72% of the ICU days of the non-survivors. CONCLUSION: patients with ruptured aortic aneurysms from treatment should not be excluded. However, a physiological scoring system after 48 h appears justifiable in order to decide on the appropriateness of continual ICU support.

Individual prediction of cognitive decline in mild cognitive impairment using support vector machine-based analysis of diffusion tensor imaging data.

Although cross-sectional diffusion tensor imaging (DTI) studies revealed significant white matter changes in mild cognitive impairment (MCI), the utility of this technique in predicting further cognitive decline is debated. Thirty-five healthy controls (HC) and 67 MCI subjects with DTI baseline data were neuropsychologically assessed at one year. Among them, there were 40 stable (sMCI; 9 single domain amnestic, 7 single domain frontal, 24 multiple domain) and 27 were progressive (pMCI; 7 single domain amnestic, 4 single domain frontal, 16 multiple domain). Fractional anisotropy (FA) and longitudinal, radial, and mean diffusivity were measured using Tract-Based Spatial Statistics. Statistics included group comparisons and individual classification of MCI cases using support vector machines (SVM). FA was significantly higher in HC compared to MCI in a distributed network including the ventral part of the corpus callosum, right temporal and frontal pathways. There were no significant group-level differences between sMCI versus pMCI or between MCI subtypes after correction for multiple comparisons. However, SVM analysis allowed for an individual classification with accuracies up to 91.4% (HC versus MCI) and 98.4% (sMCI versus pMCI). When considering the MCI subgroups separately, the minimum SVM classification accuracy for stable versus progressive cognitive decline was 97.5% in the multiple domain MCI group. SVM analysis of DTI data provided highly accurate individual classification of stable versus progressive MCI regardless of MCI subtype, indicating that this method may become an easily applicable tool for early individual detection of MCI subjects evolving to dementia.

Improvement of cardiovascular risk prediction: time to review current knowledge, debates, and fundamentals on how to assess test characteristics.

Cardiovascular risk assessment might be improved with the addition of emerging, new tests derived from atherosclerosis imaging, laboratory tests or functional tests. This article reviews relative risk, odds ratios, receiver-operating curves, posttest risk calculations based on likelihood ratios, the net reclassification improvement and integrated discrimination. This serves to determine whether a new test has an added clinical value on top of conventional risk testing and how this can be verified statistically. Two clinically meaningful examples serve to illustrate novel approaches. This work serves as a review and basic work for the development of new guidelines on cardiovascular risk prediction, taking into account emerging tests, to be proposed by members of the 'Taskforce on Vascular Risk Prediction' under the auspices of the Working Group 'Swiss Atherosclerosis' of the Swiss Society of Cardiology in the future.

Stochastic models and numerical algorithms for a class of regulatory gene networks.

Regulatory gene networks contain generic modules, like those involving feedback loops, which are essential for the regulation of many biological functions (Guido et al. in Nature 439:856-860, 2006). We consider a class of self-regulated genes which are the building blocks of many regulatory gene networks, and study the steady-state distribution of the associated Gillespie algorithm by providing efficient numerical algorithms. We also study a regulatory gene network of interest in gene therapy, using mean-field models with time delays. Convergence of the related time-nonhomogeneous Markov chain is established for a class of linear catalytic networks with feedback loops.