Straightforward Bibliography Management in R with the RefManageR Package

This work introduces the R package RefManageR, which provides tools for importing and working with bibliographic references. It extends the bibentry class in R in a number of useful ways, including providing R with previously unavailable support for BibLaTeX. BibLaTeX provides a superset of the functionality of BibTeX, including full Unicode support, no memory limitations, additional fields and entry types, and more sophisticated sorting of references. RefManageR provides functions for citing and generating a bibliography with hyperlinks for documents prepared with RMarkdown or RHTML. Existing .bib files can be read into R and converted from BibTeX to BibLaTeX and vice versa. References can also be imported via queries to NCBI's Entrez, Zotero libraries, Google Scholar, and CrossRef. Additionally, references can be created by reading PDFs stored on the user's machine with the help of Poppler. Entries stored in the reference manager can be easily searched by any field, by date ranges, and by various formats for name lists (author by last names, translator by full names, etc.). Entries can also be updated, combined, sorted, printed in a number of styles, and exported

Functional generalized additive models

We introduce the functional generalized additive model (FGAM), a novel regression model for association studies between a scalar response and a functional predictor. We model the link-transformed mean response as the integral with respect to t of F{X(t), t} where F(·, ·) is an unknown regression function and X(t) is a functional covariate. Rather than having an additive model in a finite number of principal components as by Müller and Yao (2008), our model incorporates the functional predictor directly and thus our model can be viewed as the natural functional extension of generalized additive models. We estimate F(·, ·) using tensor-product B-splines with roughness penalties. A pointwise quantile transformation of the functional predictor is also considered to ensure each tensor-product B-spline has observed data on its support. The methods are evaluated using simulated data and their predictive performance is compared with other competing scalar-on-function regression alternatives. We illustrate the usefulness of our approach through an application to brain tractography, where X(t) is a signal from diffusion tensor imaging at position, t, along a tract in the brain. In one example, the response is disease-status (case or control) and in a second example, it is the score on a cognitive test. The FGAM is implemented in R in the refund package. There are additional supplementary materials available online. © 2013 American Statistical Association, Institute of Mathematical Statistics, and Interface Foundation of North America

Bayesian Functional Generalized Additive Models with Sparsely Observed Covariates

The functional generalized additive model (FGAM) was recently proposed in McLean et al. (2013) as a more flexible alternative to the common functional linear model (FLM) for regressing a scalar on functional covariates. In this paper, we develop a Bayesian version of FGAM for the case of Gaussian errors with identity link function. Our approach allows the functional covariates to be sparsely observed and measured with error, whereas the estimation procedure of McLean et al. (2013) required that they be noiselessly observed on a regular grid. We consider both Monte Carlo and variational Bayes methods for fitting the FGAM with sparsely observed covariates. Due to the complicated form of the model posterior distribution and full conditional distributions, standard Monte Carlo and variational Bayes algorithms cannot be used. The strategies we use to handle the updating of parameters without closed-form full conditionals should be of independent interest to applied Bayesian statisticians working with nonconjugate models. Our numerical studies demonstrate the benefits of our algorithms over a two-step approach of first recovering the complete trajectories using standard techniques and then fitting a functional regression model. In a real data analysis, our methods are applied to forecasting closing price for items up for auction on the online auction website eBay

Generation and quality control of lipidomics data for the alzheimers disease neuroimaging initiative cohort.

Alzheimers disease (AD) is a major public health priority with a large socioeconomic burden and complex etiology. The Alzheimer Disease Metabolomics Consortium (ADMC) and the Alzheimer Disease Neuroimaging Initiative (ADNI) aim to gain new biological insights in the disease etiology. We report here an untargeted lipidomics of serum specimens of 806 subjects within the ADNI1 cohort (188 AD, 392 mild cognitive impairment and 226 cognitively normal subjects) along with 83 quality control samples. Lipids were detected and measured using an ultra-high-performance liquid chromatography quadruple/time-of-flight mass spectrometry (UHPLC-QTOF MS) instrument operated in both negative and positive electrospray ionization modes. The dataset includes a total 513 unique lipid species out of which 341 are known lipids. For over 95% of the detected lipids, a relative standard deviation of better than 20% was achieved in the quality control samples, indicating high technical reproducibility. Association modeling of this dataset and available clinical, metabolomics and drug-use data will provide novel insights into the AD etiology. These datasets are available at the ADNI repository at http://adni.loni.usc.edu/

