Matrix stiffening induces a pathogenic QKI-miR-7-SRSF1 signaling axis in pulmonary arterial endothelial cells.

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.

Systems Analysis of the Human Pulmonary Arterial Hypertension Lung Transcriptome.

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.

Factors Associated with Heritable Pulmonary Arterial Hypertension Exert Convergent Actions on the miR-130/301-Vascular Matrix Feedback Loop.

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γ-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.

Single Nucleotide Polymorphism rs9277336 Controls the Nuclear Alpha Actinin 4-Human Leukocyte Antigen-DPA1 Axis and Pulmonary Endothelial Pathophenotypes in Pulmonary Arterial Hypertension.

Background Pulmonary arterial hypertension (PAH) is a complex, fatal disease where disease severity has been associated with the single nucleotide polymorphism (SNP) rs2856830, located near the human leukocyte antigen DPA1 (HLA-DPA1) gene. We aimed to define the genetic architecture of functional variants associated with PAH disease severity by identifying allele-specific binding transcription factors and downstream targets that control endothelial pathophenotypes and PAH. Methods and Results Electrophoretic mobility shift assays of oligonucleotides containing SNP rs2856830 and 8 SNPs in linkage disequilibrium revealed functional SNPs via allele-imbalanced binding to human pulmonary arterial endothelial cell nuclear proteins. DNA pulldown proteomics identified SNP-binding proteins. SNP genotyping and clinical correlation analysis were performed in 84 patients with PAH at University of Pittsburgh Medical Center and in 679 patients with PAH in the All of Us database. SNP rs9277336 was identified as a functional SNP in linkage disequilibrium (r2>0.8) defined by rs2856830, and the minor allele was associated with decreased hospitalizations and improved cardiac output in patients with PAH, an index of disease severity. SNP pulldown proteomics showed allele-specific binding of nuclear ACTN4 (alpha actinin 4) protein to rs9277336 minor allele. Both ACTN4 and HLA-DPA1 were downregulated in pulmonary endothelium in human patients and rodent models of PAH. Via transcriptomic and phenotypic analyses, knockdown of HLA-DPA1 phenocopied knockdown of ACTN4, both similarly controlling cell structure pathways, immune pathways, and endothelial dysfunction. Conclusions We defined the pathogenic activity of functional SNP rs9277336, entailing the allele-specific binding of ACTN4 and controlling expression of the neighboring HLA-DPA1 gene. Through inflammatory or genetic means, downregulation of this ACTN4-HLA-DPA1 regulatory axis promotes endothelial pathophenotypes, providing a mechanistic explanation for the association between this SNP and PAH outcomes.

Improved hospitalization rates in a specialty center for heart failure with preserved ejection fraction and pulmonary hypertension.

Heart failure with preserved ejection fraction can be complicated by pulmonary hypertension. We designed a retrospective study to provide supporting evidence for referral to specialty care centers. Specialty care centers improved hospitalizations but not mortality-in part due to more aggressive medication management and guideline-directed monitoring.

miRNA/mRNA co-profiling identifies the miR-200 family as a central regulator of SMC quiescence.

miRNAs are versatile regulators of smooth muscle cell (SMC) fate and behavior in vascular development and disease. Targeted loss-of-function studies have established the relevance of specific miRNAs in controlling SMC differentiation or mediating phenotypic modulation. Our goal was to characterize SMC miRNAome and its contribution to transcriptome changes during phenotypic modulation. Small RNA sequencing revealed that dedifferentiation led to the differential expression of over 50 miRNAs in cultured SMC. miRNA/mRNA comparison predicted that over a third of SMC transcript expression was regulated by differentially expressed miRNAs. Our screen identified the miR-200 cluster as highly downregulated during dedifferentiation. miR-200 maintains SMC quiescence and represses proliferation, migration, and neointima formation, in part by targeting Quaking, a central SMC phenotypic switching mediator. Our study unraveled the substantial contribution of miRNAs in regulating the SMC transcriptome and identified the miR-200 cluster as a pro-quiescence mechanism and a potential inhibitor of vascular restenosis.

Single-cell RNA sequencing profiling of mouse endothelial cells in response to pulmonary arterial hypertension.

AIMS: Endothelial 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 RESULTS: We 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. CONCLUSION: This study uncovers PAH-induced EC transcriptomic changes at a high resolution, revealing novel targets for potential therapeutic candidate development.

Increased Mortality in Patients With Preoperative and Persistent Postoperative Pulmonary Hypertension Undergoing Mitral Valve Surgery for Mitral Regurgitation: A Cohort Study.

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.

Frataxin deficiency promotes endothelial senescence in pulmonary hypertension.

The dynamic regulation of endothelial pathophenotypes in pulmonary hypertension (PH) remains undefined. Cellular senescence is linked to PH with intracardiac shunts; however, its regulation across PH subtypes is unknown. Since endothelial deficiency of iron-sulfur (Fe-S) clusters is pathogenic in PH, we hypothesized that a Fe-S biogenesis protein, frataxin (FXN), controls endothelial senescence. An endothelial subpopulation in rodent and patient lungs across PH subtypes exhibited reduced FXN and elevated senescence. In vitro, hypoxic and inflammatory FXN deficiency abrogated activity of endothelial Fe-S-containing polymerases, promoting replication stress, DNA damage response, and senescence. This was also observed in stem cell-derived endothelial cells from Friedreich's ataxia (FRDA), a genetic disease of FXN deficiency, ataxia, and cardiomyopathy, often with PH. In vivo, FXN deficiency-dependent senescence drove vessel inflammation, remodeling, and PH, whereas pharmacologic removal of senescent cells in Fxn-deficient rodents ameliorated PH. These data offer a model of endothelial biology in PH, where FXN deficiency generates a senescent endothelial subpopulation, promoting vascular inflammatory and proliferative signals in other cells to drive disease. These findings also establish an endothelial etiology for PH in FRDA and left heart disease and support therapeutic development of senolytic drugs, reversing effects of Fe-S deficiency across PH subtypes.

