Discovery of complex oxides via automated experiments and data science.

The quest to identify materials with tailored properties is increasingly expanding into high-order composition spaces, with a corresponding combinatorial explosion in the number of candidate materials. A key challenge is to discover regions in composition space where materials have novel properties. Traditional predictive models for material properties are not accurate enough to guide the search. Herein, we use high-throughput measurements of optical properties to identify novel regions in three-cation metal oxide composition spaces by identifying compositions whose optical trends cannot be explained by simple phase mixtures. We screen 376,752 distinct compositions from 108 three-cation oxide systems based on the cation elements Mg, Fe, Co, Ni, Cu, Y, In, Sn, Ce, and Ta. Data models for candidate phase diagrams and three-cation compositions with emergent optical properties guide the discovery of materials with complex phase-dependent properties, as demonstrated by the discovery of a Co-Ta-Sn substitutional alloy oxide with tunable transparency, catalytic activity, and stability in strong acid electrolytes. These results required close coupling of data validation to experiment design to generate a reliable end-to-end high-throughput workflow for accelerating scientific discovery.

Learning causal networks using inducible transcription factors and transcriptome-wide time series.

We present IDEA (the Induction Dynamics gene Expression Atlas), a dataset constructed by independently inducing hundreds of transcription factors (TFs) and measuring timecourses of the resulting gene expression responses in budding yeast. Each experiment captures a regulatory cascade connecting a single induced regulator to the genes it causally regulates. We discuss the regulatory cascade of a single TF, Aft1, in detail; however, IDEA contains > 200 TF induction experiments with 20 million individual observations and 100,000 signal-containing dynamic responses. As an application of IDEA, we integrate all timecourses into a whole-cell transcriptional model, which is used to predict and validate multiple new and underappreciated transcriptional regulators. We also find that the magnitudes of coefficients in this model are predictive of genetic interaction profile similarities. In addition to being a resource for exploring regulatory connectivity between TFs and their target genes, our modeling approach shows that combining rapid perturbations of individual genes with genome-scale time-series measurements is an effective strategy for elucidating gene regulatory networks.

Leveraging Multi-ethnic Evidence for Risk Assessment of Quantitative Traits in Minority Populations.

An essential component of precision medicine is the ability to predict an individual's risk of disease based on genetic and non-genetic factors. For complex traits and diseases, assessing the risk due to genetic factors is challenging because it requires knowledge of both the identity of variants that influence the trait and their corresponding allelic effects. Although the set of risk variants and their allelic effects may vary between populations, a large proportion of these variants were identified based on studies in populations of European descent. Heterogeneity in genetic architecture underlying complex traits and diseases, while broadly acknowledged, remains poorly characterized. Ignoring such heterogeneity likely reduces predictive accuracy for minority individuals. In this study, we propose an approach, called XP-BLUP, which ameliorates this ethnic disparity by combining trans-ethnic and ethnic-specific information. We build a polygenic model for complex traits that distinguishes candidate trait-relevant variants from the rest of the genome. The set of candidate variants are selected based on studies in any human population, yet the allelic effects are evaluated in a population-specific fashion. Simulation studies and real data analyses demonstrate that XP-BLUP adaptively utilizes trans-ethnic information and can substantially improve predictive accuracy in minority populations. At the same time, our study highlights the importance of the continued expansion of minority cohorts.

Leveraging Multi-ethnic Evidence for Mapping Complex Traits in Minority Populations: An Empirical Bayes Approach.

