Joint Modeling and Registration of Cell Populations in Cohorts of High-Dimensional Flow Cytometric Data

In biomedical applications, an experimenter encounters different potential sources of variation in data such as individual samples, multiple experimental conditions, and multivariate responses of a panel of markers such as from a signaling network. In multiparametric cytometry, which is often used for analyzing patient samples, such issues are critical. While computational methods can identify cell populations in individual samples, without the ability to automatically match them across samples, it is difficult to compare and characterize the populations in typical experiments, such as those responding to various stimulations or distinctive of particular patients or time-points, especially when there are many samples. Joint Clustering and Matching (JCM) is a multi-level framework for simultaneous modeling and registration of populations across a cohort. JCM models every population with a robust multivariate probability distribution. Simultaneously, JCM fits a random-effects model to construct an overall batch template – used for registering populations across samples, and classifying new samples. By tackling systems-level variation, JCM supports practical biomedical applications involving large cohorts. Software for fitting the JCM models have been implemented in an R package EMMIX-JCM, available from http://www.maths.uq.edu.au/~gjm/mix_soft/EMMIX-JCM/.     

Transcervical Ultrasonography Is Feasible to Visualize and Evaluate Base of Tongue Cancers

Background

Base of tongue (BOT) is a difficult subsite to examine clinically and radiographically. Yet, anatomic delineation of the primary tumor site, its extension to adjacent sites or across midline, and endophytic vs. exophytic extent are important characteristics for staging and treatment planning. We hypothesized that ultrasound could be used to visualize and describe BOT tumors.

Methods

Transcervical ultrasound was performed using a standardized protocol in cases and controls. Cases had suspected or confirmed BOT malignancy. Controls were healthy individuals without known malignancy.

Results

100% of BOT tumors were visualized. On ultrasound BOT tumors were hypoechoic (90.9%) with irregular margins (95.5%). Ultrasound could be used to characterize adjacent site involvement, midline extent, and endophytic extent, and visualize the lingual artery. No tumors were suspected for controls.

Conclusions

Ultrasonography can be used to transcervically visualize BOT tumors and provides clinically relevant characteristics that may not otherwise be appreciable.     

Population Connectivity Shifts at High Frequency within an Open-Coast Marine Protected Area Network

A complete understanding of population connectivity via larval dispersal is of great value to the effective design and management of marine protected areas (MPA). However empirical estimates of larval dispersal distance, self-recruitment, and within season variability of population connectivity patterns and their influence on metapopulation structure remain rare. We used high-resolution otolith microchemistry data from the temperate reef fish Hypsypops rubicundus to explore biweekly, seasonal, and annual connectivity patterns in an open-coast MPA network. The three MPAs, spanning 46 km along the southern California coastline were connected by larval dispersal, but the magnitude and direction of connections reversed between 2008 and 2009. Self-recruitment, i.e. spawning, dispersal, and settlement to the same location, was observed at two locations, one of which is a MPA. Self-recruitment to this MPA ranged from 50–84%; within the entire 60 km study region, self-recruitment accounted for 45% of all individuals settling to study reefs. On biweekly time scales we observed directional variability in alongshore current data and larval dispersal trajectories; if viewed in isolation these data suggest the system behaves as a source-sink metapopulation. However aggregate biweekly data over two years reveal a reef network in which H. rubicundus behaves more like a well-mixed metapopulation. As one of the few empirical studies of population connectivity within a temperate open coast reef network, this work can inform the MPA design process, implementation of ecosystem based management plans, and facilitate conservation decisions.     

Finding Near-Optimal Groups of Epidemic Spreaders in a Complex Network

In this paper, we present algorithms to find near-optimal sets of epidemic spreaders in complex networks. We extend the notion of local-centrality, a centrality measure previously shown to correspond with a node''s ability to spread an epidemic, to sets of nodes by introducing combinatorial local centrality. Though we prove that finding a set of nodes that maximizes this new measure is NP-hard, good approximations are available. We show that a strictly greedy approach obtains the best approximation ratio unless P = NP and then formulate a modified version of this approach that leverages qualities of the network to achieve a faster runtime while maintaining this theoretical guarantee. We perform an experimental evaluation on samples from several different network structures which demonstrate that our algorithm maximizes combinatorial local centrality and consistently chooses the most effective set of nodes to spread infection under the SIR model, relative to selecting the top nodes using many common centrality measures. We also demonstrate that the optimized algorithm we develop scales effectively.     

Missense Mutation Lys18Asn in Dystrophin that Triggers X-Linked Dilated Cardiomyopathy Decreases Protein Stability,Increases Protein Unfolding,and Perturbs Protein Structure,but Does Not Affect Protein Function

