A bedr way of genomic interval processing

Background

Next-generation sequencing is making it critical to robustly and rapidly handle genomic ranges within standard pipelines. Standard use-cases include annotating sequence ranges with gene or other genomic annotation, merging multiple experiments together and subsequently quantifying and visualizing the overlap. The most widely-used tools for these tasks work at the command-line (e.g. BEDTools) and the small number of available R packages are either slow or have distinct semantics and features from command-line interfaces.

Results

To provide a robust R-based interface to standard command-line tools for genomic coordinate manipulation, we created bedr. This open-source R package can use either BEDTools or BEDOPS as a back-end and performs data-manipulation extremely quickly, creating R data structures that can be readily interfaced with existing computational pipelines. It includes data-visualization capabilities and a number of data-access functions that interface with standard databases like UCSC and COSMIC.

Conclusions

bedr package provides an open source solution to enable genomic interval data manipulation and restructuring in R programming language which is commonly used in bioinformatics, and therefore would be useful to bioinformaticians and genomic researchers.

     

Intrauterine Growth Restriction Alters Mouse Intestinal Architecture during Development

Infants with intrauterine growth restriction (IUGR) are at increased risk for neonatal and lifelong morbidities affecting multiple organ systems including the intestinal tract. The underlying mechanisms for the risk to the intestine remain poorly understood. In this study, we tested the hypothesis that IUGR affects the development of goblet and Paneth cell lineages, thus compromising the innate immunity and barrier functions of the epithelium. Using a mouse model of maternal thromboxane A₂-analog infusion to elicit maternal hypertension and resultant IUGR, we tested whether IUGR alters ileal maturation and specifically disrupts mucus-producing goblet and antimicrobial-secreting Paneth cell development. We measured body weights, ileal weights and ileal lengths from birth to postnatal day (P) 56. We also determined the abundance of goblet and Paneth cells and their mRNA products, localization of cellular tight junctions, cell proliferation, and apoptosis to interrogate cellular homeostasis. Comparison of the murine findings with human IUGR ileum allowed us to verify observed changes in the mouse were relevant to clinical IUGR. At P14 IUGR mice had decreased ileal lengths, fewer goblet and Paneth cells, reductions in Paneth cell specific mRNAs, and decreased cell proliferation. These findings positively correlated with severity of IUGR. Furthermore, the decrease in murine Paneth cells was also seen in human IUGR ileum. IUGR disrupts the normal trajectory of ileal development, particularly affecting the composition and secretory products of the epithelial surface of the intestine. We speculate that this abnormal intestinal development may constitute an inherent “first hit”, rendering IUGR intestine susceptible to further injury, infection, or inflammation.     

Calcium dependent plasticity applied to repetitive transcranial magnetic stimulation with a neural field model

The calcium dependent plasticity (CaDP) approach to the modeling of synaptic weight change is applied using a neural field approach to realistic repetitive transcranial magnetic stimulation (rTMS) protocols. A spatially-symmetric nonlinear neural field model consisting of populations of excitatory and inhibitory neurons is used. The plasticity between excitatory cell populations is then evaluated using a CaDP approach that incorporates metaplasticity. The direction and size of the plasticity (potentiation or depression) depends on both the amplitude of stimulation and duration of the protocol. The breaks in the inhibitory theta-burst stimulation protocol are crucial to ensuring that the stimulation bursts are potentiating in nature. Tuning the parameters of a spike-timing dependent plasticity (STDP) window with a Monte Carlo approach to maximize agreement between STDP predictions and the CaDP results reproduces a realistically-shaped window with two regions of depression in agreement with the existing literature. Developing understanding of how TMS interacts with cells at a network level may be important for future investigation.     

SOX2-Dependent Transcription in Clock Neurons Promotes the Robustness of the Central Circadian Pacemaker

1. Download high-res image (196KB)
2. Download full-size image

     

Non-genomic effects of aldosterone on intracellular ion regulation and cell volume in rat ventricular myocytes

