Unsupervised Bayesian linear unmixing of gene expression microarrays

Background: This paper introduces a new constrained model and the corresponding algorithm, called unsupervised Bayesian linear unmixing (uBLU), to identify biological signatures from high dimensional assays like gene expression microarrays. The basis for uBLU is a Bayesian model for the data samples which are represented as an additive mixture of random positive gene signatures, called factors, with random positive mixing coefficients, called factor scores, that specify the relative contribution of each signature to a specific sample. The particularity of the proposed method is that uBLU constrains the factor loadings to be non-negative and the factor scores to be probability distributions over the factors. Furthermore, it also provides estimates of the number of factors. A Gibbs sampling strategy is adopted here to generate random samples according to the posterior distribution of the factors, factor scores, and number of factors. These samples are then used to estimate all the unknown parameters. Results: Firstly, the proposed uBLU method is applied to several simulated datasets with known ground truth and compared with previous factor decomposition methods, such as principal component analysis (PCA), non negative matrix factorization (NMF), Bayesian factor regression modeling (BFRM), and the gradient-based algorithm for general matrix factorization (GB-GMF). Secondly, we illustrate the application of uBLU on a real time-evolving gene expression dataset from a recent viral challenge study in which individuals have been inoculated with influenza A/H3N2/Wisconsin. We show that the uBLU method significantly outperforms the other methods on the simulated and real data sets considered here. Conclusions: The results obtained on synthetic and real data illustrate the accuracy of the proposed uBLU method when compared to other factor decomposition methods from the literature (PCA, NMF, BFRM, and GB-GMF). The uBLU method identifies an inflammatory component closely associated with clinical symptom scores collected during the study. Using a constrained model allows recovery of all the inflammatory genes in a single factor

Immunosuppressive niche engineering at the onset of human colorectal cancer

The evolutionary dynamics of tumor initiation remain undetermined, and the interplay between neoplastic cells and the immune system is hypothesized to be critical in transformation. Colorectal cancer (CRC) presents a unique opportunity to study the transition to malignancy as pre-cancers (adenomas) and early-stage cancers are frequently resected. Here, we examine tumor-immune eco-evolutionary dynamics from pre-cancer to carcinoma using a computational model, ecological analysis of digital pathology data, and neoantigen prediction in 62 patient samples. Modeling predicted recruitment of immunosuppressive cells would be the most common driver of transformation. As predicted, ecological analysis reveals that progressed adenomas co-localized with immunosuppressive cells and cytokines, while benign adenomas co-localized with a mixed immune response. Carcinomas converge to a common immune “cold” ecology, relaxing selection against immunogenicity and high neoantigen burdens, with little evidence for PD-L1 overexpression driving tumor initiation. These findings suggest re-engineering the immunosuppressive niche may prove an effective immunotherapy in CRC

Current controversies in TNM for the radiological staging of rectal cancer and how to deal with them: results of a global online survey and multidisciplinary expert consensus

Objectives: To identify the main problem areas in the applicability of the current TNM staging system (8th ed.) for the radiological staging and reporting of rectal cancer and provide practice recommendations on how to handle them. Methods: A global case-based online survey was conducted including 41 image-based rectal cancer cases focusing on various items included in the TNM system. Cases reaching < 80% agreement among survey respondents were identified as problem areas and discussed among an international expert panel, including 5 radiologists, 6 colorectal surgeons, 4 radiation oncologists, and 3 pathologists. Results: Three hundred twenty-one respondents (from 32 countries) completed the survey. Sixteen problem areas were identified, related to cT staging in low-rectal cancers, definitions for cT4b and cM1a disease, definitions for mesorectal fascia (MRF) involvement, evaluation of lymph nodes versus tumor deposits, and staging of lateral lymph nodes. The expert panel recommended strategies on how to handle these, including advice on cT-stage categorization in case of involvement of different layers of the anal canal, specifications on which structures to include in the definition of cT4b disease, how to define MRF involvement by the primary tumor and other tumor-bearing structures, how to differentiate and report lymph nodes and tumor deposits on MRI, and how to anatomically localize and stage lateral lymph nodes. Conclusions: The recommendations derived from this global survey and expert panel discussion may serve as a practice guide and support tool for radiologists (and other clinicians) involved in the staging of rectal cancer and may contribute to improved consistency in radiological staging and reporting. Key Points: • Via a case-based online survey (incl. 321 respondents from 32 countries), we identified 16 problem areas related to the applicability of the TNM staging system for the radiological staging and reporting of rectal cancer. • A multidisciplinary panel of experts recommended strategies on how to handle these problem areas, including advice on cT-stage categorization in case of involvement of different layers of the anal canal, specifications on which structures to include in the definition of cT4b disease, how to define mesorectal fascia involvement by the primary tumor and other tumor-bearing structures, how to differentiate and report lymph nodes and tumor deposits on MRI, and how to anatomically localize and stage lateral lymph nodes. • These recommendations may serve as a practice guide and support tool for radiologists (and other clinicians) involved in the staging of rectal cancer and may contribute to improved consistency in radiological staging and reporting

Characterization of complex networks: A survey of measurements

