Molecular Features Research Articles

Study on the Mass Spectrometry Cleavage Pattern of Quinolone Antibiotics.

Quinolone antibiotics are extensively used clinically for human treatment and in agriculture. However, improper and excessive use can lead to the persistence of quinolone residues in animal tissues, potentially accumulating in the human body and posing health risks. Investigating the correlation between mass spectrometry cleavage patterns and molecular structural features enhances the analytical framework for detecting trace or unknown impurities in quinolones. To collect data, we employed triple quadrupole linear ion trap mass spectrometry in electrospray positive ion mode. Primary mass spectrometry scanning was utilized to confirm parent ions, while secondary mass spectrometry scanning enabled the observation of fragment ions. The cleavage characteristics and pathways of the compounds were inferred from accurate mass-to-charge ratios obtained from both primary and secondary mass spectrometry. Under soft ionization conditions, the compounds generally exhibited characteristic fragment ions of [M+H-H2O]+, [M+H-CO]+, and [M+H-H2O-CO]+. Additionally, subtle variations were observed in each compound due to differences in modifying groups. For instance, upon deacidification, the piperazine ring structure underwent breakage and rearrangement, yielding fragment ion peaks devoid of neutral molecules such as C2H5N, C3H7N, or C4H8N. Notably, compounds featuring a cyclopropyl substituent group at the N-1 position typically exhibited characteristic fragments resulting from the loss of the cyclopropyl radical (⋅C3H5). Moreover, substituents at the N-1 and C-8 positions, when linked to form a six-membered carbocyclic ring, were prone to cleavage, releasing the neutral C3H6 molecule. Quinolone antibiotics share structural similarities in their parent nuclei, leading to partially similar cleavage pathways. Nevertheless, distinct cleavage patterns emerge due to variations in functional groups. According to the difference of mass spectrometry cleavage patterns, it can provide an identification basis for the measured detection of antibiotics.

Polatuzumab vedotin in diffuse large B-cell lymphoma therapy: Literature review and experience of the Oncohematology Division of the Lapino Clinical Hospital

Study of molecular and genetic features of diffuse large B-cell lymphoma (DLBCL) suggests differential approach to treatment of this aggressive B-cell disorder characterized by short remission after each subsequent therapy line. Therefore, more effective regimens should be used at earlier stages. Identification of biological markers and their integration into prognostic scales can help to personalize therapy, especially in patients with high risk of quick progression and increase survival.Our own experience of using immunoconjugate polatuzumab vedotin in therapy of patients with DLBCL. The use of Pola-BR combination in a patient with refractory non-GCB DLBCL allowed to achieve full remission alongside satisfactory tolerability and absence of significant adverse events. An evaluation of the Pola-R-CHP regimen in treatment of a patient with newly diagnosed generalized GCB DLBCL is presented.

Molecular mechanisms of circular RNA translation.

Circular RNAs (circRNAs) are covalently closed single-stranded RNAs without a 5' cap structure and a 3' poly(A) tail typically present in linear mRNAs of eukaryotic cells. CircRNAs are predominantly generated through a back-splicing process within the nucleus. CircRNAs have long been considered non-coding RNAs seemingly devoid of protein-coding potential. However, many recent studies have challenged this idea and have provided substantial evidence that a subset of circRNAs can associate with polysomes and indeed be translated. Therefore, in this review, we primarily highlight the 5' cap-independent internal initiation of translation that occurs on circular RNAs. Several molecular features of circRNAs, including the internal ribosome entry site, N6-methyladenosine modification, and the exon junction complex deposited around the back-splicing junction after back-splicing event, play pivotal roles in their efficient internal translation. We also propose a possible relationship between the translatability of circRNAs and their stability, with a focus on nonsense-mediated mRNA decay and nonstop decay, both of which are well-characterized mRNA surveillance mechanisms. An in-depth understanding of circRNA translation will reshape and expand our current knowledge of proteomics.

Revolutionary multi-omics analysis revealing prognostic signature of thyroid cancer and subsequent in vitro validation of SNAI1 in mediating thyroid cancer progression through EMT

