NS-kNN: a modified <Emphasis Type="Italic">k</Emphasis>-nearest neighbors approach for imputing metabolomics data

Introduction

A common problem in metabolomics data analysis is the existence of a substantial number of missing values, which can complicate, bias, or even prevent certain downstream analyses. One of the most widely-used solutions to this problem is imputation of missing values using a k-nearest neighbors (kNN) algorithm to estimate missing metabolite abundances. kNN implicitly assumes that missing values are uniformly distributed at random in the dataset, but this is typically not true in metabolomics, where many values are missing because they are below the limit of detection of the analytical instrumentation.

Objectives

Here, we explore the impact of nonuniformly distributed missing values (missing not at random, or MNAR) on imputation performance. We present a new model for generating synthetic missing data and a new algorithm, No-Skip kNN (NS-kNN), that accounts for MNAR values to provide more accurate imputations.

Methods

We compare the imputation errors of the original kNN algorithm using two distance metrics, NS-kNN, and a recently developed algorithm KNN-TN, when applied to multiple experimental datasets with different types and levels of missing data.

Results

Our results show that NS-kNN typically outperforms kNN when at least 20–30% of missing values in a dataset are MNAR. NS-kNN also has lower imputation errors than KNN-TN on realistic datasets when at least 50% of missing values are MNAR.

Conclusion

Accounting for the nonuniform distribution of missing values in metabolomics data can significantly improve the results of imputation algorithms. The NS-kNN method imputes missing metabolomics data more accurately than existing kNN-based approaches when used on realistic datasets.

     

In Vitro Studies on the Putative Function of N-acetylaspartate as an Osmoregulator

Efflux and tissue content of N-acetylaspartate (NAA) and amino acids were evaluated from cultured and acutely prepared hippocampal slices in response to changes in osmolarity. The osmoregulator taurine, but not NAA, was lost from both types of slices after moderate reductions in extracellular osmolarity (−60 mOsm) for 10–48 h. Hypoosmotic shock (−166 mOsm) for 5 min resulted in unselective efflux of several amino acids from acutely prepared slices. Notably, the efflux of taurine, but not NAA, was prominent also after the shock. Efflux of NAA was markedly enhanced by NMDA and high K⁺, in particular after the stimulation period. The high K⁺-mediated efflux was decreased by high extracellular osmolarity and a NMDA-receptor antagonist. The results indicate that NAA efflux can be induced by a sudden non-physiological decrease in extracellular osmolarity but not by prolonged more moderate changes in osmolarity. The mechanisms behind the efflux of NAA by high K⁺ are complex and may involve both swelling and activation of NMDA-receptors.     

Characterization of Upstream Sequences of the <Emphasis Type=Italic>LOJ</Emphasis> Gene Leads to Identification of a Novel Enhancer Element Conferring Lateral Organ Junction-Specific Expression in <Emphasis Type=Italic>Arabidopsis thaliana</Emphasis>

Isolation and characterization of promoters are important in understanding gene regulation and genetic engineering of crop plants. Earlier, a pentatricopeptide repeat protein (PPR) encoding gene (At2g39230), designated as Lateral Organ Junction (LOJ) gene, was identified through T-DNA promoter trapping in Arabidopsis thaliana. The upstream sequence of the LOJ gene conferred on the reporter gene a novel LOJ-specific expression. The present study was aimed at identifying and characterizing the cis-regulatory motifs responsible for tissue-specific expression in the −673 and +90 bases upstream of the LOJ gene recognized as LOJ promoter. In silico analysis of the LOJ promoter revealed the presence of a few relevant regulatory motifs and a unique feature like AT-rich inverted repeat. Deletion analysis of the LOJ promoter confirmed the presence of an enhancer-like element in the distal region (−673/−214), which stimulates a minimal promoter-like sequence in the −424/−214 region in a position and orientation autonomous manner. The −136/+90 region of the LOJ promoter was efficient in driving reporter gene expression in tissues like developing anthers and seeds of Arabidopsis. A positive regulation for the seed- and anther-specific expression module was contemplated within the 5′ untranslated region of the LOJ gene. However, this function was repressed in the native context by the lateral organ junction-specific expression. The present study has led to the identification of a novel lateral organ junction-specific element and an enhancer sequence in Arabidopsis with potential applications in plant genetic engineering.     

Germination responses of three grassland species differ between native and invasive origins

