DLG1 functions upstream of SDCCAG3 and IFT20 to control ciliary targeting of polycystin-2.

Polarized vesicular trafficking directs specific receptors and ion channels to cilia, but the underlying mechanisms are poorly understood. Here we describe a role for DLG1, a core component of the Scribble polarity complex, in regulating ciliary protein trafficking in kidney epithelial cells. Conditional knockout of Dlg1 in mouse kidney causes ciliary elongation and cystogenesis, and cell-based proximity labeling proteomics and fluorescence microscopy show alterations in the ciliary proteome upon loss of DLG1. Specifically, the retromer-associated protein SDCCAG3, IFT20, and polycystin-2 (PC2) are reduced in the cilia of DLG1-deficient cells compared to control cells. This phenotype is recapitulated in vivo and rescuable by re-expression of wild-type DLG1, but not a Congenital Anomalies of the Kidney and Urinary Tract (CAKUT)-associated DLG1 variant, p.T489R. Finally, biochemical approaches and Alpha Fold modelling suggest that SDCCAG3 and IFT20 form a complex that associates, at least indirectly, with DLG1. Our work identifies a key role for DLG1 in regulating ciliary protein composition and suggests that ciliary dysfunction of the p.T489R DLG1 variant may contribute to CAKUT.

WeSA: a web server for improving analysis of affinity proteomics data.

Protein-protein interaction experiments still yield many false positive interactions. The socioaffinity metric can distinguish true protein-protein interactions from noise based on available data. Here, we present WeSA (Weighted SocioAffinity), which considers large datasets of interaction proteomics data (IntAct, BioGRID, the BioPlex) to score human protein interactions and, in a statistically robust way, flag those (even from a single experiment) that are likely to be false positives. ROC analysis (using CORUM-PDB positives and Negatome negatives) shows that WeSA improves over other measures of interaction confidence. WeSA shows consistently good results over all datasets (up to: AUC = 0.93 and at best threshold: TPR = 0.84, FPR = 0.11, Precision = 0.98). WeSA is freely available without login (wesa.russelllab.org). Users can submit their own data or look for organized information on human protein interactions using the web server. Users can either retrieve available information for a list of proteins of interest or calculate scores for new experiments. The server outputs either pre-computed or updated WeSA scores for the input enriched with information from databases. The summary is presented as a table and a network-based visualization allowing the user to remove those nodes/edges that the method considers spurious.

Reduced Protein Stability of 11 Pathogenic Missense STXBP1/MUNC18-1 Variants and Improved Disease Prediction.

BACKGROUND: Pathogenic variants in STXBP1/MUNC18-1 cause severe encephalopathies that are among the most common in genetic neurodevelopmental disorders. Different molecular disease mechanisms have been proposed, and pathogenicity prediction is limited. In this study, we aimed to define a generalized disease concept for STXBP1-related disorders and improve prediction. METHODS: A cohort of 11 disease-associated and 5 neutral variants (detected in healthy individuals) were tested in 3 cell-free assays and in heterologous cells and primary neurons. Protein aggregation was tested using gel filtration and Triton X-100 insolubility. PRESR (predicting STXBP1-related disorder), a machine learning algorithm that uses both sequence- and 3-dimensional structure-based features, was developed to improve pathogenicity prediction using 231 known disease-associated variants and comparison to our experimental data. RESULTS: Disease-associated variants, but none of the neutral variants, produced reduced protein levels. Cell-free assays demonstrated directly that disease-associated variants have reduced thermostability, with most variants denaturing around body temperature. In addition, most disease-associated variants impaired SNARE-mediated membrane fusion in a reconstituted assay. Aggregation/insolubility was observed for none of the variants in vitro or in neurons. PRESR outperformed existing tools substantially: Matthews correlation coefficient = 0.71 versus <0.55. CONCLUSIONS: These data establish intrinsic protein instability as the generalizable, primary cause for STXBP1-related disorders and show that protein-specific ortholog and 3-dimensional information improve disease prediction. PRESR is a publicly available diagnostic tool.

Multi-omic and functional analysis for classification and treatment of sarcomas with FUS-TFCP2 or EWSR1-TFCP2 fusions.

