Search

Cardiac Magnetic Resonance Imaging (MRI) allows quantifying myocardial tissue deformation and strain based on the tagging principle. In this work, we investigate accuracy and precision of strain quantification from synthetic 3D tagged MRI using equilibrated warping. To this end, synthetic biomechanical left-ventricular tagged MRI data with varying tag distance, spatial resolution and signal-to-noise ratio (SNR) were generated and processed to quantify errors in radial, circumferential and longitudinal strains relative to ground truth. Results reveal that radial strain is more sensitive to image resolution and noise than the other strain components. The study also shows robustness of quantifying circumferential and longitudinal strain in the presence of geometrical inconsistencies of 3D tagged data. In conclusion, our study points to the need for higher-resolution 3D tagged MRI than currently available in practice in order to achieve sufficient accuracy of radial strain quantification., PLoS ONE, 16 (11), ISSN:1932-6203

This work aims at investigating the interactions between the flow of fluids in the eyes and the brain and their potential implications in the development of visual impairment in astronauts, a condition also known as spaceflight associated neuro-ocular syndrome (SANS). To this end, we propose a reduced (0-dimensional) mathematical model of fluid flow in the eyes and brain, which is embedded into a simplified whole-body circulation model. In particular, the model accounts for: (i) the flows of blood and aqueous humor in the eyes; (ii) the flows of blood, cerebrospinal fluid and interstitial fluid in the brain; and (iii) their interactions. The model is used to simulate variations in intraocular pressure, intracranial pressure and blood flow due to microgravity conditions, which are thought to be critical factors in SANS. Specifically, the model predicts that both intracranial and intraocular pressures increase in microgravity, even though their respective trends may be different. In such conditions, ocular blood flow is predicted to decrease in the choroid and ciliary body circulations, whereas retinal circulation is found to be less susceptible to microgravity-induced alterations, owing to a purely mechanical component in perfusion control associated with the venous segments. These findings indicate that the particular anatomical architecture of venous drainage in the retina may be one of the reasons why most of the SANS alterations are not observed in the retina but, rather, in other vascular beds, particularly the choroid. Thus, clinical assessment of ocular venous function may be considered as a determinant SANS factor, for which astronauts could be screened on earth and in-flight.

Objective To demonstrate the incremental benefit of using free text data in addition to vital sign and demographic data to identify patients with suspected infection in the emergency department. Methods This was a retrospective, observational cohort study performed at a tertiary academic teaching hospital. All consecutive ED patient visits between 12/17/08 and 2/17/13 were included. No patients were excluded. The primary outcome measure was infection diagnosed in the emergency department defined as a patient having an infection related ED ICD-9-CM discharge diagnosis. Patients were randomly allocated to train (64%), validate (20%), and test (16%) data sets. After preprocessing the free text using bigram and negation detection, we built four models to predict infection, incrementally adding vital signs, chief complaint, and free text nursing assessment. We used two different methods to represent free text: a bag of words model and a topic model. We then used a support vector machine to build the prediction model. We calculated the area under the receiver operating characteristic curve to compare the discriminatory power of each model. Results A total of 230,936 patient visits were included in the study. Approximately 14% of patients had the primary outcome of diagnosed infection. The area under the ROC curve (AUC) for the vitals model, which used only vital signs and demographic data, was 0.67 for the training data set, 0.67 for the validation data set, and 0.67 (95% CI 0.65–0.69) for the test data set. The AUC for the chief complaint model which also included demographic and vital sign data was 0.84 for the training data set, 0.83 for the validation data set, and 0.83 (95% CI 0.81–0.84) for the test data set. The best performing methods made use of all of the free text. In particular, the AUC for the bag-of-words model was 0.89 for training data set, 0.86 for the validation data set, and 0.86 (95% CI 0.85–0.87) for the test data set. The AUC for the topic model was 0.86 for the training data set, 0.86 for the validation data set, and 0.85 (95% CI 0.84–0.86) for the test data set. Conclusion Compared to previous work that only used structured data such as vital signs and demographic information, utilizing free text drastically improves the discriminatory ability (increase in AUC from 0.67 to 0.86) of identifying infection.

This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Small molecule screens are widely used to prioritize pharmaceutical development. However, determining the pathways targeted by these molecules is challenging, since the compounds are often promiscuous. We present a network strategy that takes into account the polypharmacology of small molecules in order to generate hypotheses for their broader mode of action. We report a screen for kinase inhibitors that increase the efficacy of gemcitabine, the first-line chemotherapy for pancreatic cancer. Eight kinase inhibitors emerge that are known to affect 201 kinases, of which only three kinases have been previously identified as modifiers of gemcitabine toxicity. In this work, we use the SAMNet algorithm to identify pathways linking these kinases and genetic modifiers of gemcitabine toxicity with transcriptional and epigenetic changes induced by gemcitabine that we measure using DNaseI-seq and RNA-seq. SAMNet uses a constrained optimization algorithm to connect genes from these complementary datasets through a small set of protein-protein and protein-DNA interactions. The resulting network recapitulates known pathways including DNA repair, cell proliferation and the epithelial-to-mesenchymal transition. We use the network to predict genes with important roles in the gemcitabine response, including six that have already been shown to modify gemcitabine efficacy in pancreatic cancer and ten novel candidates. Our work reveals the important role of polypharmacology in the activity of these chemosensitiz-ing agents., National Institutes of Health (U.S.) (Grant R01-GM089903), National Institutes of Health (U.S.) (Grant U01-CA184898)

