Bayesian nonparametric discovery of isoforms and individual specific quantification

Most human protein-coding genes can be transcribed into multiple distinct mRNA isoforms. These alternative splicing patterns encourage molecular diversity, and dysregulation of isoform expression plays an important role in disease etiology. However, isoforms are difficult to characterize from short-read RNA-seq data because they share identical subsequences and occur in different frequencies across tissues and samples. Here, we develop BIISQ, a Bayesian nonparametric model for isoform discovery and individual specific quantification from short-read RNA-seq data. BIISQ does not require isoform reference sequences but instead estimates an isoform catalog shared across samples. We use stochastic variational inference for efficient posterior estimates and demonstrate superior precision and recall for simulations compared to state-of-the-art isoform reconstruction methods. BIISQ shows the most gains for low abundance isoforms, with 36% more isoforms correctly inferred at low coverage versus a multi-sample method and 170% more versus single-sample methods. We estimate isoforms in the GEUVADIS RNA-seq data and validate inferred isoforms by associating genetic variants with isoform ratios

Reverse Engineering Gene Networks with ANN: Variability in Network Inference Algorithms

Motivation :Reconstructing the topology of a gene regulatory network is one of the key tasks in systems biology. Despite of the wide variety of proposed methods, very little work has been dedicated to the assessment of their stability properties. Here we present a methodical comparison of the performance of a novel method (RegnANN) for gene network inference based on multilayer perceptrons with three reference algorithms (ARACNE, CLR, KELLER), focussing our analysis on the prediction variability induced by both the network intrinsic structure and the available data. Results: The extensive evaluation on both synthetic data and a selection of gene modules of "Escherichia coli" indicates that all the algorithms suffer of instability and variability issues with regards to the reconstruction of the topology of the network. This instability makes objectively very hard the task of establishing which method performs best. Nevertheless, RegnANN shows MCC scores that compare very favorably with all the other inference methods tested. Availability: The software for the RegnANN inference algorithm is distributed under GPL3 and it is available at the corresponding author home page (http://mpba.fbk.eu/grimaldi/regnann-supmat

HIV-1 Nef increases astrocyte sensitivity towards exogenous hydrogen peroxide

<p>Abstract</p> <p>Background</p> <p>HIV-1 infected individuals are under chronic exposure to reactive oxygen species (ROS) considered to be instrumental in the progression of AIDS and the development of HIV-1 associated dementia (HAD). Astrocytes support neuronal function and protect them against cytotoxic substances including ROS. The protein HIV-1 Nef, a progression factor in AIDS pathology is abundantly expressed in astrocytes in patients with HAD, and thus may influence its functions.</p> <p>Results</p> <p>Endogenous expressed HIV-1 Nef leads to increased sensitivity of human astrocytes towards exogenous hydrogen peroxide but not towards TNF-alpha. Cell death of <it>nef</it>-expressing astrocytes exposed to 10 μM hydrogen peroxide for 30 min occurred within 4 h.</p> <p>Conclusion</p> <p>HIV-1 Nef may contribute to neuronal dysfunction and the development of HAD by causing death of astrocytes through decreasing their tolerance for hydrogen peroxide.</p

Scuba:Scalable kernel-based gene prioritization

Abstract Background The uncovering of genes linked to human diseases is a pressing challenge in molecular biology and precision medicine. This task is often hindered by the large number of candidate genes and by the heterogeneity of the available information. Computational methods for the prioritization of candidate genes can help to cope with these problems. In particular, kernel-based methods are a powerful resource for the integration of heterogeneous biological knowledge, however, their practical implementation is often precluded by their limited scalability. Results We propose Scuba, a scalable kernel-based method for gene prioritization. It implements a novel multiple kernel learning approach, based on a semi-supervised perspective and on the optimization of the margin distribution. Scuba is optimized to cope with strongly unbalanced settings where known disease genes are few and large scale predictions are required. Importantly, it is able to efficiently deal both with a large amount of candidate genes and with an arbitrary number of data sources. As a direct consequence of scalability, Scuba integrates also a new efficient strategy to select optimal kernel parameters for each data source. We performed cross-validation experiments and simulated a realistic usage setting, showing that Scuba outperforms a wide range of state-of-the-art methods. Conclusions Scuba achieves state-of-the-art performance and has enhanced scalability compared to existing kernel-based approaches for genomic data. This method can be useful to prioritize candidate genes, particularly when their number is large or when input data is highly heterogeneous. The code is freely available at https://github.com/gzampieri/Scuba

Impacts of savanna trees on forage quality for a large African herbivore

Recently, cover of large trees in African savannas has rapidly declined due to elephant pressure, frequent fires and charcoal production. The reduction in large trees could have consequences for large herbivores through a change in forage quality. In Tarangire National Park, in Northern Tanzania, we studied the impact of large savanna trees on forage quality for wildebeest by collecting samples of dominant grass species in open grassland and under and around large Acacia tortilis trees. Grasses growing under trees had a much higher forage quality than grasses from the open field indicated by a more favourable leaf/stem ratio and higher protein and lower fibre concentrations. Analysing the grass leaf data with a linear programming model indicated that large savanna trees could be essential for the survival of wildebeest, the dominant herbivore in Tarangire. Due to the high fibre content and low nutrient and protein concentrations of grasses from the open field, maximum fibre intake is reached before nutrient requirements are satisfied. All requirements can only be satisfied by combining forage from open grassland with either forage from under or around tree canopies. Forage quality was also higher around dead trees than in the open field. So forage quality does not reduce immediately after trees die which explains why negative effects of reduced tree numbers probably go initially unnoticed. In conclusion our results suggest that continued destruction of large trees could affect future numbers of large herbivores in African savannas and better protection of large trees is probably necessary to sustain high animal densities in these ecosystems

Incorporating Existing Network Information into Gene Network Inference

One methodology that has met success to infer gene networks from gene expression data is based upon ordinary differential equations (ODE). However new types of data continue to be produced, so it is worthwhile to investigate how to integrate these new data types into the inference procedure. One such data is physical interactions between transcription factors and the genes they regulate as measured by ChIP-chip or ChIP-seq experiments. These interactions can be incorporated into the gene network inference procedure as a priori network information. In this article, we extend the ODE methodology into a general optimization framework that incorporates existing network information in combination with regularization parameters that encourage network sparsity. We provide theoretical results proving convergence of the estimator for our method and show the corresponding probabilistic interpretation also converges. We demonstrate our method on simulated network data and show that existing network information improves performance, overcomes the lack of observations, and performs well even when some of the existing network information is incorrect. We further apply our method to the core regulatory network of embryonic stem cells utilizing predicted interactions from two studies as existing network information. We show that including the prior network information constructs a more closely representative regulatory network versus when no information is provided

Facilitation or Competition? Tree Effects on Grass Biomass across a Precipitation Gradient

Savanna ecosystems are dominated by two distinct plant life forms, grasses and trees, but the interactions between them are poorly understood. Here, we quantified the effects of isolated savanna trees on grass biomass as a function of distance from the base of the tree and tree height, across a precipitation gradient in the Kruger National Park, South Africa. Our results suggest that mean annual precipitation (MAP) mediates the nature of tree-grass interactions in these ecosystems, with the impact of trees on grass biomass shifting qualitatively between 550 and 737 mm MAP. Tree effects on grass biomass were facilitative in drier sites (MAP≤550 mm), with higher grass biomass observed beneath tree canopies than outside. In contrast, at the wettest site (MAP = 737 mm), grass biomass did not differ significantly beneath and outside tree canopies. Within this overall precipitation-driven pattern, tree height had positive effect on sub-canopy grass biomass at some sites, but these effects were weak and not consistent across the rainfall gradient. For a more synthetic understanding of tree-grass interactions in savannas, future studies should focus on isolating the different mechanisms by which trees influence grass biomass, both positively and negatively, and elucidate how their relative strengths change over broad environmental gradients. © 2013 Moustakas et al

Chemokine Transfer by Liver Sinusoidal Endothelial Cells Contributes to the Recruitment of CD4+ T Cells into the Murine Liver

Leukocyte adhesion and transmigration are central features governing immune surveillance and inflammatory reactions in body tissues. Within the liver sinusoids, chemokines initiate the first crucial step of T-cell migration into the hepatic tissue. We studied molecular mechanisms involved in endothelial chemokine supply during hepatic immune surveillance and liver inflammation and their impact on the recruitment of CD4+ T cells into the liver. In the murine model of Concanavalin A-induced T cell-mediated hepatitis, we showed that hepatic expression of the inflammatory CXC chemokine ligands (CXCL)9 and CXCL10 strongly increased whereas homeostatic CXCL12 significantly decreased. Consistently, CD4+ T cells expressing the CXC chemokine receptor (CXCR)3 accumulated within the inflamed liver tissue. In histology, CXCL9 was associated with liver sinusoidal endothelial cells (LSEC) which represent the first contact site for T-cell immigration into the liver. LSEC actively transferred basolaterally internalized CXCL12, CXCL9 and CXCL10 via clathrin- coated vesicles to CD4+ T cells leading to enhanced transmigration of CXCR4+ total CD4+ T cells and CXCR3+ effector/memory CD4+ T cells, respectively in vitro. LSEC-expressed CXCR4 mediated CXCL12 transport and blockage of endothelial CXCR4 inhibited CXCL12-dependent CD4+ T-cell transmigration. In contrast, CXCR3 was not involved in the endothelial transport of its ligands CXCL9 and CXCL10. The clathrin-specific inhibitor chlorpromazine blocked endothelial chemokine internalization and CD4+ T-cell transmigration in vitro as well as migration of CD4+ T cells into the inflamed liver in vivo. Moreover, hepatic accumulation of CXCR3+ CD4+ T cells during T cell-mediated hepatitis was strongly reduced after administration of chlorpromazine. These data demonstrate that LSEC actively provide perivascularly expressed homeostatic and inflammatory chemokines by CXCR4- and clathrin-dependent intracellular transport mechanisms thereby contributing to the hepatic recruitment of CD4+ T-cell populations during immune surveillance and liver inflammation