Linking clinical multi-omics with mechanistic studies may improve the understanding of rare cancers. We leverage two precision oncology programs to investigate rhabdomyosarcoma with FUS/EWSR1-TFCP2 fusions, an orphan malignancy without effective therapies. All tumors exhibit outlier ALK expression, partly accompanied by intragenic deletions and aberrant splicing resulting in ALK variants that are oncogenic and sensitive to ALK inhibitors. Additionally, recurrent CKDN2A/MTAP co-deletions provide a rationale for PRMT5-targeted therapies. Functional studies show that FUS-TFCP2 blocks myogenic differentiation, induces transcription of ALK and truncated TERT, and inhibits DNA repair. Unlike other fusion-driven sarcomas, TFCP2-rearranged tumors exhibit genomic instability and signs of defective homologous recombination. DNA methylation profiling demonstrates a close relationship with undifferentiated sarcomas. In two patients, sarcoma was preceded by benign lesions carrying FUS-TFCP2, indicating stepwise sarcomagenesis. This study illustrates the potential of linking precision oncology with preclinical research to gain insight into the classification, pathogenesis, and therapeutic vulnerabilities of rare cancers.

DLG1 functions upstream of SDCCAG3 and IFT20 to control ciliary targeting of polycystin-2.

Polarized vesicular trafficking directs specific receptors and ion channels to cilia, but the underlying mechanisms are poorly understood. Here we describe a role for DLG1, a core component of the Scribble polarity complex, in regulating ciliary protein trafficking in kidney epithelial cells. Conditional knockout of Dlg1 in mouse kidney caused ciliary elongation and cystogenesis, and cell-based proximity labelling proteomics and fluorescence microscopy showed alterations in the ciliary proteome upon loss of DLG1. Specifically, the retromer-associated protein SDCCAG3, IFT20 and polycystin-2 (PC2) were reduced in cilia of DLG1 deficient cells compared to control cells. This phenotype was recapitulated in vivo and rescuable by re-expression of wildtype DLG1, but not a Congenital Anomalies of the Kidney and Urinary Tract (CAKUT)-associated DLG1 variant, p.T489R. Finally, biochemical approaches and Alpha Fold modelling suggested that SDCCAG3 and IFT20 form a complex that associates, at least indirectly, with DLG1. Our work identifies a key role for DLG1 in regulating ciliary protein composition and suggests that ciliary dysfunction of the p.T489R DLG1 variant may contribute to CAKUT.

De-Suppression of Mesenchymal Cell Identities and Variable Phenotypic Outcomes Associated with Knockout of Bbs1.

Bardet-Biedl syndrome (BBS) is an archetypal ciliopathy caused by dysfunction of primary cilia. BBS affects multiple tissues, including the kidney, eye and hypothalamic satiety response. Understanding pan-tissue mechanisms of pathogenesis versus those which are tissue-specific, as well as gauging their associated inter-individual variation owing to genetic background and stochastic processes, is of paramount importance in syndromology. The BBSome is a membrane-trafficking and intraflagellar transport (IFT) adaptor protein complex formed by eight BBS proteins, including BBS1, which is the most commonly mutated gene in BBS. To investigate disease pathogenesis, we generated a series of clonal renal collecting duct IMCD3 cell lines carrying defined biallelic nonsense or frameshift mutations in Bbs1, as well as a panel of matching wild-type CRISPR control clones. Using a phenotypic screen and an unbiased multi-omics approach, we note significant clonal variability for all assays, emphasising the importance of analysing panels of genetically defined clones. Our results suggest that BBS1 is required for the suppression of mesenchymal cell identities as the IMCD3 cell passage number increases. This was associated with a failure to express epithelial cell markers and tight junction formation, which was variable amongst clones. Transcriptomic analysis of hypothalamic preparations from BBS mutant mice, as well as BBS patient fibroblasts, suggested that dysregulation of epithelial-to-mesenchymal transition (EMT) genes is a general predisposing feature of BBS across tissues. Collectively, this work suggests that the dynamic stability of the BBSome is essential for the suppression of mesenchymal cell identities as epithelial cells differentiate.

Emerging Trends in the Field of Inflammation and Proteinopathy in ALS/FTD Spectrum Disorder.