This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. We describe DNase-capture, an assay that increases the analytical resolution of DNase-seq by focusing its sequencing phase on selected genomic regions. We introduce a new method to compensate for capture bias called BaseNormal that allows for accurate recovery of transcription factor protection profiles from DNase-capture data. We show that these normalized data allow for nuanced detection of transcription factor binding heterogeneity with as few as dozens of sites., National Institute of General Medical Sciences (U.S.) (Grant 1U01HG007037)

Nonhuman primates are commonly used for cognitive neuroscience research and often surgically implanted with cephalic recording chambers for electrophysiological recording. Aerobic bacterial cultures from 25 macaques identified 72 bacterial isolates, including 15 Enterococcus faecalis isolates. The E. faecalis isolates displayed multi-drug resistant phenotypes, with resistance to ciprofloxacin, enrofloxacin, trimethoprim-sulfamethoxazole, tetracycline, chloramphenicol, bacitracin, and erythromycin, as well as high-level aminoglycoside resistance. Multi-locus sequence typing showed that most belonged to two E. faecalis sequence types (ST): ST 4 and ST 55. The genomes of three representative isolates were sequenced to identify genes encoding antimicrobial resistances and other traits. Antimicrobial resistance genes identified included aac(6’)-aph(2”), aph(3’)-III, str, ant(6)-Ia, tetM, tetS, tetL, ermB, bcrABR, cat, and dfrG, and polymorphisms in parC (S80I) and gyrA (S83I) were observed. These isolates also harbored virulence factors including the cytolysin toxin genes in ST 4 isolates, as well as multiple biofilm-associated genes (esp, agg, ace, SrtA, gelE, ebpABC), hyaluronidases (hylA, hylB), and other survival genes (ElrA, tpx). Crystal violet biofilm assays confirmed that ST 4 isolates produced more biofilm than ST 55 isolates. The abundance of antimicrobial resistance and virulence factor genes in the ST 4 isolates likely relates to the loss of CRISPR-cas. This macaque colony represents a unique model for studying E. faecalis infection associated with indwelling devices, and provides an opportunity to understand the basis of persistence of this pathogen in a healthcare setting., United States. National Institutes of Health (T32 OD010978), United States. National Institutes of Health (P30 ES02109)

We derive a physically realistic model for the generation of virtual transillumination, white light microscopy images using epi-fluorescence measurements from thick, unsectioned tissue. We demonstrate this technique by generating virtual transillumination H&E images of unsectioned human breast tissue from epi-fluorescence multiphoton microscopy data. The virtual transillumination algorithm is shown to enable improved contrast and color accuracy compared with previous color mapping methods. Finally, we present an open source implementation of the algorithm in OpenGL, enabling real-time GPU-based generation of virtual transillumination microscopy images using conventional fluorescence microscopy systems.

