Mapping of long interspersed element-1 (L1) insertions by TIPseq provides information about sub chromosomal genetic variation in human embryos.

PURPOSE: Retrotransposons play important roles during early development when they are transiently de-repressed during epigenetic reprogramming. Long interspersed element-1 (L1), the only autonomous retrotransposon in humans, comprises 17% of the human genome. We applied the Single Cell Transposon Insertion Profiling by Sequencing (scTIPseq) to characterize and map L1 insertions in human embryos. METHODS: Sixteen cryopreserved, genetically tested, human blastocysts, were accessed from consenting couples undergoing IVF at NYU Langone Fertility Center. Additionally, four trios (father, mother, and embryos) were also evaluated. scTIPseq was applied to map L1 insertions in all samples, using L1 locations reported in the 1000 Genomes as controls. RESULTS: Twenty-nine unknown and unique insertions were observed in the sixteen embryos. Most were intergenic; no insertions were located in exons or immediately upstream of genes. The location or number of unknown insertions did not differ between euploid and aneuploid embryos, suggesting they are not merely markers of aneuploidy. Rather, scTIPseq provides novel information about sub-chromosomal structural variation in human embryos. Trio analyses showed a parental origin of all L1 insertions in embryos. CONCLUSION: Several studies have measured L1 expression at different stages of development in mice, but this study for the first time reports unknown insertions in human embryos that were inherited from one parent, confirming no de novo L1 insertions occurred in parental germline or during embryogenesis. Since one-third of euploid embryo transfers fail, future studies would be useful for understanding whether these sub-chromosomal genetic variants or de novo L1 insertions affect embryo developmental potential.

Predicting immune checkpoint therapy response in three independent metastatic melanoma cohorts.

Introduction: While Immune checkpoint inhibition (ICI) therapy shows significant efficacy in metastatic melanoma, only about 50% respond, lacking reliable predictive methods. We introduce a panel of six proteins aimed at predicting response to ICI therapy. Methods: Evaluating previously reported proteins in two untreated melanoma cohorts, we used a published predictive model (EaSIeR score) to identify potential proteins distinguishing responders and non-responders. Results: Six proteins initially identified in the ICI cohort correlated with predicted response in the untreated cohort. Additionally, three proteins correlated with patient survival, both at the protein, and at the transcript levels, in an independent immunotherapy treated cohort. Discussion: Our study identifies predictive biomarkers across three melanoma cohorts, suggesting their use in therapeutic decision-making.

Dormant origin firing promotes head-on transcription-replication conflicts at transcription termination sites in response to BRCA2 deficiency.

BRCA2 is a tumor suppressor protein responsible for safeguarding the cellular genome from replication stress and genotoxicity, but the specific mechanism(s) by which this is achieved to prevent early oncogenesis remains unclear. Here, we provide evidence that BRCA2 acts as a critical suppressor of head-on transcription-replication conflicts (HO-TRCs). Using Okazaki-fragment sequencing (Ok-seq) and computational analysis, we identified origins (dormant origins) that are activated near the transcription termination sites (TTS) of highly expressed, long genes in response to replication stress. Dormant origins are a source for HO-TRCs, and drug treatments that inhibit dormant origin firing led to a reduction in HO-TRCs, R-loop formation, and DNA damage. Using super-resolution microscopy, we showed that HO-TRC events track with elongating RNA polymerase II, but not with transcription initiation. Importantly, RNase H2 is recruited to sites of HO-TRCs in a BRCA2-dependent manner to help alleviate toxic R-loops associated with HO-TRCs. Collectively, our results provide a mechanistic basis for how BRCA2 shields against genomic instability by preventing HO-TRCs through both direct and indirect means occurring at predetermined genomic sites based on the pre-cancer transcriptome.

Modulation of GPR133 (ADGRD1) signaling by its intracellular interaction partner extended synaptotagmin 1.