Genetic mutations in a vital muscle protein dystrophin trigger X-linked dilated cardiomyopathy (XLDCM). However, disease mechanisms at the fundamental protein level are not understood. Such molecular knowledge is essential for developing therapies for XLDCM. Our main objective is to understand the effect of disease-causing mutations on the structure and function of dystrophin. This study is on a missense mutation K18N. The K18N mutation occurs in the N-terminal actin binding domain (N-ABD). We created and expressed the wild-type (WT) N-ABD and its K18N mutant, and purified to homogeneity. Reversible folding experiments demonstrated that both mutant and WT did not aggregate upon refolding. Mutation did not affect the protein''s overall secondary structure, as indicated by no changes in circular dichroism of the protein. However, the mutant is thermodynamically less stable than the WT (denaturant melts), and unfolds faster than the WT (stopped-flow kinetics). Despite having global secondary structure similar to that of the WT, mutant showed significant local structural changes at many amino acids when compared with the WT (heteronuclear NMR experiments). These structural changes indicate that the effect of mutation is propagated over long distances in the protein structure. Contrary to these structural and stability changes, the mutant had no significant effect on the actin-binding function as evident from co-sedimentation and depolymerization assays. These results summarize that the K18N mutation decreases thermodynamic stability, accelerates unfolding, perturbs protein structure, but does not affect the function. Therefore, K18N is a stability defect rather than a functional defect. Decrease in stability and increase in unfolding decrease the net population of dystrophin molecules available for function, which might trigger XLDCM. Consistently, XLDCM patients have decreased levels of dystrophin in cardiac muscle.     

Continuous Stroke Volume Estimation from Aortic Pressure Using Zero Dimensional Cardiovascular Model: Proof of Concept Study from Porcine Experiments

Introduction

Accurate, continuous, left ventricular stroke volume (SV) measurements can convey large amounts of information about patient hemodynamic status and response to therapy. However, direct measurements are highly invasive in clinical practice, and current procedures for estimating SV require specialized devices and significant approximation.

Method

This study investigates the accuracy of a three element Windkessel model combined with an aortic pressure waveform to estimate SV. Aortic pressure is separated into two components capturing; 1) resistance and compliance, 2) characteristic impedance. This separation provides model-element relationships enabling SV to be estimated while requiring only one of the three element values to be known or estimated. Beat-to-beat SV estimation was performed using population-representative optimal values for each model element. This method was validated using measured SV data from porcine experiments (N = 3 female Pietrain pigs, 29–37 kg) in which both ventricular volume and aortic pressure waveforms were measured simultaneously.

Results

The median difference between measured SV from left ventricle (LV) output and estimated SV was 0.6 ml with a 90% range (5^th–95^th percentile) −12.4 ml–14.3 ml. During periods when changes in SV were induced, cross correlations in between estimated and measured SV were above R = 0.65 for all cases.

Conclusion

The method presented demonstrates that the magnitude and trends of SV can be accurately estimated from pressure waveforms alone, without the need for identification of complex physiological metrics where strength of correlations may vary significantly from patient to patient.     

Unique Antibody Responses to Malondialdehyde-Acetaldehyde (MAA)-Protein Adducts Predict Coronary Artery Disease

Malondialdehyde-acetaldehyde adducts (MAA) have been implicated in atherosclerosis. The purpose of this study was to investigate the role of MAA in atherosclerotic disease. Serum samples from controls (n = 82) and patients with; non-obstructive coronary artery disease (CAD), (n = 40), acute myocardial infarction (AMI) (n = 42), or coronary artery bypass graft (CABG) surgery due to obstructive multi-vessel CAD (n = 72), were collected and tested for antibody isotypes to MAA-modifed human serum albumin (MAA-HSA). CAD patients had elevated relative levels of IgG and IgA anti-MAA, compared to control patients (p<0.001). AMI patients had a significantly increased relative levels of circulating IgG anti-MAA-HSA antibodies as compared to stable angina (p<0.03) or CABG patients (p<0.003). CABG patients had significantly increased relative levels of circulating IgA anti-MAA-HSA antibodies as compared to non-obstructive CAD (p<0.001) and AMI patients (p<0.001). Additionally, MAA-modified proteins were detected in the tissue of human AMI lesions. In conclusion, the IgM, IgG and IgA anti-MAA-HSA antibody isotypes are differentially and significantly associated with non-obstructive CAD, AMI, or obstructive multi-vessel CAD and may serve as biomarkers of atherosclerotic disease.     

The Microtubule Plus-End Tracking Protein ARMADILLO-REPEAT KINESIN1 Promotes Microtubule Catastrophe in Arabidopsis

Microtubule dynamics are critically important for plant cell development. Here, we show that Arabidopsis thaliana ARMADILLO-REPEAT KINESIN1 (ARK1) plays a key role in root hair tip growth by promoting microtubule catastrophe events. This destabilizing activity appears to maintain adequate free tubulin concentrations in order to permit rapid microtubule growth, which in turn is correlated with uniform tip growth. Microtubules in ark1-1 root hairs exhibited reduced catastrophe frequency and slower growth velocities, both of which were restored by low concentrations of the microtubule-destabilizing drug oryzalin. An ARK1-GFP (green fluorescent protein) fusion protein expressed under its endogenous promoter localized to growing microtubule plus ends and rescued the ark1-1 root hair phenotype. Transient overexpression of ARK1-RFP (red fluorescent protein) increased microtubule catastrophe frequency. ARK1-fusion protein constructs lacking the N-terminal motor domain still labeled microtubules, suggesting the existence of a second microtubule binding domain at the C terminus of ARK1. ARK1-GFP was broadly expressed in seedlings, but mutant phenotypes were restricted to root hairs, indicating that ARK1’s function is redundant in cells other than those forming root hairs.