Thyroid carcinoma (TC), the most commonly diagnosed malignancy of the endocrine system, has witnessed a significant rise in incidence over the past few decades. The integration of scRNA-seq with other sequencing approaches offers researchers a distinct perspective to explore mechanisms underlying TC progression. Therefore, it is crucial to develop a prognostic model for TC patients by utilizing a multi-omics approach. We acquired and processed transcriptomic data from the TCGA-THCA dataset, including mRNA expression profiles, lncRNA expression profiles, miRNA expression profiles, methylation chip data, gene mutation data, and clinical data. We constructed a tumor-related risk model using machine learning methods and developed a consensus machine learning-driven signature (CMLS) for accurate and stable prediction of TC patient outcomes. 2 strains of undifferentiated TC cell lines and 1 strain of PTC cell line were utilized for in vitro validation. mRNA, protein levels of hub genes, epithelial-mesenchymal transition (EMT)-associated phenotypes were detected by a series of in vitro experiments. We identified 3 molecular subtypes of TC based on integrated multi-omics clustering algorithms, which were associated with overall survival and displayed distinct molecular features. We developed a CMLS based on 28 hub genes to predict patient outcomes, and demonstrated that CMLS outperformed other prognostic models. TC patients of relatively lower CMLS score had significantly higher levels of T cells, B cells, and macrophages, indicating an immune-activated state. Fibroblasts were predominantly enriched in the high CMLS group, along with markers associated with immune suppression and evasion. We identified several drugs that could be suitable for patients with high CMLS, including Staurosporine_1034, Rapamycin_1084, gemcitabine, and topotecan. SNAI1 was elevated in both undifferentiated TC cell lines, comparing to PTC cells. Knockdown of SNAI1 reduced the cell proliferation and EMT phenotypes of undifferentiated TC cells. Our findings highlight the importance of multi-omics analysis in understanding the molecular subtypes and immune characteristics of TC, and provide a novel prognostic model and potential therapeutic targets for this disease. Moreover, we identified SNAI1 in mediating TC progression through EMT in vitro.

Models for the No-Observed-Effect Concentration (NOEC) and Maximal Half-Effective Concentration (EC50).

Typical in silico models for ecotoxicology focus on a few endpoints, but there is a need to increase the diversity of these models. This study proposes models using the NOEC for the harlequin fly (Chironomus riparius) and EC50 for swollen duckweed (Lemna gibba) for the first time. The data were derived from the EFSA OpenFoodTox database. The models were based on the correlation weights of molecular features used to calculate the 2D descriptor in CORAL software. The Monte Carlo method was used to calculate the correlation weights of the algorithms. The determination coefficients of the best models for the external validation set were 0.74 (NOAEC) and 0.85 (EC50).

TRIM35 Monoubiquitinates H2B in Cardiac Cells, Implications for Heart Failure.

The tumor suppressor and proapoptotic transcription factor P53 is induced (and activated) in several forms of heart failure, including cardiotoxicity and dilated cardiomyopathy; however, the precise mechanism that coordinates its induction with accessibility to its transcriptional promoter sites remains unresolved, especially in the setting of mature terminally differentiated (nonreplicative) cardiomyocytes. Male and female control or TRIM35 (tripartite motif containing 35) overexpression adolescent (aged 1-3 months) and adult (aged 4-6 months) transgenic mice were used for all in vivo experiments. Primary adolescent or adult mouse cardiomyocytes were isolated from control or TRIM35 overexpression transgenic mice for all in vitro experiments. Adenovirus or small-interfering RNA was used for all molecular experiments to overexpress or knockdown, respectively, target genes in primary mouse cardiomyocytes. Patient dilated cardiomyopathy or nonfailing left ventricle samples were used for translational and mechanistic insight. Chromatin immunoprecipitation and DNA sequencing or quantitative real-time polymerase chain reaction (qPCR) was used to assess P53 binding to its transcriptional promoter targets, and RNA sequencing was used to identify disease-specific signaling pathways. Here, we show that E3-ubiquitin ligase TRIM35 can directly monoubiquitinate lysine-120 (K120) on histone 2B in postnatal mature cardiomyocytes. This epigenetic modification was sufficient to promote chromatin remodeling, accessibility of P53 to its transcriptional promoter targets, and elongation of its transcribed mRNA. We found that increased P53 transcriptional activity (in cardiomyocyte-specific Trim35 overexpression transgenic mice) was sufficient to initiate heart failure and these molecular findings were recapitulated in nonischemic human LV dilated cardiomyopathy samples. These findings suggest that TRIM35 and the K120Ub-histone 2B epigenetic modification are molecular features of cardiomyocytes that can collectively predict dilated cardiomyopathy pathogenesis.

PARP-1 selectively impairs KRAS-driven phenotypic and molecular features in intrahepatic cholangiocarcinoma