GPR133 (ADGRD1) is an adhesion G-protein-coupled receptor that signals through Gαs/cyclic AMP (cAMP) and is required for the growth of glioblastoma (GBM), an aggressive brain malignancy. The regulation of GPR133 signaling is incompletely understood. Here, we use proximity biotinylation proteomics to identify ESYT1, a Ca2+-dependent mediator of endoplasmic reticulum-plasma membrane bridge formation, as an intracellular interactor of GPR133. ESYT1 knockdown or knockout increases GPR133 signaling, while its overexpression has the opposite effect, without altering GPR133 levels in the plasma membrane. The GPR133-ESYT1 interaction requires the Ca2+-sensing C2C domain of ESYT1. Thapsigargin-mediated increases in cytosolic Ca2+ relieve signaling-suppressive effects of ESYT1 by promoting ESYT1-GPR133 dissociation. ESYT1 knockdown or knockout in GBM slows tumor growth, suggesting tumorigenic functions of ESYT1. Our findings demonstrate a mechanism for the modulation of GPR133 signaling by increased cytosolic Ca2+, which reduces the signaling-suppressive interaction between GPR133 and ESYT1 to raise cAMP levels.

Predicting immune checkpoint therapy response in three independent metastatic melanoma cohorts.

While Immune checkpoint inhibition (ICI) therapy shows significant efficacy in metastatic melanoma, only about 50% respond, lacking reliable predictive methods. We introduce a panel of six proteins aimed at predicting response to ICI therapy. Evaluating previously reported proteins in two untreated melanoma cohorts, we used a published predictive model (EaSIeR score) to identify potential proteins distinguishing responders and non-responders. Six proteins initially identified in the ICI cohort correlated with predicted response in the untreated cohort. Additionally, three proteins correlated with patient survival, both at the protein, and at the transcript levels, in an independent immunotherapy treated cohort. Our study identifies predictive biomarkers across three melanoma cohorts, suggesting their use in therapeutic decision-making.

The MTIST platform: a microbiome time series inference standardized test.

The human gut microbiome is a promising therapeutic target, but interventions are hampered by our limited understanding of microbial ecosystems. Here, we present a platform to develop, evaluate, and score approaches to learn ecological interactions from microbiome time series data. The microbiome time series inference standardized test (MTIST) comprises: a simulation framework for the in silico generation of microbiome study data akin to what is obtained with quantitative next-generation sequencing approaches, a compilation of a large curated data set generated by the simulation framework representing 648 simulated microbiome studies containing 18,360 time series, with a total of 2,182,800 species abundance measurements, and a scoring method to rank ecological inference algorithms. We use the MTIST platform to rank five implementations of microbiome inference approaches, revealing that while all algorithms performed well on ecosystems with few species (3 and 10), all algorithms failed to infer most interaction in a large ecosystem with 100 member species. However, we do find that the strongest interactions within a large ecosystem are inferred with higher success by all algorithms. Finally, we use the MTIST platform to compare different microbiome study designs, characterizing tradeoffs between samples per subject and number of subjects. Interestingly, we find that when only few samples can be collected per subject, ecological inference is most successful when these samples are collected with highest feasible temporal frequency. Taken together, we provide a computational tool to aid the development of better microbiome ecosystem inference approaches, which will be crucial towards the development of reliable and predictable therapeutic approaches that target the microbiome ecosystem.

Few-shot genes selection: subset of PAM50 genes for breast cancer subtypes classification.

BACKGROUND: In recent years, researchers have made significant strides in understanding the heterogeneity of breast cancer and its various subtypes. However, the wealth of genomic and proteomic data available today necessitates efficient frameworks, instruments, and computational tools for meaningful analysis. Despite its success as a prognostic tool, the PAM50 gene signature's reliance on many genes presents challenges in terms of cost and complexity. Consequently, there is a need for more efficient methods to classify breast cancer subtypes using a reduced gene set accurately. RESULTS: This study explores the potential of achieving precise breast cancer subtype categorization using a reduced gene set derived from the PAM50 gene signature. By employing a "Few-Shot Genes Selection" method, we randomly select smaller subsets from PAM50 and evaluate their performance using metrics and a linear model, specifically the Support Vector Machine (SVM) classifier. In addition, we aim to assess whether a more compact gene set can maintain performance while simplifying the classification process. Our findings demonstrate that certain reduced gene subsets can perform comparable or superior to the full PAM50 gene signature. CONCLUSIONS: The identified gene subsets, with 36 genes, have the potential to contribute to the development of more cost-effective and streamlined diagnostic tools in breast cancer research and clinical settings.

Identifying Ferroptosis Inducers, HDAC, and RTK Inhibitor Sensitivity in Melanoma Subtypes through Unbiased Drug Target Prediction.

