PhosphoDisco: A Toolkit for Co-regulated Phosphorylation Module Discovery in Phosphoproteomic Data.

Kinases are key players in cancer-relevant pathways and are the targets of many successful precision cancer therapies. Phosphoproteomics is a powerful approach to study kinase activity and has been used increasingly for the characterization of tumor samples leading to the identification of novel chemotherapeutic targets and biomarkers. Finding co-regulated phosphorylation sites which represent potential kinase-substrate sets or members of the same signaling pathway allows us to harness these data to identify clinically relevant and targetable alterations in signaling cascades. Unfortunately, studies have found that databases of co-regulated phosphorylation sites are only experimentally supported in a small number of substrate sets. To address the inherent challenge of defining co-regulated phosphorylation modules relevant to a given dataset, we developed PhosphoDisco, a toolkit for determining co-regulated phosphorylation modules. We applied this approach to tandem mass spectrometry based phosphoproteomic data for breast and non-small cell lung cancer and identified canonical as well as putative new phosphorylation site modules. Our analysis identified several interesting modules in each cohort. Among these was a new cell cycle checkpoint module enriched in basal breast cancer samples and a module of PRKC isozymes putatively co-regulated by CDK12 in lung cancer. We demonstrate that modules defined by PhosphoDisco can be used to further personalized cancer treatment strategies by establishing active signaling pathways in a given patient tumor or set of tumors, and in providing new ways to classify tumors based on signaling activity.

Rapid TCR:Epitope Ranker (RAPTER): a primary human T cell reactivity screening assay pairing epitope and TCR at single cell resolution.

Identifying epitopes that T cells respond to is critical for understanding T cell-mediated immunity. Traditional multimer and other single cell assays often require large blood volumes and/or expensive HLA-specific reagents and provide limited phenotypic and functional information. Here, we present the Rapid TCR:Epitope Ranker (RAPTER) assay, a single cell RNA sequencing (scRNA-SEQ) method that uses primary human T cells and antigen presenting cells (APCs) to assess functional T cell reactivity. Using hash-tag oligonucleotide (HTO) coding and T cell activation-induced markers (AIM), RAPTER defines paired epitope specificity and TCR sequence and can include RNA- and protein-level T cell phenotype information. We demonstrate that RAPTER identified specific reactivities to viral and tumor antigens at sensitivities as low as 0.15% of total CD8+ T cells, and deconvoluted low-frequency circulating HPV16-specific T cell clones from a cervical cancer patient. The specificities of TCRs identified by RAPTER for MART1, EBV, and influenza epitopes were functionally confirmed in vitro. In summary, RAPTER identifies low-frequency T cell reactivities using primary cells from low blood volumes, and the resulting paired TCR:ligand information can directly enable immunogenic antigen selection from limited patient samples for vaccine epitope inclusion, antigen-specific TCR tracking, and TCR cloning for further therapeutic development.

Melanoma-Secreted Amyloid Beta Suppresses Neuroinflammation and Promotes Brain Metastasis.

Brain metastasis is a significant cause of morbidity and mortality in multiple cancer types and represents an unmet clinical need. The mechanisms that mediate metastatic cancer growth in the brain parenchyma are largely unknown. Melanoma, which has the highest rate of brain metastasis among common cancer types, is an ideal model to study how cancer cells adapt to the brain parenchyma. Our unbiased proteomics analysis of melanoma short-term cultures revealed that proteins implicated in neurodegenerative pathologies are differentially expressed in melanoma cells explanted from brain metastases compared with those derived from extracranial metastases. We showed that melanoma cells require amyloid beta (Aß) for growth and survival in the brain parenchyma. Melanoma-secreted Aß activates surrounding astrocytes to a prometastatic, anti-inflammatory phenotype and prevents phagocytosis of melanoma by microglia. Finally, we demonstrate that pharmacologic inhibition of Aß decreases brain metastatic burden. SIGNIFICANCE: Our results reveal a novel mechanistic connection between brain metastasis and Alzheimer's disease, two previously unrelated pathologies; establish Aß as a promising therapeutic target for brain metastasis; and demonstrate suppression of neuroinflammation as a critical feature of metastatic adaptation to the brain parenchyma. This article is highlighted in the In This Issue feature, p. 1171.

