Systematic evaluation of spliced alignment programs for RNA-seq data

High-throughput RNA sequencing is an increasingly accessible method for studying gene structure and activity on a genome-wide scale. A critical step in RNA-seq data analysis is the alignment of partial transcript reads to a reference genome sequence. To assess the performance of current mapping software, we invited developers of RNA-seq aligners to process four large human and mouse RNA-seq data sets. In total, we compared 26 mapping protocols based on 11 programs and pipelines and found major performance differences between methods on numerous benchmarks, including alignment yield, basewise accuracy, mismatch and gap placement, exon junction discovery and suitability of alignments for transcript reconstruction. We observed concordant results on real and simulated RNA-seq data, confirming the relevance of the metrics employed. Future developments in RNA-seq alignment methods would benefit from improved placement of multimapped reads, balanced utilization of existing gene annotation and a reduced false discovery rate for splice junctions

GENCODE: producing a reference annotation for ENCODE

BACKGROUND: The GENCODE consortium was formed to identify and map all protein-coding genes within the ENCODE regions. This was achieved by a combination of initial manual annotation by the HAVANA team, experimental validation by the GENCODE consortium and a refinement of the annotation based on these experimental results. RESULTS: The GENCODE gene features are divided into eight different categories of which only the first two (known and novel coding sequence) are confidently predicted to be protein-coding genes. 5' rapid amplification of cDNA ends (RACE) and RT-PCR were used to experimentally verify the initial annotation. Of the 420 coding loci tested, 229 RACE products have been sequenced. They supported 5' extensions of 30 loci and new splice variants in 50 loci. In addition, 46 loci without evidence for a coding sequence were validated, consisting of 31 novel and 15 putative transcripts. We assessed the comprehensiveness of the GENCODE annotation by attempting to validate all the predicted exon boundaries outside the GENCODE annotation. Out of 1,215 tested in a subset of the ENCODE regions, 14 novel exon pairs were validated, only two of them in intergenic regions. CONCLUSION: In total, 487 loci, of which 434 are coding, have been annotated as part of the GENCODE reference set available from the UCSC browser. Comparison of GENCODE annotation with RefSeq and ENSEMBL show only 40% of GENCODE exons are contained within the two sets, which is a reflection of the high number of alternative splice forms with unique exons annotated. Over 50% of coding loci have been experimentally verified by 5' RACE for EGASP and the GENCODE collaboration is continuing to refine its annotation of 1% human genome with the aid of experimental validation

Geophysical monitoring and reactive transport modeling of ureolytically-driven calcium carbonate precipitation

Ureolytically-driven calcium carbonate precipitation is the basis for a promising in-situ remediation method for sequestration of divalent radionuclide and trace metal ions. It has also been proposed for use in geotechnical engineering for soil strengthening applications. Monitoring the occurrence, spatial distribution, and temporal evolution of calcium carbonate precipitation in the subsurface is critical for evaluating the performance of this technology and for developing the predictive models needed for engineering application. In this study, we conducted laboratory column experiments using natural sediment and groundwater to evaluate the utility of geophysical (complex resistivity and seismic) sensing methods, dynamic synchrotron x-ray computed tomography (micro-CT), and reactive transport modeling for tracking ureolytically-driven calcium carbonate precipitation processes under site relevant conditions. Reactive transport modeling with TOUGHREACT successfully simulated the changes of the major chemical components during urea hydrolysis. Even at the relatively low level of urea hydrolysis observed in the experiments, the simulations predicted an enhanced calcium carbonate precipitation rate that was 3-4 times greater than the baseline level. Reactive transport modeling results, geophysical monitoring data and micro-CT imaging correlated well with reaction processes validated by geochemical data. In particular, increases in ionic strength of the pore fluid during urea hydrolysis predicted by geochemical modeling were successfully captured by electrical conductivity measurements and confirmed by geochemical data. The low level of urea hydrolysis and calcium carbonate precipitation suggested by the model and geochemical data was corroborated by minor changes in seismic P-wave velocity measurements and micro-CT imaging; the latter provided direct evidence of sparsely distributed calcium carbonate precipitation. Ion exchange processes promoted through NH4+ production during urea hydrolysis were incorporated in the model and captured critical changes in the major metal species. The electrical phase increases were potentially due to ion exchange processes that modified charge structure at mineral/water interfaces. Our study revealed the potential of geophysical monitoring for geochemical changes during urea hydrolysis and the advantages of combining multiple approaches to understand complex biogeochemical processes in the subsurface

Roles of Major Facilitator Superfamily Transporters in Phosphate Response in Drosophila

The major facilitator superfamily (MFS) transporter Pho84 and the type III transporter Pho89 are responsible for metabolic effects of inorganic phosphate in yeast. While the Pho89 ortholog Pit1 was also shown to be involved in phosphate-activated MAPK in mammalian cells, it is currently unknown, whether orthologs of Pho84 have a role in phosphate-sensing in metazoan species. We show here that the activation of MAPK by phosphate observed in mammals is conserved in Drosophila cells, and used this assay to characterize the roles of putative phosphate transporters. Surprisingly, while we found that RNAi-mediated knockdown of the fly Pho89 ortholog dPit had little effect on the activation of MAPK in Drosophila S2R+ cells by phosphate, two Pho84/SLC17A1–9 MFS orthologs (MFS10 and MFS13) specifically inhibited this response. Further, using a Xenopus oocyte assay, we show that MSF13 mediates uptake of [³³P]-orthophosphate in a sodium-dependent fashion. Consistent with a role in phosphate physiology, MSF13 is expressed highest in the Drosophila crop, midgut, Malpighian tubule, and hindgut. Altogether, our findings provide the first evidence that Pho84 orthologs mediate cellular effects of phosphate in metazoan cells. Finally, while phosphate is essential for Drosophila larval development, loss of MFS13 activity is compatible with viability indicating redundancy at the levels of the transporters.National Institutes of Health (U.S.) (NIDDK 5K08DK078361)Harvard Catalys

