Rp3: Ribosome profiling-assisted proteogenomics improves coverage and confidence during microprotein discovery.

There has been a dramatic increase in the identification of non-canonical translation and a significant expansion of the protein-coding genome. Among the strategies used to identify unannotated small Open Reading Frames (smORFs) that encode microproteins, Ribosome profiling (Ribo-Seq) is the gold standard for the annotation of novel coding sequences by reporting on smORF translation. In Ribo-Seq, ribosome-protected footprints (RPFs) that map to multiple genomic sites are removed since they cannot be unambiguously assigned to a specific genomic location. Furthermore, RPFs necessarily result in short (25-34 nucleotides) reads, increasing the chance of multi-mapping alignments, such that smORFs residing in these regions cannot be identified by Ribo-Seq. Moreover, it has been challenging to identify protein evidence for Ribo-Seq. To solve this, we developed Rp3, a pipeline that integrates proteogenomics and Ribosome profiling to provide unambiguous evidence for a subset of microproteins missed by current Ribo-Seq pipelines. Here, we show that Rp3 maximizes proteomics detection and confidence of microprotein-encoding smORFs.

Proteogenomic network analysis reveals dysregulated mechanisms and potential mediators in Parkinson's disease.

Parkinson's disease is highly heterogeneous across disease symptoms, clinical manifestations and progression trajectories, hampering the identification of therapeutic targets. Despite knowledge gleaned from genetics analysis, dysregulated proteome mechanisms stemming from genetic aberrations remain underexplored. In this study, we develop a three-phase system-level proteogenomic analytical framework to characterize disease-associated proteins and dysregulated mechanisms. Proteogenomic analysis identified 577 proteins that enrich for Parkinson's disease-related pathways, such as cytokine receptor interactions and lysosomal function. Converging lines of evidence identified nine proteins, including LGALS3, CSNK2A1, SMPD3, STX4, APOA2, PAFAH1B3, LDLR, HSPB1, BRK1, with potential roles in disease pathogenesis. This study leverages the largest population-scale proteomics dataset, the UK Biobank Pharma Proteomics Project, to characterize genetically-driven protein disturbances associated with Parkinson's disease. Taken together, our work contributes to better understanding of genome-proteome dynamics in Parkinson's disease and sets a paradigm to identify potential indirect mediators connected to GWAS signals for complex neurodegenerative disorders.

Exploring the Alternative Proteome with OpenProt and Mass Spectrometry.

Proteogenomics has revealed the translation of unannotated open reading frames (ORFs) present in mRNAs and in noncoding RNAs (ncRNAs). OpenProt annotates all ORFs with a minimum of 30 codons in the transcriptome of several species and displays many functional features associated with the corresponding proteins. Two types of proteins are annotated: reference or canonical proteins which are proteins already annotated in UniProt, RefSeq, or Ensembl and noncanonical proteins. Noncanonical proteins form two groups: predicted novel isoforms that display a significant level of homology with a reference protein and alternative proteins that are new proteins with no significant homology to known proteins. This chapter describes how to check whether a gene and/or transcript contains multiple open reading frames and how to use OpenProt databases for the detection of alternative proteins and novel isoforms by mass spectrometry-based proteomics.

Identification of Novel Bacterial Microproteins Encoded by Small Open Reading Frames Using a Computational Proteogenomics Workflow.

Genome annotation has historically ignored small open reading frames (smORFs), which encode a class of proteins shorter than 100 amino acids, collectively referred to as microproteins. This cutoff was established to avoid thousands of false positives due to limitations of pure genomics pipelines. Proteogenomics, a computational approach that combines genomics, transcriptomics, and proteomics, makes it possible to accurately identify these short sequences by overlaying different levels of omics evidence. In this chapter, we showcase the use of µProteInS, a bioinformatics pipeline developed for the identification of unannotated microproteins encoded by smORFs in bacteria. The workflow covers all the steps from quality control and transcriptome assembly to the scoring and post-processing of mass spectrometry data. Additionally, we provide an example on how to apply the pipeline's machine learning method to identify high-confidence spectra and pinpoint the most reliable identifications from large datasets.

A Microdissection Protocol for Proteogenomic Analysis of Histological Sections to Advance Drug Development.