Elucidating the genetic basis of complex traits and diseases in non-European populations is particularly challenging because US minority populations have been under-represented in genetic association studies. We developed an empirical Bayes approach named XPEB (cross-population empirical Bayes), designed to improve the power for mapping complex-trait-associated loci in a minority population by exploiting information from genome-wide association studies (GWASs) from another ethnic population. Taking as input summary statistics from two GWASs-a target GWAS from an ethnic minority population of primary interest and an auxiliary base GWAS (such as a larger GWAS in Europeans)-our XPEB approach reprioritizes SNPs in the target population to compute local false-discovery rates. We demonstrated, through simulations, that whenever the base GWAS harbors relevant information, XPEB gains efficiency. Moreover, XPEB has the ability to discard irrelevant auxiliary information, providing a safeguard against inflated false-discovery rates due to genetic heterogeneity between populations. Applied to a blood-lipids study in African Americans, XPEB more than quadrupled the discoveries from the conventional approach, which used a target GWAS alone, bringing the number of significant loci from 14 to 65. Thus, XPEB offers a flexible framework for mapping complex traits in minority populations.

Genome-wide characterization of shared and distinct genetic components that influence blood lipid levels in ethnically diverse human populations.

Blood lipid concentrations are heritable risk factors associated with atherosclerosis and cardiovascular diseases. Lipid traits exhibit considerable variation among populations of distinct ancestral origin as well as between individuals within a population. We performed association analyses to identify genetic loci influencing lipid concentrations in African American and Hispanic American women in the Women's Health Initiative SNP Health Association Resource. We validated one African-specific high-density lipoprotein cholesterol locus at CD36 as well as 14 known lipid loci that have been previously implicated in studies of European populations. Moreover, we demonstrate striking similarities in genetic architecture (loci influencing the trait, direction and magnitude of genetic effects, and proportions of phenotypic variation explained) of lipid traits across populations. In particular, we found that a disproportionate fraction of lipid variation in African Americans and Hispanic Americans can be attributed to genomic loci exhibiting statistical evidence of association in Europeans, even though the precise genes and variants remain unknown. At the same time, we found substantial allelic heterogeneity within shared loci, characterized both by population-specific rare variants and variants shared among multiple populations that occur at disparate frequencies. The allelic heterogeneity emphasizes the importance of including diverse populations in future genetic association studies of complex traits such as lipids; furthermore, the overlap in lipid loci across populations of diverse ancestral origin argues that additional knowledge can be gleaned from multiple populations.

Genetic architecture of skin and eye color in an African-European admixed population.

Variation in human skin and eye color is substantial and especially apparent in admixed populations, yet the underlying genetic architecture is poorly understood because most genome-wide studies are based on individuals of European ancestry. We study pigmentary variation in 699 individuals from Cape Verde, where extensive West African/European admixture has given rise to a broad range in trait values and genomic ancestry proportions. We develop and apply a new approach for measuring eye color, and identify two major loci (HERC2[OCA2] P = 2.3 × 10(-62), SLC24A5 P = 9.6 × 10(-9)) that account for both blue versus brown eye color and varying intensities of brown eye color. We identify four major loci (SLC24A5 P = 5.4 × 10(-27), TYR P = 1.1 × 10(-9), APBA2[OCA2] P = 1.5 × 10(-8), SLC45A2 P = 6 × 10(-9)) for skin color that together account for 35% of the total variance, but the genetic component with the largest effect (~44%) is average genomic ancestry. Our results suggest that adjacent cis-acting regulatory loci for OCA2 explain the relationship between skin and eye color, and point to an underlying genetic architecture in which several genes of moderate effect act together with many genes of small effect to explain ~70% of the estimated heritability.

Development and Validation of a Deep Learning Algorithm for Detection of Diabetic Retinopathy in Retinal Fundus Photographs.

