A wrapper generation tool for the creation of scriptable scientific applications

Journal ArticleIn recent years, there has been considerable interest in the use of scripting languages as a mechanism for controlling and developing scientific software. Scripting languages allow scientific applications to be encapsulated in an interpreted environment similar to that found in commercial scientific packages such as MATLAB, Mathematica, and IDL. This improves the usability of scientific software by providing a powerful meachanism for specifyling and controlling cimplex problems as well as giving users an interactive and exploratory problem solving environment. Scripting languages also provide a framework for building and integrating software components that allows tools be used in a more efficient manner. This streamlines the problem solving process and enable scientists to be more productive

SWIG users manual (version 1.1)

technical reportSWIG is a tool for solving problems. More specifically, SWIG is a simple tool for building interactive C, C++, or Objective-C programs with common scripting languages such as Tel, Perl, and Python. Of course, more importantly, SWIG is a tool for making C programming more enjoyable and promoting laziness (an essential feature). SWIG is not part of an overgrown software engineering project, an attempt to build some sort of monolithic programming environment, or an attempt to force everyone to rewrite all of their code (ie. code reuse). In fact, none of these things have ever been a priority. SWIG was originally developed in the Theoretical Physics Division at Los Alamos National Laboratory for building interfaces to large materials science research simulations being run on the Connection Machine 5 supercomputer

Training Big Random Forests with Little Resources

Without access to large compute clusters, building random forests on large datasets is still a challenging problem. This is, in particular, the case if fully-grown trees are desired. We propose a simple yet effective framework that allows to efficiently construct ensembles of huge trees for hundreds of millions or even billions of training instances using a cheap desktop computer with commodity hardware. The basic idea is to consider a multi-level construction scheme, which builds top trees for small random subsets of the available data and which subsequently distributes all training instances to the top trees' leaves for further processing. While being conceptually simple, the overall efficiency crucially depends on the particular implementation of the different phases. The practical merits of our approach are demonstrated using dense datasets with hundreds of millions of training instances.Comment: 9 pages, 9 Figure

Lightweight computational steering of very large scale molecular dynamics simulations

We present a computational steering approach for controlling, analyzing, and visualizing very large scale molecular dynamics simulations involving tens to hundreds of millions of atoms. Our approach relies on extensible scripting languages and an easy to use tool for building extensions and modules. The system is extremely easy to modify, works with existing C code, is memory efficient, and can be used from inexpensive workstations and networks. We demonstrate how we have used this system to manipulate data from production MD simulations involving as many as 104 million atoms running on the CM-5 and Cray T3D. We also show how this approach can be used to build systems that integrate common scripting languages (including Tcl/Tk, Perl, and Python), simulation code, user extensions, and commercial data analysis packages

Sensory neuronal sensitisation occurs through HMGB-1/RAGE and TRPV1 in high glucose conditions

Many potential causes for painful diabetic neuropathy have been proposed including actions of cytokines and growth factors. High mobility group protein B1 (HMGB1) is a RAGE agonist, increased in diabetes, that contributes to pain by modulating peripheral inflammatory responses. HMGB1 enhances nociceptive behaviour in naïve animals through an unknown mechanism. We tested the hypothesis that HMGB1 causes pain through direct neuronal activation of RAGE and alteration of nociceptive neuronal responsiveness. HMGB1 and RAGE expression were increased in skin and primary sensory (DRG) neurons of diabetic rats at times when pain behaviour was enhanced. Agonist-evoked TRPV1-mediated calcium responses increased in cultured DRG neurons from diabetic rats and in neurons from naïve rats exposed to high glucose concentrations. HMGB1-mediated increases in TRPV1-evoked calcium responses in DRG neurons were RAGE and PKC-dependent, and this was blocked by co-administration of the growth factor splice variant, VEGF-A165b. Pain behaviour and DRG RAGE expression increases were blocked by VEGF-A165b treatment of diabetic rats in vivo. HMGB-1-RAGE activation sensitizes DRG neurons in vitro. VEGF-A165b blocks HMGB-1/RAGE DRG activation, which may contribute to its analgesic properties in vivo

Vascular endothelial growth factor-A165b prevents diabetic neuropathic pain and sensory neuronal degeneration

Diabetic peripheral neuropathy affects up to half of diabetic patients. This neuronal damage leads to sensory disturbances, including allodynia and hyperalgesia. Many growth factors have been suggested as useful treatments for prevention of neurodegeneration, including the vascular endothelial growth factor (VEGF) family. VEGF-A is generated as two alternative splice variant families. The most widely studied isoform, VEGF-A165a is both pro-angiogenic and neuroprotective, but pro-nociceptive and increases vascular permeability in animal models. Streptozotocin (STZ)-induced diabetic rats develop both hyperglycaemia and many of the resulting diabetic complications seen in patients, including peripheral neuropathy. In the present study, we show that the anti-angiogenic VEGF-A splice variant, VEGF-A165b, is also a potential therapeutic for diabetic neuropathy. Seven weeks of VEGF-A165b treatment in diabetic rats reversed enhanced pain behaviour in multiple behavioural paradigms and was neuroprotective, reducing hyperglycaemia-induced activated caspase 3 (AC3) levels in sensory neuronal subsets, epidermal sensory nerve fibre loss and aberrant sciatic nerve morphology. Furthermore, VEGF-A165b inhibited a STZ-induced increase in Evans Blue extravasation in dorsal root ganglia (DRG), saphenous nerve and plantar skin of the hind paw. Increased transient receptor potential ankyrin 1 (TRPA1) channel activity is associated with the onset of diabetic neuropathy. VEGF-A165b also prevented hyperglycaemia-enhanced TRPA1 activity in an in vitro sensory neuronal cell line indicating a novel direct neuronal mechanism that could underlie the anti-nociceptive effect observed in vivo. These results demonstrate that in a model of Type I diabetes VEGF-A165b attenuates altered pain behaviour and prevents neuronal stress, possibly through an effect on TRPA1 activity