Proteinopathy and neuroinflammation are two main hallmarks of neurodegenerative diseases. They also represent rare common events in an exceptionally broad landscape of genetic, environmental, neuropathologic, and clinical heterogeneity present in patients. Here, we aim to recount the emerging trends in amyotrophic lateral sclerosis (ALS) and frontotemporal degeneration (FTD) spectrum disorder. Our review will predominantly focus on neuroinflammation and systemic immune imbalance in ALS and FTD, which have recently been highlighted as novel therapeutic targets. A common mechanism of most ALS and ~50% of FTD patients is dysregulation of TAR DNA-binding protein 43 (TDP-43), an RNA/DNA-binding protein, which becomes depleted from the nucleus and forms cytoplasmic aggregates in neurons and glia. This, in turn, via both gain and loss of function events, alters a variety of TDP-43-mediated cellular events. Experimental attempts to target TDP-43 aggregates or manipulate crosstalk in the context of inflammation will be discussed. Targeting inflammation, and the immune system in general, is of particular interest because of the high plasticity of immune cells compared to neurons.

Genomic landscape of follicular lymphoma across a wide spectrum of clinical behaviors.

While some follicular lymphoma (FL) patients do not require treatment or experience prolonged responses, others relapse early, and little is known about genetic alterations specific to patients with a particular clinical behavior. We selected 56 grade 1-3A FL patients according to their need of treatment or timing of relapse: never treated (n = 7), non-relapsed (19), late relapse (14), early relapse or POD24 (11), and primary refractory (5). We analyzed 56 diagnostic and 12 paired relapse lymphoid tissue biopsies and performed copy number alteration (CNA) analysis and next generation sequencing (NGS). We identified six focal driver losses (1p36.32, 6p21.32, 6q14.1, 6q23.3, 9p21.3, 10q23.33) and 1p36.33 copy-neutral loss of heterozygosity (CN-LOH). By integrating CNA and NGS results, the most frequently altered genes/regions were KMT2D (79%), CREBBP (67%), TNFRSF14 (46%) and BCL2 (40%). Although we found that mutations in PIM1, FOXO1 and TMEM30A were associated with an adverse clinical behavior, definitive conclusions cannot be drawn, due to the small sample size. We identified common precursor cells harboring early oncogenic alterations of the KMT2D, CREBBP, TNFRSF14 and EP300 genes and 16p13.3-p13.2 CN-LOH. Finally, we established the functional consequences of mutations by means of protein modeling (CD79B, PLCG2, PIM1, MCL1 and IRF8). These data expand the knowledge on the genomics behind the heterogeneous FL population and, upon replication in larger cohorts, could contribute to risk stratification and the development of targeted therapies.

A targeted multi-proteomics approach generates a blueprint of the ciliary ubiquitinome.

Establishment and maintenance of the primary cilium as a signaling-competent organelle requires a high degree of fine tuning, which is at least in part achieved by a variety of post-translational modifications. One such modification is ubiquitination. The small and highly conserved ubiquitin protein possesses a unique versatility in regulating protein function via its ability to build mono and polyubiquitin chains onto target proteins. We aimed to take an unbiased approach to generate a comprehensive blueprint of the ciliary ubiquitinome by deploying a multi-proteomics approach using both ciliary-targeted ubiquitin affinity proteomics, as well as ubiquitin-binding domain-based proximity labelling in two different mammalian cell lines. This resulted in the identification of several key proteins involved in signaling, cytoskeletal remodeling and membrane and protein trafficking. Interestingly, using two different approaches in IMCD3 and RPE1 cells, respectively, we uncovered several novel mechanisms that regulate cilia function. In our IMCD3 proximity labeling cell line model, we found a highly enriched group of ESCRT-dependent clathrin-mediated endocytosis-related proteins, suggesting an important and novel role for this pathway in the regulation of ciliary homeostasis and function. In contrast, in RPE1 cells we found that several structural components of caveolae (CAV1, CAVIN1, and EHD2) were highly enriched in our cilia affinity proteomics screen. Consistently, the presence of caveolae at the ciliary pocket and ubiquitination of CAV1 specifically, were found likely to play a role in the regulation of ciliary length in these cells. Cilia length measurements demonstrated increased ciliary length in RPE1 cells stably expressing a ubiquitination impaired CAV1 mutant protein. Furthermore, live cell imaging in the same cells revealed decreased CAV1 protein turnover at the cilium as the possible cause for this phenotype. In conclusion, we have generated a comprehensive list of cilia-specific proteins that are subject to regulation via ubiquitination which can serve to further our understanding of cilia biology in health and disease.