Importance: Deep learning is a family of computational methods that allow an algorithm to program itself by learning from a large set of examples that demonstrate the desired behavior, removing the need to specify rules explicitly. Application of these methods to medical imaging requires further assessment and validation. Objective: To apply deep learning to create an algorithm for automated detection of diabetic retinopathy and diabetic macular edema in retinal fundus photographs. Design and Setting: A specific type of neural network optimized for image classification called a deep convolutional neural network was trained using a retrospective development data set of 128 175 retinal images, which were graded 3 to 7 times for diabetic retinopathy, diabetic macular edema, and image gradability by a panel of 54 US licensed ophthalmologists and ophthalmology senior residents between May and December 2015. The resultant algorithm was validated in January and February 2016 using 2 separate data sets, both graded by at least 7 US board-certified ophthalmologists with high intragrader consistency. Exposure: Deep learning-trained algorithm. Main Outcomes and Measures: The sensitivity and specificity of the algorithm for detecting referable diabetic retinopathy (RDR), defined as moderate and worse diabetic retinopathy, referable diabetic macular edema, or both, were generated based on the reference standard of the majority decision of the ophthalmologist panel. The algorithm was evaluated at 2 operating points selected from the development set, one selected for high specificity and another for high sensitivity. Results: The EyePACS-1 data set consisted of 9963 images from 4997 patients (mean age, 54.4 years; 62.2% women; prevalence of RDR, 683/8878 fully gradable images [7.8%]); the Messidor-2 data set had 1748 images from 874 patients (mean age, 57.6 years; 42.6% women; prevalence of RDR, 254/1745 fully gradable images [14.6%]). For detecting RDR, the algorithm had an area under the receiver operating curve of 0.991 (95% CI, 0.988-0.993) for EyePACS-1 and 0.990 (95% CI, 0.986-0.995) for Messidor-2. Using the first operating cut point with high specificity, for EyePACS-1, the sensitivity was 90.3% (95% CI, 87.5%-92.7%) and the specificity was 98.1% (95% CI, 97.8%-98.5%). For Messidor-2, the sensitivity was 87.0% (95% CI, 81.1%-91.0%) and the specificity was 98.5% (95% CI, 97.7%-99.1%). Using a second operating point with high sensitivity in the development set, for EyePACS-1 the sensitivity was 97.5% and specificity was 93.4% and for Messidor-2 the sensitivity was 96.1% and specificity was 93.9%. Conclusions and Relevance: In this evaluation of retinal fundus photographs from adults with diabetes, an algorithm based on deep machine learning had high sensitivity and specificity for detecting referable diabetic retinopathy. Further research is necessary to determine the feasibility of applying this algorithm in the clinical setting and to determine whether use of the algorithm could lead to improved care and outcomes compared with current ophthalmologic assessment.

Patient-derived xenografts of triple-negative breast cancer reproduce molecular features of patient tumors and respond to mTOR inhibition.

INTRODUCTION: Triple-negative breast cancer (TNBC) is aggressive and lacks targeted therapies. Phosphatidylinositide 3-kinase (PI3K)/mammalian target of rapamycin (mTOR) pathways are frequently activated in TNBC patient tumors at the genome, gene expression and protein levels, and mTOR inhibitors have been shown to inhibit growth in TNBC cell lines. We describe a panel of patient-derived xenografts representing multiple TNBC subtypes and use them to test preclinical drug efficacy of two mTOR inhibitors, sirolimus (rapamycin) and temsirolimus (CCI-779). METHODS: We generated a panel of seven patient-derived orthotopic xenografts from six primary TNBC tumors and one metastasis. Patient tumors and corresponding xenografts were compared by histology, immunohistochemistry, array comparative genomic hybridization (aCGH) and phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit alpha (PIK3CA) sequencing; TNBC subtypes were determined. Using a previously published logistic regression approach, we generated a rapamycin response signature from Connectivity Map gene expression data and used it to predict rapamycin sensitivity in 1,401 human breast cancers of different intrinsic subtypes, prompting in vivo testing of mTOR inhibitors and doxorubicin in our TNBC xenografts. RESULTS: Patient-derived xenografts recapitulated histology, biomarker expression and global genomic features of patient tumors. Two primary tumors had PIK3CA coding mutations, and five of six primary tumors showed flanking intron single nucleotide polymorphisms (SNPs) with conservation of sequence variations between primary tumors and xenografts, even on subsequent xenograft passages. Gene expression profiling showed that our models represent at least four of six TNBC subtypes. The rapamycin response signature predicted sensitivity for 94% of basal-like breast cancers in a large dataset. Drug testing of mTOR inhibitors in our xenografts showed 77 to 99% growth inhibition, significantly more than doxorubicin; protein phosphorylation studies indicated constitutive activation of the mTOR pathway that decreased with treatment. However, no tumor was completely eradicated. CONCLUSIONS: A panel of patient-derived xenograft models covering a spectrum of TNBC subtypes was generated that histologically and genomically matched original patient tumors. Consistent with in silico predictions, mTOR inhibitor testing in our TNBC xenografts showed significant tumor growth inhibition in all, suggesting that mTOR inhibitors can be effective in TNBC, but will require use with additional therapies, warranting investigation of optimal drug combinations.

