ABSTRACT
AIMS: BMP9 is a high affinity ligand of ALK1 and endoglin receptors that are mutated in the rare genetic vascular disorder hereditary hemorrhagic telangiectasia (HHT). We have previously shown that loss of Bmp9 in the 129/Ola genetic background leads to spontaneous liver fibrosis via capillarization of liver sinusoidal endothelial cells (LSEC) and kidney lesions. We aimed to decipher the molecular mechanisms downstream of BMP9 to better characterize its role in vascular homeostasis in different organs. METHODS AND RESULTS: For this, we performed an RNA-seq analysis on LSEC from adult WT and Bmp9-KO mice and identified over 2000 differentially expressed genes. Gene ontology analysis showed that Bmp9 deletion led to a decrease in BMP and Notch signalling, but also LSEC capillary identity while increasing their cell cycle. The gene ontology term 'glomerulus development' was also negatively enriched in Bmp9-KO mice vs. WT supporting a role for BMP9 in kidney vascularization. Through different imaging approaches (electron microscopy, immunostainings), we found that loss of Bmp9 led to vascular enlargement of the glomeruli capillaries associated with alteration of podocytes. Importantly, we also showed for the first time that the loss of Bmp9 led to spontaneous arteriovenous malformations (AVMs) in the liver, gastrointestinal tract, and uterus. CONCLUSION: Altogether, these results demonstrate that BMP9 plays an important role in vascular quiescence both locally in the liver by regulating endothelial capillary differentiation markers and cell cycle but also at distance in many organs via its presence in the circulation. It also reveals that loss of Bmp9 is sufficient to induce spontaneous AVMs, supporting a key role for BMP9 in the pathogenesis of HHT.
Subject(s)
Arteriovenous Malformations , Endothelial Cells , Growth Differentiation Factor 2 , Mice, Knockout , Signal Transduction , Animals , Growth Differentiation Factor 2/metabolism , Growth Differentiation Factor 2/genetics , Endothelial Cells/metabolism , Endothelial Cells/pathology , Arteriovenous Malformations/metabolism , Arteriovenous Malformations/genetics , Arteriovenous Malformations/pathology , Disease Models, Animal , Mice, 129 Strain , Liver/metabolism , Liver/pathology , Liver/blood supply , Phenotype , RNA-Seq , Receptors, Notch/metabolism , Receptors, Notch/genetics , MaleABSTRACT
BACKGROUND: BMP9 and BMP10 are two major regulators of vascular homeostasis. These two ligands bind with high affinity to the endothelial type I kinase receptor ALK1, together with a type II receptor, leading to the direct phosphorylation of the SMAD transcription factors. Apart from this canonical pathway, little is known. Interestingly, mutations in this signaling pathway have been identified in two rare cardiovascular diseases, hereditary hemorrhagic telangiectasia and pulmonary arterial hypertension. METHODS: To get an overview of the signaling pathways modulated by BMP9 and BMP10 stimulation in endothelial cells, we employed an unbiased phosphoproteomic-based strategy. Identified phosphosites were validated by western blot analysis and regulated targets by RT-qPCR. Cell cycle analysis was analyzed by flow cytometry. RESULTS: Large-scale phosphoproteomics revealed that BMP9 and BMP10 treatment induced a very similar phosphoproteomic profile. These BMPs activated a non-canonical transcriptional SMAD-dependent MAPK pathway (MEKK4/P38). We were able to validate this signaling pathway and demonstrated that this activation required the expression of the protein GADD45ß. In turn, activated P38 phosphorylated the heat shock protein HSP27 and the endocytosis protein Eps15 (EGF receptor pathway substrate), and regulated the expression of specific genes (E-selectin, hyaluronan synthase 2 and cyclooxygenase 2). This study also highlighted the modulation in phosphorylation of proteins involved in transcriptional regulation (phosphorylation of the endothelial transcription factor ERG) and cell cycle inhibition (CDK4/6 pathway). Accordingly, we found that BMP10 induced a G1 cell cycle arrest and inhibited the mRNA expression of E2F2, cyclinD1 and cyclinA1. CONCLUSIONS: Overall, our phosphoproteomic screen identified numerous proteins whose phosphorylation state is impacted by BMP9 and BMP10 treatment, paving the way for a better understanding of the molecular mechanisms regulated by BMP signaling in vascular diseases.
Subject(s)
Bone Morphogenetic Proteins , Endothelial Cells , Cell Cycle Checkpoints , Phosphorylation , G1 Phase Cell Cycle CheckpointsABSTRACT
Clear-cell renal cell carcinoma (ccRCC) accounts for 75% of kidney cancers. Due to the high recurrence rate and treatment options that come with high costs and potential side effects, a correct prognosis of patient survival is essential for the successful and effective treatment of patients. Novel biomarkers could play an important role in the assessment of the overall survival of patients. COL7A1 encodes for collagen type VII, a constituent of the basal membrane. COL7A1 is associated with survival in many cancers; however, the prognostic value of COL7A1 expression as a standalone biomarker in ccRCC has not been investigated. With five publicly available independent cohorts, we used Kaplan-Meier curves and the Cox proportional hazards model to investigate the prognostic value of COL7A1, as well as gene set enrichment analysis to investigate genes co-expressed with COL7A1. COL7A1 expression stratifies patients in terms of aggressiveness, where the 5-year survival probability of each of the four groups was 72.4%, 59.1%, 34.15%, and 8.6% in order of increasing expression. Additionally, COL7A1 expression was successfully used to further divide patients of each stage and histological grade into groups of high and low risk. Similar results were obtained in independent cohorts. In vitro knockdown of COL7A1 expression significantly affected ccRCC cells' ability to migrate, leading to the hypothesis that COL7A1 may have a role in cancer aggressiveness. To conclude, we identified COL7A1 as a new prognosis marker that can stratify ccRCC patients.
ABSTRACT
(1) Background: tumor profiling enables patient survival prediction. The two essential parameters to be calibrated when designing a study based on tumor profiles from a cohort are the sequencing depth of RNA-seq technology and the number of patients. This calibration is carried out under cost constraints, and a compromise has to be found. In the context of survival data, the goal of this work is to benchmark the impact of the number of patients and of the sequencing depth of miRNA-seq and mRNA-seq on the predictive capabilities for both the Cox model with elastic net penalty and random survival forest. (2) Results: we first show that the Cox model and random survival forest provide comparable prediction capabilities, with significant differences for some cancers. Second, we demonstrate that miRNA and/or mRNA data improve prediction over clinical data alone. mRNA-seq data leads to slightly better prediction than miRNA-seq, with the notable exception of lung adenocarcinoma for which the tumor miRNA profile shows higher predictive power. Third, we demonstrate that the sequencing depth of RNA-seq data can be reduced for most of the investigated cancers without degrading the prediction abilities, allowing the creation of independent validation sets at a lower cost. Finally, we show that the number of patients in the training dataset can be reduced for the Cox model and random survival forest, allowing the use of different models on different patient subgroups.
Subject(s)
Lung Neoplasms , MicroRNAs , Humans , Proportional Hazards Models , Random Forest , Gene Expression Profiling , MicroRNAs/genetics , Lung Neoplasms/genetics , RNA, Messenger/geneticsABSTRACT
BACKGROUND: Prediction of patient survival from tumor molecular '-omics' data is a key step toward personalized medicine. Cox models performed on RNA profiling datasets are popular for clinical outcome predictions. But these models are applied in the context of "high dimension", as the number p of covariates (gene expressions) greatly exceeds the number n of patients and e of events. Thus, pre-screening together with penalization methods are widely used for dimensional reduction. METHODS: In the present paper, (i) we benchmark the performance of the lasso penalization and three variants (i.e., ridge, elastic net, adaptive elastic net) on 16 cancers from TCGA after pre-screening, (ii) we propose a bi-dimensional pre-screening procedure based on both gene variability and p-values from single variable Cox models to predict survival, and (iii) we compare our results with iterative sure independence screening (ISIS). RESULTS: First, we show that integration of mRNA-seq data with clinical data improves predictions over clinical data alone. Second, our bi-dimensional pre-screening procedure can only improve, in moderation, the C-index and/or the integrated Brier score, while excluding irrelevant genes for prediction. We demonstrate that the different penalization methods reached comparable prediction performances, with slight differences among datasets. Finally, we provide advice in the case of multi-omics data integration. CONCLUSIONS: Tumor profiles convey more prognostic information than clinical variables such as stage for many cancer subtypes. Lasso and Ridge penalizations perform similarly than Elastic Net penalizations for Cox models in high-dimension. Pre-screening of the top 200 genes in term of single variable Cox model p-values is a practical way to reduce dimension, which may be particularly useful when integrating multi-omics.