Focal structural variants revealed by whole genome sequencing disrupt the histone demethylase KDM4C in B-cell lymphomas.

Histone methylation-modifiers, such as EZH2 and KMT2D, are recurrently altered in B-cell lymphomas. To comprehensively describe the landscape of alterations affecting genes encoding histone methylation-modifiers in lymphomagenesis we investigated whole genome and transcriptome data of 186 mature B-cell lymphomas sequenced in the ICGC MMML-Seq project. Besides confirming common alterations of KMT2D (47% of cases), EZH2 (17%), SETD1B (5%), PRDM9 (4%), KMT2C (4%), and SETD2 (4%), also identified by prior exome or RNA-sequencing studies, we here found recurrent alterations to KDM4C in chromosome 9p24, encoding a histone demethylase. Focal structural variation was the main mechanism of KDM4C alterations, and was independent from 9p24 amplification. We also identified KDM4C alterations in lymphoma cell lines including a focal homozygous deletion in a classical Hodgkin lymphoma cell line. By integrating RNA-sequencing and genome sequencing data we predict that KDM4C structural variants result in loss-offunction. By functional reconstitution studies in cell lines, we provide evidence that KDM4C can act as a tumor suppressor. Thus, we show that identification of structural variants in whole genome sequencing data adds to the comprehensive description of the mutational landscape of lymphomas and, moreover, establish KDM4C as a putative tumor suppressive gene recurrently altered in subsets of B-cell derived lymphomas.

The human signal peptidase complex acts as a quality control enzyme for membrane proteins.

Cells need to detect and degrade faulty membrane proteins to maintain homeostasis. In this study, we identify a previously unknown function of the human signal peptidase complex (SPC)-the enzyme that removes endoplasmic reticulum (ER) signal peptides-as a membrane protein quality control factor. We show that the SPC cleaves membrane proteins that fail to correctly fold or assemble into their native complexes at otherwise hidden cleavage sites, which our study reveals to be abundant in the human membrane proteome. This posttranslocational cleavage synergizes with ER-associated degradation to sustain membrane protein homeostasis and contributes to cellular fitness. Cryptic SPC cleavage sites thus serve as predetermined breaking points that, when exposed, help to target misfolded or surplus proteins for degradation, thereby maintaining a healthy membrane proteome.

Biochemically validated structural model of the 15-subunit intraflagellar transport complex IFT-B.

Cilia are ubiquitous eukaryotic organelles impotant for cellular motility, signaling, and sensory reception. Cilium formation requires intraflagellar transport of structural and signaling components and involves 22 different proteins organized into intraflagellar transport (IFT) complexes IFT-A and IFT-B that are transported by molecular motors. The IFT-B complex constitutes the backbone of polymeric IFT trains carrying cargo between the cilium and the cell body. Currently, high-resolution structures are only available for smaller IFT-B subcomplexes leaving > 50% structurally uncharacterized. Here, we used Alphafold to structurally model the 15-subunit IFT-B complex. The model was validated using cross-linking/mass-spectrometry data on reconstituted IFT-B complexes, X-ray scattering in solution, diffraction from crystals as well as site-directed mutagenesis and protein-binding assays. The IFT-B structure reveals an elongated and highly flexible complex consistent with cryo-electron tomographic reconstructions of IFT trains. The IFT-B complex organizes into IFT-B1 and IFT-B2 parts with binding sites for ciliary cargo and the inactive IFT dynein motor, respectively. Interestingly, our results are consistent with two different binding sites for IFT81/74 on IFT88/70/52/46 suggesting the possibility of different structural architectures for the IFT-B1 complex. Our data present a structural framework to understand IFT-B complex assembly, function, and ciliopathy variants.

Never-homozygous genetic variants in healthy populations are potential recessive disease candidates.