BlackSheep: A Bioconductor and Bioconda Package for Differential Extreme Value Analysis.

Unbiased assays such as shotgun proteomics and RNA-seq provide high-resolution molecular characterization of tumors. These assays measure molecules with highly varied distributions, making interpretation and hypothesis testing challenging. Samples with the most extreme measurements for a molecule can reveal the most interesting biological insights yet are often excluded from analysis. Furthermore, rare disease subtypes are, by definition, underrepresented in cancer cohorts. To provide a strategy for identifying molecules aberrantly enriched in small sample cohorts, we present BlackSheep, a package for nonparametric description and differential analysis of genome-wide data, available from Bioconductor (https://www.bioconductor.org/packages/release/bioc/html/blacksheepr.html) and Bioconda (https://bioconda.github.io/recipes/blksheep/README.html). BlackSheep is a complementary tool to other differential expression analysis methods, which is particularly useful when analyzing small subgroups in a larger cohort.

Proteogenomic Landscape of Breast Cancer Tumorigenesis and Targeted Therapy.

The integration of mass spectrometry-based proteomics with next-generation DNA and RNA sequencing profiles tumors more comprehensively. Here this "proteogenomics" approach was applied to 122 treatment-naive primary breast cancers accrued to preserve post-translational modifications, including protein phosphorylation and acetylation. Proteogenomics challenged standard breast cancer diagnoses, provided detailed analysis of the ERBB2 amplicon, defined tumor subsets that could benefit from immune checkpoint therapy, and allowed more accurate assessment of Rb status for prediction of CDK4/6 inhibitor responsiveness. Phosphoproteomics profiles uncovered novel associations between tumor suppressor loss and targetable kinases. Acetylproteome analysis highlighted acetylation on key nuclear proteins involved in the DNA damage response and revealed cross-talk between cytoplasmic and mitochondrial acetylation and metabolism. Our results underscore the potential of proteogenomics for clinical investigation of breast cancer through more accurate annotation of targetable pathways and biological features of this remarkably heterogeneous malignancy.

Hyperactive CDK2 Activity in Basal-like Breast Cancer Imposes a Genome Integrity Liability that Can Be Exploited by Targeting DNA Polymerase ε.

Knowledge of fundamental differences between breast cancer subtypes has driven therapeutic advances; however, basal-like breast cancer (BLBC) remains clinically intractable. Because BLBC exhibits alterations in DNA repair enzymes and cell-cycle checkpoints, elucidation of factors enabling the genomic instability present in this subtype has the potential to reveal novel anti-cancer strategies. Here, we demonstrate that BLBC is especially sensitive to suppression of iron-sulfur cluster (ISC) biosynthesis and identify DNA polymerase epsilon (POLE) as an ISC-containing protein that underlies this phenotype. In BLBC cells, POLE suppression leads to replication fork stalling, DNA damage, and a senescence-like state or cell death. In contrast, luminal breast cancer and non-transformed mammary cells maintain viability upon POLE suppression but become dependent upon an ATR/CHK1/CDC25A/CDK2 DNA damage response axis. We find that CDK1/2 targets exhibit hyperphosphorylation selectively in BLBC tumors, indicating that CDK2 hyperactivity is a genome integrity vulnerability exploitable by targeting POLE.

Hypercluster: a flexible tool for parallelized unsupervised clustering optimization.

BACKGROUND: Unsupervised clustering is a common and exceptionally useful tool for large biological datasets. However, clustering requires upfront algorithm and hyperparameter selection, which can introduce bias into the final clustering labels. It is therefore advisable to obtain a range of clustering results from multiple models and hyperparameters, which can be cumbersome and slow. RESULTS: We present hypercluster, a python package and SnakeMake pipeline for flexible and parallelized clustering evaluation and selection. Users can efficiently evaluate a huge range of clustering results from multiple models and hyperparameters to identify an optimal model. CONCLUSIONS: Hypercluster improves ease of use, robustness and reproducibility for unsupervised clustering application for high throughput biology. Hypercluster is available on pip and bioconda; installation, documentation and example workflows can be found at: https://github.com/ruggleslab/hypercluster .

Proteogenomic Characterization Reveals Therapeutic Vulnerabilities in Lung Adenocarcinoma.

To explore the biology of lung adenocarcinoma (LUAD) and identify new therapeutic opportunities, we performed comprehensive proteogenomic characterization of 110 tumors and 101 matched normal adjacent tissues (NATs) incorporating genomics, epigenomics, deep-scale proteomics, phosphoproteomics, and acetylproteomics. Multi-omics clustering revealed four subgroups defined by key driver mutations, country, and gender. Proteomic and phosphoproteomic data illuminated biology downstream of copy number aberrations, somatic mutations, and fusions and identified therapeutic vulnerabilities associated with driver events involving KRAS, EGFR, and ALK. Immune subtyping revealed a complex landscape, reinforced the association of STK11 with immune-cold behavior, and underscored a potential immunosuppressive role of neutrophil degranulation. Smoking-associated LUADs showed correlation with other environmental exposure signatures and a field effect in NATs. Matched NATs allowed identification of differentially expressed proteins with potential diagnostic and therapeutic utility. This proteogenomics dataset represents a unique public resource for researchers and clinicians seeking to better understand and treat lung adenocarcinomas.

Proteogenomic Characterization of Endometrial Carcinoma.

We undertook a comprehensive proteogenomic characterization of 95 prospectively collected endometrial carcinomas, comprising 83 endometrioid and 12 serous tumors. This analysis revealed possible new consequences of perturbations to the p53 and Wnt/ß-catenin pathways, identified a potential role for circRNAs in the epithelial-mesenchymal transition, and provided new information about proteomic markers of clinical and genomic tumor subgroups, including relationships to known druggable pathways. An extensive genome-wide acetylation survey yielded insights into regulatory mechanisms linking Wnt signaling and histone acetylation. We also characterized aspects of the tumor immune landscape, including immunogenic alterations, neoantigens, common cancer/testis antigens, and the immune microenvironment, all of which can inform immunotherapy decisions. Collectively, our multi-omic analyses provide a valuable resource for researchers and clinicians, identify new molecular associations of potential mechanistic significance in the development of endometrial cancers, and suggest novel approaches for identifying potential therapeutic targets.

Impaired Expression of Rearranged Immunoglobulin Genes and Premature p53 Activation Block B Cell Development in BMI1 Null Mice.

B cell development is a highly regulated process that requires stepwise rearrangement of immunoglobulin genes to generate a functional B cell receptor (BCR). The polycomb group protein BMI1 is required for B cell development, but its function in developing B cells remains poorly defined. We demonstrate that BMI1 functions in a cell-autonomous manner at two stages during early B cell development. First, loss of BMI1 results in a differentiation block at the pro-B cell to pre-B cell transition due to the inability of BMI1-deficient cells to transcribe newly rearranged Igh genes. Accordingly, introduction of a pre-rearranged Igh allele partially restored B cell development in Bmi1-/- mice. In addition, BMI1 is required to prevent premature p53 signaling, and as a consequence, Bmi1-/- large pre-B cells fail to properly proliferate. Altogether, our results clarify the role of BMI1 in early B cell development and uncover an unexpected function of BMI1 during VDJ recombination.

Capturing the Onset of PRC2-Mediated Repressive Domain Formation.

Polycomb repressive complex 2 (PRC2) maintains gene silencing by catalyzing methylation of histone H3 at lysine 27 (H3K27me2/3) within chromatin. By designing a system whereby PRC2-mediated repressive domains were collapsed and then reconstructed in an inducible fashion in vivo, a two-step mechanism of H3K27me2/3 domain formation became evident. First, PRC2 is stably recruited by the actions of JARID2 and MTF2 to a limited number of spatially interacting "nucleation sites," creating H3K27me3-forming Polycomb foci within the nucleus. Second, PRC2 is allosterically activated via its binding to H3K27me3 and rapidly spreads H3K27me2/3 both in cis and in far-cis via long-range contacts. As PRC2 proceeds further from the nucleation sites, its stability on chromatin decreases such that domains of H3K27me3 remain proximal, and those of H3K27me2 distal, to the nucleation sites. This study demonstrates the principles of de novo establishment of PRC2-mediated repressive domains across the genome.

The Conserved ATM Kinase RAG2-S365 Phosphorylation Site Limits Cleavage Events in Individual Cells Independent of Any Repair Defect.

Many DNA lesions associated with lymphoid malignancies are linked to off-target cleavage by the RAG1/2 recombinase. However, off-target cleavage has mostly been analyzed in the context of DNA repair defects, confounding any mechanistic understanding of cleavage deregulation. We identified a conserved SQ phosphorylation site on RAG2 365 to 366 that is involved in feedback control of RAG cleavage. Mutation of serine 365 to a non-phosphorylatable alanine permits bi-allelic and bi-locus RAG-mediated breaks in the same cell, leading to reciprocal translocations. This phenomenon is analogous to the phenotype we described for ATM kinase inactivation. Here, we establish deregulated cleavage itself as a driver of chromosomal instability without the associated repair defect. Intriguingly, a RAG2-S365E phosphomimetic rescues the deregulated cleavage of ATM inactivation, reducing the incidence of reciprocal translocations. These data support a model in which feedback control of cleavage and maintenance of genome stability involves ATM-mediated phosphorylation of RAG2.

miRNAs Are Essential for the Regulation of the PI3K/AKT/FOXO Pathway and Receptor Editing during B Cell Maturation.

B cell development is a tightly regulated process dependent on sequential rearrangements of immunoglobulin loci that encode the antigen receptor. To elucidate the role of microRNAs (miRNAs) in the orchestration of B cell development, we ablated all miRNAs at the earliest stage of B cell development by conditionally targeting the enzymes critical for RNAi in early B cell precursors. Absence of any one of these enzymes led to a block at the pro- to pre-B cell transition due to increased apoptosis and a failure of pre-B cells to proliferate. Expression of a Bcl2 transgene allowed for partial rescue of B cell development, however, the majority of the rescued B cells had low surface immunoglobulin expression with evidence of ongoing light chain editing. Our analysis revealed that miRNAs are critical for the regulation of the PTEN-AKT-FOXO1 pathway that in turn controls Rag expression during B cell development.

RAG Off-Target Activity Is in the Loop.

The 'off-target' activity of RAG recombinases contributes to mutations and cancer. Recent studies show that the influence of DNA regulatory elements is largely constrained by the formation of chromatin loops and interaction frequencies. Now, the Alt lab has identified major RAG off-target activity bearing similar limitations, with joining events restricted to convergent paired RSS elements in loop domains.

Phylogenetic analysis of algal symbionts associated with four North American amphibian egg masses.

Egg masses of the yellow-spotted salamander Ambystoma maculatum form an association with the green alga "Oophila amblystomatis" (Lambert ex Wille), which, in addition to growing within individual egg capsules, has recently been reported to invade embryonic tissues and cells. The binomial O. amblystomatis refers to the algae that occur in A. maculatum egg capsules, but it is unknown whether this population of symbionts constitutes one or several different algal taxa. Moreover, it is unknown whether egg masses across the geographic range of A. maculatum, or other amphibians, associate with one or multiple algal taxa. To address these questions, we conducted a phylogeographic study of algae sampled from egg capsules of A. maculatum, its allopatric congener A. gracile, and two frogs: Lithobates sylvatica and L. aurora. All of these North American amphibians form associations with algae in their egg capsules. We sampled algae from egg capsules of these four amphibians from localities across North America, established representative algal cultures, and amplified and sequenced a region of 18S rDNA for phylogenetic analysis. Our combined analysis shows that symbiotic algae found in egg masses of four North American amphibians are closely related to each other, and form a well-supported clade that also contains three strains of free-living chlamydomonads. We designate this group as the 'Oophila' clade, within which the symbiotic algae are further divided into four distinct subclades. Phylogenies of the host amphibians and their algal symbionts are only partially congruent, suggesting that host-switching and co-speciation both play roles in their associations. We also established conditions for isolating and rearing algal symbionts from amphibian egg capsules, which should facilitate further study of these egg mass specialist algae.