Faculty Bibliography

Advances in mass-spectrometry have generated increasingly large-scale proteomics datasets containing tens of thousands of phosphorylation sites (phosphosites) that require prioritization. We develop a bioinformatics tool called HotPho and systematically discover 3D co-clustering of phosphosites and cancer mutations on protein structures. HotPho identifies 474 such hybrid clusters containing 1255 co-clustering phosphosites, including RET p.S904/Y928, the conserved HRAS/KRAS p.Y96, and IDH1 p.Y139/IDH2 p.Y179 that are adjacent to recurrent mutations on protein structures not found by linear proximity approaches. Hybrid clusters, enriched in histone and kinase domains, frequently include expression-associated mutations experimentally shown as activating and conferring genetic dependency. Approximately 300 co-clustering phosphosites are verified in patient samples of 5 cancer types or previously implicated in cancer, including CTNNB1 p.S29/Y30, EGFR p.S720, MAPK1 p.S142, and PTPN12 p.S275. In summary, systematic 3D clustering analysis highlights nearly 3,000 likely functional mutations and over 1000 cancer phosphosites for downstream investigation and evaluation of potential clinical relevance.

D-type cyclins are central regulators of the cell division cycle and are among the most frequently deregulated therapeutic targets in human cancer¹, but the mechanisms that regulate their turnover are still being debated^2,3. Here, by combining biochemical and genetics studies in somatic cells, we identify CRL4^AMBRA1 (also known as CRL4^DCAF3) as the ubiquitin ligase that targets all three D-type cyclins for degradation. During development, loss of Ambra1 induces the accumulation of D-type cyclins and retinoblastoma (RB) hyperphosphorylation and hyperproliferation, and results in defects of the nervous system that are reduced by treating pregnant mice with the FDA-approved CDK4 and CDK6 (CDK4/6) inhibitor abemaciclib. Moreover, AMBRA1 acts as a tumour suppressor in mouse models and low AMBRA1 mRNA levels are predictive of poor survival in cancer patients. Cancer hotspot mutations in D-type cyclins abrogate their binding to AMBRA1 and induce their stabilization. Finally, a whole-genome, CRISPR-Cas9 screen identified AMBRA1 as a regulator of the response to CDK4/6 inhibition. Loss of AMBRA1 reduces sensitivity to CDK4/6 inhibitors by promoting the formation of complexes of D-type cyclins with CDK2. Collectively, our results reveal the molecular mechanism that controls the stability of D-type cyclins during cell-cycle progression, in development and in human cancer, and implicate AMBRA1 as a critical regulator of the RB pathway.

Recent studies have revealed diverse amino acid, post-translational, and noncanonical modifications of proteins in diverse organisms and tissues. However, their unbiased detection and analysis remain hindered by technical limitations. Here, we present a spectral alignment method for the identification of protein modifications using high-resolution mass spectrometry proteomics. Termed SAMPEI for spectral alignment-based modified peptide identification, this open-source algorithm is designed for the discovery of functional protein and peptide signaling modifications, without prior knowledge of their identities. Using synthetic standards and controlled chemical labeling experiments, we demonstrate its high specificity and sensitivity for the discovery of substoichiometric protein modifications in complex cellular extracts. SAMPEI mapping of mouse macrophage differentiation revealed diverse post-translational protein modifications, including distinct forms of cysteine itaconatylation. SAMPEI's robust parametrization and versatility are expected to facilitate the discovery of biological modifications of diverse macromolecules. SAMPEI is implemented as a Python package and is available open-source from BioConda and GitHub (https://github.com/FenyoLab/SAMPEI).

BACKGROUND:Long INterspersed Element-1 (LINE-1) is an autonomous retroelement able to "copy-and-paste" itself into new loci of the host genome through a process called retrotransposition. The LINE-1 bicistronic mRNA codes for two proteins, ORF1p, a nucleic acid chaperone, and ORF2p, a protein with endonuclease and reverse transcriptase activity. Both proteins bind LINE-1 mRNA in cis and are necessary for retrotransposition. While LINE-1 transcription is usually repressed in most healthy somatic cells through a plethora of mechanisms, ORF1p expression has been observed in nearly 50% of tumors, and new LINE-1 insertions have been documented in a similar fraction of tumors, including prostate cancer. RESULTS:Here, we utilized RNA ImmunoPrecipitation (RIP) and the L1EM analysis software to identify ORF1p bound RNA in prostate cancer cells. We identified LINE-1 loci that were expressed in parental androgen sensitive and androgen independent clonal derivatives. In all androgen independent cells, we found higher levels of LINE-1 RNA, as well as unique expression patterns of LINE-1 loci. Interestingly, we observed that ORF1p bound many non-LINE-1 mRNA in all prostate cancer cell lines evaluated, and polyA RNA, and RNA localized in p-bodies were especially enriched. Furthermore, the expression levels of RNAs identified in our ORF1p RIP correlated with RNAs expressed in LINE-1 positive tumors from The Cancer Genome Atlas (TCGA). CONCLUSION/CONCLUSIONS:Our results show a significant remodeling of LINE-1 loci expression in androgen independent cell lines when compared to parental androgen dependent cells. Additionally, we found that ORF1p bound a significant amount of non-LINE-1 mRNA, and that the enriched ORF1p bound mRNAs are also amplified in LINE-1 expressing TCGA prostate tumors, indicating the biological relevance of our findings to prostate cancer.