The rapid pace with which genetic variants are now being determined means there is a pressing need to understand how they affect biological systems. Variants from healthy individuals have previously been used to study blood groups or HLA diversity and to identify genes that can apparently be nonfunctional in healthy people. These studies and others have observed a lower than expected frequency of homozygous individuals for potentially deleterious alleles, which would suggest that several of these alleles can lead to recessive disorders. Here we exploited this principle to hunt for potential disease variants in genomes from healthy people. We identified at least 108 exclusively heterozygous variants with evidence for an impact on biological function. We discuss several examples of candidate variants/genes including CCDC8, PANK3, RHD and NLRP12. Overall, the results suggest there are many, comparatively frequent, potentially lethal or disease-causing variants lurking in healthy human populations.

The S-palmitoylome and DHHC-PAT interactome of Drosophila melanogaster S2R+ cells indicate a high degree of conservation to mammalian palmitoylomes.

Protein S-palmitoylation, the addition of a long-chain fatty acid to target proteins, is among the most frequent reversible protein modifications in Metazoa, affecting subcellular protein localization, trafficking and protein-protein interactions. S-palmitoylated proteins are abundant in the neuronal system and are associated with neuronal diseases and cancer. Despite the importance of this post-translational modification, it has not been thoroughly studied in the model organism Drosophila melanogaster. Here we present the palmitoylome of Drosophila S2R+ cells, comprising 198 proteins, an estimated 3.5% of expressed genes in these cells. Comparison of orthologs between mammals and Drosophila suggests that S-palmitoylated proteins are more conserved between these distant phyla than non-S-palmitoylated proteins. To identify putative client proteins and interaction partners of the DHHC family of protein acyl-transferases (PATs) we established DHHC-BioID, a proximity biotinylation-based method. In S2R+ cells, ectopic expression of the DHHC-PAT dHip14-BioID in combination with Snap24 or an interaction-deficient Snap24-mutant as a negative control, resulted in biotinylation of Snap24 but not the Snap24-mutant. DHHC-BioID in S2R+ cells using 10 different DHHC-PATs as bait identified 520 putative DHHC-PAT interaction partners of which 48 were S-palmitoylated and are therefore putative DHHC-PAT client proteins. Comparison of putative client protein/DHHC-PAT combinations indicates that CG8314, CG5196, CG5880 and Patsas have a preference for transmembrane proteins, while S-palmitoylated proteins with the Hip14-interaction motif are most enriched by DHHC-BioID variants of approximated and dHip14. Finally, we show that BioID is active in larval and adult Drosophila and that dHip14-BioID rescues dHip14 mutant flies, indicating that DHHC-BioID is non-toxic. In summary we provide the first systematic analysis of a Drosophila palmitoylome. We show that DHHC-BioID is sensitive and specific enough to identify DHHC-PAT client proteins and provide DHHC-PAT assignment for ca. 25% of the S2R+ cell palmitoylome, providing a valuable resource. In addition, we establish DHHC-BioID as a useful concept for the identification of tissue-specific DHHC-PAT interactomes in Drosophila.

PRECOGx: exploring GPCR signaling mechanisms with deep protein representations.

