Neoantigen-based immunotherapy has emerged as a transformative approach in cancer treatment, offering precision medicine strategies that target tumor-specific antigens derived from genetic, transcriptomic, and proteomic alterations unique to cancer cells. These neoantigens serve as highly specific targets for personalized therapies, promising more effective and tailored treatments. The aim of this article is to explore the advances in neoantigen-based therapies, highlighting successful treatments such as vaccines, tumor-infiltrating lymphocyte (TIL) therapy, T-cell receptor-engineered T cells therapy (TCR-T), and chimeric antigen receptor T cells therapy (CAR-T), particularly in cancer types like glioblastoma (GBM). Advances in technologies such as next-generation sequencing, RNA-based platforms, and CRISPR gene editing have accelerated the identification and validation of neoantigens, moving them closer to clinical application. Despite promising results, challenges such as tumor heterogeneity, immune evasion, and resistance mechanisms persist. The integration of AI-driven tools and multi-omic data has refined neoantigen discovery, while combination therapies are being developed to address issues like immune suppression and scalability. Additionally, the article discusses the ongoing development of personalized immunotherapies targeting tumor mutations, emphasizing the need for continued collaboration between computational and experimental approaches. Ultimately, the integration of cutting-edge technologies in neoantigen research holds the potential to revolutionize cancer care, offering hope for more effective and targeted treatments.

Gastric cancer (GC) remains among the top five global health problems. Therefore, comprehending the tumor energetic behavior is critical to understanding its progression. This study aimed to investigate mitochondrial DNA (mtDNA) alterations in GC cancer cell lines in an animal model.

Four mitochondrial genes (COI, ATP8, ND1, and ND3) were analyzed in GC (AGP01, ACP02, ACP03, and PG100) and control (Walker 256 carcinosarcoma) cell lines inoculated in Sapajus apella, exposed and not exposed to N-methyl-N-nitrosourea.

Two synonymous alterations were identified in ND1. In ND3, a non-synonymous alteration (A10398G ➔ Thr114Ala) may decrease the respiratory chain Complex I efficiency, enhancing cellular reactive oxygen species and contributing to mtDNA damage. As alterations in ND1 and ND3 were observed in highly aggressive cell lines, our results suggest these genes may play crucial roles in energetic efficiency and gastric carcinogenesis.

Alternative splicing (AS) is one of the most potent mechanisms for expanding the diversity of proteomes. During infection, human oncogenic viruses may exploit the AS to facilitate their replication cycle. Moreover, persistently infecting viruses can target key genes involved in classical signaling pathways to promote viral persistence and tumor progression. Here, we highlight how oncogenic viruses hijack AS system to manipulate host biological processes, and the host's AS system in turn modulates viral infection and replication. In addition, we have summarized the relatively underexplored involvement of noncoding RNAs in AS following tumor virus infection. This bidirectional interaction provides novel insights into interaction of virus-host and opens new avenues for therapeutic strategies targeting oncogenic viral infections.

Triple-negative breast cancer (TNBC) is an aggressive subtype with a poor prognosis. Although circular RNAs (circRNAs) have been implicated in cancer progression, their roles in TNBC remain poorly understood. In this study, we aimed to develop a prognostic model for TNBC by constructing a competing endogenous RNA (ceRNA) network. This network integrates circRNAs, long noncoding RNAs (lncRNAs), microRNAs (miRNAs), and messenger RNAs (mRNAs) to identify potential biomarkers and therapeutic targets for improving clinical outcomes.

Differentially expressed circRNAs, lncRNAs, and mRNAs were identified from GEO datasets (144 samples: 94 TNBC and 50 normal tissues). A ceRNA network was constructed, and key genes were validated using The Cancer Genome Atlas (TCGA) dataset (115 TNBC and 113 para-cancer tissues). Multivariate Cox regression analysis was performed to develop a prognostic model, and Gene Set Enrichment Analysis (GSEA) was performed to identify associated pathways.

Nine genes (SH3BGRL2, CA12, LRP8, NAV3, GFRA1, DCDC2, CDC7, ABAT, NPTX1) were identified as key factors in the prognostic model, which demonstrated an area under the curve (AUC) of 0.90. Patients classified as high-risk patients exhibited significantly shorter overall survival (median OS: 8.12 years vs. 9.51 years, P < 0.01). The mitogen-activated protein kinase (MAPK) signaling pathway was identified as a key regulatory pathway, with circRNAs (hsa_circ_0005455, hsa_circ_000632, hsa_circ_0001666, and hsa_circ_0000069) regulating CA12, GFRA1, and NPTX1 expression.

This study developed a novel prognostic model based on a ceRNA network analysis, highlighting the critical role of circRNAs and the MAPK signaling pathway in TNBC progression. These findings offer valuable insights into potential biomarkers for TNBC prognosis and reveal promising therapeutic targets for improving patient outcomes.

Nonstop extension or stop-loss mutations lead to the extension of a protein at its carboxyl terminus. Recently, nonstop mutations in the tumor suppressor SMAD Family Member 4 (SMAD4) have been discovered to lead to proteasomal SMAD4 degradation. However, this mutation type has not been studied in other cancer genes. Here, we explore somatic nonstop mutations in the tumor suppressor genes BRCA1 Associated Protein 1 (BAP1) and Von Hippel-Lindau (VHL) enriched in renal cell carcinoma. For BAP1, nonstop mutations generate an extremely long extension. Instead of proteasomal degradation, the extension decreases translation and depletes BAP1 messenger RNA from heavy polysomes. For VHL, the short extension leads to proteasomal degradation. Unexpectedly, the mutation alters the selection of the translational start site shifting VHL isoforms. We identify germline VHL nonstop mutations in patients leading to the early onset of severe disease manifestations. In summary, nonstop extension mutations inhibit the expression of renal tumor suppressor genes with pleiotropic effects on translation and protein stability.