A minimum variance method for genome-wide data-driven normalization of quantitative real-time polymerase chain reaction expression data.

Advances in multiplex qRT-PCR have enabled increasingly accurate and robust quantification of RNA, even at lower concentrations, facilitating RNA expression profiling in clinical and environmental samples. Here we describe a data-driven qRT-PCR normalization method, the minimum variance method, and evaluate it on clinically derived Mycobacterium tuberculosis samples with variable transcript detection percentages. For moderate to significant amounts of nondetection (∼50%), our minimum variance method consistently produces the lowest false discovery rates compared to commonly used data-driven normalization methods.

Ancestral components of admixed genomes in a Mexican cohort.

For most of the world, human genome structure at a population level is shaped by interplay between ancient geographic isolation and more recent demographic shifts, factors that are captured by the concepts of biogeographic ancestry and admixture, respectively. The ancestry of non-admixed individuals can often be traced to a specific population in a precise region, but current approaches for studying admixed individuals generally yield coarse information in which genome ancestry proportions are identified according to continent of origin. Here we introduce a new analytic strategy for this problem that allows fine-grained characterization of admixed individuals with respect to both geographic and genomic coordinates. Ancestry segments from different continents, identified with a probabilistic model, are used to construct and study "virtual genomes" of admixed individuals. We apply this approach to a cohort of 492 parent-offspring trios from Mexico City. The relative contributions from the three continental-level ancestral populations-Africa, Europe, and America-vary substantially between individuals, and the distribution of haplotype block length suggests an admixing time of 10-15 generations. The European and Indigenous American virtual genomes of each Mexican individual can be traced to precise regions within each continent, and they reveal a gradient of Amerindian ancestry between indigenous people of southwestern Mexico and Mayans of the Yucatan Peninsula. This contrasts sharply with the African roots of African Americans, which have been characterized by a uniform mixing of multiple West African populations. We also use the virtual European and Indigenous American genomes to search for the signatures of selection in the ancestral populations, and we identify previously known targets of selection in other populations, as well as new candidate loci. The ability to infer precise ancestral components of admixed genomes will facilitate studies of disease-related phenotypes and will allow new insight into the adaptive and demographic history of indigenous people.

Transcript levels of androgen receptor variant AR-V1 or AR-V7 do not predict recurrence in patients with prostate cancer at indeterminate risk for progression.

PURPOSE: AR-V7, a ligand independent splice variant of androgen receptor, may support the growth of castration resistant prostate cancer and have prognostic value. Another variant, AR-V1, interferes with AR-V7 activity. We investigated whether AR-V7 or V1 expression would predict biochemical recurrence in men at indeterminate (about 50%) risk for progression following radical prostatectomy. MATERIALS AND METHODS: AR-V7 and V1 transcripts in a mixed grade cohort of 53 men in whom cancer contained 30% to 70% Gleason grade 4/5 and in a grade 3 only cohort of 52 were measured using a branched chain DNA assay. Spearman rank correlations of the transcripts, and histomorphological and clinical variables were determined. AR-V7 and V1 levels were assessed as determinants of recurrence in the mixed grade cohort by logistic regression and survival analysis. The impact of TMPRSS2-ERG gene fusion on prognosis was also evaluated. RESULTS: Neither AR-V7 nor V1 levels in grade 3 or 4/5 cancer in the mixed grade cohort were associated with recurrence or time to recurrence. However, AR-V7 and V1 inversely correlated with serum prostate specific antigen and positively correlated with age. The AR-V1 level in grade 3 cancer in the grade 3 only cohort was higher than in grade 3 or grade 4/5 components of mixed grade cancer. TMPRSS2-ERG fusion was not associated with AR-V7, AR-V1 or recurrence but it was associated with the percent of grade 4/5 cancer. CONCLUSIONS: The AR-V1 or V7 transcript level does not predict recurrence in patients with high grade prostate cancer at indeterminate risk for progression. Grade 3 cancer in mixed grade tumors may differ from 100% grade 3 cancer, at least in AR-V1 expression.