Combining proteogenomics with laser capture microdissection (LCM) in cancer research offers a targeted way to explore the intricate interactions between tumor cells and the different microenvironment components. This is especially important for immuno-oncology (IO) research where improvements in the predictability of IO-based drugs are sorely needed, and depends on a better understanding of the spatial relationships involving the tumor, blood supply, and immune cell interactions, in the context of their associated microenvironments. LCM is used to isolate and obtain distinct histological cell types, which may be routinely performed on complex and heterogeneous solid tumor specimens. Once cells have been captured, nucleic acids and proteins may be extracted for in-depth multimodality molecular profiling assays. Optimizing the minute tissue quantities from LCM captured cells is challenging. Following the isolation of nucleic acids, RNA-seq may be performed for gene expression and DNA sequencing performed for the discovery and analysis of actionable mutations, copy number variation, methylation profiles, etc. However, there remains a need for highly sensitive proteomic methods targeting small-sized samples. A significant part of this protocol is an enhanced liquid chromatography mass spectrometry (LC-MS) analysis of micro-scale and/or nano-scale tissue sections. This is achieved with a silver-stained one-dimensional sodium dodecyl sulfate polyacrylamide gel electrophoresis (1D-SDS-PAGE) approach developed for LC-MS analysis of fresh-frozen tissue specimens obtained via LCM. Included is a detailed in-gel digestion method adjusted and specifically designed to maximize the proteome coverage from amount-limited LCM samples to better facilitate in-depth molecular profiling. Described is a proteogenomic approach leveraged from microdissected fresh frozen tissue. The protocols may also be applicable to other types of specimens having limited nucleic acids, protein quantity, and/or sample volume.

Proteogenomic Workflow for Characterization of Microphthalmia Transcription Factor (MiT) Family Translocation Renal Cell Carcinoma.

Microphthalmia transcription factor (MiT) family translocation renal cell carcinoma (tRCC) is a rare, aggressive, and heterogeneous subtype of kidney cancer, which is not well characterized. Since genetic alterations are always associated with carcinogenesis, and proteins are the major executors of biological features, multi-omics studies can reveal the systematic tRCC biological process comprehensively. Here, we describe the proteogenomic workflow for characterization of tRCC in detail to provide the knowledge foundation for integrated proteogenomic analysis of tRCC and other malignant tumors in the future.

Orthogonal proteogenomic analysis identifies the druggable PA2G4-MYC axis in 3q26 AML.

The overexpression of the ecotropic viral integration site-1 gene (EVI1/MECOM) marks the most lethal acute myeloid leukemia (AML) subgroup carrying chromosome 3q26 abnormalities. By taking advantage of the intersectionality of high-throughput cell-based and gene expression screens selective and pan-histone deacetylase inhibitors (HDACis) emerge as potent repressors of EVI1. To understand the mechanism driving on-target anti-leukemia activity of this compound class, here we dissect the expression dynamics of the bone marrow leukemia cells of patients treated with HDACi and reconstitute the EVI1 chromatin-associated co-transcriptional complex merging on the role of proliferation-associated 2G4 (PA2G4) protein. PA2G4 overexpression rescues AML cells from the inhibitory effects of HDACis, while genetic and small molecule inhibition of PA2G4 abrogates EVI1 in 3q26 AML cells, including in patient-derived leukemia xenografts. This study positions PA2G4 at the crosstalk of the EVI1 leukemogenic signal for developing new therapeutics and urges the use of HDACis-based combination therapies in patients with 3q26 AML.

A Proteogenomic Approach to Unravel New Proteins Encoded in the Leishmania donovani (HU3) Genome.

The high-throughput proteomics data generated by increasingly more sensible mass spectrometers greatly contribute to our better understanding of molecular and cellular mechanisms operating in live beings. Nevertheless, proteomics analyses are based on accurate genomic and protein annotations, and some information may be lost if these resources are incomplete. Here, we show that most proteomics data may be recovered by interconnecting genomics and proteomics approaches (i.e., following a proteogenomic strategy), resulting, in turn, in an improvement of gene/protein models. In this study, we generated proteomics data from Leishmania donovani (HU3 strain) promastigotes that allowed us to detect 1908 proteins in this developmental stage on the basis of the currently annotated proteins available in public databases. However, when the proteomics data were searched against all possible open reading frames existing in the L. donovani genome, twenty new protein-coding genes could be annotated. Additionally, 43 previously annotated proteins were extended at their N-terminal ends to accommodate peptides detected in the proteomics data. Also, different post-translational modifications (phosphorylation, acetylation, methylation, among others) were found to occur in a large number of Leishmania proteins. Finally, a detailed comparative analysis of the L. donovani and Leishmania major experimental proteomes served to illustrate how inaccurate conclusions can be raised if proteomes are compared solely on the basis of the listed proteins identified in each proteome. Finally, we have created data entries (based on freely available repositories) to provide and maintain updated gene/protein models. Raw data are available via ProteomeXchange with the identifier PXD051920.