The utilization of PD1 and CTLA4 inhibitors has revolutionized the treatment of malignant melanoma (MM). However, resistance to targeted and immune-checkpoint-based therapies still poses a significant problem. Here we mine large scale MM proteogenomic data integrating it with MM cell line dependency screen, and drug sensitivity data to identify druggable targets and forecast treatment efficacy and resistance. Leveraging protein profiles from established MM subtypes and molecular structures of 82 cancer treatment drugs, we identified nine candidate hub proteins, mTOR, FYN, PIK3CB, EGFR, MAPK3, MAP4K1, MAP2K1, SRC and AKT1, across five distinct MM subtypes. These proteins serve as potential drug targets applicable to one or multiple MM subtypes. By analyzing transcriptomic data from 48 publicly accessible melanoma cell lines sourced from Achilles and CRISPR dependency screens, we forecasted 162 potentially targetable genes. We also identified genetic resistance in 260 genes across at least one melanoma subtype. In addition, we employed publicly available compound sensitivity data (Cancer Therapeutics Response Portal, CTRPv2) on the cell lines to assess the correlation of compound effectiveness within each subtype. We have identified 20 compounds exhibiting potential drug impact in at least one melanoma subtype. Remarkably, employing this unbiased approach, we have uncovered compounds targeting ferroptosis, that demonstrate a striking 30x fold difference in sensitivity among different subtypes. This implies that the proteogenomic classification of melanoma has the potential to predict sensitivity to ferroptosis compounds. Our results suggest innovative and novel therapeutic strategies by stratifying melanoma samples through proteomic profiling, offering a spectrum of novel therapeutic interventions and prospects for combination therapy. Highlights: (1) Proteogenomic subtype classification can define the landscape of genetic dependencies in melanoma (2) Nine proteins from molecular subtypes were identified as potential drug targets for specified MM patients (3) 20 compounds identified that show potential effectiveness in at least one melanoma subtype (4) Proteogenomics can predict specific ferroptosis inducers, HDAC, and RTK Inhibitor sensitivity in melanoma subtypes.

CDK-independent role of D-type cyclins in regulating DNA mismatch repair.

Although mismatch repair (MMR) is essential for correcting DNA replication errors, it can also recognize other lesions, such as oxidized bases. In G0 and G1, MMR is kept in check through unknown mechanisms as it is error-prone during these cell cycle phases. We show that in mammalian cells, D-type cyclins are recruited to sites of oxidative DNA damage in a PCNA- and p21-dependent manner. D-type cyclins inhibit the proteasomal degradation of p21, which competes with MMR proteins for binding to PCNA, thereby inhibiting MMR. The ability of D-type cyclins to limit MMR is CDK4- and CDK6-independent and is conserved in G0 and G1. At the G1/S transition, the timely, cullin-RING ubiquitin ligase (CRL)-dependent degradation of D-type cyclins and p21 enables MMR activity to efficiently repair DNA replication errors. Persistent expression of D-type cyclins during S-phase inhibits the binding of MMR proteins to PCNA, increases the mutational burden, and promotes microsatellite instability.

D-type cyclins regulate DNA mismatch repair in the G1 and S phases of the cell cycle, maintaining genome stability.

The large majority of oxidative DNA lesions occurring in the G1 phase of the cell cycle are repaired by base excision repair (BER) rather than mismatch repair (MMR) to avoid long resections that can lead to genomic instability and cell death. However, the molecular mechanisms dictating pathway choice between MMR and BER have remained unknown. Here, we show that, during G1, D-type cyclins are recruited to sites of oxidative DNA damage in a PCNA- and p21-dependent manner. D-type cyclins shield p21 from its two ubiquitin ligases CRL1SKP2 and CRL4CDT2 in a CDK4/6-independent manner. In turn, p21 competes through its PCNA-interacting protein degron with MMR components for their binding to PCNA. This inhibits MMR while not affecting BER. At the G1/S transition, the CRL4AMBRA1-dependent degradation of D-type cyclins renders p21 susceptible to proteolysis. These timely degradation events allow the proper binding of MMR proteins to PCNA, enabling the repair of DNA replication errors. Persistent expression of cyclin D1 during S-phase increases the mutational burden and promotes microsatellite instability. Thus, the expression of D-type cyclins inhibits MMR in G1, whereas their degradation is necessary for proper MMR function in S.

