rpartOrdinal: An R Package for Deriving a Classification Tree for Predicting an Ordinal Response

This paper describes an R package, rpartOrdinal, that implements alternative splitting functions for fitting a classification tree when interest lies in predicting an ordinal response. This includes the generalized Gini impurity function, which was introduced as a method for predicting an ordinal response by including costs of misclassification into the impurity function, as well as an alternative ordinal impurity function due to Piccarreta (2008) that does not require the assignment of misclassification costs. The ordered twoing splitting method, which is not defined as a decrease in node impurity, is also included in the package. Since, in the ordinal response setting, misclassifying observations to adjacent categories is a less egregious error than misclassifying observations to distant categories, this package also includes a function for estimating an ordinal measure of association, the gamma statistic

Graphical technique for identifying a monotonic variance stabilizing transformation for absolute gene intensity signals

BACKGROUND: The usefulness of log(2 )transformation for cDNA microarray data has led to its widespread application to Affymetrix data. For Affymetrix data, where absolute intensities are indicative of number of transcripts, there is a systematic relationship between variance and magnitude of measurements. Application of the log(2 )transformation expands the scale of genes with low intensities while compressing the scale of genes with higher intensities thus reversing the mean by variance relationship. The usefulness of these transformations needs to be examined. RESULTS: Using an Affymetrix GeneChip(® )dataset, problems associated with applying the log(2 )transformation to absolute intensity data are demonstrated. Use of the spread-versus-level plot to identify an appropriate variance stabilizing transformation is presented. For the data presented, the spread-versus-level plot identified a power transformation that successfully stabilized the variance of probe set summaries. CONCLUSION: The spread-versus-level plot is helpful to identify transformations for variance stabilization. This is robust against outliers and avoids assumption of models and maximizations

Empirical validation of the S-Score algorithm in the analysis of gene expression data

BACKGROUND: Current methods of analyzing Affymetrix GeneChip(® )microarray data require the estimation of probe set expression summaries, followed by application of statistical tests to determine which genes are differentially expressed. The S-Score algorithm described by Zhang and colleagues is an alternative method that allows tests of hypotheses directly from probe level data. It is based on an error model in which the detected signal is proportional to the probe pair signal for highly expressed genes, but approaches a background level (rather than 0) for genes with low levels of expression. This model is used to calculate relative change in probe pair intensities that converts probe signals into multiple measurements with equalized errors, which are summed over a probe set to form the S-Score. Assuming no expression differences between chips, the S-Score follows a standard normal distribution, allowing direct tests of hypotheses to be made. Using spike-in and dilution datasets, we validated the S-Score method against comparisons of gene expression utilizing the more recently developed methods RMA, dChip, and MAS5. RESULTS: The S-score showed excellent sensitivity and specificity in detecting low-level gene expression changes. Rank ordering of S-Score values more accurately reflected known fold-change values compared to other algorithms. CONCLUSION: The S-score method, utilizing probe level data directly, offers significant advantages over comparisons using only probe set expression summaries

Penalized Ordinal Regression Methods for Predicting Stage of Cancer in High-Dimensional Covariate Spaces

The pathological description of the stage of a tumor is an important clinical designation and is considered, like many other forms of biomedical data, an ordinal outcome. Currently, statistical methods for predicting an ordinal outcome using clinical, demographic, and high-dimensional correlated features are lacking. In this paper, we propose a method that fits an ordinal response model to predict an ordinal outcome for high-dimensional covariate spaces. Our method penalizes some covariates (high-throughput genomic features) without penalizing others (such as demographic and/or clinical covariates). We demonstrate the application of our method to predict the stage of breast cancer. In our model, breast cancer subtype is a nonpenalized predictor, and CpG site methylation values from the Illumina Human Methylation 450K assay are penalized predictors. The method has been made available in the ordinalgmifs package in the R programming environment

Epigenetic Alterations and an Increased Frequency of Micronuclei in Women with Fibromyalgia

Fibromyalgia (FM), characterized by chronic widespread pain, fatigue, and cognitive/mood disturbances, leads to reduced workplace productivity and increased healthcare expenses. To determine if acquired epigenetic/genetic changes are associated with FM, we compared the frequency of spontaneously occurring micronuclei (MN) and genome-wide methylation patterns in women with FM () to those seen in comparably aged healthy controls ( (MN); (methylation)). The mean (sd) MN frequency of women with FM (51.4 (21.9)) was significantly higher than that of controls (15.8 (8.5)) (; df = 1; ). Significant differences ( sites) in methylation patterns were observed between cases and controls considering a 5% false discovery rate. The majority of differentially methylated (DM) sites (91%) were attributable to increased values in the women with FM. The DM sites included significant biological clusters involved in neuron differentiation/nervous system development, skeletal/organ system development, and chromatin compaction. Genes associated with DM sites whose function has particular relevance to FM included BDNF, NAT15, HDAC4, PRKCA, RTN1, and PRKG1. Results support the need for future research to further examine the potential role of epigenetic and acquired chromosomal alterations as a possible biological mechanism underlying FM