Glioblastoma (GBM) is the most aggressive nervous system cancer. Understanding its molecular pathogenesis is crucial to improving diagnosis and treatment. Integrated analysis of genomic, proteomic, post-translational modification and metabolomic data on 99 treatment-naive GBMs provides insights to GBM biology. We identify key phosphorylation events (e.g., phosphorylated PTPN11 and PLCG1) as potential switches mediating oncogenic pathway activation, as well as potential targets for EGFR-, TP53-, and RB1-altered tumors. Immune subtypes with distinct immune cell types are discovered using bulk omics methodologies, validated by snRNA-seq, and correlated with specific expression and histone acetylation patterns. Histone H2B acetylation in classical-like and immune-low GBM is driven largely by BRDs, CREBBP, and EP300. Integrated metabolomic and proteomic data identify specific lipid distributions across subtypes and distinct global metabolic changes in IDH-mutated tumors. This work highlights biological relationships that could contribute to stratification of GBM patients for more effective treatment.

The ability to monitor DNA replication fork directionality at the genome-wide scale is paramount for a greater understanding of how genetic and environmental perturbations can impact replication dynamics in human cells. Here we describe a detailed protocol for isolating and sequencing Okazaki fragments from asynchronously growing mammalian cells, termed Okazaki fragment sequencing (Ok-seq), for the purpose of quantitatively determining replication initiation and termination frequencies around specific genomic loci by meta-analyses. Briefly, cells are pulsed with 5-ethynyl-2'-deoxyuridine (EdU) to label newly synthesized DNA, and collected for DNA extraction. After size fractionation on a sucrose gradient, Okazaki fragments are concentrated and purified before click chemistry is used to tag the EdU label with a biotin conjugate that is cleavable under mild conditions. Biotinylated Okazaki fragments are then captured on streptavidin beads and ligated to Illumina adapters before library preparation for Illumina sequencing. The use of Ok-seq to interrogate genome-wide replication fork initiation and termination efficiencies can be applied to all unperturbed, asynchronously growing mammalian cells or under conditions of replication stress, and the assay can be performed in less than 2 weeks.

Rationale: The cardiac sodium channel NaV1.5 has a fundamental role in excitability and conduction. Previous studies have shown that sodium channels cluster together in specific cellular subdomains. Their association with intracellular organelles in defined regions of the myocytes, and the functional consequences of that association, remain to be defined. Objective: To characterize a subcellular domain formed by sodium channel clusters in the crest region of the myocytes, and the subjacent subsarcolemmal mitochondria (SSM).Methods and Results: Through a combination of imaging approaches including super-resolution microscopy and electron microscopy we identified, in adult cardiac myocytes, a NaV1.5 subpopulation in close proximity to SSM; we further found that SSM preferentially host the mitochondrial Na+/Ca2+ exchanger (NCLX). This anatomical proximity led us to investigate functional changes in mitochondria resulting from sodium channel activity. Upon TTX exposure, mitochondria near NaV1.5 channels accumulated more Ca2+ and showed increased ROS production when compared to interfibrillar mitochondria. Finally, crosstalk between NaV1.5 channels and mitochondria was analyzed at a transcriptional level. We found that SCN5A and SLC8B1 (which encode NaV1.5 and NCLX, respectively) are negatively correlated both in a human transcriptome dataset (GTEx) and in human-induced pluripotent stem cell-derived cardiac myocytes deficient in SCN5A. Conclusions: We describe an anatomical hub (a couplon) formed by sodium channel clusters and SSM. Preferential localization of NCLX to this domain allows for functional coupling where the extrusion of Ca2+ from the mitochondria is powered, at least in part, by the entry of sodium through NaV1.5 channels. These results provide a novel entry-point into a mechanistic understanding of the intersection between electrical and structural functions of the heart.

Cellular processes are largely carried out by macromolecular assemblies, most of which are dynamic, having components that are in constant flux. One such assembly is the nuclear pore complex (NPC), an ∼50 MDa assembly comprised of ∼30 different proteins called Nups that mediates selective macromolecular transport between the nucleus and cytoplasm. We developed a proteomics method to provide a comprehensive picture of the yeast NPC component dynamics. We discovered that, although all Nups display uniformly slow turnover, their exchange rates vary considerably. Surprisingly, this exchange rate was relatively unrelated to each Nup's position, accessibility, or role in transport but correlated with its structural role; scaffold-forming Nups exchange slowly, whereas flexible connector Nups threading throughout the NPC architecture exchange more rapidly. Targeted perturbations in the NPC structure revealed a dynamic resilience to damage. Our approach opens a new window into macromolecular assembly dynamics.