In this study we show that protein language models can encode structural and functional information of GPCR sequences that can be used to predict their signaling and functional repertoire. We used the ESM1b protein embeddings as features and the binding information known from publicly available studies to develop PRECOGx, a machine learning predictor to explore GPCR interactions with G protein and ß-arrestin, which we made available through a new webserver (https://precogx.bioinfolab.sns.it/). PRECOGx outperformed its predecessor (e.g. PRECOG) in predicting GPCR-transducer couplings, being also able to consider all GPCR classes. The webserver also provides new functionalities, such as the projection of input sequences on a low-dimensional space describing essential features of the human GPCRome, which is used as a reference to track GPCR variants. Additionally, it allows inspection of the sequence and structural determinants responsible for coupling via the analysis of the most important attention maps used by the models as well as through predicted intramolecular contacts. We demonstrate applications of PRECOGx by predicting the impact of disease variants (ClinVar) and alternative splice forms from healthy tissues (GTEX) of human GPCRs, revealing the power to dissect system biasing mechanisms in both health and disease.

Next Generation Protein Structure Predictions and Genetic Variant Interpretation.

The need to make sense of the thousands of genetic variants uncovered every day in terms of pathology or biological mechanism is acute. Many insights into how genetic changes impact protein function can be gleaned if three-dimensional structures of the associated proteins are available. The availability of a highly accurate method of predicting structures from amino acid sequences (e.g. Alphafold2) is thus potentially a great boost to those wanting to understand genetic changes. In this paper we discuss the current state of protein structures known for the human and other proteomes and how Alphafold2 might impact on variant interpretation efforts. For the human proteome in particular, the state of the available structural data suggests that the impact on variant interpretation might be less than anticipated. We also discuss additional efforts in structure prediction that could further aid the understanding of genetic variants.

Mechnetor: a web server for exploring protein mechanism and the functional context of genetic variants.

Advances in DNA sequencing and proteomics mean that researchers must now regularly interrogate thousands of positional gene/protein changes in order to find those relevant for potential clinical application or biological insights. The abundance of already known information on protein interactions, mechanism, and tertiary structure provides the possible means to understand these changes rapidly, though a careful and systematic integration of these diverse datasets is first needed. For this purpose, we developed Mechnetor, a tool that allows users to quickly explore and visualize integrated mechanistic data for proteins or interactions of interest. Central to the system is a careful cataloguing of diverse sources of protein interaction mechanism, and an efficient means to visualize interactions between relevant and/or known protein regions. The result is a finer resolution interaction network that provides more immediate clues as to points of intervention or mechanistic understanding. Users can import protein, interactions, genetic variants or post-translational modifications and see these data in the best known mechanistic context. We demonstrate the tool with topical examples in human genetic diseases and cancer genomics. The tool is freely available at: mechnetor.russelllab.org.

RhoGAP19D inhibits Cdc42 laterally to control epithelial cell shape and prevent invasion.

Cdc42-GTP is required for apical domain formation in epithelial cells, where it recruits and activates the Par-6-aPKC polarity complex, but how the activity of Cdc42 itself is restricted apically is unclear. We used sequence analysis and 3D structural modeling to determine which Drosophila GTPase-activating proteins (GAPs) are likely to interact with Cdc42 and identified RhoGAP19D as the only high-probability Cdc42GAP required for polarity in the follicular epithelium. RhoGAP19D is recruited by α-catenin to lateral E-cadherin adhesion complexes, resulting in exclusion of active Cdc42 from the lateral domain. rhogap19d mutants therefore lead to lateral Cdc42 activity, which expands the apical domain through increased Par-6/aPKC activity and stimulates lateral contractility through the myosin light chain kinase, Genghis khan (MRCK). This causes buckling of the epithelium and invasion into the adjacent tissue, a phenotype resembling that of precancerous breast lesions. Thus, RhoGAP19D couples lateral cadherin adhesion to the apical localization of active Cdc42, thereby suppressing epithelial invasion.

Leveraging biochemical reactions to unravel functional impacts of cancer somatic variants affecting protein interaction interfaces.

Background: Considering protein mutations in their biological context is essential for understanding their functional impact, interpretation of high-dimensional datasets and development of effective targeted therapies in personalized medicine. Methods: We combined the curated knowledge of biochemical reactions from Reactome with the analysis of interaction-mediating 3D interfaces from Mechismo. In addition, we provided a software tool for users to explore and browse the analysis results in a multi-scale perspective starting from pathways and reactions to protein-protein interactions and protein 3D structures. Results: We analyzed somatic mutations from TCGA, revealing several significantly impacted reactions and pathways in specific cancer types. We found examples of genes not yet listed as oncodrivers, whose rare mutations were predicted to affect cancer processes similarly to known oncodrivers. Some identified processes lack any known oncodrivers, which suggests potentially new cancer-related processes (e.g. complement cascade reactions). Furthermore, we found that mutations perturbing certain processes are significantly associated with distinct phenotypes (i.e. survival time) in specific cancer types (e.g. PIK3CA centered pathways in LGG and UCEC cancer types), suggesting the translational potential of our approach for patient stratification. Our analysis also uncovered several druggable processes (e.g. GPCR signalling pathways) containing enriched reactions, providing support for new off-label therapeutic options. Conclusions: In summary, we have established a multi-scale approach to study genetic variants based on protein-protein interaction 3D structures. Our approach is different from previously published studies in its focus on biochemical reactions and can be applied to other data types (e.g. post-translational modifications) collected for many types of disease.