The expression profile and tumorigenic mechanisms of CD97 (ADGRE5) in glioblastoma render it a targetable vulnerability.

Glioblastoma (GBM) is the most common and aggressive primary brain malignancy. Adhesion G protein-coupled receptors (aGPCRs) have attracted interest for their potential as treatment targets. Here, we show that CD97 (ADGRE5) is the most promising aGPCR target in GBM, by virtue of its de novo expression compared to healthy brain tissue. CD97 knockdown or knockout significantly reduces the tumor initiation capacity of patient-derived GBM cultures (PDGCs) in vitro and in vivo. We find that CD97 promotes glycolytic metabolism via the mitogen-activated protein kinase (MAPK) pathway, which depends on phosphorylation of its C terminus and recruitment of ß-arrestin. We also demonstrate that THY1/CD90 is a likely CD97 ligand in GBM. Lastly, we show that an anti-CD97 antibody-drug conjugate selectively kills tumor cells in vitro. Our studies identify CD97 as a regulator of tumor metabolism, elucidate mechanisms of receptor activation and signaling, and provide strong scientific rationale for developing biologics to target it therapeutically in GBM.

Engineered transcription-associated Cas9 targeting in eukaryotic cells.

DNA targeting Class 2 CRISPR-Cas effector nucleases, including the well-studied Cas9 proteins, evolved protospacer-adjacent motif (PAM) and guide RNA interactions that sequentially license their binding and cleavage activities at protospacer target sites. Both interactions are nucleic acid sequence specific but function constitutively; thus, they provide intrinsic spatial control over DNA targeting activities but naturally lack temporal control. Here we show that engineered Cas9 fusion proteins which bind to nascent RNAs near a protospacer can facilitate spatiotemporal coupling between transcription and DNA targeting at that protospacer: Transcription-associated Cas9 Targeting (TraCT). Engineered TraCT is enabled when suboptimal PAM interactions limit basal activity in vivo and when one or more nascent RNA substrates are still tethered to the actively transcribing target DNA in cis. We further show that this phenomenon can be exploited for selective editing at one of two identical targets in distinct gene loci, or, in diploid allelic loci that are differentially transcribed. Our work demonstrates that temporal control over Cas9's targeting activity at specific DNA sites may be engineered without modifying Cas9's core domains and guide RNA components or their expression levels. More broadly, it establishes RNA binding in cis as a mechanism that can conditionally stimulate CRISPR-Cas DNA targeting in eukaryotes.

Proteogenomics Reveal the Overexpression of HLA-I in Cancer.

An accurate quantification of HLA class I gene expression is important in understanding the interplay with the tumor microenvironment of antitumor cytotoxic T cell activities. Because HLA-I sequences are highly variable, standard RNAseq and mass spectrometry-based quantification workflows using common genome and protein sequence references do not provide HLA-I allele specific quantifications. Here, we used personalized HLA-I nucleotide and protein reference sequences based on the subjects' HLA-I genotypes and surveyed tumor and adjacent normal samples from patients across nine cancer types. Mass spectrometry using data dependent acquisition data was validated to be sufficient to estimate HLA-A protein expression at the allele level. We found that HLA-I proteins were present in significantly higher levels in tumors compared to adjacent normal tissues from 41 to 63% of head and neck squamous cell carcinoma, uterine corpus endometrial carcinoma, and clear cell renal cell carcinoma patients, and this was driven by increased levels of HLA-I gene transcripts. Most immune cell types are universally enriched in HLA-I high tumors, while endothelial and neuronal cells showed divergent relationships with HLA-I. Pathway analysis revealed that tumor senescence and autophagy activity influence the level of HLA-I proteins in glioblastoma. Genes correlated to HLA-I protein expression are mostly the ones directly involved in HLA-I function in immune response and cell death, while glycosylation genes are exclusively co-expressed with HLA-I at the protein level.

MED19 encodes two unique protein isoforms that confer prostate cancer growth under low androgen through distinct gene expression programs.