Antibodies specifically target AML antigen NuSAP1 after allogeneic bone marrow transplantation.

Identifying the targets of immune response after allogeneic hematopoietic cell transplantation (HCT) promises to provide relevant immune therapy candidate proteins. We used protein microarrays to serologically identify nucleolar and spindle-associated protein 1 (NuSAP1) and chromatin assembly factor 1, subunit B (p60; CHAF1b) as targets of new antibody responses that developed after allogeneic HCT. Western blots and enzyme-linked immunosorbent assays (ELISA) validated their post-HCT recognition and enabled ELISA testing of 120 other patients with various malignancies who underwent allo-HCT. CHAF1b-specific antibodies were predominantly detected in patients with acute myeloid leukemia (AML), whereas NuSAP1-specific antibodies were exclusively detected in patients with AML 1 year after transplantation (P < .001). Complete genomic exon sequencing failed to identify a nonsynonymous single nucleotide polymorphism (SNP) for NuSAP1 and CHAF1b between the donor and recipient cells. Expression profiles and reverse transcriptase-polymerase chain reaction (RT-PCR) showed NuSAP1 was predominately expressed in the bone marrow CD34(+)CD90(+) hematopoietic stem cells, leukemic cell lines, and B lymphoblasts compared with other tissues or cells. Thus, NuSAP1 is recognized as an immunogenic antigen in 65% of patients with AML following allogeneic HCT and suggests a tumor antigen role.

Morphological Characterization of Antibiotic Combinations.

Combination therapies are common in many therapeutic contexts, including infectious diseases and cancer. A common approach for evaluating combinations in vitro is to assess effects on cell growth as synergistic, antagonistic, or neutral using "checkerboard" experiments to systematically sample combinations of agents in multiple doses. To further understand the effects of antibiotic combinations, we employed high-content imaging to study the morphological changes caused by combination treatments in checkerboard experiments. Using an automated, unsupervised image analysis approach to group morphologies, and an "expert-in-the-loop" to annotate them, we attributed the heterogeneous morphological changes ofEscherichia coli cells to varying doses of both single-agent and combination treatments. We identified patterns of morphological change, including morphological potentiation, competition, and the emergence of unexpected morphologies. We found these frequently did not correlate with synergistic or antagonistic effects on viability, suggesting morphological approaches may provide a distinctive signature of the biological interaction between compounds over a range of conditions. Among the unexpected morphologies we observed, there were transitional changes associated with intermediate doses of compounds and uncharacterized phenotypes associated with well-studied antibiotics. Our approach exemplifies how unsupervised image analysis and expert knowledge can be combined to reckon with complex phenotypic changes arising from combination screening, dose titrations, or polypharmacology. In this way, quantification of morphological diversity across treatment conditions could aid in analysis and prioritization of complementary pairings of antibiotic agents by more closely capturing the true signature of the integrated cellular response.

ProtSeq: Toward high-throughput, single-molecule protein sequencing via amino acid conversion into DNA barcodes.

We demonstrate early progress toward constructing a high-throughput, single-molecule protein sequencing technology utilizing barcoded DNA aptamers (binders) to recognize terminal amino acids of peptides (targets) tethered on a next-generation sequencing chip. DNA binders deposit unique, amino acid-identifying barcodes on the chip. The end goal is that, over multiple binding cycles, a sequential chain of DNA barcodes will identify the amino acid sequence of a peptide. Toward this, we demonstrate successful target identification with two sets of target-binder pairs: DNA-DNA and Peptide-Protein. For DNA-DNA binding, we show assembly and sequencing of DNA barcodes over six consecutive binding cycles. Intriguingly, our computational simulation predicts that a small set of semi-selective DNA binders offers significant coverage of the human proteome. Toward this end, we introduce a binder discovery pipeline that ultimately could merge with the chip assay into a technology called ProtSeq, for future high-throughput, single-molecule protein sequencing.