ObjectiveIntrahepatic cholangiocarcinoma (iCCA) is the second most common primary liver cancer with limited therapeutic options. KRAS mutations are among the most abundant genetic alterations in iCCA associated with poor clinical...

Identification of a gene expression signature of vascular invasion and recurrence in stage I lung adenocarcinoma via bulk and spatial transcriptomics.

Microscopic vascular invasion (VI) is predictive of recurrence and benefit from lobectomy in stage I lung adenocarcinoma (LUAD) but is difficult to assess in resection specimens and cannot be accurately predicted prior to surgery. Thus, new biomarkers are needed to identify this aggressive subset of stage I LUAD tumors. To assess molecular and microenvironment features associated with angioinvasive LUAD we profiled 162 resected stage I tumors with and without VI by RNA-seq and explored spatial patterns of gene expression in a subset of 15 samples by high-resolution spatial transcriptomics (stRNA-seq). Despite the small size of invaded blood vessels, we identified a gene expression signature of VI from the bulk RNA-seq discovery cohort (n=103) and found that it was associated with VI foci, desmoplastic stroma, and high-grade patterns in our stRNA-seq data. We observed a stronger association with high-grade patterns from VI + compared with VI - tumors. Using the discovery cohort, we developed a transcriptomic predictor of VI, that in an independent validation cohort (n=60) was associated with VI (AUROC=0.86; p=5.42×10 -6 ) and predictive of recurrence-free survival (HR=1.98; p=0.024), even in VI - LUAD (HR=2.76; p=0.003). To determine our VI predictor's robustness to intra-tumor heterogeneity we used RNA-seq data from multi-region sampling of stage I LUAD cases in TRACERx, where the predictor scores showed high correlation (R=0.87, p<2.2×10 -16 ) between two randomly sampled regions of the same tumor. Our study suggests that VI-associated gene expression changes are detectable beyond the site of intravasation and can be used to predict the presence of VI. This may enable the prediction of angioinvasive LUAD from biopsy specimens, allowing for more tailored medical and surgical management of stage I LUAD.

Exploration of prognostic models for chronic rhinosinusitis with nasal polyps based on machine learning

Objective: To analysis the molecular characteristics of chronic rhinosinusitis with nasal polyps (CRSwNP), to unravel its pathophysiological mechanisms, and to develop a prognostic model capable of effectively predicting postoperative recurrence. Methods: The data from three datasets (GSE198950, GSE179265, and GSE136825) were integrated, comprising 39 control cases, 16 cases of chronic rhinosinusitis without nasal polyps, and 89 cases of CRSwNP. Differential expression genes (DEGs) were identified based on adjusted P<0.05 and Log2FC>1. KEGG and GO enrichment analyses, as well as STRING node scoring, were conducted. Variable selection was performed using random forest and least absolute shrinkage and selection operator regression (LASSO), with key nodes identified through intersection analysis. Mann-Whitney U test was applied, and variables with P<0.05 were included in the model. A prognostic model for CRSwNP was constructed using logistic regression, externally validated using RNA-seq data, and evaluated with receiver operating characteristic (ROC) curve analysis to calculate the area under the curve (AUC). Results: This research illustrated both upregulated and downregulated DEGs in CRSwNP, activating pathways like neuroactive ligand-receptor interaction and IL-17 signaling, while inhibiting calcium signaling and gap junctions. Key nodes identified through random forest and LASSO, including G protein subunit γ4 (U=3.00 P=0.028), Cholecystokinin (U=0.50, P=0.006), Epidermal growth factor (U=1.00 P=0.008), and Neurexin-1 (U=0.00, P=0.004), showing statistical significance in external validation. The prognostic model, visualized in a line graph, exhibited high reliability (C-index=0.875,AUC=0.866). The ROC curve in external validation indicated its effectiveness in predicting postoperative recurrence (AUC=0.859). Conclusions: This study integrates multiple datasets on CRSwNP to provide a comprehensive description of its molecular features. The prognostic model, built upon key nodes identified through random forest and LASSO analyses, demonstrates high accuracy in both internal and external validations, thus providing robust support for the development of personalized treatment strategies for CRSwNP.

Does the accounting of the local symmetry fragments in quasi-SMILES improve the predictive potential of the QSAR models of toxicity toward tadpoles?

