Search

Metabolites significantly changed by p62 knockdown in Tsc2-/- MEFs

p62/sequestosome-1 (SQSTM1) is a multifunctional adaptor protein and autophagic substrate that accumulates in cells with hyperactive mTORC1, such as kidney cells with mutations in the tumor suppressor genes tuberous sclerosis complex (TSC)1 or TSC2. Here we report that p62 is a critical mediator of TSC2-driven tumorigenesis, as Tsc2+/− and Tsc2f/f CAGGCreERT2+ mice crossed to p62−/− mice were protected from renal tumor development. Metabolic profiling revealed that depletion of p62 in Tsc2-null cells decreased intracellular glutamine, glutamate, and glutathione (GSH). p62 positively regulated the glutamine transporter Slc1a5 and increased glutamine uptake in Tsc2-null cells. We also observed p62-dependent changes in Gcl, Gsr, Nqo1, and Srxn1, which were decreased by p62 attenuation and implicated in GSH production and utilization. p62 attenuation altered mitochondrial morphology, reduced mitochondrial membrane polarization and maximal respiration, and increased mitochondrial reactive oxygen species and mitophagy marker PINK1. These mitochondrial phenotypes were rescued by addition of exogenous GSH and overexpression of Sod2, which suppressed indices of mitochondrial damage and promoted growth of Tsc2-null cells. Finally, p62 depletion sensitized Tsc2-null cells to both oxidative stress and direct inhibition of GSH biosynthesis by buthionine sulfoximine. Our findings show how p62 helps maintain intracellular pools of GSH needed to limit mitochondrial dysfunction in tumor cells with elevated mTORC1, highlighting p62 and redox homeostasis as nodal vulnerabilities for therapeutic targeting in these tumors. Cancer Res; 77(12); 3255–67. ©2017 AACR.

Sod2 overexpression reduces the percentage of MitoSOX positive cells assessed by flow cytometry

Supplementary Materials and Methods

Raw data and bioinformatic analysis of metabolomics and ion torrent expression profiling

Utilizing publicly available ribonucleic acid sequencing data, we identified SCUBE1 as a BMPR2-related gene differentially expressed between induced pluripotent stem cell-endothelial cells derived from pulmonary arterial hypertension (PAH) patients carrying pathogenic BMPR2 mutations and control patients without mutations. Endothelial SCUBE1 expression was decreased by known triggers of PAH, and its down-regulation recapitulated known BMPR2-associated endothelial pathophenotypes in vitro. Meanwhile, SCUBE1 concentrations were reduced in plasma obtained from PAH rodent models and patients with PAH, whereas plasma concentrations were tightly correlated with hemodynamic markers of disease severity. Taken together, these data implicate SCUBE1 as a novel contributor to PAH pathogenesis with potential therapeutic, diagnostic, and prognostic applications.

Description, A network dependency platform was generated to define a landscape of cancer drug mechanisms in pulmonary hypertension., Cancer therapies are being considered for treating rare noncancerous diseases like pulmonary hypertension (PH), but effective computational screening is lacking. Via transcriptomic differential dependency analyses leveraging parallels between cancer and PH, we mapped a landscape of cancer drug functions dependent upon rewiring of PH gene clusters. Bromodomain and extra-terminal motif (BET) protein inhibitors were predicted to rely upon several gene clusters inclusive of galectin-8 (LGALS8). Correspondingly, LGALS8 was found to mediate the BET inhibitor–dependent control of endothelial apoptosis, an essential role for PH in vivo. Separately, a piperlongumine analog’s actions were predicted to depend upon the iron-sulfur biogenesis gene ISCU. Correspondingly, the analog was found to inhibit ISCU glutathionylation, rescuing oxidative metabolism, decreasing endothelial apoptosis, and improving PH. Thus, we identified crucial drug-gene axes central to endothelial dysfunction and therapeutic priorities for PH. These results establish a wide-ranging, network dependency platform to redefine cancer drugs for use in noncancerous conditions.