Integrating deep learning and unbiased automated high-content screening to identify complex disease signatures in human fibroblasts.

Drug discovery for diseases such as Parkinson's disease are impeded by the lack of screenable cellular phenotypes. We present an unbiased phenotypic profiling platform that combines automated cell culture, high-content imaging, Cell Painting, and deep learning. We applied this platform to primary fibroblasts from 91 Parkinson's disease patients and matched healthy controls, creating the largest publicly available Cell Painting image dataset to date at 48 terabytes. We use fixed weights from a convolutional deep neural network trained on ImageNet to generate deep embeddings from each image and train machine learning models to detect morphological disease phenotypes. Our platform's robustness and sensitivity allow the detection of individual-specific variation with high fidelity across batches and plate layouts. Lastly, our models confidently separate LRRK2 and sporadic Parkinson's disease lines from healthy controls (receiver operating characteristic area under curve 0.79 (0.08 standard deviation)), supporting the capacity of this platform for complex disease modeling and drug screening applications.

Performance of a Deep-Learning Algorithm vs Manual Grading for Detecting Diabetic Retinopathy in India.

IMPORTANCE: More than 60 million people in India have diabetes and are at risk for diabetic retinopathy (DR), a vision-threatening disease. Automated interpretation of retinal fundus photographs can help support and scale a robust screening program to detect DR. OBJECTIVE: To prospectively validate the performance of an automated DR system across 2 sites in India. DESIGN, SETTING, AND PARTICIPANTS: This prospective observational study was conducted at 2 eye care centers in India (Aravind Eye Hospital and Sankara Nethralaya) and included 3049 patients with diabetes. Data collection and patient enrollment took place between April 2016 and July 2016 at Aravind and May 2016 and April 2017 at Sankara Nethralaya. The model was trained and fixed in March 2016. INTERVENTIONS: Automated DR grading system compared with manual grading by 1 trained grader and 1 retina specialist from each site. Adjudication by a panel of 3 retinal specialists served as the reference standard in the cases of disagreement. MAIN OUTCOMES AND MEASURES: Sensitivity and specificity for moderate or worse DR or referable diabetic macula edema. RESULTS: Of 3049 patients, 1091 (35.8%) were women and the mean (SD) age for patients at Aravind and Sankara Nethralaya was 56.6 (9.0) years and 56.0 (10.0) years, respectively. For moderate or worse DR, the sensitivity and specificity for manual grading by individual nonadjudicator graders ranged from 73.4% to 89.8% and from 83.5% to 98.7%, respectively. The automated DR system's performance was equal to or exceeded manual grading, with an 88.9% sensitivity (95% CI, 85.8-91.5), 92.2% specificity (95% CI, 90.3-93.8), and an area under the curve of 0.963 on the data set from Aravind Eye Hospital and 92.1% sensitivity (95% CI, 90.1-93.8), 95.2% specificity (95% CI, 94.2-96.1), and an area under the curve of 0.980 on the data set from Sankara Nethralaya. CONCLUSIONS AND RELEVANCE: This study shows that the automated DR system generalizes to this population of Indian patients in a prospective setting and demonstrates the feasibility of using an automated DR grading system to expand screening programs.

H-Y antibody development associates with acute rejection in female patients with male kidney transplants.