Models of toxicity to tadpoles have been developed as single parameters based on special descriptors which are sums of correlation weights, molecular features, and experimental conditions. This information is presented by quasi-SMILES. Fragments of local symmetry (FLS) are involved in the development of the model and the use of FLS correlation weights improves their predictive potential. In addition, the index of ideality correlation (IIC) and correlation intensity index (CII) are compared. These two potential predictive criteria were tested in models built through Monte Carlo optimization. The CII was more effective than IIC for the models considered here.

Clinicopathological analysis of EWSR1/FUS::NFATC2 rearranged sarcoma in the left forearm: A case report

BACKGROUND We present a case of an EWSR1/FUS::NFATC2 rearranged sarcoma in the left forearm and analyze its clinicopathological and molecular features. CASE SUMMARY The patient is a 23-year-old woman. Microscopically, the tumor cells were medium-sized round cells arranged in small nests. The cytoplasm was clear, nuclei were relatively uniform, chromatin was dense, nucleoli were visible, and mitotic figures were rare. Immunohistochemically, the tumor cells were positive for Vimentin, INI-1, CD99, NKX2.2, CyclinD1, friend leukaemia virus integration 1, and NKX3.1. Next-generation sequencing revealed the presence of the EWSR1-NFATC2 fusion gene. EWSR1/FUS::NFATC2 rearranged sarcomas are rare and can easily be misdiagnosed. CONCLUSION Clinical imaging, immunohistochemistry, and molecular pathology should be considered to confirm the diagnosis.

Intrinsic epigenetic state of primary osteosarcoma drives metastasis.

Osteosarcoma (OS) is the most common primary malignant bone tumor affecting the pediatric population with high potential to metastasize. However, insights into the molecular features enabling its metastatic potential are limited. We mapped the active chromatin landscapes of OS tumors by integrating histone H3 lysine acetylated chromatin state (n=13), chromatin accessibility profiles (n=11) and gene expression (n=13) to understand the differences in their active chromatin profiles and its impact on molecular mechanisms driving the malignant phenotypes. Primary OS tumors from patients with metastasis (primary met) have a distinct active chromatin landscape compared to those without metastasis (localized). This difference shapes the transcriptional profile of OS. We identified novel candidate genes, including PPP1R1B, PREX1 and IGF2BP1, which exhibit increased chromatin activity in primary met. Loss of PREX1 in primary met OS cells significantly diminishes OS proliferation, invasion, migration, and colony formation capacity. Differential chromatin activity in primary met is associated with genes regulating cytoskeleton organization, cellular adhesion, and extracellular matrix suggestive of their role in facilitating OS metastasis. Chromatin profiling of tumors from metastatic lung lesions shows increased chromatin activity in genes involved in cell migration and Wnt pathway. This data demonstrates that metastatic potential is intrinsically present in primary metastatic tumors, with cellular chromatin profiles further adapting for successful dissemination, migration, and colonization at the distal site. Implications: Our study demonstrates that metastatic potential is intrinsic to primary metastatic osteosarcoma tumors, with chromatin profiles further adapting for successful dissemination, migration and colonization at distal metastatic site.

Free Salivary Amino Acid Profile in Breast Cancer: Clinicopathological and Molecular Biological Features.

The study of salivary amino acid profiles has attracted the attention of researchers, since amino acids are actively involved in most metabolic processes, including breast cancer. In this study, we analyzed the amino acid profile of saliva in a sample including all molecular biological subtypes of breast cancer to obtain a more complete picture and evaluate the potential utility of individual amino acids or their combinations for diagnostic purposes. This study included 116 patients with breast cancer, 24 patients with benign breast disease, and 25 healthy controls. From all patients, strictly before the start of treatment, saliva samples were collected, and the quantitative content of 26 amino acids was determined. Statistically significant differences between the three groups are shown in the content of Asp, Gly, Leu + Ile, Orn, Phe, Pro, Thr, and Tyr. To differentiate the three groups from each other, a decision tree was built. To construct it, we selected those amino acids for which the change in concentrations in the subgroups was multidirectional (GABA, Hyl, Arg, His, Pro, and Car). For the first time, it is shown that the amino acid profile of saliva depends on the molecular biological subtype of breast cancer. The most significant differences are shown for the luminal B HER2-positive and TNBC subgroups. In our opinion, it is critically important to consider the molecular biological subtype of breast cancer when searching for potential diagnostic markers.

Clinical and molecular features of patients with IDH1 wild-type primary glioblastoma presenting unexpected short-term survival after gross total resection.