Reduced Expression of Inflammatory Genes in Deceased Donor Kidneys Undergoing Pulsatile Pump Preservation

Background The use of expanded criteria donor kidneys (ECD) had been associated with worse outcomes. Whole gene expression of pre-implantation allograft biopsies from deceased donor kidneys (DDKs) was evaluated to compare the effect of pulsatile pump preservation (PPP) vs. cold storage preservation (CSP) on standard and ECD kidneys. Methodology/Principal Findings 99 pre-implantation DDK biopsies were studied using gene expression with GeneChips. Kidneys transplant recipients were followed post transplantation for 35.8 months (range = 24–62). The PPP group included 60 biopsies (cold ischemia time (CIT) = 1,367+/−509 minutes) and the CSP group included 39 biopsies (CIT = 1,022+/−485 minutes) (P Conclusions/Significance Inflammation was the most important up-regulated pattern associated with pre-implantation biopsies undergoing CSP even when the PPP group has a larger number of ECD kidneys. No significant difference was observed in delayed graft function incidence and graft function post-transplantation. These findings support the use of PPP in ECD donor kidneys

A MEIG1/PACRG complex in the manchette is essential for building the sperm flagella

A key event in the process of spermiogenesis is the formation of the flagella, which enables sperm to reach eggs for fertilization. Yeast two-hybrid studies revealed that meiosis-expressed gene 1 (MEIG1) and Parkin co-regulated gene (PACRG) interact, and that sperm-associated antigen 16, which encodes an axoneme central apparatus protein, is also a binding partner of MEIG1. In spermatocytes of wild-type mice, MEIG1 is expressed in the whole germ cell bodies, but the protein migrates to the manchette, a unique structure at the base of elongating spermatid that directs formation of the flagella. In the elongating spermatids of wild-type mice, PACRG colocalizes with α-tubulin, a marker for the manchette, whereas this localization was not changed in the few remaining elongating spermatids of Meig1-deficient mice. In addition, MEIG1 no longer localizes to the manchette in the remaining elongating spermatids of Pacrg-deficient mice, indicating that PACRG recruits MEIG1 to the manchette. PACRG is not stable in mammalian cells, but can be stabilized by MEIG1 or by inhibition of proteasome function. SPAG16L is present in the spermatocyte cytoplasm of wild-type mice, and in the manchette of elongating spermatids, but in the Meig1 or Pacrg-deficient mice, SPAG16L no longer localizes to the manchette. By contrast, MEIG1 and PACRG are still present in the manchette of Spag16L-deficient mice, indicating that SPAG16L is a downstream partner of these two proteins. Together, our studies demonstrate that MEIG1/PACRG forms a complex in the manchette and that this complex is necessary to transport cargos, such as SPAG16L, to build the sperm flagella

Biomarkers of disease differentiation: HCV recurrence versus acute cellular rejection

The wound-healing process induced by chronic hepatitis C virus (HCV) infection triggers liver damage characterized by fibrosis development and finally cirrhosis. Liver Transplantation (LT) is the optimal surgical treatment for HCV-cirrhotic patients at end-stage liver disease. However, acute cellular rejection (ACR) and HCV recurrence disease represent two devastating complications post-LT. The accurate differential diagnosis between both conditions is critical for treatment choice, and similar histological features represent a challenge for pathologists. Moreover, the HCV recurrence disease severity is highly variable post-LT. HCV recurrence disease progression is characterized by an accelerated fibrogenesis process, and almost 30% of those patients develop cirrhosis at 5-years of follow-up. Whole-genome gene expression (WGE) analyses through well-defined oligonucleotide microarray platforms represent a powerful tool for the molecular characterization of biological process. In the present manuscript, the utility of microarray technology is applied for the ACR and HCV-recurrence biological characterization in post-LT liver biopsy samples. Moreover, WGE analysis was performed to identify predictive biomarkers of HCV recurrence severity in formalin-fixed paraffin-embedded liver biopsies prospectively collected