MED19, a component of the mediator complex and a co-regulator of the androgen receptor (AR), is pivotal in prostate cancer cell proliferation. MED19 has two isoforms: a full-length "canonical" and a shorter "alternative" variant. Specific antibodies were developed to investigate these isoforms. Both exhibit similar expression in normal prostate development and adult prostate tissue, but the canonical isoform is elevated in prostate adenocarcinomas. Overexpression of canonical MED19 in LNCaP cells promotes growth under conditions of androgen deprivation in vitro and in vivo, mirroring earlier findings with alternative MED19-overexpressing LNCaP cells. Interestingly, alternative MED19 cells displayed strong colony formation in clonogenic assays under conditions of androgen deprivation, while canonical MED19 cells did not, suggesting distinct functional roles. These isoforms also modulated gene expression differently. Canonical MED19 triggered genes related to extracellular matrix remodeling while suppressing those involved in androgen-inactivating glucuronidation. In contrast, alternative MED19 elevated genes tied to cell movement and reduced those associated with cell adhesion and differentiation. The ratio of MED19 isoform expression in prostate cancers shifts with the disease stage. Early-stage cancers exhibit higher canonical MED19 expression than alternative MED19, consistent with canonical MED19's ability to promote cell proliferation under androgen deprivation. Conversely, alternative MED19 levels were higher in later-stage metastatic prostate cancer than in canonical MED19, reflecting alternative MED19's capability to enhance cell migration and autonomous cell growth. Our findings suggest that MED19 isoforms play unique roles in prostate cancer progression and highlights MED19 as a potential therapeutic target for both early and late-stage prostate cancer.

Proteogenomic data and resources for pan-cancer analysis.

The National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium (CPTAC) investigates tumors from a proteogenomic perspective, creating rich multi-omics datasets connecting genomic aberrations to cancer phenotypes. To facilitate pan-cancer investigations, we have generated harmonized genomic, transcriptomic, proteomic, and clinical data for >1000 tumors in 10 cohorts to create a cohesive and powerful dataset for scientific discovery. We outline efforts by the CPTAC pan-cancer working group in data harmonization, data dissemination, and computational resources for aiding biological discoveries. We also discuss challenges for multi-omics data integration and analysis, specifically the unique challenges of working with both nucleotide sequencing and mass spectrometry proteomics data.

Pan-cancer proteogenomics connects oncogenic drivers to functional states.

Cancer driver events refer to key genetic aberrations that drive oncogenesis; however, their exact molecular mechanisms remain insufficiently understood. Here, our multi-omics pan-cancer analysis uncovers insights into the impacts of cancer drivers by identifying their significant cis-effects and distal trans-effects quantified at the RNA, protein, and phosphoprotein levels. Salient observations include the association of point mutations and copy-number alterations with the rewiring of protein interaction networks, and notably, most cancer genes converge toward similar molecular states denoted by sequence-based kinase activity profiles. A correlation between predicted neoantigen burden and measured T cell infiltration suggests potential vulnerabilities for immunotherapies. Patterns of cancer hallmarks vary by polygenic protein abundance ranging from uniform to heterogeneous. Overall, our work demonstrates the value of comprehensive proteogenomics in understanding the functional states of oncogenic drivers and their links to cancer development, surpassing the limitations of studying individual cancer types.

Integrative multi-omic cancer profiling reveals DNA methylation patterns associated with therapeutic vulnerability and cell-of-origin.

DNA methylation plays a critical role in establishing and maintaining cellular identity. However, it is frequently dysregulated during tumor development and is closely intertwined with other genetic alterations. Here, we leveraged multi-omic profiling of 687 tumors and matched non-involved adjacent tissues from the kidney, brain, pancreas, lung, head and neck, and endometrium to identify aberrant methylation associated with RNA and protein abundance changes and build a Pan-Cancer catalog. We uncovered lineage-specific epigenetic drivers including hypomethylated FGFR2 in endometrial cancer. We showed that hypermethylated STAT5A is associated with pervasive regulon downregulation and immune cell depletion, suggesting that epigenetic regulation of STAT5A expression constitutes a molecular switch for immunosuppression in squamous tumors. We further demonstrated that methylation subtype-enrichment information can explain cell-of-origin, intra-tumor heterogeneity, and tumor phenotypes. Overall, we identified cis-acting DNA methylation events that drive transcriptional and translational changes, shedding light on the tumor's epigenetic landscape and the role of its cell-of-origin.