Purpose To characterize qualitative and quantitative features of the choroid in normal eyes using enface swept-source optical coherence tomography (SS-OCT). Methods Fifty-two eyes of 26 consecutive normal subjects were prospectively recruited to obtain multiple three-dimensional 12x12mm volumetric scans using a long-wavelength high-speed SS-OCT prototype. A motion-correction algorithm merged multiple SS-OCT volumes to improve signal. Retinal pigment epithelium (RPE) was segmented as the reference and enface images were extracted at varying depths every 4.13μm intervals. Systematic analysis of the choroid at different depths was performed to qualitatively assess the morphology of the choroid and quantify the absolute thicknesses as well as the relative thicknesses of the choroidal vascular layers including the choroidal microvasculature (choriocapillaris, terminal arterioles and venules; CC) and choroidal vessels (CV) with respect to the subfoveal total choroidal thickness (TC). Subjects were divided into two age groups: younger (, National Institutes of Health (U.S.) (R01-EY011289-27), National Institutes of Health (U.S.) (R44EY022864-01), National Institutes of Health (U.S.) (R01-CA075289-16), United States. Air Force Office of Scientific Research (FA9550-10-1-0551), United States. Air Force Office of Scientific Research (FA9550-10-1-0063)

The pulsatile release of cortisol from the adrenal glands is controlled by a hierarchical system that involves corticotropin releasing hormone (CRH) from the hypothalamus, adrenocorticotropin hormone (ACTH) from the pituitary, and cortisol from the adrenal glands. Determining the number, timing, and amplitude of the cortisol secretory events and recovering the infusion and clearance rates from serial measurements of serum cortisol levels is a challenging problem. Despite many years of work on this problem, a complete satisfactory solution has been elusive. We formulate this question as a non-convex optimization problem, and solve it using a coordinate descent algorithm that has a principled combination of (i) compressed sensing for recovering the amplitude and timing of the secretory events, and (ii) generalized cross validation for choosing the regularization parameter. Using only the observed serum cortisol levels, we model cortisol secretion from the adrenal glands using a second-order linear differential equation with pulsatile inputs that represent cortisol pulses released in response to pulses of ACTH. Using our algorithm and the assumption that the number of pulses is between 15 to 22 pulses over 24 hours, we successfully deconvolve both simulated datasets and actual 24-hr serum cortisol datasets sampled every 10 minutes from 10 healthy women. Assuming a one-minute resolution for the secretory events, we obtain physiologically plausible timings and amplitudes of each cortisol secretory event with R[superscript 2] above 0.92. Identification of the amplitude and timing of pulsatile hormone release allows (i) quantifying of normal and abnormal secretion patterns towards the goal of understanding pathological neuroendocrine states, and (ii) potentially designing optimal approaches for treating hormonal disorders., National Science Foundation (U.S.). Graduate Research Fellowship Program, National Institutes of Health (U.S.) (NIH DP1 OD003646), National Science Foundation (U.S.) (0836720), National Science Foundation (U.S.). Office of Emerging Frontiers in Research and Innovation (EFRI-0735956)

Background: Circulating tumor cells (CTCs) are cancer cells that can be isolated via liquid biopsy from blood and can be phenotypically and genetically characterized to provide critical information for guiding cancer treatment. Current analysis of CTCs is hindered by the throughput, selectivity and specificity of devices or assays used in CTC detection and isolation. Methodology/Principal Findings: Here, we enriched and characterized putative CTCs from blood samples of patients with both advanced stage metastatic breast and lung cancers using a novel multiplexed spiral microfluidic chip. This system detected putative CTCs under high sensitivity (100%, n = 56) (Breast cancer samples: 12–1275 CTCs/ml; Lung cancer samples: 10–1535 CTCs/ml) rapidly from clinically relevant blood volumes (7.5 ml under 5 min). Blood samples were completely separated into plasma, CTCs and PBMCs components and each fraction were characterized with immunophenotyping (Pan-cytokeratin/CD45, CD44/CD24, EpCAM), fluorescence in-situ hybridization (FISH) (EML4-ALK) or targeted somatic mutation analysis. We used an ultra-sensitive mass spectrometry based system to highlight the presence of an EGFR-activating mutation in both isolated CTCs and plasma cell-free DNA (cf-DNA), and demonstrate concordance with the original tumor-biopsy samples. Conclusions/Significance: We have clinically validated our multiplexed microfluidic chip for the ultra high-throughput, low-cost and label-free enrichment of CTCs. Retrieved cells were unlabeled and viable, enabling potential propagation and real-time downstream analysis using next generation sequencing (NGS) or proteomic analysis., Singapore-MIT Alliance for Research and Technology

AC electrokinetics is a versatile tool for contact-less manipulation or characterization of cells and has been widely used for separation based on genotype translation to electrical phenotypes. Cells responses to an AC electric field result in a complex combination of electrokinetic phenomena, mainly dielectrophoresis and electrohydrodynamic forces. Human cells behaviors to AC electrokinetics remain unclear over a large frequency spectrum as illustrated by the self-rotation effect observed recently. We here report and analyze human cells behaviors in different conditions of medium conductivity, electric field frequency and magnitude. We also observe the self-rotation of human cells, in the absence of a rotational electric field. Based on an analytical competitive model of electrokinetic forces, we propose an explanation of the cell self-rotation. These experimental results, coupled with our model, lead to the exploitation of the cell behaviors to measure the intrinsic dielectric properties of JURKAT, HEK and PC3 human cell lines.

Boolean implications (if-then rules) provide a conceptually simple, uniform and highly scalable way to find associations between pairs of random variables. In this paper, we propose to use Boolean implications to find relationships between variables of different data types (mutation, copy number alteration, DNA methylation and gene expression) from the glioblastoma (GBM) and ovarian serous cystadenoma (OV) data sets from The Cancer Genome Atlas (TCGA). We find hundreds of thousands of Boolean implications from these data sets. A direct comparison of the relationships found by Boolean implications and those found by commonly used methods for mining associations show that existing methods would miss relationships found by Boolean implications. Furthermore, many relationships exposed by Boolean implications reflect important aspects of cancer biology. Examples of our findings include cis relationships between copy number alteration, DNA methylation and expression of genes, a new hierarchy of mutations and recurrent copy number alterations, loss-of-heterozygosity of well-known tumor suppressors, and the hypermethylation phenotype associated with IDH1 mutations in GBM. The Boolean implication results used in the paper can be accessed at http://crookneck.stanford.edu/microarray​/TCGANetworks/., Stanford University (Chinese Undergraduate Visiting Research Program)

Objective: To optimally leverage the scalability and unique features of the electronic health records (EHR) for research that would ultimately improve patient care, we need to accurately identify patients and extract clinically meaningful measures. Using multiple sclerosis (MS) as a proof of principle, we showcased how to leverage routinely collected EHR data to identify patients with a complex neurological disorder and derive an important surrogate measure of disease severity heretofore only available in research settings. Methods: In a cross-sectional observational study, 5,495 MS patients were identified from the EHR systems of two major referral hospitals using an algorithm that includes codified and narrative information extracted using natural language processing. In the subset of patients who receive neurological care at a MS Center where disease measures have been collected, we used routinely collected EHR data to extract two aggregate indicators of MS severity of clinical relevance multiple sclerosis severity score (MSSS) and brain parenchymal fraction (BPF, a measure of whole brain volume). Results: The EHR algorithm that identifies MS patients has an area under the curve of 0.958, 83% sensitivity, 92% positive predictive value, and 89% negative predictive value when a 95% specificity threshold is used. The correlation between EHR-derived and true MSSS has a mean R[superscript 2] = 0.38±0.05, and that between EHR-derived and true BPF has a mean R[superscript 2] = 0.22±0.08. To illustrate its clinical relevance, derived MSSS captures the expected difference in disease severity between relapsing-remitting and progressive MS patients after adjusting for sex, age of symptom onset and disease duration (p = 1.56×10[superscript −12]). Conclusion: Incorporation of sophisticated codified and narrative EHR data accurately identifies MS patients and provides estimation of a well-accepted indicator of MS severity that is widely used in research settings but not part of the routine medical records. Similar approaches could be applied to other complex neurological disorders., National Institute of General Medical Sciences (U.S.) (NIH U54-LM008748)

IGF signaling is involved in cell proliferation, differentiation and apoptosis in a wide range of tissues, both normal and diseased, and so IGF-IR has been the focus of intense interest as a promising drug target. In this computational study on cartilage, we focus on two questions: (i) what are the key factors influencing IGF-IR complex formation, and (ii) how might cells regulate IGF-IR complex formation? We develop a reaction-diffusion computational model of the IGF system involving twenty three parameters. A series of parametric and sensitivity studies are used to identify the key factors influencing IGF signaling. From the model we predict the free IGF and IGF-IR complex concentrations throughout the tissue. We estimate the degradation half-lives of free IGF-I and IGFBPs in normal cartilage to be 20 and 100 mins respectively, and conclude that regulation of the IGF half-life, either directly or indirectly via extracellular matrix IGF-BP protease concentrations, are two critical factors governing the IGF-IR complex formation in the cartilage. Further we find that cellular regulation of IGF-II production, the IGF-IIR concentration and its clearance rate, all significantly influence IGF signaling. It is likely that negative feedback processes via regulation of these factors tune IGF signaling within a tissue, which may help explain the recent failures of single target drug therapies aimed at modifying IGF signaling., National Health and Medical Research Council (Australia) (APP1051455)

As neurons develop, several immature processes (i.e., neurites) grow out of the cell body. Over time, each neuron breaks symmetry when only one of its neurites grows much longer than the rest, becoming an axon. This symmetry breaking is an important step in neurodevelopment, and aberrant symmetry breaking is associated with several neuropsychiatric diseases, including schizophrenia and autism. However, the effects of neurite count in neuronal symmetry breaking have never been studied. Existing models for neuronal polarization disagree: some predict that neurons with more neurites polarize up to several days later than neurons with fewer neurites, while others predict that neurons with different neurite counts polarize synchronously. We experimentally find that neurons with different neurite counts polarize synchronously. We also show that despite the significant differences among the previously proposed models, they all agree with our experimental findings when the expression levels of the proteins responsible for symmetry breaking increase with neurite count. Consistent with these results, we observe that the expression levels of two of these proteins, HRas and shootin1, significantly correlate with neurite count. This coordinated symmetry breaking we observed among neurons with different neurite counts may be important for synchronized polarization of neurons in developing organisms., PLoS ONE, 8 (2), ISSN:1932-6203

Synthetic DNA has great propensity for efficiently and stably storing non-biological information. With DNA writing and reading technologies rapidly advancing, new applications for synthetic DNA are emerging in data storage and communication. Traditionally, DNA communication has focused on the encoding and transfer of complete sets of information. Here, we explore the use of DNA for the communication of short messages that are fragmented across multiple distinct DNA molecules. We identified three pivotal points in a communication-data encoding, data transfer & data extraction—and developed novel tools to enable communication via molecules of DNA. To address data encoding, we designed DNA-based individualized keyboards (iKeys) to convert plaintext into DNA, while reducing the occurrence of DNA homopolymers to improve synthesis and sequencing processes. To address data transfer, we implemented a secret-sharing system-Multiplexed Sequence Encoding (MuSE)-that conceals messages between multiple distinct DNA molecules, requiring a combination key to reveal messages. To address data extraction, we achieved the first instance of chromatogram patterning through multiplexed sequencing, thereby enabling a new method for data extraction. We envision these approaches will enable more widespread communication of information via DNA., National Institutes of Health (U.S.) (NIH Grant 1R01EB017755), National Institutes of Health (U.S.) (NIH Grant 1DP2OD008435), National Institutes of Health (U.S.) (NIH Grant 1P50GM098792), United States. Defense Advanced Research Projects Agency (Air Force Contract #FA8721-05-C-0002)

Using a nuclease-dead Cas9 mutant, we show that Cas9 reproducibly induces chromatin accessibility at previously inaccessible genomic loci. Cas9 chromatin opening is sufficient to enable adjacent binding and transcriptional activation by the settler transcription factor retinoic acid receptor at previously unbound motifs. Thus, we demonstrate a new use for Cas9 in increasing surrounding chromatin accessibility to alter local transcription factor binding., National Institutes of Health (U.S.) (5UL1DE019581), National Institutes of Health (U.S.) (RL1DE019021), National Institutes of Health (U.S.) (1K01DK101684-01), National Institutes of Health (U.S.) (1U01HG007037), National Institutes of Health (U.S.) (5P01NS055923)

Reverse-engineering of biological networks is a central problem in systems biology. The use of intervention data, such as gene knockouts or knockdowns, is typically used for teasing apart causal relationships among genes. Under time or resource constraints, one needs to carefully choose which intervention experiments to carry out. Previous approaches for selecting most informative interventions have largely been focused on discrete Bayesian networks. However, continuous Bayesian networks are of great practical interest, especially in the study of complex biological systems and their quantitative properties. In this work, we present an efficient, information-theoretic active learning algorithm for Gaussian Bayesian networks (GBNs), which serve as important models for gene regulatory networks. In addition to providing linear-algebraic insights unique to GBNs, leading to significant runtime improvements, we demonstrate the effectiveness of our method on data simulated with GBNs and the DREAM4 network inference challenge data sets. Our method generally leads to faster recovery of underlying network structure and faster convergence to final distribution of confidence scores over candidate graph structures using the full data, in comparison to random selection of intervention experiments.

International audience; Fifty years ago, FitzHugh introduced a phase portrait that became famous for a twofold reason: it captured in a physiological way the qualitative behavior of Hodgkin-Huxley model and it revealed the power of simple dynamical models to unfold complex firing patterns. To date, in spite of the enormous progresses in qualitative and quantitative neural modeling, this phase portrait has remained a core picture of neuronal excitability. Yet, a major difference between the neurophysiology of 1961 and of 2011 is the recognition of the prominent role of calcium channels in firing mechanisms. We show that including this extra current in Hodgkin-Huxley dynamics leads to a revision of FitzHugh-Nagumo phase portrait that affects in a fundamental way the reduced modeling of neural excitability. The revisited model considerably enlarges the modeling power of the original one. In particular, it captures essential electrophysiological signatures that otherwise require non-physiological alteration or considerable complexification of the classical model. As a basic illustration, the new model is shown to highlight a core dynamical mechanism by which calcium channels control the two distinct firing modes of thalamocortical neurons.

We have developed an implantable fuel cell that generates power through glucose oxidation, producing 3:4 mW cm{2steady-state power and up to 180 mW cm{2 peak power. The fuel cell is manufactured using a novel approach, employing semiconductor fabrication techniques, and is therefore well suited for manufacture together with integrated circuits on a single silicon wafer. Thus, it can help enable implantable microelectronic systems with long-lifetime power sources that harvest energy from their surrounds. The fuel reactions are mediated by robust, solid state catalysts. Glucose is oxidized at the nanostructured surface of an activated platinum anode. Oxygen is reduced to water at the surface of a self-assembled network of single-walled carbon nanotubes, embedded in a Nafion film that forms the cathode and is exposed to the biological environment. The catalytic electrodes are separated by a Nafion membrane. The availability of fuel cell reactants, oxygen and glucose, only as a mixture in the physiologic environment, has traditionally posed a design challenge: Net current production requires oxidation and reduction to occur separately and selectively at the anode and cathode, respectively, to prevent electrochemical short circuits. Our fuel cell is configured in a half-open geometry that shields the anode while exposing the cathode, resulting in an oxygen gradient that strongly favors oxygen reduction at the cathode. Glucose reaches the shielded anode by diffusing through the nanotube mesh, which does not catalyze glucose oxidation, and the Nafion layers, which are permeable to small neutral and cationic species. We demonstrate computationally that the natural recirculation of cerebrospinal fluid around the human brain theoretically permits glucose energy harvesting at a rate on the order of at least 1 mW with no adverse physiologic effects. Low-power brain–machine interfaces can thus potentially benefit from having their implanted units powered or recharged by glucose fuel cells., National Institutes of Health (U.S.) (grant NS-056140), United States. Office of Naval Research (grant N00014-09-1- 1015), United States. Defense Advanced Research Projects Agency (Research Assistantship to MIT Lincoln Laboratory)

Despite the known relevance of genomic structural variants to pathogen behavior, cancer, development, and evolution, certain repeat based structural variants may evade detection by existing high-throughput techniques. Here, we present ruler arrays, a technique to detect genomic structural variants including insertions and deletions (indels), duplications, and translocations. A ruler array exploits DNA polymerase’s processivity to detect physical distances between defined genomic sequences regardless of the intervening sequence. The method combines a sample preparation protocol, tiling genomic microarrays, and a new computational analysis. The analysis of ruler array data from two genomic samples enables the identification of structural variation between the samples. In an empirical test between two closely related haploid strains of yeast ruler arrays detected 78% of the structural variants larger than 100 bp., United States. National Institutes of Health (Grant R01GM069676)

ERK activation is enhanced by the scaffolding proteins KSR and MP1, localized near the cell membrane and late endosomes respectively, but little is known about their dynamic interplay. We develop here a mathematical model with ordinary differential equations to describe the dynamic activation of EGFR-ERK signaling under a conventional pathway without scaffolds, a KSR-scaffolded pathway, and an MP1-scaffolded pathway, and their impacts were examined under the influence of the endosomal regulators, Cbl-CIN85 and Endophilin A1. This new integrated model, validated against experimental results and computational constraints, shows that changes of ERK activation and EGFR endocytosis in response to EGF concentrations (i.e ligand sensitivity) depend on these scaffold proteins and regulators. The KSR-scaffolded and the conventional pathways act synergistically and are sensitive to EGF stimulation. When the KSR level is high, the sensitivity of ERK activation from this combined pathway remains low when Cbl-CIN85 level is low. But, such sensitivity can be increased with increasing levels of Endophilin if Cbl-CIN85 level becomes high. However, reduced KSR levels already present high sensitivity independent of Endophilin levels. In contrast, ERK activation by MP1 is additive to that of KSR but it shows little ligand-sensitivity under high levels of EGF. This can be partly reversed by increasing level of Endophilin while keeping Cbl-CIN85 level low. Further analyses showed that high levels of KSR affect ligand-sensitivity of EGFR endocytosis whereas MP1 ensures the robustness of endosomal ERK activation. These simulations constitute a multi-dimensional exploration of how EGF-dependent EGFR endocytosis and ERK activation are dynamically affected by scaffolds KSR and MP1, co-regulated by Cbl-CIN85 and Endophilin A1. Together, these results provide a detailed and quantitative demonstration of how regulators and scaffolds can collaborate to fine-tune the ligand-dependent sensitivity of EGFR endocytosis and ERK activation which could underlie differences during normal physiology, disease states and drug responses., Singapore-MIT Alliance for Research and Technology, National University of Singapore (Academic Research Fund (R-148-000-081-112/101)), Singapore. National Research Foundation, Singapore. Ministry of Education

Intermediate filaments (IFs) assembled in vitro from recombinantly expressed proteins have a diameter of 8–12 nm and can reach several micrometers in length. IFs assemble from a soluble pool of subunits, tetramers in the case of vimentin. Upon salt addition, the subunits form first unit length filaments (ULFs) within seconds and then assembly proceeds further by end-to-end fusion of ULFs and short filaments. So far, IF subunits have mainly been observed by electron microscopy of glycerol sprayed and rotary metal shadowed specimens. Due to the shear forces during spraying the IF subunits appear generally as straight thin rods. In this study, we used atomic force microscopy (AFM), cryo-electron microscopy (cryo-EM) combined with molecular modeling to investigate the conformation of the subunits of vimentin, desmin and keratin K5/K14 IFs in various conditions. Due to their anisotropic shape the subunits are difficult to image at high resolution by cryo-EM. In order to enhance contrast we used a cryo-negative staining approach. The subunits were clearly identified as thin, slightly curved rods. However the staining agent also forced the subunits to aggregate into two-dimensional networks of dot-like structures. To test this conformational change further, we imaged dried unfixed subunits on mica by AFM revealing a mixture of extended and dot-like conformations. The use of divalent ions such as calcium and magnesium, as well as glutaraldehyde exposure favored compact conformations over elongated ones. These experimental results as well as coarse-grained molecular dynamics simulations of a vimentin tetramer highlight the plasticity of IF subunits., Natural Sciences and Engineering Research Council of Canada, United States. Air Force Office of Scientific Research (Grant No. FA9550081-0321)

Background Generating autologous pluripotent stem cells for therapeutic applications will require the development of efficient DNA-free reprogramming techniques. Transfecting cells with in vitro-transcribed, protein-encoding RNA is a straightforward method of directly expressing high levels of reprogramming proteins without genetic modification. However, long-RNA transfection triggers a potent innate immune response characterized by growth inhibition and the production of inflammatory cytokines. As a result, repeated transfection with protein-encoding RNA causes cell death. Methodology/Principal Findings RNA viruses have evolved methods of disrupting innate immune signaling by destroying or inhibiting specific proteins to enable persistent infection. Starting from a list of known viral targets, we performed a combinatorial screen to identify siRNA cocktails that could desensitize cells to exogenous RNA. We show that combined knockdown of interferon-β (Ifnb1), Eif2ak2, and Stat2 rescues cells from the innate immune response triggered by frequent long-RNA transfection. Using this technique, we were able to transfect primary human fibroblasts every 24 hours with RNA encoding the reprogramming proteins Oct4, Sox2, Klf4, and Utf1. We provide evidence that the encoded protein is active, and we show that expression can be maintained for many days, through multiple rounds of cell division. Conclusions/Significance Our results demonstrate that suppressing innate immunity enables frequent transfection with protein-encoding RNA. This technique represents a versatile tool for investigating expression dynamics and protein interactions by enabling precise control over levels and timing of protein expression. Our finding also opens the door for the development of reprogramming and directed-differentiation methods based on long-RNA transfection., PLoS ONE, 5 (7), ISSN:1932-6203

Background Typically, algorithms to classify phenotypes using electronic medical record (EMR) data were developed to perform well in a specific patient population. There is increasing interest in analyses which can allow study of a specific outcome across different diseases. Such a study in the EMR would require an algorithm that can be applied across different patient populations. Our objectives were: (1) to develop an algorithm that would enable the study of coronary artery disease (CAD) across diverse patient populations; (2) to study the impact of adding narrative data extracted using natural language processing (NLP) in the algorithm. Additionally, we demonstrate how to implement CAD algorithm to compare risk across 3 chronic diseases in a preliminary study. Methods and Results We studied 3 established EMR based patient cohorts: diabetes mellitus (DM, n = 65,099), inflammatory bowel disease (IBD, n = 10,974), and rheumatoid arthritis (RA, n = 4,453) from two large academic centers. We developed a CAD algorithm using NLP in addition to structured data (e.g. ICD9 codes) in the RA cohort and validated it in the DM and IBD cohorts. The CAD algorithm using NLP in addition to structured data achieved specificity >95% with a positive predictive value (PPV) 90% in the training (RA) and validation sets (IBD and DM). The addition of NLP data improved the sensitivity for all cohorts, classifying an additional 17% of CAD subjects in IBD and 10% in DM while maintaining PPV of 90%. The algorithm classified 16,488 DM (26.1%), 457 IBD (4.2%), and 245 RA (5.0%) with CAD. In a cross-sectional analysis, CAD risk was 63% lower in RA and 68% lower in IBD compared to DM (p, National Institutes of Health (U.S.). Informatics for Integrating Biology and the Bedside Project (U54LM008748)

Intracellular delivery of biomolecules, such as proteins and siRNAs, into primary immune cells, especially resting lymphocytes, is a challenge. Here we describe the design and testing of microfluidic intracellular delivery systems that cause temporary membrane disruption by rapid mechanical deformation of human and mouse immune cells. Dextran, antibody and siRNA delivery performance is measured in multiple immune cell types and the approach’s potential to engineer cell function is demonstrated in HIV infection studies., National Institutes of Health (U.S.) (Grant R01GM101420-01A1), Kathy and Curt Marble Cancer Research Fund, Ragon Institute of MGH, MIT and Harvard, National Cancer Institute (U.S.) (Cancer Center Support (Core) Grant P30-CA14051)

Despite many compelling applications in economics, sociobiology, and evolutionary psychology, group selection is still one of the most hotly contested ideas in evolutionary biology. Here we propose a simple evolutionary model of behavior and show that what appears to be group selection may, in fact, simply be the consequence of natural selection occurring in stochastic environments with reproductive risks that are correlated across individuals. Those individuals with highly correlated risks will appear to form “groups”, even if their actions are, in fact, totally autonomous, mindless, and, prior to selection, uniformly randomly distributed in the population. This framework implies that a separate theory of group selection is not strictly necessary to explain observed phenomena such as altruism and cooperation. At the same time, it shows that the notion of group selection does captures a unique aspect of evolution—selection with correlated reproductive risk–that may be sufficiently widespread to warrant a separate term for the phenomenon., MIT Laboratory for Financial Engineering

Purpose To determine the reproducibility of automated segmentation of the three-dimensional (3D) lamina cribrosa (LC) microarchitecture scanned in-vivo using optical coherence tomography (OCT). Methods Thirty-nine eyes (8 healthy, 19 glaucoma suspects and 12 glaucoma) from 49 subjects were scanned twice using swept-source (SS−) OCT in a 3.5×3.5×3.64 mm (400×400×896 pixels) volume centered on the optic nerve head, with the focus readjusted after each scan. The LC was automatically segmented and analyzed for microarchitectural parameters, including pore diameter, pore diameter standard deviation (SD), pore aspect ratio, pore area, beam thickness, beam thickness SD, and beam thickness to pore diameter ratio. Reproducibility of the parameters was assessed by computing the imprecision of the parameters between the scans. Results The automated segmentation demonstrated excellent reproducibility. All LC microarchitecture parameters had an imprecision of less or equal to 4.2%. There was little variability in imprecision with respect to diagnostic category, although the method tends to show higher imprecision amongst healthy subjects. Conclusion The proposed automated segmentation of the LC demonstrated high reproducibility for 3D LC parameters. This segmentation analysis tool will be useful for in-vivo studies of the LC., Eye & Ear Foundation of Pittsburgh, Research to Prevent Blindness, Inc. (United States), National Institutes of Health (U.S.) (contract R01-EY013178), National Institutes of Health (U.S.) (contract R01-EY011289), National Institutes of Health (U.S.) (contract P30-EY008098)

We demonstrate in vivo choriocapillaris and choroidal microvasculature imaging in normal human subjects using optical coherence tomography (OCT). An ultrahigh speed swept source OCT prototype at 1060 nm wavelengths with a 400 kHz A-scan rate is developed for three-dimensional ultrahigh speed imaging of the posterior eye. OCT angiography is used to image three-dimensional vascular structure without the need for exogenous fluorophores by detecting erythrocyte motion contrast between OCT intensity cross-sectional images acquired rapidly and repeatedly from the same location on the retina. En face OCT angiograms of the choriocapillaris and choroidal vasculature are visualized by acquiring cross-sectional OCT angiograms volumetrically via raster scanning and segmenting the three-dimensional angiographic data at multiple depths below the retinal pigment epithelium (RPE). Fine microvasculature of the choriocapillaris, as well as tightly packed networks of feeding arterioles and draining venules, can be visualized at different en face depths. Panoramic ultra-wide field stitched OCT angiograms of the choriocapillaris spanning ~32 mm on the retina show distinct vascular structures at different fundus locations. Isolated smaller fields at the central fovea and ~6 mm nasal to the fovea at the depths of the choriocapillaris and Sattler's layer show vasculature structures consistent with established architectural morphology from histological and electron micrograph corrosion casting studies. Choriocapillaris imaging was performed in eight healthy volunteers with OCT angiograms successfully acquired from all subjects. These results demonstrate the feasibility of ultrahigh speed OCT for in vivo dye-free choriocapillaris and choroidal vasculature imaging, in addition to conventional structural imaging., National Institutes of Health (U.S.) (NIH R01-EY011289-27), National Institutes of Health (U.S.) (NIH R01-EY013178-12), National Institutes of Health (U.S.) (NIH R44-EY022864-01), National Institutes of Health (U.S.) (NIH R01-CA075289-16), United States. Air Force Office of Scientific Research (AFOSR FA9550-10-1-0551), United States. Air Force Office of Scientific Research (AFOSR FA9550-12-1-0499)

Real-time brain-machine interfaces (BMI) have focused on either estimating the continuous movement trajectory or target intent. However, natural movement often incorporates both. Additionally, BMIs can be modeled as a feedback control system in which the subject modulates the neural activity to move the prosthetic device towards a desired target while receiving real-time sensory feedback of the state of the movement. We develop a novel real-time BMI using an optimal feedback control design that jointly estimates the movement target and trajectory of monkeys in two stages. First, the target is decoded from neural spiking activity before movement initiation. Second, the trajectory is decoded by combining the decoded target with the peri-movement spiking activity using an optimal feedback control design. This design exploits a recursive Bayesian decoder that uses an optimal feedback control model of the sensorimotor system to take into account the intended target location and the sensory feedback in its trajectory estimation from spiking activity. The real-time BMI processes the spiking activity directly using point process modeling. We implement the BMI in experiments consisting of an instructed-delay center-out task in which monkeys are presented with a target location on the screen during a delay period and then have to move a cursor to it without touching the incorrect targets. We show that the two-stage BMI performs more accurately than either stage alone. Correct target prediction can compensate for inaccurate trajectory estimation and vice versa. The optimal feedback control design also results in trajectories that are smoother and have lower estimation error. The two-stage decoder also performs better than linear regression approaches in offline cross-validation analyses. Our results demonstrate the advantage of a BMI design that jointly estimates the target and trajectory of movement and more closely mimics the sensorimotor control system., National Institutes of Health (U.S.) (NIH grant No.DP1-0D003646-01), National Institutes of Health (U.S.) (NIH grant R01-EB006385)

Recent analysis of genomic signatures in mammals, flies, and worms indicates that functional translational stop codon readthrough is considerably more abundant in metazoa than previously recognized, but this analysis provides only limited clues about the function or mechanism of readthrough. If an mRNA known to be read through in one species is also read through in another, perhaps these questions can be studied in a simpler setting. With this end in mind, we have investigated whether some of the readthrough genes in human, fly, and worm also exhibit readthrough when expressed in S. cerevisiae. We found that readthrough was highest in a gene with a post-stop hexamer known to trigger readthrough, while other metazoan readthrough genes exhibit borderline readthrough in S. cerevisiae., National Institutes of Health (U.S.) (5U54HG004555-03)

The ability to decode neural activity into meaningful control signals for prosthetic devices is critical to the development of clinically useful brain– machine interfaces (BMIs). Such systems require input from tens to hundreds of brain-implanted recording electrodes in order to deliver robust and accurate performance; in serving that primary function they should also minimize power dissipation in order to avoid damaging neural tissue; and they should transmit data wirelessly in order to minimize the risk of infection associated with chronic, transcutaneous implants. Electronic architectures for brain– machine interfaces must therefore minimize size and power consumption, while maximizing the ability to compress data to be transmitted over limited-bandwidth wireless channels. Here we present a system of extremely low computational complexity, designed for real-time decoding of neural signals, and suited for highly scalable implantable systems. Our programmable architecture is an explicit implementation of a universal computing machine emulating the dynamics of a network of integrate-and-fire neurons; it requires no arithmetic operations except for counting, and decodes neural signals using only computationally inexpensive logic operations. The simplicity of this architecture does not compromise its ability to compress raw neural data by factors greater than . We describe a set of decoding algorithms based on this computational architecture, one designed to operate within an implanted system, minimizing its power consumption and data transmission bandwidth; and a complementary set of algorithms for learning, programming the decoder, and postprocessing the decoded output, designed to operate in an external, nonimplanted unit. The implementation of the implantable portion is estimated to require fewer than 5000 operations per second. A proof-of-concept, 32-channel field-programmable gate array (FPGA) implementation of this portion is consequently energy efficient. We validate the performance of our overall system by decoding electrophysiologic data from a behaving rodent., United States. National Institutes of Health (Grant NS056140)

It is well known that damage to the peripheral auditory system causes deficits in tone detection as well as pitch and loudness perception across a wide range of frequencies. However, the extent to which to which the auditory cortex plays a critical role in these basic aspects of spectral processing, especially with regard to speech, music, and environmental sound perception, remains unclear. Recent experiments indicate that primary auditory cortex is necessary for the normally-high perceptual acuity exhibited by humans in pure-tone frequency discrimination. The present study assessed whether the auditory cortex plays a similar role in the intensity domain and contrasted its contribution to sensory versus discriminative aspects of intensity processing. We measured intensity thresholds for pure-tone detection and pure-tone loudness discrimination in a population of healthy adults and a middle-aged man with complete or near-complete lesions of the auditory cortex bilaterally. Detection thresholds in his left and right ears were 16 and 7 dB HL, respectively, within clinically-defined normal limits. In contrast, the intensity threshold for monaural loudness discrimination at 1 kHz was 6.5±2.1 dB in the left ear and 6.5±1.9 dB in the right ear at 40 dB sensation level, well above the means of the control population (left ear: 1.6±0.22 dB; right ear: 1.7±0.19 dB). The results indicate that auditory cortex lowers just-noticeable differences for loudness discrimination by approximately 5 dB but is not necessary for tone detection in quiet. Previous human and Old-world monkey experiments employing lesion-effect, neurophysiology, and neuroimaging methods to investigate the role of auditory cortex in intensity processing are reviewed., United States. National Institutes of Health (DC03328), United States. National Institutes of Health (DC006353), United States. National Institutes of Health (DC00117), United States. National Institutes of Health (T32-DC00038)

Background Autocrine & paracrine signaling are widespread both in vivo and in vitro, and are particularly important in embryonic stem cell (ESC) pluripotency and lineage commitment. Although autocrine signaling via fibroblast growth factor-4 (FGF4) is known to be required in mouse ESC (mESC) neuroectodermal specification, the question of whether FGF4 autocrine signaling is sufficient, or whether other soluble ligands are also involved in fate specification, is unknown. The spatially confined and closed-loop nature of diffusible signaling makes its experimental control challenging; current experimental approaches typically require prior knowledge of the factor/receptor in order to modulate the loop. A new approach explored in this work is to leverage transport phenomena at cellular resolution to downregulate overall diffusible signaling through the physical removal of cell-secreted ligands. Methodology/Principal Findings We develop a multiplex microfluidic platform to continuously remove cell-secreted (autocrine\paracrine) factors to downregulate diffusible signaling. By comparing cell growth and differentiation in side-by-side chambers with or without added cell-secreted factors, we isolate the effects of diffusible signaling from artifacts such as shear, nutrient depletion, and microsystem effects, and find that cell-secreted growth factor(s) are required during neuroectodermal specification. Then we induce FGF4 signaling in minimal chemically defined medium (N2B27) and inhibit FGF signaling in fully supplemented differentiation medium with cell-secreted factors to determine that the non-FGF cell-secreted factors are required to promote growth of differentiating mESCs. Conclusions/Significance Our results demonstrate for the first time that flow can downregulate autocrine\paracrine signaling and examine sufficiency of extracellular factors. We show that autocrine\paracrine signaling drives neuroectodermal commitment of mESCs through both FGF4-dependent and -independent pathways. Overall, by uncovering autocrine\paracrine processes previously hidden in conventional culture systems, our results establish microfluidic perfusion as a technique to study and manipulate diffusible signaling in cell systems., National Institutes of Health (U.S.) (NIH grant No. EB007278), Swiss National Science Foundation