Vascular endothelial growth factor-A 165 b ameliorates outer-retinal barrier and vascular dysfunction in the diabetic retina

Diabetic retinopathy (DR) is one of the leading causes of blindness in the developed world. Characteristic features of DR are retinal neurodegeneration, pathological angiogenesis and breakdown of both the inner and outer retinal barriers of the retinal vasculature and retinal pigmented epithelial (RPE)–choroid respectively. Vascular endothelial growth factor (VEGF-A), a key regulator of angiogenesis and permeability, is the target of most pharmacological interventions of DR. VEGF-A can be alternatively spliced at exon 8 to form two families of isoforms, pro- and anti-angiogenic. VEGF-A165a is the most abundant pro-angiogenic isoform, is pro-inflammatory and a potent inducer of permeability. VEGF-A165b is anti-angiogenic, anti-inflammatory, cytoprotective and neuroprotective. In the diabetic eye, pro-angiogenic VEGF-A isoforms are up-regulated such that they overpower VEGF-A165b. We hypothesized that this imbalance may contribute to increased breakdown of the retinal barriers and by redressing this imbalance, the pathological angiogenesis, fluid extravasation and retinal neurodegeneration could be ameliorated. VEGF-A165b prevented VEGF-A165a and hyperglycaemia-induced tight junction (TJ) breakdown and subsequent increase in solute flux in RPE cells. In streptozotocin (STZ)-induced diabetes, there was an increase in Evans Blue extravasation after both 1 and 8 weeks of diabetes, which was reduced upon intravitreal and systemic delivery of recombinant human (rh)VEGF-A165b. Eight-week diabetic rats also showed an increase in retinal vessel density, which was prevented by VEGF-A165b. These results show rhVEGF-A165b reduces DR-associated blood–retina barrier (BRB) dysfunction, angiogenesis and neurodegeneration and may be a suitable therapeutic in treating DR

Patterns of Cis Regulatory Variation in Diverse Human Populations

The genetic basis of gene expression variation has long been studied with the aim to understand the landscape of regulatory variants, but also more recently to assist in the interpretation and elucidation of disease signals. To date, many studies have looked in specific tissues and population-based samples, but there has been limited assessment of the degree of inter-population variability in regulatory variation. We analyzed genome-wide gene expression in lymphoblastoid cell lines from a total of 726 individuals from 8 global populations from the HapMap3 project and correlated gene expression levels with HapMap3 SNPs located in cis to the genes. We describe the influence of ancestry on gene expression levels within and between these diverse human populations and uncover a non-negligible impact on global patterns of gene expression. We further dissect the specific functional pathways differentiated between populations. We also identify 5,691 expression quantitative trait loci (eQTLs) after controlling for both non-genetic factors and population admixture and observe that half of the cis-eQTLs are replicated in one or more of the populations. We highlight patterns of eQTL-sharing between populations, which are partially determined by population genetic relatedness, and discover significant sharing of eQTL effects between Asians, European-admixed, and African subpopulations. Specifically, we observe that both the effect size and the direction of effect for eQTLs are highly conserved across populations. We observe an increasing proximity of eQTLs toward the transcription start site as sharing of eQTLs among populations increases, highlighting that variants close to TSS have stronger effects and therefore are more likely to be detected across a wider panel of populations. Together these results offer a unique picture and resource of the degree of differentiation among human populations in functional regulatory variation and provide an estimate for the transferability of complex trait variants across populations

Genetic mechanisms of critical illness in COVID-19.

Host-mediated lung inflammation is present1, and drives mortality2, in the critical illness caused by coronavirus disease 2019 (COVID-19). Host genetic variants associated with critical illness may identify mechanistic targets for therapeutic development3. Here we report the results of the GenOMICC (Genetics Of Mortality In Critical Care) genome-wide association study in 2,244 critically ill patients with COVID-19 from 208 UK intensive care units. We have identified and replicated the following new genome-wide significant associations: on chromosome 12q24.13 (rs10735079, P = 1.65 × 10-8) in a gene cluster that encodes antiviral restriction enzyme activators (OAS1, OAS2 and OAS3); on chromosome 19p13.2 (rs74956615, P = 2.3 × 10-8) near the gene that encodes tyrosine kinase 2 (TYK2); on chromosome 19p13.3 (rs2109069, P = 3.98 ×  10-12) within the gene that encodes dipeptidyl peptidase 9 (DPP9); and on chromosome 21q22.1 (rs2236757, P = 4.99 × 10-8) in the interferon receptor gene IFNAR2. We identified potential targets for repurposing of licensed medications: using Mendelian randomization, we found evidence that low expression of IFNAR2, or high expression of TYK2, are associated with life-threatening disease; and transcriptome-wide association in lung tissue revealed that high expression of the monocyte-macrophage chemotactic receptor CCR2 is associated with severe COVID-19. Our results identify robust genetic signals relating to key host antiviral defence mechanisms and mediators of inflammatory organ damage in COVID-19. Both mechanisms may be amenable to targeted treatment with existing drugs. However, large-scale randomized clinical trials will be essential before any change to clinical practice