The excision of introns from pre-mRNA is a crucial process in the expression of the majority of genes. Alternative splicing allows a single gene to generate diverse mRNA and protein products. Aberrant RNA splicing is recognized as a molecular characteristic present in almost all types of tumors. Therefore, identifying cancer-specific subtypes from aberrant processing offers new opportunities for therapeutic development. Numerous splicing modulators, each utilizing different mechanisms, have been developed as promising anticancer therapies, some of which are in clinical trials. In this review, we summarize the splice-altered signatures of cancer cell transcriptomes and the contributions of splicing aberrations to tumorigenesis and progression. Especially, we discuss current and emerging RNA splicing-targeted strategies for cancer therapy, including pharmacological approaches and splice-switching antisense oligonucleotides (ASOs). Finally, we address the challenges and opportunities in translating these findings into clinical practice.

Transforming growth factor-β (TGF-β) superfamily plays an important role in tumor progression and metastasis. Activin A receptor type 1C (ACVR1C) is a TGF-β type I receptor that is involved in tumorigenesis through binding to different ligands.

To evaluate the correlation between single nucleotide polymorphisms (SNPs) of ACVR1C and susceptibility to esophageal squamous cell carcinoma (ESCC) in Chinese Han population.

In this hospital-based cohort study, 1043 ESCC patients and 1143 healthy controls were enrolled. Five SNPs (rs4664229, rs4556933, rs77886248, rs77263459, rs6734630) of ACVR1C were assessed by the ligation detection reaction method. Hardy-Weinberg equilibrium test, genetic model analysis, stratified analysis, linkage disequilibrium test, and haplotype analysis were conducted.

Participants carrying ACVR1C rs4556933 GA mutant had significantly decreased risk of ESCC, and those with rs77886248 TA mutant were related with higher risk, especially in older male smokers. In the haplotype analysis, ACVR1C Trs4664229Ars4556933Trs77886248Crs77263459Ars6734630 increased risk of ESCC, while Trs4664229Grs4556933Trs77886248Crs77263459Ars6734630 was associated with lower susceptibility to ESCC.

ACVR1C rs4556933 and rs77886248 SNPs were associated with the susceptibility to ESCC, which could provide a potential target for early diagnosis and treatment of ESCC in Chinese Han population.

The choroid plexus (ChP) is the secretory epithelial structure located in the brain ventricles. Choroid plexus tumors (CPTs) are rare neoplasms predominantly occurring in young patients with intensified malignancy in children. CPT treatment is hindered by insufficient knowledge of tumor pathology and the limited availability of valid models.

Genomic and transcriptomic data from CPT patients were analyzed to identify the putative pathological pathway. Cellular and molecular techniques were employed to validate bioinformatic results in CPT patient samples. Pharmacologic inhibition of Wnt/β-catenin signaling was assessed in CPT cells. Cell-based assays of ChP cell lines were performed following CRISPR-Cas9-derived knockout and overexpression of Wnt/β-catenin pathway genes. A 3D CPT model was generated through CRISPR-Cas9-derived knockout of APC.

We discovered that Wnt/β-catenin signaling is activated in human CPTs, likely as a consequence of large-scale chromosomal instability events of the CPT genomes. We demonstrated that CPT-derived cells depend on autocrine Wnt/β-catenin signaling for survival. Constitutive Wnt/β-catenin pathway activation, either through knockout of the negative regulator APC or overexpression of the ligand WNT3A, induced tumorigenic properties in ChP 2D in vitro models. Increased activation of the Wnt/β-catenin pathway in ChP organoids, through treatment with a potent GSK3β inhibitor, reduced the differentiation of mature ChP epithelial cells. Remarkably, the depletion of APC was sufficient to induce the oncogenic transformation of ChP organoids.

Our research identifies Wnt/β-catenin signaling as a critical driver of CPT tumorigenesis and provides the first 3D in vitro model for future pathological and therapeutic studies of CPT.

Despite the widespread use of high-throughput genotyping strategies, certain mutation types remain understudied. We provide an overview of these often overlooked mutation types, with representative examples from common hereditary cancer syndromes.

We conducted a comprehensive review of the literature and locus-specific variant databases to summarize the germline pathogenic variants discovered through non-routine genotyping methods. We evaluated appropriate detection and analysis methods tailored for these specific genetic aberrations. Additionally, we performed in silico splice predictions on deep intronic variants registered in the ClinVar database.

Our study suggests that, aside from founder mutations, most cases are sporadic. However, we anticipate a relatively high likelihood of splice effects for deep intronic variants. The findings underscore the significant clinical utility of genome sequencing techniques and the importance of applying relevant analysis methods.

Colorectal cancer (CRC) ranks among the leading causes of cancer-related deaths worldwide. Enhancing CRC diagnosis and prognosis requires the development of improved biomarkers and therapeutic targets. Emerging evidence suggests that the unfolded protein response (UPR) plays a pivotal role in CRC progression, presenting new opportunities for diagnosis, treatment, and prevention. This study hypothesizes that genetic variants in endoplasmic reticulum (ER) stress response genes influence CRC susceptibility. We examined the frequencies of SNPs in PERK (rs13045) and GRP78/BiP (rs430397) within a South Iranian cohort. We mapped the cellular and molecular features of PERK and GRP78 genes in colorectal cancer, observing their differential expressions in tumor and metastatic tissues. We constructed co-expression and protein-protein interaction networks and performed gene set enrichment analysis, highlighting autophagy as a significant pathway through KEGG. Furthermore, the study included 64 CRC patients and 60 control subjects. DNA extraction and genotyping were conducted using high-resolution melting (HRM) analysis. Significant differences in PERK and GRP78 expressions were observed between CRC tissues and controls. Variations in PERK and GRP78 genotypes were significantly correlated with CRC risk. Utilizing a Multi-Target Directed Ligands approach, a dual PERK/GRP78 inhibitor was designed and subjected to molecular modeling studies. Docking experiments indicated high-affinity binding between the proposed inhibitor and both genes, PERK and GRP78, suggesting a novel therapy for CRC. These findings highlight the importance of understanding genetic backgrounds in different populations to assess CRC risk. Polymorphisms in UPR signaling pathway elements may serve as potential markers for predicting CRC susceptibility, paving the way for personalized therapeutic strategies.

