Methods to study splicing from high-throughput RNA Sequencing data

The development of novel high-throughput sequencing (HTS) methods for RNA (RNA-Seq) has provided a very powerful mean to study splicing under multiple conditions at unprecedented depth. However, the complexity of the information to be analyzed has turned this into a challenging task. In the last few years, a plethora of tools have been developed, allowing researchers to process RNA-Seq data to study the expression of isoforms and splicing events, and their relative changes under different conditions. We provide an overview of the methods available to study splicing from short RNA-Seq data. We group the methods according to the different questions they address: 1) Assignment of the sequencing reads to their likely gene of origin. This is addressed by methods that map reads to the genome and/or to the available gene annotations. 2) Recovering the sequence of splicing events and isoforms. This is addressed by transcript reconstruction and de novo assembly methods. 3) Quantification of events and isoforms. Either after reconstructing transcripts or using an annotation, many methods estimate the expression level or the relative usage of isoforms and/or events. 4) Providing an isoform or event view of differential splicing or expression. These include methods that compare relative event/isoform abundance or isoform expression across two or more conditions. 5) Visualizing splicing regulation. Various tools facilitate the visualization of the RNA-Seq data in the context of alternative splicing. In this review, we do not describe the specific mathematical models behind each method. Our aim is rather to provide an overview that could serve as an entry point for users who need to decide on a suitable tool for a specific analysis. We also attempt to propose a classification of the tools according to the operations they do, to facilitate the comparison and choice of methods.Comment: 31 pages, 1 figure, 9 tables. Small corrections adde

Patient preference and acceptability of calcium plus vitamin D3 supplementation: a randomised, open, cross-over trial

Preference for a drug formulation is important in adherence to long-term medication for chronic illnesses such as osteoporosis. We investigated the preference for and acceptability of chewable tablet containing calcium and vitamin D (Calci Chew D3, Nycomed) compared to that of a sachet containing calcium and vitamin D3 (Cad, Will-Pharma). This open, randomised, cross-over trial was set up to compare the preference and acceptability of two calcium plus vitamin D3 formulations (both with 500 mg calcium and 400/440 IU vitamin D3), given twice a day in patients with osteoporosis. Preference and acceptability were assessed by means of questionnaires. Preference was determined by asking the question, which treatment the patient preferred, and acceptability was measured by scoring five variables, using rating scales. Of the 102 patients indicating a preference for a trial medication, 67% preferred the chewable tablet, 19% the sachet with calcium and vitamin D3, and 15% stated no preference. The significant preference for Calci Chew D3 (p < 0.0001) was associated with higher scores for all five acceptability variables. The two formulations were tolerated equally well. A significant greater number of patients considered the chewable tablet as preferable and acceptable to the sachet, containing calcium and vitamin D3. Trial registration: Current Controlled Trials ISRCTN18822358

Search for Gravitational Waves from Primordial Black Hole Binary Coalescences in the Galactic Halo

We use data from the second science run of the LIGO gravitational-wave detectors to search for the gravitational waves from primordial black hole (PBH) binary coalescence with component masses in the range 0.2--$1.0 M_\odot$. The analysis requires a signal to be found in the data from both LIGO observatories, according to a set of coincidence criteria. No inspiral signals were found. Assuming a spherical halo with core radius 5 kpc extending to 50 kpc containing non-spinning black holes with masses in the range 0.2--$1.0 M_\odot$, we place an observational upper limit on the rate of PBH coalescence of 63 per year per Milky Way halo (MWH) with 90% confidence.Comment: 7 pages, 4 figures, to be submitted to Phys. Rev.

HCMV Targets the Metabolic Stress Response through Activation of AMPK Whose Activity Is Important for Viral Replication

Human Cytomegalovirus (HCMV) infection induces several metabolic activities that have been found to be important for viral replication. The cellular AMP-activated protein kinase (AMPK) is a metabolic stress response kinase that regulates both energy-producing catabolic processes and energy-consuming anabolic processes. Here we explore the role AMPK plays in generating an environment conducive to HCMV replication. We find that HCMV infection induces AMPK activity, resulting in the phosphorylation and increased abundance of several targets downstream of activated AMPK. Pharmacological and RNA-based inhibition of AMPK blocked the glycolytic activation induced by HCMV-infection, but had little impact on the glycolytic pathway of uninfected cells. Furthermore, inhibition of AMPK severely attenuated HCMV replication suggesting that AMPK is an important cellular factor for HCMV replication. Inhibition of AMPK attenuated early and late gene expression as well as viral DNA synthesis, but had no detectable impact on immediate-early gene expression, suggesting that AMPK activity is important at the immediate early to early transition of viral gene expression. Lastly, we find that inhibition of the Ca2+-calmodulin-dependent kinase kinase (CaMKK), a kinase known to activate AMPK, blocks HCMV-mediated AMPK activation. The combined data suggest a model in which HCMV activates AMPK through CaMKK, and depends on their activation for high titer replication, likely through induction of a metabolic environment conducive to viral replication

Fast Coding of Orientation in Primary Visual Cortex

Understanding how populations of neurons encode sensory information is a major goal of systems neuroscience. Attempts to answer this question have focused on responses measured over several hundred milliseconds, a duration much longer than that frequently used by animals to make decisions about the environment. How reliably sensory information is encoded on briefer time scales, and how best to extract this information, is unknown. Although it has been proposed that neuronal response latency provides a major cue for fast decisions in the visual system, this hypothesis has not been tested systematically and in a quantitative manner. Here we use a simple ‘race to threshold’ readout mechanism to quantify the information content of spike time latency of primary visual (V1) cortical cells to stimulus orientation. We find that many V1 cells show pronounced tuning of their spike latency to stimulus orientation and that almost as much information can be extracted from spike latencies as from firing rates measured over much longer durations. To extract this information, stimulus onset must be estimated accurately. We show that the responses of cells with weak tuning of spike latency can provide a reliable onset detector. We find that spike latency information can be pooled from a large neuronal population, provided that the decision threshold is scaled linearly with the population size, yielding a processing time of the order of a few tens of milliseconds. Our results provide a novel mechanism for extracting information from neuronal populations over the very brief time scales in which behavioral judgments must sometimes be made