Pan-cancer proteogenomics expands the landscape of therapeutic targets.

Fewer than 200 proteins are targeted by cancer drugs approved by the Food and Drug Administration (FDA). We integrate Clinical Proteomic Tumor Analysis Consortium (CPTAC) proteogenomics data from 1,043 patients across 10 cancer types with additional public datasets to identify potential therapeutic targets. Pan-cancer analysis of 2,863 druggable proteins reveals a wide abundance range and identifies biological factors that affect mRNA-protein correlation. Integration of proteomic data from tumors and genetic screen data from cell lines identifies protein overexpression- or hyperactivation-driven druggable dependencies, enabling accurate predictions of effective drug targets. Proteogenomic identification of synthetic lethality provides a strategy to target tumor suppressor gene loss. Combining proteogenomic analysis and MHC binding prediction prioritizes mutant KRAS peptides as promising public neoantigens. Computational identification of shared tumor-associated antigens followed by experimental confirmation nominates peptides as immunotherapy targets. These analyses, summarized at https://targets.linkedomics.org, form a comprehensive landscape of protein and peptide targets for companion diagnostics, drug repurposing, and therapy development.

Proteogenomic characterization of difficult-to-treat breast cancer with tumor cells enriched through laser microdissection.

BACKGROUND: Breast cancer (BC) is the most commonly diagnosed cancer and the leading cause of cancer death among women globally. Despite advances, there is considerable variation in clinical outcomes for patients with non-luminal A tumors, classified as difficult-to-treat breast cancers (DTBC). This study aims to delineate the proteogenomic landscape of DTBC tumors compared to luminal A (LumA) tumors. METHODS: We retrospectively collected a total of 117 untreated primary breast tumor specimens, focusing on DTBC subtypes. Breast tumors were processed by laser microdissection (LMD) to enrich tumor cells. DNA, RNA, and protein were simultaneously extracted from each tumor preparation, followed by whole genome sequencing, paired-end RNA sequencing, global proteomics and phosphoproteomics. Differential feature analysis, pathway analysis and survival analysis were performed to better understand DTBC and investigate biomarkers. RESULTS: We observed distinct variations in gene mutations, structural variations, and chromosomal alterations between DTBC and LumA breast tumors. DTBC tumors predominantly had more mutations in TP53, PLXNB3, Zinc finger genes, and fewer mutations in SDC2, CDH1, PIK3CA, SVIL, and PTEN. Notably, Cytoband 1q21, which contains numerous cell proliferation-related genes, was significantly amplified in the DTBC tumors. LMD successfully minimized stromal components and increased RNA-protein concordance, as evidenced by stromal score comparisons and proteomic analysis. Distinct DTBC and LumA-enriched clusters were observed by proteomic and phosphoproteomic clustering analysis, some with survival differences. Phosphoproteomics identified two distinct phosphoproteomic profiles for high relapse-risk and low relapse-risk basal-like tumors, involving several genes known to be associated with breast cancer oncogenesis and progression, including KIAA1522, DCK, FOXO3, MYO9B, ARID1A, EPRS, ZC3HAV1, and RBM14. Lastly, an integrated pathway analysis of multi-omics data highlighted a robust enrichment of proliferation pathways in DTBC tumors. CONCLUSIONS: This study provides an integrated proteogenomic characterization of DTBC vs LumA with tumor cells enriched through laser microdissection. We identified many common features of DTBC tumors and the phosphopeptides that could serve as potential biomarkers for high/low relapse-risk basal-like BC and possibly guide treatment selections.

Unraveling the adaptive strategies of Mycoplasma hominis through proteogenomic profiling of clinical isolates.