We analyzed the impact of nine previously identified missense PKD1 variants from our studies, including c.6928G > A p.G2310R, c.8809G > A p.E2937K, c.2899 T > C p.W967R, c.6284A > G p.D2095G, c.6644G > A p.R2215Q, c.7810G > A p.D2604N, c.11249G > C p.R3750P, c.1001C > T p.T334M, and c.3101A > G p.N1034S on RNA structures and PC1 protein structure dynamics utilizing computational tools. RNA structure analysis was done using short RNA snippets of 41 nucleotides with the variant position at the 21st nucleotide, ensuring 20 bases on both sides. The secondary structures of these RNA snippets were predicted using RNAstructure. Structural changes of the mutants compared to the wild type were analyzed using the MutaRNA webserver. Molecular dynamics (MD) simulation of PC1 wild-type and mutant protein regions were performed using GROMACS 2018 (GROMOS96 54a7 force field). Findings revealed that five variants including c.8809G > A (p.E2937K), c.11249G > C (p.R3750P), c.3101A > G (p.N1034S), c.6928G > A (p.G2310R), c.6644G > A (p.R2215Q) exhibited major alterations in RNA structures and thereby their interactions with other proteins or RNAs affecting protein structure dynamics. While certain variants have minimal impact on RNA conformations, their observed alterations in MD simulations indicate impact on protein structure dynamics highlighting the importance of evaluating the functional consequences of genetic variants by considering both RNA and protein levels. The study also emphasizes that each missense variant exerts a unique impact on RNA stability, and protein structure dynamics, potentially contributing to the heterogeneous clinical manifestations and progression observed in Autosomal Dominant Polycystic Kidney Disease (ADPKD) patients offering a novel perspective in this direction. Thus, the utility of studying the structure dynamics through computational tools can help in prioritizing the variants for their functional implications, understanding the molecular mechanisms underlying variability in ADPKD presentation and developing targeted therapeutic interventions.

The online version contains supplementary material available at 10.1007/s13205-024-04057-9.

Mathematical models of biochemical reaction networks are an important and emerging tool for the study of cell signaling networks involved in disease processes. One promising potential application of such mathematical models is the study of how disease-causing mutations promote the signaling phenotype that contributes to the disease. It is commonly assumed that one must have a thorough characterization of the network readily available for mathematical modeling to be useful, but we hypothesized that mathematical modeling could be useful when there is incomplete knowledge and that it could be a tool for discovery that opens new areas for further exploration. In the present study, we first develop a mechanistic mathematical model of a G-protein coupled receptor signaling network that is mutated in almost all cases of uveal melanoma and use model-driven explorations to uncover and explore multiple new areas for investigating this disease. Modeling the two major, mutually-exclusive, oncogenic mutations (Gαq/11 and CysLT2R) revealed the potential for previously unknown qualitative differences between seemingly interchangeable disease-promoting mutations, and our experiments confirmed oncogenic CysLT2R was impaired at activating the FAK/YAP/TAZ pathway relative to Gαq/11. This led us to hypothesize that CYSLTR2 mutations in UM must co-occur with other mutations to activate FAK/YAP/TAZ signaling, and our bioinformatic analysis uncovers a role for co-occurring mutations involving the plexin/semaphorin pathway, which has been shown capable of activating this pathway. Overall, this work highlights the power of mechanism-based computational systems biology as a discovery tool that can leverage available information to open new research areas.

RNA splicing is pivotal in post-transcriptional gene regulation, yet the exponential expansion of intron length in humans poses a challenge for accurate splicing. Here, we identify hnRNPM as an essential RNA-binding protein that suppresses cryptic splicing through binding to deep introns, maintaining human transcriptome integrity. Long interspersed nuclear elements (LINEs) in introns harbor numerous pseudo splice sites. hnRNPM preferentially binds at intronic LINEs to repress pseudo splice site usage for cryptic splicing. Remarkably, cryptic exons can generate long dsRNAs through base-pairing of inverted ALU transposable elements interspersed among LINEs and consequently trigger an interferon response, a well-known antiviral defense mechanism. Significantly, hnRNPM-deficient tumors show upregulated interferon-associated pathways and elevated immune cell infiltration. These findings unveil hnRNPM as a guardian of transcriptome integrity by repressing cryptic splicing and suggest that targeting hnRNPM in tumors may be used to trigger an inflammatory immune response, thereby boosting cancer surveillance.

The translation of RNA by ribosomes represents a central biological process and one of the most dysregulated processes in cancer. While translation is traditionally thought to occur exclusively in the protein-coding regions of messenger RNAs (mRNAs), recent transcriptome-wide approaches have shown abundant ribosome activity across diverse stretches of RNA transcripts. The most common type of this kind of ribosome activity occurs in gene leader sequences, also known as 5' untranslated regions (UTRs) of the mRNA, that precede the main coding sequence. Translation of these upstream open reading frames (uORFs) is now known to occur in upwards of 25% of all protein-coding genes. With diverse functions from RNA regulation to microprotein generation, uORFs are rapidly igniting a new arena of cancer biology, where they are linked to cancer genetics, cancer signaling, and tumor-immune interactions. This review focuses on the contributions of uORFs and their associated 5'UTR sequences to cancer biology.

Understanding the function of rare non-coding variants represents a significant challenge. Using MapUTR, a screening method, we studied the function of rare 3' UTR variants affecting mRNA abundance post-transcriptionally. Among 17,301 rare gnomAD variants, an average of 24.5% were functional, with 70% in cancer-related genes, many in critical cancer pathways. This observation motivated an interrogation of 11,929 somatic mutations, uncovering 3928 (33%) functional mutations in 155 cancer driver genes. Functional MapUTR variants were enriched in microRNA- or protein-binding sites and may underlie outlier gene expression in tumors. Further, we introduce untranslated tumor mutational burden (uTMB), a metric reflecting the amount of somatic functional MapUTR variants of a tumor and show its potential in predicting patient survival. Through prime editing, we characterized three variants in cancer-relevant genes (MFN2, FOSL2, and IRAK1), demonstrating their cancer-driving potential. Our study elucidates the function of tens of thousands of non-coding variants, nominates non-coding cancer driver mutations, and demonstrates their potential contributions to cancer.