The germination stage is critical in plant life-history and is also a key process during the expansion of species’ ranges into new environments. In this study we investigated the germination patterns of three plant species (Achillea millefolium, Hieracium pilosella and Hypericum perforatum) that are invasive to New Zealand (NZ) and native to Central Europe. We asked whether the species show differences in germination temperature requirements, germination speed and maximum germination rates, and thus, whether they display evidence of adaptation to different conditions in the invasive range. Seeds from three populations per species and region were subjected to three different temperature regimes to compare germination rates among origins and across temperature conditions. For Achillea millefolium and Hypericum perforatum, germination rates were significantly higher for invasive NZ provenances than for native German ones. Seeds from invasive populations of all three species displayed increased maximum germination at medium temperature conditions when compared to native populations, which indicates altered germination strategies in the invaded range. Changes in temporal development patterns were most conspicuous for invasive Hieracium pilosella and Hypericum perforatum populations. These findings imply that adaptation in germination patterns towards different climatic conditions in invasive populations has occurred. Our study emphasises the importance of the germination stage during plant invasion and its role in explaining range expansion of these species.     

A landscape of shifting-mosaic steady state in Lassen Volcanic National Park,California

It is generally perceived that landscape patterns or textures in a given protected area are spatially stationary. The findings of this study suggest that this common perception is only partially correct. Over the course of 52 years, equilibrium in landscape shifting was detected using digital data for the Lassen Volcanic National Park (USA). Vertical aerial photographs taken of the park in 1941 were geo-referenced with the digital orthophoto quarter quadrangle (DOQQ) images of the same area from 1993 to identify landscape compositions and to measure change. Spatial analysis was used to observe pattern changes over time. The results suggested that landscape development maintained equilibrium while patches were in various stages of a successional sequence. The total area of each landscape component held steady, although over time patches throughout the landscape changed—a shifting-mosaic steady state (SMSS). These findings reflect the limitations of contemporary environmental conservation theory. They also suggest the importance of considering landscape change in policies that currently govern park planning and management.     

Long-term population dynamics in a Mediterranean aquatic snake

A review of several long-term studies has recently suggested that snakes might be declining in large parts of the world. Additional data from other long-term studies are therefore urgently needed in order to assess the generalities of such suggested declines. Based on a 20-year study, we analyzed demographic data on adult dice snakes (Natrix tessellata) studied in central Italy between 1985 and 2004. Both male and female dice snakes were relatively long-lived, with no significant differences in longevity between the sexes. Individual males and females were observed over a maximum of 10 and 14 years, respectively. However, the among-year recapture rates between the year the snakes were initially captured and the subsequent year (i.e., year 1 to year 2) was significantly lower (45%) than the among-year recapture rates during subsequent years (74%; i.e., year 2 to year 3), suggesting that a large proportion of the snakes at first capture were in fact not resident within our study area, and hence many snakes were migrating in and out of our 2-km stream study site, with no inter-sexual difference in dispersal rates. Sex ratio was virtually equal if we consider the study period as a whole. Significant annual fluctuations were, however, observed through the study. In 1985–1990, 1993–1995, 1998 and 1999 the sex ratio was male-biased, whereas in 2000–2004 female-biased. In terms of both survival and recapture probabilities, model selection showed that Akaike’s information criterion favored the model incorporating body size, with the model incorporating year having an intermediate likelihood, and the model with sex included being the most disfavored. Total population number estimates suggest an average 86 adult individuals along the 2 km of stream with only minor annual variations. However, a significant decrease in the number of males occurred during the last 6 years of our study. Thus, further monitoring of this population is warranted in light of the decline of snake populations reported recently.     

Inflammation-associated S100 proteins: new mechanisms that regulate function

This review focuses on new aspects of extracellular roles of the calgranulins. S100A8, S100A9 and S100A12 are constitutively expressed in neutrophils and induced in several cell types. The S100A8 and S100A9 genes are regulated by pro- and anti-inflammatory mediators and their functions may depend on cell type, mediators within a particular inflammatory milieu, receptors involved in their recognition and their post-translational modification. The S100A8 gene induction in macrophages is dependent on IL-10 and potentiated by immunosuppressive agents. S100A8 and S100A9 are oxidized by peroxide, hypochlorite and nitric oxide (NO). HOCl generates intra-chain sulfinamide bonds; stronger oxidation promotes cross-linked forms that are seen in human atheroma. S100A8 is >200-fold more sensitive to oxidative cross-linking than low-density lipoprotein and may reduce oxidative damage. S100A8 and S100A9 can be S-nitrosylated. S100A8–SNO suppresses mast cell activation and inflammation in the microcirculation and may act as an NO transporter to regulate vessel tone in inflammatory lesions. S100A12 activates mast cells and is a monocyte and mast cell chemoattractant; a G-protein-coupled mechanism may be involved. Structure–function studies are discussed in relation to conservation and divergence of functions in S100A8. S100A12 induces cytokines in mast cells, but not monocytes/macrophages. It forms complexes with Zn²⁺ and, by chelating Zn²⁺, S100A12 significantly inhibits MMPs. Zn²⁺ in S100A12 complexes co-localize with MMP-9 in foam cells in atheroma. In summary, S100A12 has pro-inflammatory properties that are likely to be stable in an oxidative environment, because it lacks Cys and Met residues. Conversely, S100A8 and S100A9 oxidation and S-nitrosylation may have important protective mechanisms in inflammation.     