Depletion of Murine Intestinal Microbiota: Effects on Gut Mucosa and Epithelial Gene Expression

Background Inappropriate cross talk between mammals and their gut microbiota may trigger intestinal inflammation and drive extra-intestinal immune-mediated diseases. Epithelial cells constitute the interface between gut microbiota and host tissue, and may regulate host responses to commensal enteric bacteria. Gnotobiotic animals represent a powerful approach to study bacterial-host interaction but are not readily accessible to the wide scientific community. We aimed at refining a protocol that in a robust manner would deplete the cultivable intestinal microbiota of conventionally raised mice and that would prove to have significant biologic validity. Methodology/Principal Findings Previously published protocols for depleting mice of their intestinal microbiota by administering broad-spectrum antibiotics in drinking water were difficult to reproduce. We show that twice daily delivery of antibiotics by gavage depleted mice of their cultivable fecal microbiota and reduced the fecal bacterial DNA load by 400 fold while ensuring the animals' health. Mice subjected to the protocol for 17 days displayed enlarged ceca, reduced Peyer's patches and small spleens. Antibiotic treatment significantly reduced the expression of antimicrobial factors to a level similar to that of germ-free mice and altered the expression of 517 genes in total in the colonic epithelium. Genes involved in cell cycle were significantly altered concomitant with reduced epithelial proliferative activity in situ assessed by Ki-67 expression, suggesting that commensal microbiota drives cellular proliferation in colonic epithelium. Conclusion We present a robust protocol for depleting conventionally raised mice of their cultivatable intestinal microbiota with antibiotics by gavage and show that the biological effect of this depletion phenocopies physiological characteristics of germ-free mice

The Bacterium Endosymbiont of Crithidia deanei Undergoes Coordinated Division with the Host Cell Nucleus

In trypanosomatids, cell division involves morphological changes and requires coordinated replication and segregation of the nucleus, kinetoplast and flagellum. In endosymbiont-containing trypanosomatids, like Crithidia deanei, this process is more complex, as each daughter cell contains only a single symbiotic bacterium, indicating that the prokaryote must replicate synchronically with the host protozoan. In this study, we used light and electron microscopy combined with three-dimensional reconstruction approaches to observe the endosymbiont shape and division during C. deanei cell cycle. We found that the bacterium replicates before the basal body and kinetoplast segregations and that the nucleus is the last organelle to divide, before cytokinesis. In addition, the endosymbiont is usually found close to the host cell nucleus, presenting different shapes during the protozoan cell cycle. Considering that the endosymbiosis in trypanosomatids is a mutualistic relationship, which resembles organelle acquisition during evolution, these findings establish an excellent model for the understanding of mechanisms related with the establishment of organelles in eukaryotic cells

Fortunella margarita Transcriptional Reprogramming Triggered by Xanthomonas citri subsp. citri

<p>Abstract</p> <p>Background</p> <p>Citrus canker disease caused by the bacterial pathogen <it>Xanthomonas citri </it>subsp. <it>citri (</it>Xcc) <it>has </it>become endemic in areas where high temperature, rain, humidity, and windy conditions provide a favourable environment for the dissemination of the bacterium. Xcc is pathogenic on many commercial citrus varieties but appears to elicit an incompatible reaction on the citrus relative <it>Fortunella margarita </it>Swing (kumquat), in the form of a very distinct delayed necrotic response. We have developed subtractive libraries enriched in sequences expressed in kumquat leaves during both early and late stages of the disease. The isolated differentially expressed transcripts were subsequently sequenced. Our results demonstrate how the use of microarray expression profiling can help assign roles to previously uncharacterized genes and elucidate plant pathogenesis-response related mechanisms. This can be considered to be a case study in a citrus relative where high throughput technologies were utilized to understand defence mechanisms in <it>Fortunella </it>and citrus at the molecular level.</p> <p>Results</p> <p><b>cDNAs from sequenced kumquat libraries (ESTs) made from subtracted RNA populations, healthy vs. infected, were used to make this microarray</b>. Of 2054 selected genes on a customized array, 317 were differentially expressed (P < 0.05) in Xcc challenged kumquat plants compared to mock-inoculated ones. This study identified components of the incompatible interaction such as reactive oxygen species (ROS) and programmed cell death (PCD). Common defence mechanisms and a number of resistance genes were also identified. In addition, there were a considerable number of differentially regulated genes that had no homologues in the databases. This could be an indication of either a specialized set of genes employed by kumquat in response to canker disease or new defence mechanisms in citrus.</p> <p>Conclusion</p> <p>Functional categorization of kumquat Xcc-responsive genes revealed an enhanced defence-related metabolism as well as a number of resistant response-specific genes in the kumquat transcriptome in response to Xcc inoculation. Gene expression profile(s) were analyzed to assemble a comprehensive and inclusive image of the molecular interaction in the kumquat/Xcc system. This was done in order to elucidate molecular mechanisms associated with the development of the hypersensitive response phenotype in kumquat leaves. These data will be used to perform comparisons among citrus species to evaluate means to enhance the host immune responses against bacterial diseases.</p