This study investigated the factors influencing short-term survivors (STS) after gross total resection (GTR) in patients with IDH1 wild-type primary glioblastoma. We analyzed five independent cohorts who underwent GTR, including 83 patients from Kitasato University (K-cohort), and four validation cohorts of 148 patients from co-investigators (V-cohort), 66 patients from the Kansai Molecular Diagnosis Network for the Central Nervous System tumors, 109 patients from the Cancer Genome Atlas, and 40 patients from the Glioma Longitudinal AnalySiS. The study defined STS as those who had an overall survival ≤ 12months after GTR with subsequent radiation therapy, and concurrent and adjuvant temozolomide (TMZ). The study included 446 patients with glioblastoma. All cohorts experienced unexpected STS after GTR, with a range of 15.0-23.9% of the cases. Molecular profiling revealed no significant difference in major genetic alterations between the STS and non-STS groups, including MGMT, TERT, EGFR, PTEN, and CDKN2A. Clinically, the STS group had a higher incidence of non-local recurrence early in their treatment course, with 60.0% of non-local recurrence in the K-cohort and 43.5% in the V-cohort. The study revealed that unexpected STS after GTR in patients with glioblastoma is not uncommon and such tumors tend to present early non-local recurrence. Interestingly, we did not find any significant genetic alterations in the STS group, indicating that such major alterations are characteristics of GB rather than being reliable predictors for recurrence patterns or development of unexpected STS.

Endometrioid Endometrial RNA Index Predicts Recurrence in Stage I Patients.

Risk prediction with genomic and transcriptomic data has the potential to improve patient outcomes by enabling clinicians to identify patients requiring adjuvant treatment approaches, while sparing low-risk patients from unnecessary interventions. Endometrioid endometrial carcinoma (EEC) is the most common cancer in women in developed countries, and rates of endometrial cancer are increasing. We collected a 105-patient case-control cohort of stage I EEC comprised of 45 patients who experienced recurrence less than 6 years after excision, and 60 FIGO grade matched controls without recurrence. We first utilized two RNA based, previously validated machine learning approaches, namely EcoTyper and Complexity Index in Sarcoma (CINSARC). We developed Endometrioid Endometrial RNA Index (EERI) which uses RNA expression data from 46 genes to generate a personalized risk score for each patient. EERI was trained on our 105-patient cohort and tested on a publicly available cohort of 263 stage I EEC patients. EERI was able to predict recurrences with 94% accuracy in the training set and 81% accuracy in the test set. In the test set, patients assigned as EERI high-risk were significantly more likely to experience recurrence (30%) than the EERI low-risk group (1%) with a hazard ratio of 9.9 (95% CI 4.1-23.8, P <0.001). Tumors with high-risk genetic features may require additional treatment or closer monitoring and are not readily identified using traditional clinicopathologic and molecular features. EERI performs with high sensitivity and modest specificity, which may benefit from further optimization and validation in larger independent cohorts.

Hyperornithinemia-hyperammonemia-homocitrullinuria: a rare neurometabolic disorder in two siblings.

Hyperornithinemia-hyperammonemia-homocitrullinuria syndrome is an extremely rare disorder of urea cycle, with few patients reported worldwide. Despite hyperammonemia control, the long-term outcome remains poor with progressive neurological deterioration. We report the clinical, biochemical, and molecular features of two Lebanese siblings diagnosed with this disorder and followed for 8 and 15years, respectively. Variable clinical manifestations and neurological outcome were observed. The patient with earlier onset of symptoms had a severe neurological deterioration while the other developed a milder form of the disease at an older age. Diagnosis was challenging in the absence of the complete biochemical triad and the non-specific clinical presentations. Whole exome sequencing revealed a homozygous variant, p.Phe188del, in the SLC25A15 gene, a French- Canadian founder mutation previously unreported in Arab patients. Hyperammonemia was controlled in both patients but hyperonithinemia persisted. Frequent hyperalaninemia spikes and lactic acidosis occured concomitantly with the onset of seizures in one of the siblings. Variable neurological deterioration and outcome were observed within the same family. This is the first report from the Arab population of the long-term outcome of this devastating neurometabolic disorder.

Characterization of a novel TFG variant causing autosomal recessive pure hereditary spastic paraplegia.