Each complex network (or class of networks) presents specific topological features which characterize its connectivity and highly influence the dynamics of processes executed on the network. The analysis, discrimination, and synthesis of complex networks therefore rely on the use of measurements capable of expressing the most relevant topological features. This article presents a survey of such measurements. It includes general considerations about complex network characterization, a brief review of the principal models, and the presentation of the main existing measurements. Important related issues covered in this work comprise the representation of the evolution of complex networks in terms of trajectories in several measurement spaces, the analysis of the correlations between some of the most traditional measurements, perturbation analysis, as well as the use of multivariate statistics for feature selection and network classification. Depending on the network and the analysis task one has in mind, a specific set of features may be chosen. It is hoped that the present survey will help the proper application and interpretation of measurements.Comment: A working manuscript with 78 pages, 32 figures. Suggestions of measurements for inclusion are welcomed by the author

Integrin β3 Crosstalk with VEGFR Accommodating Tyrosine Phosphorylation as a Regulatory Switch

Integrins mediate cell adhesion, migration, and survival by connecting intracellular machinery with the surrounding extracellular matrix. Previous studies demonstrated the importance of the interaction between β3 integrin and VEGF type 2 receptor (VEGFR2) in VEGF-induced angiogenesis. Here we present in vitro evidence of the direct association between the cytoplasmic tails (CTs) of β3 and VEGFR2. Specifically, the membrane-proximal motif around 801YLSI in VEGFR2 mediates its binding to non-phosphorylated β3CT, accommodating an α-helical turn in integrin bound conformation. We also show that Y747 phosphorylation of β3 enhances the above interaction. To demonstrate the importance of β3 phosphorylation in endothelial cell functions, we synthesized β3CT-mimicking Y747 phosphorylated and unphosphorylated membrane permeable peptides. We show that a peptide containing phospho-Y747 but not F747 significantly inhibits VEGF-induced signaling and angiogenesis. Moreover, phospho-Y747 peptide exhibits inhibitory effect only in WT but not in β3 integrin knock-out or β3 integrin knock-in cells expressing β3 with two tyrosines substituted for phenylalanines, demonstrating its specificity. Importantly, these peptides have no effect on fibroblast growth factor receptor signaling. Collectively these data provide novel mechanistic insights into phosphorylation dependent cross-talk between integrin and VEGFR2

Biomarker discovery in heterogeneous tissue samples -taking the in-silico deconfounding approach

<p>Abstract</p> <p>Background</p> <p>For heterogeneous tissues, such as blood, measurements of gene expression are confounded by relative proportions of cell types involved. Conclusions have to rely on estimation of gene expression signals for homogeneous cell populations, e.g. by applying micro-dissection, fluorescence activated cell sorting, or <it>in-silico </it>deconfounding. We studied feasibility and validity of a non-negative matrix decomposition algorithm using experimental gene expression data for blood and sorted cells from the same donor samples. Our objective was to optimize the algorithm regarding detection of differentially expressed genes and to enable its use for classification in the difficult scenario of reversely regulated genes. This would be of importance for the identification of candidate biomarkers in heterogeneous tissues.</p> <p>Results</p> <p>Experimental data and simulation studies involving noise parameters estimated from these data revealed that for valid detection of differential gene expression, quantile normalization and use of non-log data are optimal. We demonstrate the feasibility of predicting proportions of constituting cell types from gene expression data of single samples, as a prerequisite for a deconfounding-based classification approach.</p> <p>Classification cross-validation errors with and without using deconfounding results are reported as well as sample-size dependencies. Implementation of the algorithm, simulation and analysis scripts are available.</p> <p>Conclusions</p> <p>The deconfounding algorithm without decorrelation using quantile normalization on non-log data is proposed for biomarkers that are difficult to detect, and for cases where confounding by varying proportions of cell types is the suspected reason. In this case, a deconfounding ranking approach can be used as a powerful alternative to, or complement of, other statistical learning approaches to define candidate biomarkers for molecular diagnosis and prediction in biomedicine, in realistically noisy conditions and with moderate sample sizes.</p

Performance of CMS muon reconstruction in pp collision events at sqrt(s) = 7 TeV

The performance of muon reconstruction, identification, and triggering in CMS has been studied using 40 inverse picobarns of data collected in pp collisions at sqrt(s) = 7 TeV at the LHC in 2010. A few benchmark sets of selection criteria covering a wide range of physics analysis needs have been examined. For all considered selections, the efficiency to reconstruct and identify a muon with a transverse momentum pT larger than a few GeV is above 95% over the whole region of pseudorapidity covered by the CMS muon system, abs(eta) < 2.4, while the probability to misidentify a hadron as a muon is well below 1%. The efficiency to trigger on single muons with pT above a few GeV is higher than 90% over the full eta range, and typically substantially better. The overall momentum scale is measured to a precision of 0.2% with muons from Z decays. The transverse momentum resolution varies from 1% to 6% depending on pseudorapidity for muons with pT below 100 GeV and, using cosmic rays, it is shown to be better than 10% in the central region up to pT = 1 TeV. Observed distributions of all quantities are well reproduced by the Monte Carlo simulation.Comment: Replaced with published version. Added journal reference and DO