AimsEndothelial cell (EC) dysfunction drives the initiation and pathogenesis of pulmonary arterial hypertension (PAH). We aimed to characterize EC dynamics in PAH at single-cell resolution.Methods and resultsWe carried out single-cell RNA sequencing (scRNA-seq) of lung ECs isolated from an EC lineage-tracing mouse model in Control and SU5416/hypoxia-induced PAH conditions. EC populations corresponding to distinct lung vessel types, including two discrete capillary populations, were identified in both Control and PAH mice. Differential gene expression analysis revealed global PAH-induced EC changes that were confirmed by bulk RNA-seq. This included upregulation of the major histocompatibility complex class II pathway, supporting a role for ECs in the inflammatory response in PAH. We also identified a PAH response specific to the second capillary EC population including upregulation of genes involved in cell death, cell motility, and angiogenesis. Interestingly, four genes with genetic variants associated with PAH were dysregulated in mouse ECs in PAH. To compare relevance across PAH models and species, we performed a detailed analysis of EC heterogeneity and response to PAH in rats and humans through whole-lung PAH scRNA-seq datasets, revealing that 51% of up-regulated mouse genes were also up-regulated in rat or human PAH. We identified promising new candidates to target endothelial dysfunction including CD74, the knockdown of which regulates EC proliferation and barrier integrity in vitro. Finally, with an in silico cell ordering approach, we identified zonation-dependent changes across the arteriovenous axis in mouse PAH and showed upregulation of the Serine/threonine-protein kinase Sgk1 at the junction between the macro- and microvasculature.ConclusionThis study uncovers PAH-induced EC transcriptomic changes at a high resolution, revealing novel targets for potential therapeutic candidate development.

Background Preoperative pulmonary hypertension (PH) is associated with excess mortality among patients with severe mitral regurgitation undergoing mitral valve surgery (MVS). However, the links between PH phenotype, pulmonary vascular remodeling, and persistent postoperative PH are not well understood. We aimed to describe the associations between components of pulmonary hemodynamics as well as postoperative residual PH with longitudinal mortality in patients with severe mitral regurgitation who received MVS. Methods and Results Patients undergoing MVS for severe mitral regurgitation from 2011 to 2016 were retrospectively identified within our health system (n=488). Mean pulmonary artery pressure and other hemodynamic variables were determined by presurgical right‐heart catheterization. Postoperative pulmonary artery systolic pressure was assessed on echocardiogram 42 to 365 days post‐MVS. Longitudinal survival over a mean 3.9 years of follow‐up was evaluated using Cox proportional hazards modeling to compare survival after adjustment for demographics, surgical characteristics, and comorbidities. Pre‐MVS prevalence of PH was high at 85%. After adjustment, each 10‐mm Hg increase in preoperative mean pulmonary artery pressure was associated with a 1.38‐fold increase in risk of death (95% CI, 1.13–1.68). Elevated preoperative pulmonary vascular resistance, transpulmonary gradient, and right atrial pressure were similarly associated with increased mortality. Among 231 patients with postoperative echocardiogram, evidence of PH on echocardiogram (pulmonary artery systolic pressure ≥35 mm Hg) was associated with increased risk of death (hazard ratio [HR], 2.02 [95% CI, 1.17–3.47]); however, this was no longer statistically significant after adjustment (HR, 1.55 [95% CI, 0.85–2.85]). Conclusions In patients undergoing MVS for mitral regurgitation, preoperative PH, and postoperative PH were associated with increased mortality.

International audience; Pulmonary arterial hypertension (PAH) refers to a set of heterogeneous vascular diseases defined by elevation of pulmonary arterial pressure (PAP) and pulmonary vascular resistance (PVR), leading to right ventricular (RV) remodeling and often death. Early increases in pulmonary artery stiffness in PAH drive pathogenic alterations of pulmonary arterial endothelial cells (PAECs), leading to vascular remodeling. Dysregulation of microRNAs can drive PAEC dysfunction. However, the role of vascular stiffness in regulating pathogenic microRNAs in PAH is incompletely understood. Here, we demonstrated that extracellular matrix (ECM) stiffening downregulated miR-7 levels in PAECs. The RNA-binding protein quaking (QKI) has been implicated in the biogenesis of miR-7. Correspondingly, we found that ECM stiffness upregulated QKI, and QKI knockdown led to increased miR-7. Downstream of the QKI-miR-7 axis, the serine and arginine-rich splicing factor 1 (SRSF1) was identified as a direct target of miR-7. Correspondingly, SRSF1 was reciprocally upregulated in PAECs exposed to stiff ECM and was negatively correlated with miR-7. Decreased miR-7 and increased QKI and SRSF1 were observed in lungs from patients with PAH and PAH rats exposed to SU5416/hypoxia. Lastly, miR-7 upregulation inhibited human PAEC migration, whereas forced SRSF1 expression reversed this phenotype, proving that miR-7 depended upon SRSF1 to control migration. In aggregate, these results define the QKI-miR-7-SRSF1 axis as a mechanosensitive mechanism linking pulmonary arterial vascular stiffness to pathogenic endothelial function. These findings emphasize implications relevant to PAH and suggest the potential benefit of developing therapies that target this miRNA-dependent axis in PAH.