Introduction: Mycoplasma hominis (M. hominis) belongs to the class Mollicutes, characterized by a very small genome size, reduction of metabolic pathways, including transcription factors, and the absence of a cell wall. Despite this, they adapt well not only to specific niches within the host organism but can also spread throughout the body, colonizing various organs and tissues. The adaptation mechanisms of M. hominis, as well as their regulatory pathways, are poorly understood. It is known that, when adapting to adverse conditions, Mycoplasmas can undergo phenotypic switches that may persist for several generations. Methods: To investigate the adaptive properties of M. hominis related to survival in the host, we conducted a comparative phenotypic and proteogenomic analysis of eight clinical isolates of M. hominis obtained from patients with urogenital infections and the laboratory strain H-34. Results: We have shown that clinical isolates differ in phenotypic features from the laboratory strain, form biofilms more effectively and show resistance to ofloxacin. The comparative proteogenomic analysis revealed that, unlike the laboratory strain, the clinical isolates possess several features related to stress survival: they switch carbon metabolism, activating the energetically least advantageous pathway of nucleoside utilization, which allows slowing down cellular processes and transitioning to a starvation state; they reconfigure the repertoire of membrane proteins; they have integrative conjugative elements in their genomes, which are key mediators of horizontal gene transfer. The upregulation of the methylating subunit of the restriction-modification (RM) system type I and the additional components of RM systems found in clinical isolates suggest that DNA methylation may play a role in regulating the adaptation mechanisms of M. hominis in the host organism. It has been shown that based on the proteogenomic profile, namely the genome sequence, protein content, composition of the RM systems and additional subunits HsdM, HsdS and HsdR, composition and number of transposable elements, as well as the sequence of the main variable antigen Vaa, we can divide clinical isolates into two phenotypes: typical colonies (TC), which have a high growth rate, and atypical (aTC) mini-colonies, which have a slow growth rate and exhibit properties similar to persisters. Discussion: We believe that the key mechanism of adaptation of M. hominis in the host is phenotypic restructuring, leading to a slowing down cellular processes and the formation of small atypical colonies. This is due to a switch in carbon metabolism and activation the pathway of nucleoside utilization. We hypothesize that DNA methylation may play a role in regulating this switch.

Comprehensive proteogenomic characterization of rare kidney tumors.

Non-clear cell renal cell carcinomas (non-ccRCCs) encompass diverse malignant and benign tumors. Refinement of differential diagnosis biomarkers, markers for early prognosis of aggressive disease, and therapeutic targets to complement immunotherapy are current clinical needs. Multi-omics analyses of 48 non-ccRCCs compared with 103 ccRCCs reveal proteogenomic, phosphorylation, glycosylation, and metabolic aberrations in RCC subtypes. RCCs with high genome instability display overexpression of IGF2BP3 and PYCR1. Integration of single-cell and bulk transcriptome data predicts diverse cell-of-origin and clarifies RCC subtype-specific proteogenomic signatures. Expression of biomarkers MAPRE3, ADGRF5, and GPNMB differentiates renal oncocytoma from chromophobe RCC, and PIGR and SOSTDC1 distinguish papillary RCC from MTSCC. This study expands our knowledge of proteogenomic signatures, biomarkers, and potential therapeutic targets in non-ccRCC.

Integrated proteogenomic analysis for inherited bone marrow failure syndrome.

Recent advances in in-depth data-independent acquisition proteomic analysis have enabled comprehensive quantitative analysis of >10,000 proteins. Herein, an integrated proteogenomic analysis for inherited bone marrow failure syndrome (IBMFS) was performed to reveal their biological features and to develop a proteomic-based diagnostic assay in the discovery cohort; dyskeratosis congenita (n = 12), Fanconi anemia (n = 11), Diamond-Blackfan anemia (DBA, n = 9), Shwachman-Diamond syndrome (SDS, n = 6), ADH5/ALDH2 deficiency (n = 4), and other IBMFS (n = 18). Unsupervised proteomic clustering identified eight independent clusters (C1-C8), with the ribosomal pathway specifically downregulated in C1 and C2, enriched for DBA and SDS, respectively. Six patients with SDS had significantly decreased SBDS protein expression, with two of these not diagnosed by DNA sequencing alone. Four patients with ADH5/ALDH2 deficiency showed significantly reduced ADH5 protein expression. To perform a large-scale rapid IBMFS screening, targeted proteomic analysis was performed on 417 samples from patients with IBMFS-related hematological disorders (n = 390) and healthy controls (n = 27). SBDS and ADH5 protein expressions were significantly reduced in SDS and ADH5/ALDH2 deficiency, respectively. The clinical application of this first integrated proteogenomic analysis would be useful for the diagnosis and screening of IBMFS, where appropriate clinical screening tests are lacking.