Molecular dynamics simulations play a pivotal role in elucidating the dynamic behaviors of RNA structures, offering a valuable complement to traditional methods such as nuclear magnetic resonance or X-ray. Despite this, the current precision of RNA force fields lags behind that of protein force fields. In this work, we systematically compared the performance of four RNA force fields (ff99bsc0χOL3, AMBERDES, ff99OL3_CMAP1, AMBERMaxEnt) across diverse RNA structures. Our findings highlight significant challenges in maintaining stability, particularly with regard to cross-strand and cross-loop hydrogen bonds. Furthermore, we observed the limitations in accurately describing the conformations of nonhelical structural motif, terminal nucleotides, and also base pairing and base stacking interactions by the tested RNA force fields. The identified deficiencies in existing RNA force fields provide valuable insights for subsequent force field development. Concurrently, these findings offer recommendations for selecting appropriate force fields in RNA simulations.

Cancer research has long relied on non-silent mutations. Yet, it has become overwhelmingly clear that silent mutations can affect gene expression and cancer cell fitness. One fundamental mechanism that apparently silent mutations can severely disrupt is alternative splicing. Here we introduce Oncosplice, a tool that scores mutations based on models of proteomes generated using aberrant splicing predictions. Oncosplice leverages a highly accurate neural network that predicts splice sites within arbitrary mRNA sequences, a greedy transcript constructor that considers alternate arrangements of splicing blueprints, and an algorithm that grades the functional divergence between proteins based on evolutionary conservation. By applying this tool to 12M somatic mutations we identify 8K deleterious variants that are significantly depleted within the healthy population; we demonstrate the tool's ability to identify clinically validated pathogenic variants with a positive predictive value of 94%; we show strong enrichment of predicted deleterious mutations across pan-cancer drivers. We also achieve improved patient survival estimation using a proposed set of novel cancer-involved genes. Ultimately, this pipeline enables accelerated insight-gathering of sequence-specific consequences for a class of understudied mutations and provides an efficient way of filtering through massive variant datasets - functionalities with immediate experimental and clinical applications.

Synonymous single nucleotide variants (sSNVs) do not alter the primary structure of a protein, thus it was previously accepted that they were neutral. Recently, several studies demonstrated their significance to a range of diseases. Still, variant prioritization strategies lack focus on sSNVs. Here, we identified 22,841 deleterious synonymous variants in 125,748 human exomes using two in silico predictors (SilVA and CADD). While 98.2% of synonymous variants are classified as neutral, 1.8% are predicted to be deleterious, yielding an average of 9.82 neutral and 0.18 deleterious sSNVs per exome. Further investigation of prediction features via Heterogeneous Ensemble Feature Selection revealed that impact on amino acid sequence and conservation carry the most weight for a deleterious prediction. Thirty nine detrimental sSNVs are not rare and are located on disease associated genes. Ten distinct putatively non-deleterious sSNVs are likely to be under positive selection in the North-Western European and East Asian populations. Taken together our analysis gives voice to the so-called silent mutations as we propose a robust framework for evaluating the deleteriousness of sSNVs in variant prioritization studies.

The TP53 3'UTR variant rs78378222 A>C has been detected in different tumor types as a somatic alteration that reduces p53 expression through modification of polyadenylation and miRNA regulation. Its prevalence is not yet known in all tumors. Herein, we examine tumor tissue prevalence of rs7837822 in Brazilian cohorts of patients from south and southeast regions diagnosed with lung adenocarcinoma (LUAD, n=586), sarcoma (SARC, n=188) and uterine leiomyoma (ULM, n=41). The minor allele (C) was identified in heterozygosity in 6/586 LUAD tumors (prevalence = 1.02 %) and none of the SARC and ULM samples. Additionally, next generation sequencing analysis revealed that all variant-positive tumors (n=4) with sample availability had additional pathogenic or likely pathogenic somatic variants in the TP53 coding regions. Among them, 3/4 (75 %) had the same pathogenic or likely pathogenic sequence variant (allele frequency <0.05 in tumor DNA) namely c.751A>C (p.Ile251Leu). Our results indicate a low somatic prevalence of rs78378222 in LUAD, ULM and SARC tumors from Brazilian patients, which suggests that no further analysis of this variant in the specific studied regions of Brazil is warranted. However, these findings should not exclude tumor molecular testing of this TP53 3'UTR functional variant for different populations.

Disrupted pre-mRNA splicing is a frequent deleterious mechanism in hereditary cancer. We aimed to functionally analyze candidate spliceogenic variants of the breast cancer susceptibility gene CHEK2 by splicing reporter minigenes.