Arterial oxygen content is precisely maintained by graded erythrocytotic responses in settings of high/normal serum iron levels, and predicts exercise capacity: an observational study of hypoxaemic patients with pulmonary arteriovenous malformations.

Oxygen, haemoglobin and cardiac output are integrated components of oxygen transport: each gram of haemoglobin transports 1.34 mls of oxygen in the blood. Low arterial partial pressure of oxygen (PaO2), and haemoglobin saturation (SaO2), are the indices used in clinical assessments, and usually result from low inspired oxygen concentrations, or alveolar/airways disease. Our objective was to examine low blood oxygen/haemoglobin relationships in chronically compensated states without concurrent hypoxic pulmonary vasoreactivity.165 consecutive unselected patients with pulmonary arteriovenous malformations were studied, in 98 cases, pre/post embolisation treatment. 159 (96%) had hereditary haemorrhagic telangiectasia. Arterial oxygen content was calculated by SaO2 x haemoglobin x 1.34/100.There was wide variation in SaO2 on air (78.5-99, median 95)% but due to secondary erythrocytosis and resultant polycythaemia, SaO2 explained only 0.1% of the variance in arterial oxygen content per unit blood volume. Secondary erythrocytosis was achievable with low iron stores, but only if serum iron was high-normal: Low serum iron levels were associated with reduced haemoglobin per erythrocyte, and overall arterial oxygen content was lower in iron deficient patients (median 16.0 [IQR 14.9, 17.4]mls/dL compared to 18.8 [IQR 17.4, 20.1]mls/dL, p<0.0001). Exercise tolerance appeared unrelated to SaO2 but was significantly worse in patients with lower oxygen content (p<0.0001). A pre-defined athletic group had higher Hb:SaO2 and serum iron:ferritin ratios than non-athletes with normal exercise capacity. PAVM embolisation increased SaO2, but arterial oxygen content was precisely restored by a subsequent fall in haemoglobin: 86 (87.8%) patients reported no change in exercise tolerance at post-embolisation follow-up.Haemoglobin and oxygen measurements in isolation do not indicate the more physiologically relevant oxygen content per unit blood volume. This can be maintained for SaO2 ≥78.5%, and resets to the same arterial oxygen content after correction of hypoxaemia. Serum iron concentrations, not ferritin, seem to predict more successful polycythaemic responses

Fin-Tail Coordination during Escape and Predatory Behavior in Larval Zebrafish

Larval zebrafish innately perform a suite of behaviors that are tightly linked to their evolutionary past, notably escape from threatening stimuli and pursuit and capture of prey. These behaviors have been carefully examined in the past, but mostly with regard to the movements of the trunk and tail of the larvae. Here, we employ kinematics analyses to describe the movements of the pectoral fins during escape and predatory behavior. In accord with previous studies, we find roles for the pectoral fins in slow swimming and immediately after striking prey. We find novel roles for the pectoral fins in long-latency, but not in short-latency C-bends. We also observe fin movements that occur during orienting J-turns and S-starts that drive high-velocity predatory strikes. Finally, we find that the use of pectoral fins following a predatory strike is scaled to the velocity of the strike, supporting a role for the fins in braking. The implications of these results for central control of coordinated movements are discussed, and we hope that these results will provide baselines for future analyses of cross-body coordination using mutants, morphants, and transgenic approaches

A Multi-Variant, Viral Dynamic Model of Genotype 1 HCV to Assess the in vivo Evolution of Protease-Inhibitor Resistant Variants