BACKGROUND: Human minor histocompatibility antigens (mHA) and clinically relevant immune responses to them have not been well defined in organ transplantation. We hypothesized that women with male kidney transplants would develop antibodies against H-Y, the mHA encoded on the Y-chromosome, in association with graft rejection. METHODS: We tested sera from 118 consecutive transplant recipients with kidney biopsies. Antibodies that specifically recognized the recombinant H-Y antigens RPS4Y1 or DDX3Y were detected by IgG enzyme-linked immunosorbent assay and western blotting. Immunogenic epitopes were further identified using overlapping H-Y antigen peptides for both the H-Y proteins. RESULTS: In the 26 female recipients of male kidneys, H-Y antibody development posttransplant (1) was more frequent (46%) than in other gender combinations (P<0.001), (2) showed strong correlation with acute rejection (P=0.00048), (3) correlated with plasma cell infiltrates in biopsied kidneys (P=0.04), and (4) did not correlate with C4d deposition or donor-specific anti-human leukocyte antigen (HLA) antibodies. Of the two H-Y antigens, RPS4Y1 was more frequently recognized (P=0.005). CONCLUSION: This first demonstration of a strong association between H-Y antibody development and acute rejection in kidney transplant recipients shows that in solid organ allografts, humoral immune responses against well defined mHA have clear clinical correlates, can be easily monitored, and warrant study for possible effects on long-term graft function.

The significance of monoamine oxidase-A expression in high grade prostate cancer.

PURPOSE: Gleason grade 4/5 prostate cancer is a determinant for recurrence following radical prostatectomy. Monoamine oxidase-A is over expressed in grade 4/5 compared to grade 3 cancer. Monoamine oxidase-A is also expressed by normal basal cells and in vitro studies suggest that its function is to repress secretory differentiation. Therefore, monoamine oxidase-A in grade 4/5 cancer might reflect dedifferentiation to a basal cell-like phenotype. We investigated whether monoamine oxidase-A expression correlates with another basal cell protein, CD44, in high grade cancer and whether either is associated with an aggressive phenotype. MATERIALS AND METHODS: A total of 133 grade 4/5 archival cancers from a cohort previously used to evaluate the prognostic significance of histomorphological variables were scored for monoamine oxidase-A and CD44 immunohistochemical labeling. Spearman rank correlations of the proteins, and histomorphological and clinical variables were determined. The univariate and multivariate value of each variable as a determinant of biochemical recurrence was assessed by logistic regression. RESULTS: Monoamine oxidase-A expression correlated with CD44. Neither was prognostic for biochemical recurrence. However, monoamine oxidase-A expression positively correlated with preoperative serum prostate specific antigen and the percent of grade 4/5 cancer. CONCLUSIONS: Concurrent expression of monoamine oxidase-A and CD44 suggests that grade 4/5 cancer may be basal cell-like in nature, despite the absence of other classic basal cell biomarkers such as cytokeratins 5 and 14, and p63. The correlation of monoamine oxidase-A expression with prostate specific antigen and the percent of grade 4/5 cancer suggests that monoamine oxidase-A may contribute to growth of high grade cancer and that antidepressant drugs that target monoamine oxidase-A may have applications in treating prostate cancer.

A suite of algorithms for the comprehensive analysis of complex protein mixtures using high-resolution LC-MS.

MOTIVATION: Comparing two or more complex protein mixtures using liquid chromatography mass spectrometry (LC-MS) requires multiple analysis steps to locate and quantitate natural peptides within a single experiment and to align and normalize findings across multiple experiments. RESULTS: We describe msInspect, an open-source application comprising algorithms and visualization tools for the analysis of multiple LC-MS experimental measurements. The platform integrates novel algorithms for detecting signatures of natural peptides within a single LC-MS measurement and combines multiple experimental measurements into a peptide array, which may then be mined using analysis tools traditionally applied to genomic array analysis. The platform supports quantitation by both label-free and isotopic labeling approaches. The software implementation has been designed so that many key components may be easily replaced, making it useful as a workbench for integrating other novel algorithms developed by a growing research community. AVAILABILITY: The msInspect software is distributed freely under an Apache 2.0 license. The software as well as a Zip file with all peptide feature files and scripts needed to generate the tables and figures in this article are available at http://proteomics.fhcrc.org/.