Aldosterone has non-genomic effects that express within minutes and modulate intracellular ion milieu and cellular function. However, it is still undefined whether aldosterone actually alters intracellular ion concentrations or cellular contractility. To clarify the non-genomic effects of aldosterone, we measured [Na+]i, Ca2+ transient (CaT), and cell volume in dye-loaded rat ventricular myocytes, and we also evaluated myocardial contractility. We found the following: (i) aldosterone increased [Na+]i at the concentrations of 100 nmol/L to 10 micromol/L; (ii) aldosterone (up to 10 micromol/L) did not alter CaT and cell shortening in isolated myocytes, developed tension in papillary muscles, or left ventricular developed pressure in Langendorff-perfused hearts; (iii) aldosterone (100 nmol/L) increased the cell volume from 47.5 +/- 3.6 pL to 49.8 +/- 3.7 pL (n=8, p<0.05); (iv) both the increases in [Na+]i and cell volume were blocked by a Na+-K+-2Cl- co-transporter (NKCCl) inhibitor, bumetanide, or by a Na+/H+ exchange (NHE) inhibitor, 5-(N-ethyl-N-isopropyl) amiloride; and (v) spironolactone by itself increased in [Na+]i and cell volume. In conclusion, aldosterone rapidly increased [Na+]i and cell volume via NKCC1 and NHE, whereas there were no changes in CaT or myocardial contractility. Hence the non-genomic effects of aldosterone may be related to cell swelling rather than the increase in contractility.     

crm-1 facilitates BMP signaling to control body size in Caenorhabditis elegans

We have identified in Caenorhabditis elegans a homologue of the vertebrate Crim1, crm-1, which encodes a putative transmembrane protein with multiple cysteine-rich (CR) domains known to have bone morphogenetic proteins (BMPs) binding activity. Using the body morphology of C. elegans as an indicator, we showed that attenuation of crm-1 activity leads to a small body phenotype reminiscent of that of BMP pathway mutants. We showed that the crm-1 loss-of-function phenotype can be rescued by constitutive supply of sma-4 activity. crm-1 can enhance BMP signaling and this activity is dependent on the presence of the DBL-1 ligand and its receptors. crm-1 is expressed in neurons at the ventral nerve cord, where the DBL-1 ligand is produced. However, ectopic expression experiments reveal that crm-1 gene products act outside the DBL-1 producing cells and function non-autonomously to facilitate dbl/sma pathway signaling to control body size.     

Correlation of gene expression with physiological functions: Examples of pulmonary blood vessel rheology, hypoxic hypertension, and tissue remodeling

Microarray gene chip technology is a powerful invention looking for applications. A general principle is proposed here to direct the power of the technology toward physiology, medicine, and pharmacology. Our principle is to match quantitative measures of gene expression with the trend of mathematical parameters that describe biological functions. Mathematical parameterization is the heart. The procedure is illustrated by lung physiology, including the hypoxic hypertension, rheological properties of the tissues, and the remodeling of the pulmonary arterial wall under hypertensive stress. We show first how to reduce the experimental results on these physiological functions into mathematical formulas, and how the parameters of these formulas describe the functional trends precisely. Then under the assumption that the microarray reveals gene activities quantitatively, we match the trends of the gene activity with the trends of the functional parameters. Genes whose trends do match are interpreted as relevant to the functions. Those that do not match are considered irrelevant to the functions. The more functions we consider, the fewer will be the number of genes that are relevant to all functions. Thus we learn about the generality and specificity of the influence of genes on physiology.     

An improved synthesis of the saponin, polyphyllin D

Polyphyllin D, namely diosgenyl alpha-L-rhamnopyranosyl-(1 --> 2)- [(alpha-L-arabinofuranosyl)-(1 --> 4)]-[beta-D-glucopyranoside, was synthesized from diosgenyl-beta-D-glucopyranoside in four steps and in 30% overall yield, taking advantage of regioselective pivaloylation and alpha-L-rhamnopyranosylation reactions.