A cross-sectional study of the prevalence of intensity of infection with Schistosoma japonicum in 50 irrigated and rain-fed villages in Samar Province, the Philippines

BACKGROUND: Few studies have described heterogeneity in Schistosoma japonicum infection intensity, and none were done in Philippines. The purpose of this report is to describe the village-to-village variation in the prevalence of two levels of infection intensity across 50 villages of Samar Province, the Philippines. METHODS: This cross-sectional study was conducted in 25 rain-fed and 25 irrigated villages endemic for S. japonicum between August 2003 and November 2004. Villages were selected based on irrigation and farming criteria. A maximum of 35 eligible households were selected per village. Each participant was asked to provide stool samples on three consecutive days. All those who provided at least one stool sample were included in the analysis. A Bayesian three category outcome hierarchical cumulative logit regression model with adjustment for age, sex, occupation and measurement error of the Kato-Katz technique was used for analysis. RESULTS: A total of 1427 households and 6917 individuals agreed to participate in the study. A total of 5624 (81.3%) participants provided at least one stool sample. The prevalences of those lightly and at least moderately infected varied from 0% (95% Bayesian credible interval (BCI): 0%–3.1%) to 45.2% (95% BCI: 36.5%–53.9%) and 0% to 23.0% (95% BCI: 16.4%–31.2%) from village-to-village, respectively. Using the 0–7 year old group as a reference category, the highest odds ratio (OR) among males and females were that of being aged 17–40-year old (OR = 8.76; 95% BCI: 6.03–12.47) and 11–16-year old (OR = 8.59; 95% BCI: 4.74–14.28), respectively. People who did not work on a rice farm had a lower prevalence of infection than those working full time on a rice farm. The OR for irrigated villages compared to rain-fed villages was 1.41 (95% BCI: 0.50–3.21). DISCUSSION: We found very important village-to-village variation in prevalence of infection intensity. This variation is probably due to village-level variables other than that explained by a crude classification of villages into the irrigated and non-irrigated categories. We are planning to capture this spatial heterogeneity by updating our initial transmission dynamics model with the data reported here combined with 1-year post-treatment follow-up of study participants

Roles of residues in the interface of transient protein-protein complexes before complexation

Transient protein-protein interactions play crucial roles in all facets of cellular physiology. Here, using an analysis on known 3-D structures of transient protein-protein complexes, their corresponding uncomplexed forms and energy calculations we seek to understand the roles of protein-protein interfacial residues in the unbound forms. We show that there are conformationally near invariant and evolutionarily conserved interfacial residues which are rigid and they account for ∼65% of the core interface. Interestingly, some of these residues contribute significantly to the stabilization of the interface structure in the uncomplexed form. Such residues have strong energetic basis to perform dual roles of stabilizing the structure of the uncomplexed form as well as the complex once formed while they maintain their rigid nature throughout. This feature is evolutionarily well conserved at both the structural and sequence levels. We believe this analysis has general bearing in the prediction of interfaces and understanding molecular recognition

Prognostic factors for recovery following acute lateral ankle ligament sprain: a systematic review.

BACKGROUND: One-third of individuals who sustain an acute lateral ankle ligament sprain suffer significant disability due to pain, functional instability, mechanical instability or recurrent sprain after recovery plateaus at 1 to 5 years post injury. The identification of early prognostic factors associated with poor recovery may provide an opportunity for early-targeted intervention and improve outcome. METHODS: We performed a comprehensive search of AMED, EMBASE, Psych Info, CINAHL, SportDiscus, PubMed, CENTRAL, PEDro, OpenGrey, abstracts and conference proceedings from inception to September 2016. Prospective studies investigating the association between baseline prognostic factors and recovery over time were included. Two independent assessors performed the study selection, data extraction and quality assessment of the studies. A narrative synthesis is presented due to inability to meta-analyse results due to clinical and statistical heterogeneity. RESULTS: The search strategy yielded 3396 titles/abstracts after duplicates were removed. Thirty-six full text articles were then assessed, nine of which met the study inclusion criteria. Six were prospective cohorts, and three were secondary analyses of randomised controlled trials. Results are presented for nine studies that presented baseline prognostic factors for recovery after an acute ankle sprain. Age, female gender, swelling, restricted range of motion, limited weight bearing ability, pain (at the medial joint line and on weight-bearing dorsi-flexion at 4 weeks, and pain at rest at 3 months), higher injury severity rating, palpation/stress score, non-inversion mechanism injury, lower self-reported recovery, re-sprain within 3 months, MRI determined number of sprained ligaments, severity and bone bruise were found to be independent predictors of poor recovery. Age was one prognostic factor that demonstrated a consistent association with outcome in three studies, however cautious interpretation is advised. CONCLUSIONS: The associations between prognostic factors and poor recovery after an acute lateral ankle sprain are largely inconclusive. At present, there is insufficient evidence to recommend any factor as an independent predictor of outcome. There is a need for well-conducted prospective cohort studies with adequate sample size and long-term follow-up to provide robust evidence on prognostic factors of recovery following an acute lateral ankle sprain. TRIAL REGISTRATION: Prospero registration: CRD42014014471