Human Alzheimer’s disease synaptic <Emphasis Type=Italic>O</Emphasis>-GlcNAc site mapping and iTRAQ expression proteomics with ion trap mass spectrometry

Neuronal synaptic functional deficits are linked to impaired learning and memory in Alzheimer’s disease (AD). We recently demonstrated that O-GlcNAc, a novel cytosolic and nuclear carbohydrate post-translational modification, is enriched at neuronal synapses and positively regulates synaptic plasticity linked to learning and memory in mice. Reduced levels of O-GlcNAc have been observed in AD, suggesting a possible link to deficits in synaptic plasticity. Using lectin enrichment and mass spectrometry, we mapped several human cortical synaptic O-GlcNAc modification sites. Overlap in patterns of O-GlcNAcation between mouse and human appears to be high, as previously mapped mouse synaptic O-GlcNAc sites in Bassoon, Piccolo, and tubulin polymerization promoting protein p25 were identified in human. Novel O-GlcNAc modification sites were identified on Mek2 and RPN13/ADRM1. Mek2 is a signaling component of the Erk 1/2 pathway involved in synaptic plasticity. RPN13 is a component of the proteasomal degradation pathway. The potential interplay of phosphorylation with mapped O-GlcNAc sites, and possible implication of those sites in synaptic plasticity in normal versus AD states is discussed. iTRAQ is a powerful differential isotopic quantitative approach in proteomics. Pulsed Q dissociation (PQD) is a recently introduced fragmentation strategy that enables detection of low mass iTRAQ reporter ions in ion trap mass spectrometry. We optimized LTQ ion trap settings for PQD-based iTRAQ quantitation and demonstrated its utility in O-GlcNAc site mapping. Using iTRAQ, abnormal synaptic expression levels of several proteins previously implicated in AD pathology were observed in addition to novel changes in synaptic specific protein expression including Synapsin II.     

Phase-resetting curve determines how BK currents affect neuronal firing

BK channels are large conductance potassium channels gated by calcium and voltage. Paradoxically, blocking these channels has been shown experimentally to increase or decrease the firing rate of neurons, depending on the neural subtype and brain region. The mechanism for how this current can alter the firing rates of different neurons remains poorly understood. Using phase-resetting curve (PRC) theory, we determine when BK channels increase or decrease the firing rates in neural models. The addition of BK currents always decreases the firing rate when the PRC has only a positive region. When the PRC has a negative region (type II), BK currents can increase the firing rate. The influence of BK channels on firing rate in the presence of other conductances, such as I _m and I _h , as well as with different amplitudes of depolarizing input, were also investigated. These results provide a formal explanation for the apparently contradictory effects of BK channel antagonists on firing rates.     

Vanadyl bisacetylacetonate protects β cells from palmitate-induced cell death through the unfolded protein response pathway

Endoplasmic reticulum (ER) stress induced by free fatty acids (FFA) is important to β-cell loss during the development of type 2 diabetes. To test whether vanadium compounds could influence ER stress and the responses in their mechanism of antidiabetic effects, we investigated the effects and the mechanism of vanadyl bisacetylacetonate [VO(acac)₂] on β cells upon treatment with palmitate, a typical saturated FFA. The experimental results showed that VO(acac)₂ could enhance FFA-induced signaling pathways of unfolded protein responses by upregulating the prosurvival chaperone immunoglobulin heavy-chain binding protein/78-kDa glucose-regulated protein and downregulating the expression of apoptotic C/EBP homologous protein, and consequently the reduction of insulin synthesis. VO(acac)₂ also ameliorated FFA-disturbed Ca²⁺ homeostasis in β cells. Overall, VO(acac)₂ enhanced stress adaption, thus protecting β cells from palmitate-induced apoptosis. This study provides some new insights into the mechanisms of antidiabetic vanadium compounds.