A total of 128 CHEK2 splice-site variants identified in the Breast Cancer After Diagnostic Gene Sequencing (BRIDGES) project (https://cordis.europa.eu/project/id/634935) were analyzed with MaxEntScan and subsetted to 52 variants predicted to impact splicing. Three CHEK2 minigenes, which span all 15 exons, were constructed and validated. The 52 selected variants were then genetically engineered into the minigenes and assayed in MCF-7 (human breast adenocarcinoma) cells.

Of 52 variants, 46 (88.5%) impaired splicing. Some of them led to complex splicing patterns with up to 11 different transcripts. Thirty-four variants induced splicing anomalies without any trace or negligible amounts of the full-length transcript. A total of 89 different transcripts were annotated, which derived from different events: single- or multi-exon skipping, alternative site-usage, mutually exclusive exon inclusion, intron retention or combinations of the abovementioned events. Fifty-nine transcripts were predicted to introduce premature termination codons, 7 kept the original open-reading frame, 5 removed the translation start codon, 6 affected the 5'UTR (Untranslated Region), and 2 included missense variations. Analysis of variant c.684-2A > G revealed the activation of a non-canonical TG-acceptor site and exon 6 sequences critical for its recognition.

Incorporation of minigene read-outs into an ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based classification scheme allowed us to classify 32 CHEK2 variants (27 pathogenic/likely pathogenic and 5 likely benign). However, 20 variants (38%) remained of uncertain significance, reflecting in part the complex splicing patterns of this gene.

Investigating the human genome is vital for identifying risk factors and devising effective therapies to combat genetic disorders and cancer. Despite the extensive knowledge of the "light genome", the poorly understood "dark genome" remains understudied. In this study, we integrated data from 20,412 protein-coding genes in Pharos and 8,395 patient-derived tumours from The Cancer Genome Atlas (TCGA) to examine the genetic and pharmacological dependencies in human cancers and their treatment implications. We discovered that dark genes exhibited high mutation rates in certain cancers, similar to light genes. By combining the drug response profiles of cancer cells with cell fitness post-CRISPR-mediated gene knockout, we identified the crucial vulnerabilities associated with both dark and light genes. Our analysis also revealed that tumours harbouring dark gene mutations displayed worse overall and disease-free survival rates than those without such mutations. Furthermore, dark gene expression levels significantly influenced patient survival outcomes. Our findings demonstrated a similar distribution of genetic and pharmacological dependencies across the light and dark genomes, suggesting that targeting the dark genome holds promise for cancer treatment. This study underscores the need for ongoing research on the dark genome to better comprehend the underlying mechanisms of cancer and develop more effective therapies.

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is an X-linked enzymopathy caused by mutations in the G6PD gene. A medical concern associated with G6PD deficiency is acute hemolytic anemia induced by certain foods, drugs, and infections. Although phenotypic tests can correctly identify hemizygous males, as well as homozygous and compound heterozygous females, heterozygous females with a wide range of G6PD activity may be misclassified as normal. This study aimed to develop multiplex high-resolution melting (HRM) analyses to enable the accurate detection of G6PD mutations, especially among females with heterozygous deficiency. Multiplex HRM assays were developed to detect six G6PD variants, i.e., G6PD Gaohe (c.95A>G), G6PD Chinese-4 (c.392G>T), G6PD Mahidol (c.487G>A), G6PD Viangchan (c.871G>A), G6PD Chinese-5 (c.1024C>T), and G6PD Union (c.1360C>T) in two reactions. The assays were validated and then applied to genotype G6PD mutations in 248 Thai females. The sensitivity of the HRM assays developed was 100% [95% confidence interval (CI): 94.40%-100%] with a specificity of 100% (95% CI: 88.78%-100%) for detecting these six mutations. The prevalence of G6PD deficiency was estimated as 3.63% (9/248) for G6PD deficiency and 31.05% (77/248) for intermediate deficiency by phenotypic assay. The developed HRM assays identified three participants with normal enzyme activity as heterozygous for G6PD Viangchan. Interestingly, a deletion in intron 5 nucleotide position 637/638 (c.486-34delT) was also detected by the developed HRM assays. G6PD genotyping revealed a total of 12 G6PD genotypes, with a high prevalence of intronic variants. Our results suggested that HRM analysis-based genotyping is a simple and reliable approach for detecting G6PD mutations, and could be used to prevent the misdiagnosis of heterozygous females by phenotypic assay. This study also sheds light on the possibility of overlooking intronic variants, which could affect G6PD expression and contribute to enzyme deficiency.

The advantage of an increasing amount of bioinformatics data on leukemias intrigued us to explore the hot-spot mutation profiles and investigate the implications of those hot-spot mutations in patient survival. We retrieved somatic mutations and their distribution in protein domains through data analysis of The Cancer Genome Atlas and cBioPortal databases. After determining differentially expressed mutant genes related to leukemia, we further conducted principal component analysis and single-factor Cox regression analyses. Moreover, survival analysis was performed for the obtained candidate genes, followed by a multi-factor Cox proportional hazard model method for the impacts of the candidate genes on the survival and prognosis of patients with leukemia. At last, the signaling pathways involved in leukemia were investigated by gene set enrichment analysis. There were 223 somatic missense mutation hot-spots identified with pertinence to leukemia, which were distributed in 41 genes. Differential expression in leukemia was witnessed in 39 genes. We found a close correlation between seven genes and the prognosis of leukemia patients, among which, three genes could significantly influence the survival rate. In addition, among these three genes, CD74 and P2RY8 were highlighted due to close pertinence with survival conditions of leukemia patients. Finally, data suggested that B cell receptor, Hedgehog, and TGF-beta signaling pathways were enriched in low-hazard patients. In conclusion, these data underline the involvement of hot-spot mutations of CD74 and P2RY8 genes in survival status of leukemia patients, highlighting their as novel therapeutic targets or prognostic indicators for leukemia patients. Summary of Graphical Abstract: We identified 223 leukemia-associated somatic missense mutation hotspots concentrated in 41 different genes from 2297 leukemia patients in the TCGA database. Differential analysis of leukemic and normal samples from the TCGA and GTEx databases revealed that 39 of these 41 genes showed significant differential expression in leukemia. These 39 genes were subjected to PCA analysis, univariate Cox analysis, survival analysis, multivariate Cox regression analysis, GSEA pathway enrichment analysis, and then the association with leukemia survival prognosis and related pathways were investigated.

Hu sheep is an excellent short fat-tailed breed in China. Fat deposition in Hu sheep tail affects carcass quality and consumes a lot of energy, leading to an increase in feed cost. The objective of this study was to analyze the effects of HMGA1 polymorphism on tail fat weight in Hu sheep. Partial coding and non-coding sequences of HMGA1 were amplified with PCR and single nucleotide polymorphisms (SNP) of HMGA1 in 1163 Hu sheep were detected using DNA sequencing and KASPar technology. RT-qPCR analysis was performed to test HMGA1 expression in different tissues. The results showed that the expression of HMGA1 was higher in the duodenum, liver, spleen, kidney, and lung than in the heart, muscle, rumen, tail fat, and lymph. A mutation, g.5312 C > T, was detected in HMGA1; g.5312 C > T was significantly associated with tail fat weight, relative weight of tail fat (body weight), and relative weight of tail fat (carcass) (p < 0.05). The tail fat weight of the TT genotype was remarkably higher than that of the CC and TC genotypes. Therefore, HMGA1 can be used as a genetic marker for marker-assisted selection of tail fat weight in Hu sheep.

A number of data indicate that the sources of different kinds of PDAC may be discovered at the transcription/transduction stage. RNA metabolism is manipulated at various steps by different RNA-binding proteins (RBPs), and the deregulation or irregular activity of RBPs is known to contribute to tumor promotion and progression. The insulin-like growth factor 2 mRNA-binding protein family (IMPs), and IMP1 in particular, has been linked with a poor prognosis in PDAC patients; however, little is known about its contribution in PDAC carcinogenesis. In this study, we investigated the function of IMP1 in PDAC. To evaluate IMP1 expression and correlation with PDAC prognosis, we utilized several public databases. Using a specific siRNA IMP1, we analyzed cell death and cell cycle progression in PDAC cell lines and 3D spheroids. The role of IMP1 was also evaluated in vivo in a Panc-1-derived tumor xenograft murine model. Public data suggest that PDAC patients with higher expression of IMP1 showed poor overall and progression-free survival. IMP1 silencing leads to reduced cell growth in PDAC cells and three-dimensional spheroids. Abrogation of IMP1 in PDAC cells showed lower levels of CDC25A, increased phosphorylation of the cyclin-dependent kinase (CDK)2, and accumulation of PDAC cells in the G1 phase. Immunoprecipitation experiments revealed that IMP1 binds CDC25A mRNA, thus controlling cell-cycle progression. Ultimately, we proved that suppression of IMP1 blocked in vivo growth of Panc-1 transferred into immunodeficient mice. Our results indicate that IMP1 drives the PDCA cell cycle and represents a novel strategy for overcoming PDCA cell proliferation.

Epigenetic reprogramming represents a series of essential events during many cellular processes including oncogenesis. The genome of Kaposi's sarcoma-associated herpesvirus (KSHV), an oncogenic herpesvirus, is predetermined for a well-orchestrated epigenetic reprogramming once it enters into the host cell. The initial epigenetic reprogramming of the KSHV genome allows restricted expression of encoded genes and helps to hide from host immune recognition. Infection with KSHV is associated with Kaposi's sarcoma, multicentric Castleman's disease, KSHV inflammatory cytokine syndrome, and primary effusion lymphoma. The major epigenetic modifications associated with KSHV can be labeled under three broad categories: DNA methylation, histone modifications, and the role of noncoding RNAs. These epigenetic modifications significantly contribute toward the latent-lytic switch of the KSHV lifecycle. This review gives a brief account of the major epigenetic modifications affiliated with the KSHV genome in infected cells and their impact on pathogenesis.

Recent advances in neoantigen research have accelerated the development of tumor immunotherapies, including adoptive cell therapies (ACTs), cancer vaccines and antibody-based therapies, particularly for solid tumors. With the development of next-generation sequencing and bioinformatics technology, the rapid identification and prediction of tumor-specific antigens (TSAs) has become possible. Compared with tumor-associated antigens (TAAs), highly immunogenic TSAs provide new targets for personalized tumor immunotherapy and can be used as prospective indicators for predicting tumor patient survival, prognosis, and immune checkpoint blockade response. Here, the identification and characterization of neoantigens and the clinical application of neoantigen-based TCR-T immunotherapy strategies are summarized, and the current status, inherent challenges, and clinical translational potential of these strategies are discussed.

Keratoconus (KTCN) is the most common corneal ectasia resulting in a conical shape of the cornea. Here, genomic variation in the corneal epithelium (CE) across the keratoconic cone surface in patients with KTCN and its relevance in the functioning of the immune system were assessed.

Samples from four unrelated adolescent patients with KTCN and two control individuals were obtained during the CXL and PRK procedures, respectively. Three topographic regions, central, middle, and peripheral, were separated towards the whole-genome sequencing (WGS) study embracing a total of 18 experimental samples. The coding and non-coding sequence variation, including structural variation, was assessed and then evaluated together with the previously reported transcriptomic outcomes for the same CE samples and full-thickness corneas.

First, pathway enrichment analysis of genes with identified coding variants pointed to "Antigen presentation" and "Interferon alpha/beta signaling" as the most overrepresented pathways, indicating the involvement of inflammatory responses in KTCN. Both coding and non-coding sequence variants were found in genes (or in their close proximity) linked to the previously revealed KTCN-specific cellular components, namely, "Actin cytoskeleton", "Extracellular matrix", "Collagen-containing extracellular matrix", "Focal adhesion", "Hippo signaling pathway", and "Wnt signaling" pathways. No genomic heterogeneity across the corneal surface was found comparing the assessed topographic regions. Thirty-five chromosomal regions enriched in both coding and non-coding KTCN-specific sequence variants were revealed, with a most representative 5q locus previously recognized as involved in KTCN.

The identified genomic features indicate the involvement of innate and adaptive immune system responses in KTCN pathogenesis.

Chimeric antigen receptor (CAR) T cells and natural killer (NK) cells are genetically engineered immune cells that can detect target antigens on the surface of target cells and eliminate them following adoptive transfer. Recent progress in CAR-based therapies has led to outstanding clinical success in certain patients with leukemias and lymphomas and offered therapeutic benefits to those resistant to conventional therapies. The universal approach to stable CAR transgene delivery into the T/NK cells is the use of viral particles. Such approaches mediate semi-random transgene insertions spanning the entire genome with a high preference for integration into sites surrounding highly-expressed genes and active loci. Regardless of the variable CAR expression level based on the integration site of the CAR transgene, foreign integrated DNA fragments may affect the neighboring endogenous genes and chromatin structure and potentially change a transduced T/NK cell behavior and function or even favor cellular transformation. In contrast, site-specific integration of CAR constructs using recent genome-editing technologies could overcome the limitations and disadvantages of universal random gene integration. Herein, we explain random and site-specific integration of CAR transgenes in CAR-T/NK cell therapies. Also, we tend to summarize the methods for site-specific integration as well as the clinical outcomes of certain gene disruptions or enhancements due to CAR transgene integration. Also, the advantages and limitations of using site-specific integration methods are discussed in this review. Ultimately, we will introduce the genomic safe harbor (GSH) standards and suggest some appropriate safety prospects for CAR integration in CAR-T/NK cell therapies.

Metastasis, particularly brain metastasis, continues to puzzle researchers to this day, and exploring its molecular basis promises to break ground in developing new strategies for combatting this deadly cancer. In recent years, the research focus has shifted toward the earliest steps in the formation of metastasis. In this regard, significant progress has been achieved in understanding how the primary tumor affects distant organ sites before the arrival of tumor cells. The term pre-metastatic niche was introduced for this concept and encompasses all influences on sites of future metastases, ranging from immunological modulation and ECM remodeling to the softening of the blood-brain barrier. The mechanisms governing the spread of metastasis to the brain remain elusive. However, we begin to understand these processes by looking at the earliest steps in the formation of metastasis. This review aims to present recent findings on the brain pre-metastatic niche and to discuss existing and emerging methods to further explore the field. We begin by giving an overview of the pre-metastatic and metastatic niches in general before focusing on their manifestations in the brain. To conclude, we reflect on the methods usually employed in this field of research and discuss novel approaches in imaging and sequencing.

Retinoblastomas form in response to biallelic RB1 mutations or MYCN amplification and progress to more aggressive and therapy-resistant phenotypes through accumulation of secondary genomic changes. Progression-related changes include recurrent somatic copy number alterations and typically non-recurrent nucleotide variants, including synonymous and non-coding variants, whose significance has been unclear. To determine if nucleotide variants recurrently affect specific biological processes, we identified altered genes and over-represented variant gene ontologies in 168 exome or whole-genome-sequenced retinoblastomas and 12 tumor-matched cell lines. In addition to RB1 mutations, MYCN amplification, and established retinoblastoma somatic copy number alterations, the analyses revealed enrichment of variant genes related to diverse biological processes including histone monoubiquitination, mRNA processing (P) body assembly, and mitotic sister chromatid segregation and cytokinesis. Importantly, non-coding and synonymous variants increased the enrichment significance of each over-represented biological process term. To assess the effects of such mutations, we examined the consequences of a 3' UTR variant of PCGF3 (a BCOR-binding component of Polycomb repressive complex I), dual 3' UTR variants of CDC14B (a regulator of sister chromatid segregation), and a synonymous variant of DYNC1H1 (a regulator of P-body assembly). One PCGF3 and one of two CDC14B 3' UTR variants impaired gene expression whereas a base-edited DYNC1H1 synonymous variant altered protease sensitivity and stability. Retinoblastoma cell lines retained only ~50% of variants detected in tumors and enriched for new variants affecting p53 signaling. These findings reveal potentially important differences in retinoblastoma cell lines and tumors and implicate synonymous and non-coding variants, along with non-synonymous variants, in retinoblastoma oncogenesis.

CRC cells evolve a variety of strategies to limit or circumvent apoptosis cell death. RNA binding proteins (RBPs) regulate many of the molecular mechanisms that underlie the development of cancer. The insulin-like growth factor II mRNA-binding proteins (IMP) family are oncofoetal RBPs, consisting of IMP1, IMP2 and IMP3, which have an important role in RNA metabolism. IMP3 is highly expressed in colorectal cancer (CRC) tissue, where its expression often correlates with poor prognosis. However, the role of IMP3 in CRC is not fully understood. IMP3 expression was analysed using a public database and by Western blotting and immunohistochemistry in human colon samples derived from patients with sporadic CRC and healthy subjects. To address whether IMP3 controls cancer cell survival, we analysed cell death pathways in in vitro and in vivo experiments after IMP3 downregulation by siRNA or an antisense oligonucleotide. IMP3 was highly expressed in CRC samples compared to normal control tissues. The knockdown of IMP3 enhanced a caspase-independent cell death in CRC cell lines. Furthermore, the treatment of CRC cells with IMP3 siRNA did not alter the expression of GSDMD, GPX-4 and the activated form of RIP3, three key molecules that govern pyroptosis, ferroptosis and necroptosis, respectively. Abrogation of IMP3 in CRC significantly reduced Bcl-2 and Bcl-xL mRNA and was associated with an altered mitochondrial membrane potential that allowed the nuclear migration of the apoptosis-inducing factor (AIF). Moreover, specific immunoprecipitation experiments on CRC human cell lines indicated that IMP3 binds Bcl-2 and Bcl-xL mRNA, suggesting that IMP3 acts as a regulator of the intrinsic apoptotic pathway through the surveillance of anti-apoptotic Bcl mRNA metabolism. Finally, we showed that IMP3 block inhibited the growth of CRC cell lines in vivo after transplantation into immunodeficient mice. Altogether, these data support a novel role for IMP3 in controlling the intrinsic caspase-independent apoptotic pathway in CRC.

Deregulation of lncRNAs has been observed in human osteosarcoma. This study explored the diagnostic and prognostic significance of EPB41L4A-AS1 and UNC5B-AS1 in osteosarcoma.

Relative levels of EPB41L4A-AS1 and UNC5B-AS1 were detected in osteosarcoma tissue samples and cells. The ability to distinguish osteosarcoma from health was assessed by receiver operating characteristic (ROC) curve construction. Kaplan-Meier (K-M) and Cox proportional-hazards analyses were performed for prognosis factors. The bioinformatics approach was used to identify targeting miRNA for EPB41L4A-AS1 and UNC5B-AS1. Kaplan-Meier survival curves and Whitney Mann U tests were conducted for validating the statistical significance. In cell culture experiments, the influence of EPB41L4A-AS1 and UNC5B-AS1 on proliferation, migration, and invasion of the osteosarcoma cell line was examined by CCK-8 and Transwell assays.

Levels of EPB41L4A-AS1 and UNC5B-AS1 were upregulated in osteosarcoma patients and cells compared with the healthy participants and normal cell lines. EPB41L4A-AS1 and UNC5B-AS1 have a potent ability to distinguish the patients with osteosarcoma from the health. EPB41L4A-AS1 and UNC5B-AS1 levels correlated with SSS stage. Patients with high levels of EPB41L4A-AS1 and UNC5B-AS1 had significantly shorter survival times. EPB41L4A-AS1 and UNC5B-AS1 were independent prognostic indexes for overall survival. miR-1306-5p was a common target for EPB41L4A-AS1 and UNC5B-AS1. A propulsive impact on cell proliferation, migration, and invasion by EPB41L4A-AS1 and UNC5B-AS1 was observed, but can be rescued by miR-1306-5p.

It was concluded that upregulations of EPB41L4A-AS1 and UNC5B-AS1 expression were diagnostic and prognostic biomarkers for human osteosarcoma. EPB41L4A-AS1 and UNC5B-AS1 contribute to the biological behavior of osteosarcoma via miR-1306-5p.

Arginine methylation is a form of posttranslational modification that regulates many cellular functions such as development, DNA damage repair, inflammatory response, splicing, and signal transduction, among others. Protein arginine methyltransferase 5 (PRMT5) is one of nine identified methyltransferases, and it can methylate both histone and non-histone targets. It has pleiotropic functions, including recruitment of repair machinery to a chromosomal DNA double strand break (DSB) and coordinating the interplay between repair and checkpoint activation. Thus, PRMT5 has been actively studied as a cancer treatment target, and small molecule inhibitors of its enzymatic activity have already been developed. In this report, we analyzed all reported PRMT5 mutations appearing in cancer cells using data from the Catalogue of Somatic Mutations in Cancers (COSMIC). Our goal is to classify mutations as either drivers or passengers to understand which ones are likely to promote cellular transformation. Using gold standard artificial intelligence algorithms, we uncovered several key driver mutations in the active site of the enzyme (D306H, L315P, and N318K). In silico protein modeling shows that these mutations may affect the affinity of PRMT5 for S-adenosylmethionine (SAM), which is required as a methyl donor. Electrostatic analysis of the enzyme active site shows that one of these mutations creates a tunnel in the vicinity of the SAM binding site, which may allow interfering molecules to enter the enzyme active site and decrease its activity. We also identified several non-coding mutations that appear to affect PRMT5 splicing. Our analyses provide insights into the role of PRMT5 mutations in cancer cells. Additionally, since PRMT5 single molecule inhibitors have already been developed, this work may uncover future directions in how mutations can affect targeted inhibition.

Cancer is among the leading causes of mortality worldwide. While considerable attention has been given to genetic and epigenetic sources of cancer-specific cellular activities, the role of alternative mRNA splicing has only recently received attention as a major contributor to cancer initiation and progression. The distribution of alternate mRNA splicing variants in cancer cells is different from their non-cancer counterparts, and cancer cells are more sensitive than non-cancer cells to drugs that target components of the splicing regulatory network. While many of the alternatively spliced mRNAs in cancer cells may represent "noise" from splicing dysregulation, certain recurring splicing variants have been shown to contribute to tumor progression. Some pathogenic splicing disruption events result from mutations in cis-acting splicing regulatory sequences in disease-associated genes, while others may result from shifts in balance among naturally occurring alternate splicing variants among mRNAs that participate in cell cycle progression and the regulation of apoptosis. This review provides examples of cancer-related alternate splicing events resulting from each step of mRNA processing and the promising therapies that may be used to address them.

Currently, colorectal cancer (CRC) represents the third most common malignancy and the second most deadly cancer worldwide, with a higher incidence in developed countries. Like other solid tumors, CRC is a heterogeneous genomic disease in which various alterations, such as point mutations, genomic rearrangements, gene fusions or chromosomal copy number alterations, can contribute to the disease development. However, because of its orderly natural history, easily accessible onset location and high lifetime incidence, CRC is ideally suited for preventive intervention, but the many screening efforts of the last decades have been compromised by performance limitations and low penetrance of the standard screening tools. The advent of next-generation sequencing (NGS) has both facilitated the identification of previously unrecognized CRC features such as its relationship with gut microbial pathogens and revolutionized the speed and throughput of cataloguing CRC-related genomic alterations. Hence, in this review, we summarized the several diagnostic tools used for CRC screening in the past and the present, focusing on recent NGS approaches and their revolutionary role in the identification of novel genomic CRC characteristics, the advancement of understanding the CRC carcinogenesis and the screening of clinically actionable targets for personalized medicine.

Synonymous variations have always been ignored while studying the underlying genetic mechanisms for most of the human diseases. However, recent studies have suggested that these silent changes in the genome can alter the protein expression and folding.

CSRP3, which is a well-known candidate gene associated with dilated cardiomyopathy (DCM) and hypertrophic cardiomyopathy (HCM), was screened for 100 idiopathic DCM cases and 100 controls. Three synonymous variations were identified viz., c.96G > A, p.K32=; c.336G > A, p.A112=; c.354G > A, p.E118=. A comprehensive in silico analysis was performed using various web based widely accepted tools, Mfold, Codon Usage, HSF3.1 and RNA22. Mfold predicted structural changes in all the variants except c.96 G > A (p.K32=), however it predicted changes in the stability of mRNA due to all the synonymous variants. Codon bias was observed as evident by the Relative Synonymous Codon Usage and Log Ratio of Codon Usage Frequencies. The Human Splicing Finder also predicted remarkable changes in the regulatory elements in the variants c.336G > A and c.354 G > A. The miRNA target prediction using varied modes available in RNA22 revealed that 70.6% of the target sites of miRNAs in CSRP3 were altered due to variant c.336G > A while 29.41% sites were completely lost.

Findings of the present study suggest that synonymous variants revealed striking deviations in the structural conformation of mRNA, stability of mRNA, relative synonymous codon usage, splicing and miRNA binding sites from the wild type suggesting their possible role in the pathogenesis of DCM, either by destabilizing the mRNA structure, or codon usage bias or else altering the cis-acting regulatory elements during splicing.