Variants resistant to compounds specifically targeting HCV are observed in clinical trials. A multi-variant viral dynamic model was developed to quantify the evolution and in vivo fitness of variants in subjects dosed with monotherapy of an HCV protease inhibitor, telaprevir. Variant fitness was estimated using a model in which variants were selected by competition for shared limited replication space. Fitness was represented in the absence of telaprevir by different variant production rate constants and in the presence of telaprevir by additional antiviral blockage by telaprevir. Model parameters, including rate constants for viral production, clearance, and effective telaprevir concentration, were estimated from 1) plasma HCV RNA levels of subjects before, during, and after dosing, 2) post-dosing prevalence of plasma variants from subjects, and 3) sensitivity of variants to telaprevir in the HCV replicon. The model provided a good fit to plasma HCV RNA levels observed both during and after telaprevir dosing, as well as to variant prevalence observed after telaprevir dosing. After an initial sharp decline in HCV RNA levels during dosing with telaprevir, HCV RNA levels increased in some subjects. The model predicted this increase to be caused by pre-existing variants with sufficient fitness to expand once available replication space increased due to rapid clearance of wild-type (WT) virus. The average replicative fitness estimates in the absence of telaprevir ranged from 1% to 68% of WT fitness. Compared to the relative fitness method, the in vivo estimates from the viral dynamic model corresponded more closely to in vitro replicon data, as well as to qualitative behaviors observed in both on-dosing and long-term post-dosing clinical data. The modeling fitness estimates were robust in sensitivity analyses in which the restoration dynamics of replication space and assumptions of HCV mutation rates were varied

Elution of gentamicin and vancomycin from polymethylmethacrylate beads and hip spacers in vivo

Background and purpose Late infections after total hip arthroplasty are still a problem. Treatment procedures include resection arthroplasty with implantation of antibiotic-loaded beads or implantation of an antibiotic-impreganted spacer. However, little is known about antibiotic elution from bone cement beyond the first 2–3 postoperative days in humans

CRISPR/Cas9 DNA cleavage at SNP-derived PAM enables both in vitro and in vivo KRT12 mutation-specific targeting

CRISPR/Cas9-based therapeutics hold the possibility for permanent treatment of genetic disease. The potency and specificity of this system has been used to target dominantly inherited conditions caused by heterozygous missense mutations through inclusion of the mutated base in the short-guide RNA (sgRNA) sequence. This research evaluates a novel approach for targeting heterozygous single-nucleotide polymorphisms (SNPs) using CRISPR/Cas9. We determined that a mutation within KRT12, which causes Meesmann's epithelial corneal dystrophy (MECD), leads to the occurrence of a novel protospacer adjacent motif (PAM). We designed an sgRNA complementary to the sequence adjacent to this SNP-derived PAM and evaluated its potency and allele specificity both in vitro and in vivo. This sgRNA was found to be highly effective at reducing the expression of mutant KRT12 mRNA and protein in vitro. To assess its activity in vivo we injected a combined Cas9/sgRNA expression construct into the corneal stroma of a humanized MECD mouse model. Sequence analysis of corneal genomic DNA revealed non-homologous end-joining repair resulting in frame-shifting deletions within the mutant KRT12 allele. This study is the first to demonstrate in vivo gene editing of a heterozygous disease-causing SNP that results in a novel PAM, further highlighting the potential for CRISPR/Cas9-based therapeutics

Search for new phenomena in final states with an energetic jet and large missing transverse momentum in pp collisions at √ s = 8 TeV with the ATLAS detector

Results of a search for new phenomena in final states with an energetic jet and large missing transverse momentum are reported. The search uses 20.3 fb−1 of √ s = 8 TeV data collected in 2012 with the ATLAS detector at the LHC. Events are required to have at least one jet with pT > 120 GeV and no leptons. Nine signal regions are considered with increasing missing transverse momentum requirements between Emiss T > 150 GeV and Emiss T > 700 GeV. Good agreement is observed between the number of events in data and Standard Model expectations. The results are translated into exclusion limits on models with either large extra spatial dimensions, pair production of weakly interacting dark matter candidates, or production of very light gravitinos in a gauge-mediated supersymmetric model. In addition, limits on the production of an invisibly decaying Higgs-like boson leading to similar topologies in the final state are presente