Proteogenomic analysis identifies neoantigens and bacterial peptides as immunotherapy targets in colorectal cancer.

Considerable progress has recently been made in cancer immunotherapy, including immune checkpoint blockade, cancer vaccine, and adoptive T cell methods. The lack of effective targets is a major cause of the low immunotherapy response rate in colorectal cancer (CRC). Here, we used a proteogenomic strategy comprising immunopeptidomics, whole exome sequencing, and 16 S ribosomal DNA sequencing analyses of 8 patients with CRC to identify neoantigens and bacterial peptides that can serve as antitumor targets. This study directly identified several personalized neoantigens and bacterial immunopeptides. Immunoassays showed that all neoantigens and 5 of 8 bacterial immunopeptides could be recognized by autologous T cells. Additionally, T cell receptor (TCR) αß sequencing revealed the TCR repertoire of epitope-reactive CD8+ T cells. Functional studies showed that T cell receptor-T (TCR-T) could be activated by epitope pulsed lymphoblastoid cells. Overall, this study comprehensively profiled the CRC immunopeptidome, revealing several neoantigens and bacterial peptides with potential to serve as immunotherapy targets in CRC.

Understanding bacterial pathogen diversity: A proteogenomic analysis and use of an array of genome assemblies to identify novel virulence factors of the honey bee bacterial pathogen Paenibacillus larvae.

Mass spectrometry proteomics data are typically evaluated against publicly available annotated sequences, but the proteogenomics approach is a useful alternative. A single genome is commonly utilized in custom proteomic and proteogenomic data analysis. We pose the question of whether utilizing numerous different genome assemblies in a search database would be beneficial. We reanalyzed raw data from the exoprotein fraction of four reference Enterobacterial Repetitive Intergenic Consensus (ERIC) I-IV genotypes of the honey bee bacterial pathogen Paenibacillus larvae and evaluated them against three reference databases (from NCBI-protein, RefSeq, and UniProt) together with an array of protein sequences generated by six-frame direct translation of 15 genome assemblies from GenBank. The wide search yielded 453 protein hits/groups, which UpSet analysis categorized into 50 groups based on the success of protein identification by the 18 database components. Nine hits that were not identified by a unique peptide were not considered for marker selection, which discarded the only protein that was not identified by the reference databases. We propose that the variability in successful identifications between genome assemblies is useful for marker mining. The results suggest that various strains of P. larvae can exhibit specific traits that set them apart from the established genotypes ERIC I-V.

Proteogenomic insights into early-onset endometrioid endometrial carcinoma: predictors for fertility-sparing therapy response.

Endometrial carcinoma remains a public health concern with a growing incidence, particularly in younger women. Preserving fertility is a crucial consideration in the management of early-onset endometrioid endometrial carcinoma (EEEC), particularly in patients under 40 who maintain both reproductive desire and capacity. To illuminate the molecular characteristics of EEEC, we undertook a large-scale multi-omics study of 215 patients with endometrial carcinoma, including 81 with EEEC. We reveal an unexpected association between exposome-related mutational signature and EEEC, characterized by specific CTNNB1 and SIGLEC10 hotspot mutations and disruption of downstream pathways. Interestingly, SIGLEC10Q144K mutation in EEECs resulted in aberrant SIGLEC-10 protein expression and promoted progestin resistance by interacting with estrogen receptor alpha. We also identified potential protein biomarkers for progestin response in fertility-sparing treatment for EEEC. Collectively, our study establishes a proteogenomic resource of EEECs, uncovering the interactions between exposome and genomic susceptibilities that contribute to the development of primary prevention and early detection strategies for EEECs.

Proteogenomic Characterization Reveals Estrogen Signaling as a Target for Never-Smoker Lung Adenocarcinoma Patients without EGFR or ALK Alterations.