TFG mutations have previously been implicated in autosomal recessive hereditary spastic paraplegia (HSP), also known as SPG57. This study aimed to investigate the clinical and molecular features of TFG mutations in a Taiwanese HSP cohort. Genetic analysis of TFG was conducted in 242 unrelated Taiwanese HSP patients using a targeted resequencing panel covering the entire coding regions of TFG. Functional assays were performed using an in vitro cell model to assess the impact of TFG variants on protein function. Additionally, other representative TFG mutant proteins were examined to understand the broader implications of TFG mutations in HSP. The study identified a novel homozygous TFG c.177A>C (p.(Lys59Asn)) variant in a family with adolescent-onset, pure form HSP. Functional analysis revealed that the Lys59Asn TFG variant, similar to other HSP-associated TFG mutants, exhibited a low affinity between TFG monomers and abnormal assembly of TFG homo-oligomers. These structural alterations led to aberrant intracellular distribution, compromising TFG's protein secretion function and resulting in decreased cellular viability. These findings confirm that the homozygous TFG c.177A>C (p.(Lys59Asn)) variant is a novel cause of SPG57. The study expands our understanding of the clinical and mutational spectrum of TFG-associated diseases, highlighting the functional defects associated with this specific TFG variant. Overall, this research contributes to the broader comprehension of the genetic and molecular mechanisms underlying HSP.

Integrated analysis of -omic landscapes in breast cancer subtypes

The subtypes of breast cancer exhibit diverse histology, molecular features, therapeutic response, aggressiveness, and patient outcomes. Multi-omics high-throughput technologies, which are widely used in cancer research, generated waste amounts of multimodal omic datasets calling for new approaches of integrated analyses to uncover patterns of transcriptomic, genomic, and epigenetic changes in breast cancer subtypes and connect them to disease clinical characteristics. Here, we applied multi-layer self-organizing map (ml-SOM) algorithms to PAM50-classified TCGA breast cancer samples to disentangle the diversity of the effects of gene expression, methylation, copy number, and somatic single nucleotide variation in the disease subtypes. Furthermore, we studied the association of perturbed gene modules with survival, prognosis, and other clinical characteristics. Our findings highlight the power of multi-omic analyses to offer a better understanding of the molecular diversity of breast cancer subtypes compared to single-omic analyses. Moreover, they highlight the complex subtype-characteristic associations between gene expression and epigenetic/genomic factors and their implications for survival and clinical outcomes.

Influence of Temperatures on Physicochemical Properties and Structural Features of Tamarind Seed Polysaccharide.

Due to the high content of impurities such as proteins in tamarind seed polysaccharide (TSP), they must be separated and purified before it can be used. TSP can disperse in cold water, but a solution can only be obtained by heating the mixture. Therefore, it is important to understand the dispersion and dissolution process of TSP at different temperatures to expand the application of TSP. In this study, pasting behavior and rheological properties as a function of temperature were characterized in comparison with potato starch (PS), and their relationship with TSP molecular features and microstructure was revealed. Pasting behavior showed that TSP had higher peak viscosity and stronger thermal stability than PS. Rheological properties exhibited that G' and G'' of TSP gradually increased with the increase in temperature, without exhibiting typical starch gelatinization behavior. The crystalline or amorphous structure of TSP and starch was disrupted under different temperature treatment conditions. The SEM results show that TSP particles directly transformed into fragments with the temperature increase, while PS granules first expanded and then broken down into fragments. Therefore, TSP and PS underwent different dispersion mechanisms during the dissolution process: As the temperature gradually increased, TSP possibly underwent a straightforward dispersion and was then dissolved in aqueous solution, while PS granules initially expanded, followed by disintegration and dispersion.

DUX4 is a common driver of immune evasion and immunotherapy failure in metastatic cancers

Cancer immune evasion contributes to checkpoint immunotherapy failure in many patients with metastatic cancers. The embryonic transcription factor DUX4 was recently characterized as a suppressor of interferon-γ signaling and antigen presentation that is aberrantly expressed in a small subset of primary tumors. Here, we report that DUX4 expression is a common feature of metastatic tumors, with ~10–50% of advanced bladder, breast, kidney, prostate, and skin cancers expressing DUX4. DUX4 expression is significantly associated with immune cell exclusion and decreased objective response to PD-L1 blockade in a large cohort of urothelial carcinoma patients. DUX4 expression is a significant predictor of survival even after accounting for tumor mutational burden and other molecular and clinical features in this cohort, with DUX4 expression associated with a median reduction in survival of over 1 year. Our data motivate future attempts to develop DUX4 as a biomarker and therapeutic target for checkpoint immunotherapy resistance.