Computational repurposing of therapeutic small molecules from cancer to 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.

Outcomes of Pulmonary Arterial Hypertension Are Improved in a Specialty Care Center.

BACKGROUND: Pulmonary arterial hypertension (PAH) is characterized by elevated pulmonary arterial pressures and is managed by vasodilator therapies. Current guidelines encourage PAH management in specialty care centers (SCCs), but evidence is sparse regarding improvement in clinical outcomes and correlation to vasodilator use with referral. RESEARCH QUESTION: Is PAH management at SCCs associated with improved clinical outcomes? STUDY DESIGNAND METHODS: A single-center, retrospective study was performed at the University of Pittsburgh Medical Center (UPMC; overseeing 40 hospitals). Patients with PAH were identified between 2008 and 2018 and classified into an SCC or non-SCC cohort. Cox proportional hazard modeling was done to compare for all-cause mortality, as was negative binomial regression modeling for hospitalizations. Vasodilator therapy was included to adjust outcomes. RESULTS: Of 580 patients with PAH at UPMC, 455 (78%) were treated at the SCC, comprising a younger (58.8 vs 64.8 years; P < .001) and more often female (68.4% vs 51.2%; P < .001) population with more comorbidities without differences in race or income. SCC patients demonstrated improved survival (hazard ratio, 0.68; P = .012) and fewer hospitalizations (incidence ratio, 0.54; P < .001), and provided more frequent disease monitoring. Early patient referral to SCC (< 6 months from time of diagnosis) was associated with improved outcomes compared with non-SCC patients. SCC patients were more frequently prescribed vasodilators (P < .001) and carried more diagnostic PAH coding (P < .001). Vasodilators were associated with improved outcomes irrespective of location but without statistical significance when comparing between locations (P > .05). INTERPRETATION: The UPMC SCC demonstrated improved outcomes in mortality and hospitalizations. The SCC benefit was multifactorial, with more frequent vasodilator therapy and disease monitoring. These findings provide robust evidence for early and regular referral of patients with PAH to SCCs.

SCUBE1 Controls BMPR2-Relevant Pulmonary Endothelial Function: Implications for Diagnostic Marker Development in Pulmonary Arterial Hypertension.

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.

Atrial arrhythmias are associated with increased mortality in pulmonary arterial hypertension.

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 < 0.001). This finding was consistent with multivariable adjustment of hemodynamic, cardiac structural, and heart rate indices as well as in sensitivity analyses of patients with paroxysmal versus non-paroxysmal arrhythmias. In these PAH/CTEPH patients, presence of AF/AFL significantly increased mortality risk. Mortality remained elevated in the absence of a high burden of uncontrolled or persistent arrhythmias, thus suggesting additional etiologies beyond rapid heart rate as an explanation. Future studies are warranted to confirm this observation and interrogate whether other therapies beyond rate and rhythm control are necessary to mitigate this risk.

Rapamycin-induced miR-21 promotes mitochondrial homeostasis and adaptation in mTORC1 activated cells.

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<0.01) and ∼75% (p<0.0001), respectively, and combined rapamycin treatment decreased soft agar growth by ∼90% (p<0.0001). miR-21 inhibition also increased sensitivity to apoptosis. Through a network biology-driven integration of RNAseq data, we discovered that miR-21 promotes mitochondrial adaptation and homeostasis in Tsc2-deficient cells. miR-21 inhibition reduced mitochondrial polarization and function in Tsc2-deficient cells, with and without co-treatment with rapamycin. Importantly, miR-21 inhibition limited Tsc2-deficient tumor growth in vivo, reducing tumor size by approximately 3-fold (p<0.0001). When combined with rapamcyin, miR-21 inhibition showed even more striking efficacy, both during treatment and after treatment cessation, with a 4-fold increase in median survival following rapamycin cessation (p=0.0008). We conclude that miR-21 promotes mTORC1-driven tumorigenesis via a mechanism that involves the mitochondria, and that miR-21 is a potential therapeutic target for TSC-associated hamartomas and other mTORC1-driven tumors, with the potential for synergistic efficacy when combined with rapalogs.

p62/SQSTM1 Cooperates with Hyperactive mTORC1 to Regulate Glutathione Production, Maintain Mitochondrial Integrity, and Promote Tumorigenesis.

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 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 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.

Is the average shortest path length of gene set a reflection of their biological relatedness?

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.

How can we improve Science, Technology, Engineering, and Math education to encourage careers in Biomedical and Pathology Informatics?

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.

Schizophrenia interactome with 504 novel protein-protein interactions.

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.

LENS: web-based lens for enrichment and network studies of human proteins.

BACKGROUND: 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. RESULTS: 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. CONCLUSION: 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.