Never-smoker lung adenocarcinoma (NSLA) is prevalent in Asian populations, particularly in women. EGFR mutations and anaplastic lymphoma kinase (ALK) fusions are major genetic alterations observed in NSLA, and NSLA with these alterations have been well studied and can be treated with targeted therapies. To provide insights into the molecular profile of NSLA without EGFR and ALK alterations (NENA), we selected 141 NSLA tissues and performed proteogenomic characterization, including whole genome sequencing (WGS), transcriptomic, methylation EPIC array, total proteomic, and phosphoproteomic analyses. Forty patients with NSLA harboring EGFR and ALK alterations and seven patients with NENA with microsatellite instability were excluded. Genome analysis revealed that TP53 (25%), KRAS (22%), and SETD2 (11%) mutations and ROS1 fusions (14%) were the most frequent genetic alterations in NENA patients. Proteogenomic impact analysis revealed that STK11 and ERBB2 somatic mutations had broad effects on cancer-associated genes in NENA. DNA copy number alteration analysis identified 22 prognostic proteins that influenced transcriptomic and proteomic changes. Gene set enrichment analysis revealed estrogen signaling as the key pathway activated in NENA. Increased estrogen signaling was associated with proteogenomic alterations, such as copy number deletions in chromosomes 14 and 21, STK11 mutation, and DNA hypomethylation of LLGL2 and ST14. Finally, saracatinib, an Src inhibitor, was identified as a potential drug for targeting activated estrogen signaling in NENA and was experimentally validated in vitro. Collectively, this study enhanced our understanding of NENA NSLA by elucidating the proteogenomic landscape and proposed saracatinib as a potential treatment for this patient population that lacks effective targeted therapies. SIGNIFICANCE: The proteogenomic landscape in never-smoker lung cancer without known driver mutations reveals prognostic proteins and enhanced estrogen signaling that can be targeted as a potential therapeutic strategy to improve patient outcomes.

Proteogenomic Gene Structure Validation in the Pineapple Genome.

MD2 pineapple (Ananas comosus) is the second most important tropical crop that preserves crassulacean acid metabolism (CAM), which has high water-use efficiency and is fast becoming the most consumed fresh fruit worldwide. Despite the significance of environmental efficiency and popularity, until very recently, its genome sequence has not been determined and a high-quality annotated proteome has not been available. Here, we have undertaken a pilot proteogenomic study, analyzing the proteome of MD2 pineapple leaves using liquid chromatography-mass spectrometry (LC-MS/MS), which validates 1781 predicted proteins in the annotated F153 (V3) genome. In addition, a further 603 peptide identifications are found that map exclusively to an independent MD2 transcriptome-derived database but are not found in the standard F153 (V3) annotated proteome. Peptide identifications derived from these MD2 transcripts are also cross-referenced to a more recent and complete MD2 genome annotation, resulting in 402 nonoverlapping peptides, which in turn support 30 high-quality gene candidates novel to both pineapple genomes. Many of the validated F153 (V3) genes are also supported by an independent proteomics data set collected for an ornamental pineapple variety. The contigs and peptides have been mapped to the current F153 genome build and are available as bed files to display a custom gene track on the Ensembl Plants region viewer. These analyses add to the knowledge of experimentally validated pineapple genes and demonstrate the utility of transcript-derived proteomics to discover both novel genes and genetic structure in a plant genome, adding value to its annotation.

Integrated proteogenomic and metabolomic characterization of papillary thyroid cancer with different recurrence risks.

Although papillary thyroid cancer (PTC) has a good prognosis, its recurrence rate is high and remains a core concern in the clinic. Molecular factors contributing to different recurrence risks (RRs) remain poorly defined. Here, we perform an integrative proteogenomic and metabolomic characterization of 102 Chinese PTC patients with different RRs. Genomic profiling reveals that mutations in MUC16 and TERT promoter as well as multiple gene fusions like NCOA4-RET are enriched by the high RR. Integrative multi-omics analyses further describe the multi-dimensional characteristics of PTC, especially in metabolism pathways, and delineate dominated molecular patterns of different RRs. Moreover, the PTC patients are clustered into four subtypes (CS1: low RR and BRAF-like; CS2: high RR and metabolism type, worst prognosis; CS3: high RR and immune type, better prognosis; CS4: high RR and BRAF-like) based on the omics data. Notably, the subtypes display significant differences considering BRAF and TERT promoter mutations, metabolism and immune pathway profiles, epithelial cell compositions, and various clinical factors (especially RRs and prognosis) as well as druggable targets. This study can provide insights into the complex molecular characteristics of PTC recurrences and help promote early diagnosis and precision treatment of recurrent PTC.