Introduction: Early increases in pulmonary artery stiffness drive pathogenic alterations of pulmonary arterial endothelial cells (PAECs), leading to vascular remodeling, a prominent feature of pulmonary hypertension (PH). However, the molecular mechanisms by which PAECs respond to altered extracellular matrix (ECM) mechanics remain unclear. Recent RNA sequencing study of PAECs exposed to stiff versus soft matrices has identified miRNA-7 (miR-7) as a mechanosensitive downregulated miRNA. The biogenesis of miR-7 is inhibited by the RNA binding protein quaking (QKI), which is shown to be correspondingly upregulated in ECM stiffness. Moreover, subsequent RNA sequencing and computational analyses have identified serine and arginine rich splicing factor 1 (SRSF1) as a miR-7 target in PAECs with ECM stiffening. Hypothesis: The QKI-miR-7-SRSF1 axis is a pivotal ECM stiffening-mediated pathogenic signal in PAECs. Methods: Human PAECs were cultured on soft (0.5 or 1 kPa) or stiff (50 kPa) matrices. Lungs from PH patients and rats (SU5416/Hypoxia- and MCT-exposed models) were used to evaluate genes and miRNA expression. Boyden chamber migration analysis examined PAEC mobility. Results: In human and rat PH lungs, miRNA-7 (miR-7) was down-regulated ( P < 0.05). To validate RNA sequencing data by RT-qPCR, QKI level was increased (2.13 ± 0.19, P P P P Conclusions: Our data demonstrate that mechanoactivation of a QKI-miR-7-SRSF1 signal axis promotes PAEC pathophenotypes driven by ECM stiffness. These insights broaden our understanding of RNA pathobiology in vascular stiffening and carry important therapeutic implications in PH.

Introduction: Pulmonary arterial hypertension (PAH) is a morbid vascular disease where mutations of bone morphogenetic protein receptor 2 (BMPR2) control pulmonary endothelial pathophenotypes. Transcriptomic screening from endothelial cells (ECs) derived from hereditary BMPR2-driven disease offers an opportunity to identify novel effectors in PAH pathogenesis. Methods: Public RNA-sequencing data were analyzed from inducible pluripotent stem (iPS) cell-derived endothelial cells with and without BMPR2 mutations. Candidate genes were assessed in cultured pulmonary arterial ECs (PAECs) and in rodent and human PAH. Correlations were assessed of human plasma expression with clinical disease indices. Results: Signal Peptide CUB-EGF-Domain Containing Protein 1 (SCUBE1), a putative binding partner to BMPR2, was differentially expressed in ECs carrying BMPR2 mutations. SCUBE1 was enriched in PAECs, dependent on hypoxia inducible factor-1α, and downregulated by PAH triggers, including BMPR2 knockdown, hypoxia, and IL-1β exposure. In vitro analyses defined SCUBE1 as a pathogenic effector activating BMPR2-associated SMAD1/5/9, thus regulating endothelial angiogenic potential, proliferation, and apoptosis. SCUBE1 was decreased specifically in plasma and lungs in rodents and patients with PAH but not in those with pulmonary hypertension due to left heart disease, or in the patients of other acute and chronic cardiopulmonary pathologies, including pneumonia, acute lung injury, chronic obstructive pulmonary disease, ischemic heart disease, and cardiomyopathy. An optimal plasma SCUBE1 cut point of 5.02 ng/mL was defined to diagnose PAH from control cohort with a sensitivity of 0.65 and a specificity of 0.82. In PAH patients, plasma SCUBE1 levels negatively correlated with pulmonary arterial pressure, pulmonary vascular resistance, and right ventricular dysfunction. Conclusions: Guided by iPSC-EC sequencing, SCUBE1 was identified as a downregulated secreted factor in PAH, controlling endothelial pathophenotypes and correlated with disease indices. Clinically, SCUBE1 is a sensitive and specific PAH diagnostic marker. It may also serve as a therapeutic target given its inherent links to controlling PAH predisposition and severity.

Pulmonary arterial hypertension (PAH) is a deadly vascular disease, characterized by increased pulmonary arterial pressures and right heart failure. Considering prior non-US studies of atrial arrhythmias in PAH, this retrospective, regional multi-center US study sought to define more completely the risk factors and impact of paroxysmal and non-paroxysmal forms of atrial fibrillation and flutter (AF/AFL) on mortality in this disease. We identified patients seen between 2010 and 2014 at UPMC (Pittsburgh) hospitals with hemodynamic and clinical criteria for PAH or chronic thromboembolic pulmonary hypertension (CTEPH) and determined those meeting electrocardiographic criteria for AF/AFL. We used Cox proportional hazards regression with time-varying covariates to analyze the association between AF/AFL occurrence and survival with adjustments for potential cofounders and hemodynamic severity. Of 297 patients with PAH/CTEPH, 79 (26.5%) suffered from AF/AFL at some point. AF/AFL was first identified after PAH diagnosis in 42 (53.2%), identified prior to PAH diagnosis in 27 (34.2%), and had unclear timing in the remainder. AF/AFL patients were older, more often male, had lower left ventricular ejection fractions, and greater left atrial volume indices and right atrial areas than patients without AF/AFL. AF/AFL (whether diagnosed before or after PAH) was associated with a 3.81-fold increase in the hazard of death (95% CI 2.64–5.52, p

// Hilaire C. Lam 1 , Heng-Jia Liu 1 , Christian V. Baglini 1 , Harilaos Filippakis 1 , Nicola Alesi 1 , Julie Nijmeh 1 , Heng Du 1 , Alicia Llorente Lope 1 , Katherine A. Cottrill 2 , Adam Handen 2 , John M. Asara 3 , David J. Kwiatkowski 1 , Issam Ben-Sahra 4 , William M. Oldham 1 , Stephen Y. Chan 2 and Elizabeth P. Henske 1 1 Department of Medicine, Pulmonary and Critical Care Medicine, Brigham and Women’s Hospital and Harvard Medical School, Boston, MA, USA 2 Department of Medicine, Division of Cardiology, Center for Pulmonary Vascular Biology and Medicine, Pittsburgh Heart, Lung, Blood, and Vascular Medicine Institute, University of Pittsburgh School of Medicine and University of Pittsburgh Medical Center, Pittsburgh, PA, USA 3 Department of Medicine, Division of Signal Transduction, Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, MA, USA 4 Department of Biochemistry and Molecular Genetics, Northwestern University, Chicago, IL, USA Correspondence to: Elizabeth P. Henske, email: // Keywords : tuberous sclerosis complex, mTORC1, rapamycin, miR-21, mitochondria Received : June 20, 2017 Accepted : June 25, 2017 Published : August 04, 2017 Abstract mTORC1 hyperactivation drives the multi-organ hamartomatous disease tuberous sclerosis complex (TSC). Rapamycin inhibits mTORC1, inducing partial tumor responses; however, the tumors regrow following treatment cessation. We discovered that the oncogenic miRNA, miR-21, is increased in Tsc2-deficient cells and, surprisingly, further increased by rapamycin. To determine the impact of miR-21 in TSC, we inhibited miR-21 in vitro . miR-21 inhibition significantly repressed the tumorigenic potential of Tsc2-deficient cells and increased apoptosis sensitivity. Tsc2-deficient cells’ clonogenic and anchorage independent growth were reduced by ~50% ( p

p62/sequestosome-1 (SQSTM1) is a multifunctional adaptor protein and autophagic substrate that accumulates in cells with hyperactive mTORC1, such as kidney cells with mutations in the tumor suppressor genes tuberous sclerosis complex (TSC)1 or TSC2. Here we report that p62 is a critical mediator of TSC2-driven tumorigenesis, as Tsc2+/− and Tsc2f/f CAGGCreERT2+ mice crossed to p62−/− mice were protected from renal tumor development. Metabolic profiling revealed that depletion of p62 in Tsc2-null cells decreased intracellular glutamine, glutamate, and glutathione (GSH). p62 positively regulated the glutamine transporter Slc1a5 and increased glutamine uptake in Tsc2-null cells. We also observed p62-dependent changes in Gcl, Gsr, Nqo1, and Srxn1, which were decreased by p62 attenuation and implicated in GSH production and utilization. p62 attenuation altered mitochondrial morphology, reduced mitochondrial membrane polarization and maximal respiration, and increased mitochondrial reactive oxygen species and mitophagy marker PINK1. These mitochondrial phenotypes were rescued by addition of exogenous GSH and overexpression of Sod2, which suppressed indices of mitochondrial damage and promoted growth of Tsc2-null cells. Finally, p62 depletion sensitized Tsc2-null cells to both oxidative stress and direct inhibition of GSH biosynthesis by buthionine sulfoximine. Our findings show how p62 helps maintain intracellular pools of GSH needed to limit mitochondrial dysfunction in tumor cells with elevated mTORC1, highlighting p62 and redox homeostasis as nodal vulnerabilities for therapeutic targeting in these tumors. Cancer Res; 77(12); 3255–67. ©2017 AACR.

Pulmonary arterial hypertension (PAH) is characterized by increased pulmonary artery pressure and vascular resistance, typically leading to right heart failure and death. Current therapies improve quality of life of the patients but have a modest effect on long-term survival. A detailed transcriptomics and systems biology view of the PAH lung is expected to provide new testable hypotheses for exploring novel treatments. We completed transcriptomics analysis of PAH and control lung tissue to develop disease-specific and clinical data/tissue pathology gene expression classifiers from expression datasets. Gene expression data were integrated into pathway analyses. Gene expression microarray data were collected from 58 PAH and 25 control lung tissues. The strength of the dataset and its derived disease classifier was validated using multiple approaches. Pathways and upstream regulators analyses was completed with standard and novel graphical approaches. The PAH lung dataset identified expression patterns specific to PAH subtypes, clinical parameters, and lung pathology variables. Pathway analyses indicate the important global role of TNF and transforming growth factor signaling pathways. In addition, novel upstream regulators and insight into the cellular and innate immune responses driving PAH were identified. Finally, WNT-signaling pathways may be a major determinant underlying the observed sex differences in PAH. This study provides a transcriptional framework for the PAH-diseased lung, supported by previously reported findings, and will be a valuable resource to the PAH research community. Our investigation revealed novel potential targets and pathways amenable to further study in a variety of experimental systems.

p62/sequestosome-1 (SQSTM1) is a multifunctional adaptor protein and autophagic substrate which accumulates in cells with hyperactive mTORC1, such as kidney cells with mutations in the tumor suppressor genes TSC1 or TSC2. Here we report that p62 is a critical mediator of TSC2-driven tumorigenesis, as Tsc2+/− and Tsc2f/f Ksp-CreERT2+ mice crossed to p62−/− mice were protected from renal tumor development. Metabolic profiling revealed that depletion of p62 in Tsc2-null cells decreased intracellular glutamine, glutamate, and glutathione (GSH). p62 positively regulated the glutamine transporter Slc1a5 and increased glutamine uptake in Tsc2-null cells. We also observed p62-dependent changes in Gcl, Gsr, Nqo1 and Srxn1 which were decreased by p62 attenuation and implicated in GSH production and utilization. p62 attenuation altered mitochondrial morphology, reduced mitochondrial membrane polarization and maximal respiration, and increased mitochondrial ROS and mitophagy marker PINK1. These mitochondrial phenotypes were rescued by addition of exogenous GSH and overexpression of Sod2, which suppressed indices of mitochondrial damage and promoted growth of Tsc2-null cells. Finally, p62 depletion sensitized Tsc2-null cells to both oxidative stress and direct inhibition of glutathione biosynthesis by buthionine sulfoximine (BSO). Our findings show how p62 helps maintain intracellular pools of glutathione needed to limit mitochondrial dysfunction in tumor cells with elevated mTORC1, highlighting p62 and redox homeostasis as nodal vulnerabilities for therapeutic targeting in these tumors.

While the implications and outcomes of pulmonary hypertension (PH) complicating mitral stenosis have been well described, PH in patients undergoing mitral valve surgery (MVS) for severe mitral regurgitation (MR) has been comparatively understudied. We aimed to describe post-operative mortality in

Pulmonary arterial hypertension (PAH) is characterized by occlusion of lung arterioles, leading to marked increases in pulmonary vascular resistance. Although heritable forms of PAH are known to be driven by genetic mutations that share some commonality of function, the extent to which these effectors converge to regulate shared processes in this disease is unknown. We have causally connected extracellular matrix (ECM) remodeling and mechanotransduction to the miR-130/301 family in a feedback loop that drives vascular activation and downstream PAH. However, the molecular interconnections between factors genetically associated with PAH and this mechano-driven feedback loop remain undefined. We performed systematic manipulation of matrix stiffness, the miR-130/301 family, and factors genetically associated with PAH in primary human pulmonary arterial cells and assessed downstream and reciprocal consequences on their expression. We found that a network of factors linked to heritable PAH converges upon the matrix stiffening-miR-130/301-PPAR&gamma, LRP8 axis in order to remodel the ECM. Furthermore, we leveraged a computational network biology approach to predict a number of additional molecular circuits functionally linking this axis to the ECM. These results demonstrate that multiple genes associated with heritable PAH converge to control the miR-130/301 circuit, triggering a self-amplifying feedback process central to pulmonary vascular stiffening and disease.

When a set of genes are identified to be related to a disease, say through gene expression analysis, it is common to examine the average distance among their protein products in the human interactome as a measure of biological relatedness of these genes. The reasoning for this is that, genes associated with a disease would tend to be functionally related, and that functionally related genes would be closely connected to each other in the interactome. Typically, average shortest path length (ASPL) of disease genes (although referred to as genes in the context of disease-associations, the interactions are among protein-products of these genes) is compared to ASPL of randomly selected genes or to ASPL in a randomly permuted network. We examined whether the ASPL of a set of genes is indeed a good measure of biological relatedness or whether it is simply a characteristic of the degree distribution of those genes. We examined the ASPL of genes sets of some disease and pathway associations and compared them to ASPL of three types of randomly selected control sets: uniform selection, from entire proteome, degree-matched selection, and random permutation of the network. We found that disease associated genes and their degree-matched random genes have comparable ASPL. In other words, ASPL is a characteristic of the degree of the genes and the network topology, and not that of functional coherence.

Sleep apnea increases pulmonary artery pressure at night, but whether sleep apnea causes persistent daytime pulmonary hypertension is less clear. Adults who underwent an overnight sleep study within our large heath care system between 3/1999 and 12/2017 were eligible for inclusion if they also

Genome-wide association studies of schizophrenia (GWAS) have revealed the role of rare and common genetic variants, but the functional effects of the risk variants remain to be understood. Protein interactome-based studies can facilitate the study of molecular mechanisms by which the risk genes relate to schizophrenia (SZ) genesis, but protein–protein interactions (PPIs) are unknown for many of the liability genes. We developed a computational model to discover PPIs, which is found to be highly accurate according to computational evaluations and experimental validations of selected PPIs. We present here, 365 novel PPIs of liability genes identified by the SZ Working Group of the Psychiatric Genomics Consortium (PGC). Seventeen genes that had no previously known interactions have 57 novel interactions by our method. Among the new interactors are 19 drug targets that are targeted by 130 drugs. In addition, we computed 147 novel PPIs of 25 candidate genes investigated in the pre-GWAS era. While there is little overlap between the GWAS genes and the pre-GWAS genes, the interactomes reveal that they largely belong to the same pathways, thus reconciling the apparent disparities between the GWAS and prior gene association studies. The interactome including 504 novel PPIs overall, could motivate other systems biology studies and trials with repurposed drugs. The PPIs are made available on a webserver, called Schizo-Pi at http://severus.dbmi.pitt.edu/schizo-pi with advanced search capabilities.

The Computer Science, Biology, and Biomedical Informatics (CoSBBI) program was initiated in 2011 to expose the critical role of informatics in biomedicine to talented high school students.[1] By involving them in Science, Technology, Engineering, and Math (STEM) training at the high school level and providing mentorship and research opportunities throughout the formative years of their education, CoSBBI creates a research infrastructure designed to develop young informaticians. Our central premise is that the trajectory necessary to be an expert in the emerging fields of biomedical informatics and pathology informatics requires accelerated learning at an early age.In our 4(th) year of CoSBBI as a part of the University of Pittsburgh Cancer Institute (UPCI) Academy (http://www.upci.upmc.edu/summeracademy/), and our 2nd year of CoSBBI as an independent informatics-based academy, we enhanced our classroom curriculum, added hands-on computer science instruction, and expanded research projects to include clinical informatics. We also conducted a qualitative evaluation of the program to identify areas that need improvement in order to achieve our goal of creating a pipeline of exceptionally well-trained applicants for both the disciplines of pathology informatics and biomedical informatics in the era of big data and personalized medicine.

Background Protein-protein interaction (PPI) networks (interactomes) of most organisms, except for some model organisms, are largely unknown. Experimental methods including high-throughput techniques are highly resource intensive. Therefore, computational discovery of PPIs can accelerate biological discovery by presenting "most-promising" pairs of proteins that are likely to interact. For many bacteria, genome sequence, and thereby genomic context of proteomes, is readily available; additionally, for some of these proteomes, localization and functional annotations are also available, but interactomes are not available. We present here a method for rapid development of computational system to predict interactome of bacterial proteomes. While other studies have presented methods to transfer interologs across species, here, we propose transfer of computational models to benefit from cross-species annotations, thereby predicting many more novel interactions even in the absence of interologs. Mycobacterium tuberculosis (Mtb) and Clostridium difficile (CD) have been used to demonstrate the work. Results We developed a random forest classifier over features derived from Gene Ontology annotations and genetic context scores provided by STRING database for predicting Mtb and CD interactions independently. The Mtb classifier gave a precision of 94% and a recall of 23% on a held out test set. The Mtb model was then run on all the 8 million protein pairs of the Mtb proteome, resulting in 708 new interactions (at 94% expected precision) or 1,595 new interactions at 80% expected precision. The CD classifier gave a precision of 90% and a recall of 16% on a held out test set. The CD model was run on all the 8 million protein pairs of the CD proteome, resulting in 143 new interactions (at 90% expected precision) or 580 new interactions (at 80% expected precision). We also compared the overlap of predictions of our method with STRING database interactions for CD and Mtb and also with interactions identified recently by a bacterial 2-hybrid system for Mtb. To demonstrate the utility of transfer of computational models, we made use of the developed Mtb model and used it to predict CD protein-pairs. The cross species model thus developed yielded a precision of 88% at a recall of 8%. To demonstrate transfer of features from other organisms in the absence of feature-based and interaction-based information, we transferred missing feature values from Mtb orthologs into the CD data. In transferring this data from orthologs (not interologs), we showed that a large number of interactions can be predicted. Conclusions Rapid discovery of (partial) bacterial interactome can be made by using existing set of GO and STRING features associated with the organisms. We can make use of cross-species interactome development, when there are not even sufficient known interactions to develop a computational prediction system. Computational model of well-studied organism(s) can be employed to make the initial interactome prediction for the target organism. We have also demonstrated successfully, that annotations can be transferred from orthologs in well-studied organisms enabling accurate predictions for organisms with no annotations. These approaches can serve as building blocks to address the challenges associated with feature coverage, missing interactions towards rapid interactome discovery for bacterial organisms. Availability The predictions for all Mtb and CD proteins are made available at: http://severus.dbmi.pitt.edu/TB and http://severus.dbmi.pitt.edu/CD respectively for browsing as well as for download.

Network analysis is a common approach for the study of genetic view of diseases and biological pathways. Typically, when a set of genes are identified to be of interest in relation to a disease, say through a genome wide association study (GWAS) or a different gene expression study, these genes are typically analyzed in the context of their protein-protein interaction (PPI) networks. Further analysis is carried out to compute the enrichment of known pathways and disease-associations in the network. Having tools for such analysis at the fingertips of biologists without the requirement for computer programming or curation of data would accelerate the characterization of genes of interest. Currently available tools do not integrate network and enrichment analysis and their visualizations, and most of them present results in formats not most conducive to human cognition. We developed the tool Lens for Enrichment and Network Studies of human proteins (LENS) that performs network and pathway and diseases enrichment analyses on genes of interest to users. The tool creates a visualization of the network, provides easy to read statistics on network connectivity, and displays Venn diagrams with statistical significance values of the network's association with drugs, diseases, pathways, and GWASs. We used the tool to analyze gene sets related to craniofacial development, autism, and schizophrenia. LENS is a web-based tool that does not require and download or plugins to use. The tool is free and does not require login for use, and is available at http://severus.dbmi.pitt.edu/LENS .