Somatic mutation of DNAH genes implicated higher chemotherapy response rate in gastric adenocarcinoma patients

Machine learning-driven estimation of mutational burden highlights DNAH5 as a prognostic marker in colorectal cancer

BACKGROUND

Tumor Mutational Burden (TMB) have emerged as pivotal predictive biomarkers in determining prognosis and response to immunotherapy in colorectal cancer (CRC) patients. While Whole Exome Sequencing (WES) stands as the gold standard for TMB assessment, carry substantial costs and demand considerable time commitments. Additionally, the heterogeneity among high-TMB patients remains poorly characterized.

METHODS

We employed eight advanced machine learning algorithms to develop gene-panel-based models for TMB estimation. To rigorously compare and validate these TMB estimation models, four external cohorts, involving 1,956 patients, were used. Furthermore, we computed the Pearson correlation coefficient between the estimated TMB and tumor neoantigen levels to elucidate their association. CD8+ tumor-infiltrating lymphocyte (TIL) density was assessed via immunohistochemistry.

RESULTS

The TMB estimation model based on the Lasso algorithm, incorporating 20 genes, exhibiting satisfactory performance across multiple independent cohorts (R2 ≥ 0.859). This 20-gene TMB model proved to be an independent prognostic indicator for the progression-free survival (PFS) of CRC patients (p = 0.001). DNAH5 mutations were associated with a more favorable prognosis in high-TMB CRC patients, and correlated strongly with tumor neoantigen levels and CD8+ TIL density.

CONCLUSIONS

The 20-gene model offers a cost-efficient approach to precisely estimating TMB, providing prognosis in patients with CRC. Incorporating DNAH5 within this model further refines the categorization of patients with elevated TMB. Utilizing the 20-gene model facilitates the stratification of patients with CRC, enabling more precise treatment planning.

Identification of competing endogenous RNA networks associated with circRNA and lncRNA in TCDD-induced cleft palate development

2,3,7,8 -tetrachlorodibenzo-p-dioxin (TCDD) is a teratogen that can induce cleft palate formation, a common birth defect. Competing endogenous RNAs (ceRNAs), including circular RNAs (circRNAs) and long non-coding RNAs (lncRNAs), indirectly regulate gene expression via sharing microRNAs (miRNAs). Nevertheless, the mechanism by which they act as ceRNAs to regulate palatal development remains to be explored in greater detail. Here, the cleft palate model of C57BL/6 N pregnant mice was constructed by gavage of TCDD (64 ug/kg) on gestation day (GD) 10.5, and the palatal shelves were taken on gestation day (GD) 14.5 for whole-transcriptome sequencing to investigate the underlying mechanisms of the roles of circRNAs and lncRNAs as ceRNAs in cleft palate. Sequencing results revealed that 293 lncRNA, 589 circRNA, 47 miRNA, and 138 messenger RNA (mRNA) were significantly dysregulated, and the cytochrome P450 (CYP) enzymes and the aryl hydrocarbon receptor (AhR) pathway play key roles in the induction of cleft palate upon exposure to TCDD. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed the function of TCDD function was mainly related to the metabolic processes of intracellular compounds, including the metabolic processes of cellular aromatic compounds and the metabolism of exogenous drugs by cytochrome P450, etc. Furthermore, quantitative reverse transcription polymerase chain reaction (qRT-PCR) indicated that the circRNA_1781/miR-30c-1-3p/PKIB and XR_380026.2/miR-1249-3p/DNAH10 ceRNA networks were hypothesized to be a hub involved in palatal development suggesting that the circRNA_1781/miR-30c-1-3p/PKIB and XR_380026.2/miR-1249-3p/DNAH10 ceRNA networks may be critical for palatogenesis, setting the foundation for the investigation of cleft palate.

CAN008 prolongs overall survival in patients with newly diagnosed GBM characterized by high tumor mutational burden

BACKGROUND

A previous phase 1 dose-escalation study in Taiwan indicated CAN008 (asunercept) with standard concurrent chemoradiotherapy (CCRT) improved progression-free survival (PFS) in newly diagnosed glioblastoma (GBM) patients. This study evaluates the efficacy of CAN008 in promoting overall survival (OS) and identifies genetic alterations associated with treatment responses.

METHODS

We compared OS of 5-year follow-ups from 9 evaluable CAN008 cohort patients (6 received high-dose and 3 received low-dose) to a historical Taiwanese GBM cohort with 164 newly diagnosed patients. CAN008 treatment response-associated genetic alterations were identified by whole-exome sequencing and comparing variant differences between response groups. Associations among patient survival, tumor mutational burden (TMB), and genetic alterations were analyzed using CAN008 cohort and TCGA-GBM dataset.

RESULTS

OS for high-dose CAN008 patients at 2 and 5 years was 83% and 67%, respectively, and 40.1% and 8.8% for the historical GBM cohort, respectively. Better OS was observed in the high-dose CAN008 cohort (without reaching the median survival) than the historical GBM cohort (median OS: 20 months; p = 0.0103). Five high-dose CAN008 patients were divided into good and poor response groups based on their PFS. A higher variant count and TMB were observed in good response patients, whereas no significant association was observed between TMB and patient survival in the newly diagnosed TCGA-GBM dataset, suggesting TMB may modulate patient CAN008 response.

CONCLUSION

CAN008 combined with standard CCRT treatment prolonged the PFS and OS of newly diagnosed GBM patients compared to standard therapy alone. Higher treatment efficacy was associated with higher TMB.

Quantitative proteomics of sperm tail in asthenozoospermic patients: exploring the molecular pathways affecting sperm motility

Asthenozoospermia, characterized by low sperm motility, is one of the most common causes of male infertility. While many intrinsic and extrinsic factors are involved in the etiology of asthenozoospermia, the molecular basis of this condition remains unclear. Since sperm motility results from a complex flagellar structure, an in-depth proteomic analysis of the sperm tail can uncover mechanisms underlying asthenozoospermia. This study quantified the proteomic profile of 40 asthenozoospermic sperm tails and 40 controls using TMT-LC-MS/MS. Overall, 2140 proteins were identified and quantified where 156 proteins have not been described earlier in sperm tail. There were 409 differentially expressed proteins (250 upregulated and 159 downregulated) in asthenozoospermia which by far is the highest number reported earlier. Further, bioinformatics analysis revealed several biological processes, including mitochondrial-related energy production, oxidative phosphorylation (OXPHOS), citric acid cycle (CAC), cytoskeleton, stress response, and protein metabolism altered in asthenozoospermic sperm tail samples. Collectively, our findings reveal the importance of mitochondrial energy production and induced stress response as potential mechanisms involved in the loss of sperm motility in asthenozoospermia.

A Novel Four Mitochondrial Respiration-Related Signature for Predicting Biochemical Recurrence of Prostate Cancer

The biochemical recurrence (BCR) of patients with prostate cancer (PCa) after radical prostatectomy is high, and mitochondrial respiration is reported to be associated with the metabolism in PCa development. This study aimed to establish a mitochondrial respiratory gene-based risk model to predict the BCR of PCa. RNA sequencing data of PCa were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, and mitochondrial respiratory-related genes (MRGs) were sourced via GeneCards. The differentially expressed mitochondrial respiratory and BCR-related genes (DE-MR-BCRGs) were acquired through overlapping BCR-related differentially expressed genes (BCR-DEGs) and differentially expressed MRGs (DE-MRGs) between PCa samples and controls. Further, univariate Cox, least absolute shrinkage and selection operator (LASSO), and multivariate Cox analyses were performed to construct a DE-MRGs-based risk model. Then, a nomogram was established by analyzing the independent prognostic factor of five clinical features and risk scores. Moreover, Gene Set Enrichment Analysis (GSEA), tumor microenvironment, and drug susceptibility analyses were employed between high- and low-risk groups of PCa patients with BCR. Finally, qRT-PCR was utilized to validate the expression of prognostic genes. We identified 11 DE-MR-BCRGs by overlapping 132 DE-MRGs and 13 BCR-DEGs and constructed a risk model consisting of 4 genes (APOE, DNAH8, EME2, and KIF5A). Furthermore, we established an accurate nomogram, including a risk score and a Gleason score, for the BCR prediction of PCa patients. The GSEA result suggested the risk model was related to the PPAR signaling pathway, the cholesterol catabolic process, the organic hydroxy compound biosynthetic process, the small molecule catabolic process, and the steroid catabolic process. Simultaneously, we found six immune cell types relevant to the risk model: resting memory CD4+ T cells, monocytes, resting mast cells, activated memory CD4+ T cells, regulatory T cells (Tregs), and macrophages M2. Moreover, the risk model could affect the IC50 of 12 cancer drugs, including Lapatinib, Bicalutamide, and Embelin. Finally, qRT-PCR showed that APOE, EME2, and DNAH8 were highly expressed in PCa, while KIF5A was downregulated in PCa. Collectively, a mitochondrial respiratory gene-based nomogram including four genes and one clinical feature was established for BCR prediction in patients with PCa, which could provide novel strategies for further studies.

Clinical detection, diagnosis and treatment of morphological abnormalities of sperm flagella: A review of literature

Sperm carries male genetic information, and flagella help move the sperm to reach oocytes. When the ultrastructure of the flagella is abnormal, the sperm is unable to reach the oocyte and achieve insemination. Multiple morphological abnormalities of sperm flagella (MMAF) is a relatively rare idiopathic condition that is mainly characterized by multiple defects in sperm flagella. In the last decade, with the development of high-throughput DNA sequencing approaches, many genes have been revealed to be related to MMAF. However, the differences in sperm phenotypes and reproductive outcomes in many cases are attributed to different pathogenic genes or different pathogenic mutations in the same gene. Here, we will review information about the various phenotypes resulting from different pathogenic genes, including sperm ultrastructure and encoding proteins with their location and functions as well as assisted reproductive technology (ART) outcomes. We will share our clinical detection and diagnosis experience to provide additional clinical views and broaden the understanding of this disease.

Cross-species identification of cancer resistance-associated genes that may mediate human cancer risk

Cancer is a predominant disease across animals. We applied a comparative genomics approach to systematically characterize genes whose conservation levels correlate positively (PC) or negatively (NC) with cancer resistance estimates across 193 vertebrates. Pathway analysis reveals that NC genes are enriched for metabolic functions and PC genes in cell cycle regulation, DNA repair, and immune response, pointing to their corresponding roles in mediating cancer risk. We find that PC genes are less tolerant to loss-of-function (LoF) mutations, are enriched in cancer driver genes, and are associated with germline mutations that increase human cancer risk. Their relevance to cancer risk is further supported via the analysis of mouse functional genomics and cancer mortality of zoo mammals' data. In sum, our study describes a cross-species genomic analysis pointing to candidate genes that may mediate human cancer risk.

Somatic Mutation of FAT Family Genes Implicated Superior Prognosis in Patients With Stomach Adenocarcinoma

FAT family genes encode protocadherin, which regulates tumor cell proliferation and migration. Although transcriptional levels of FAT family members had been reported in multiple malignant tumors, the association between mutation and prognosis of the FAT family in stomach adenocarcinoma (STAD) has not been investigated. Herein, we performed a multi-omics integrative bioinformatics analysis using genomic and mRNA expression data to explore the role of gene mutations across the FAT family on clinical outcomes of STAD. The results showed that FAT mutations occurred in 174 of 435 (40%) of the samples. Patients with FAT mutations possessed significantly better progression-free survival (P = 0.019) and overall survival (P = 0.034) than those with non-FAT mutations, and FAT mutations exhibited significantly higher tumor mutational burden (TMB) and microsatellite instability. Notably, FAT mutations had a greater effect on somatic single-nucleotide variation than copy number variation and resulted in more abundant DNA damage repair (DDR) mutations. Further investigation demonstrated that FAT mutations contributed to an inflammatory tumor microenvironment (TME), as indicated by significantly increased numbers of activated CD4 and CD8 T cells, and significantly decreased numbers of mast cell, plasmacytoid dendritic cell, type 2 T helper cell, and high expression of immune-promoting genes. Moreover, biological process antigen processing and presentation, DNA replication, and DDR-related pathways were significantly upregulated in patients with FAT mutations. Collectively, FAT mutations significantly improved the survival of patients with STAD by enhancing tumor immunogenicity (e.g., TMB and DDR mutations) and an inflamed TME, indicating that the FAT family might be a potential prognostic and therapeutic biomarker for STAD.

Impact of neoadjuvant chemotherapy on somatic mutation status in high-grade serous ovarian carcinoma

BACKGROUND

Patients treated with neoadjuvant chemotherapy (NACT) for advanced high-grade serous ovarian carcinoma (HGSC) have a higher rate and shorter time to platinum-resistant recurrence compared to patients treated with primary cytoreductive surgery (PCS) and adjuvant chemotherapy. The purpose of this study is to determine the impact of NACT on somatic mutation status in platinum-sensitive and resistant HGSC. Patients with advanced HGSC who had a documented response to platinum-based NACT, a banked blood sample, and a banked tumor sample before and after NACT were identified. Whole exome and/or targeted deep sequencing was performed in matched normal and pre/post-NACT tumor samples from 3 platinum-resistant and 2 platinum-sensitive patients to identify somatic non-synonymous mutations at each time point.

RESULTS

When comparing exonic non-synonymous mutations in pre-NACT and post-NACT samples from the same patient, an average of 41% (1-68%) of genes were mutated at both time points. There were no trends detected in the mutational burden following exposure to NACT in platinum-resistant vs. platinum-sensitive cases. The majority of mutated genes were unique to each case. We identified several genes that were commonly mutated in pre-NACT samples specific to platinum-resistant (CSPG4, SLC35G5, TUBA3D) or sensitive (CYP2D6, NUTM1, DNAH5) cases. Four mutated genes emerged exclusively in the platinum-resistant cases (ADGRV1, MUC17, MUC20, PAK2) following NACT.

CONCLUSIONS

Patients with advanced HGSC present with significant intra-tumor heterogeneity. NACT significantly impacts the somatic mutation status irrespective of the time to recurrence. The mutated genes detected in chemo-naive pre-NACT tumor samples from either resistant or sensitive cases could potentially have a role in the prediction of chemotherapy response in patients scheduled to receive NACT; larger studies are required to further validate these genes.

The regulatory landscape of retinoblastoma: a pathway analysis perspective

Retinoblastoma (Rb) is a rare intraocular tumour in early childhood, with an approximate incidence of 1 in 18 000 live births. Experimental studies for Rb are complex due to the challenges associated with obtaining a normal retina to contrast with diseased tissue. In this work, we reanalyse a dataset that contains normal retina samples. We identified the individual genes whose expression is different in Rb in contrast with normal tissue, determined the pathways whose global expression pattern is more distant from the global expression observed in normal tissue, and finally, we identified which transcription factors regulate the highest number of differentially expressed genes (DEGs) and proposed as transcriptional master regulators (TMRs). The enrichment of DEGs in the phototransduction and retrograde endocannabinoid signalling pathways could be associated with abnormal behaviour of the processes leading to cellular differentiation and cellular proliferation. On the other hand, the TMRs nuclear receptor subfamily 5 group A member 2 and hepatocyte nuclear factor 4 gamma are involved in hepatocyte differentiation. Therefore, the enrichment of aberrant expression in these transcription factors could suggest an abnormal retina development that could be involved in Rb origin and progression.

Machine Learning-Based Comparative Analysis of Pan-Cancer and Pan-Normal Tissues Identifies Pan-Cancer Tissue-Enriched circRNAs Related to Cancer Mutations as Potential Exosomal Biomarkers

A growing body of evidence has shown that circular RNA (circRNA) is a promising exosomal cancer biomarker candidate. However, global circRNA alterations in cancer and the underlying mechanism, essential for identification of ideal circRNA cancer biomarkers, remain under investigation. We comparatively analyzed the circRNA landscape in pan-cancer and pan-normal tissues. Using co-expression and LASSO regularization analyses, as well as a support vector machine, we analyzed 265 pan-cancer and 319 pan-normal tissues in order to identify the circRNAs with the highest ability to distinguish between pan-cancer and pan-normal tissues. We further studied their expression in plasma exosomes from patients with cancer and their relation with cancer mutations and tumor microenvironment landscape. We discovered that circRNA expression was globally reduced in pan-cancer tissues and plasma exosomes from cancer patients than in pan-normal tissues and plasma exosomes from healthy controls. We identified dynein axonemal heavy chain 14 (DNAH14), the top back-spliced gene exclusive to pan-cancer tissues, as the host gene of three pan-cancer tissue-enriched circRNAs. Among these three circRNAs, chr1_224952669_224968874_+ was significantly elevated in plasma exosomes from hepatocellular carcinoma and colorectal cancer patients. It was also related to the cancer mutation chr1:224952669: G>A, a splice acceptor variant, and was increasingly transcription-driven in cancer tissues. Moreover, pan-cancer tissue-enriched and pan-normal tissue-enriched circRNAs were associated with distinct tumor microenvironment patterns. Our machine learning-based analysis provides insights into the aberrant landscape and biogenesis of circRNAs in cancer and highlights cancer mutation-related and DNAH14-derived circRNA, chr1_224952669_224968874_+, as a potential cancer biomarker.

Changes in Protein Structural Motifs upon Post-Translational Modification in Kidney Cancer

Post-translational modification (PTM) leads to conformational changes in protein structure, modulates the biological function of proteins, and, consequently, changes the signature of metabolic transformations and the immune response in the body. Common PTMs are reversible and serve as a mechanism for modulating metabolic trans-formations in cells. It is likely that dysregulation of post-translational cellular signaling leads to abnormal proliferation and oncogenesis. We examined protein PTMs in the blood samples from patients with kidney cancer. Conformational changes in proteins after modification were analyzed. The proteins were analyzed using ultra-high resolution HPLC-MS/MS and structural analysis was performed with the AMBER and GROMACS software packages. Fifteen proteins containing PTMs were identified in blood samples from patients with kidney cancer. For proteins with PDB structures, a comparative analysis of the structural changes accompanying the modifications was performed. Results revealed that PTMs are localized in stable and compact space protein globule motifs that are exposed to a solvent. The phenomenon of modification is accompanied, as a rule, by an increase in the area available for the solvent of the modified amino acid residue and its active environment.

Whole genome sequencing identifies rare germline variants enriched in cancer related genes in first degree relatives of familial pancreatic cancer patients

First-degree relatives (FDRs) of familial pancreatic cancer (FPC) patients have increased risk of developing pancreatic ductal adenocarcinoma (PDAC). Investigating and understanding the genetic basis for PDAC susceptibility in FPC predisposed families may contribute toward future risk-assessment and management of high-risk individuals. Using a Danish cohort of 27 FPC families, we performed whole-genome sequencing of 61 FDRs of FPC patients focusing on rare genetic variants that may contribute to familial aggregation of PDAC. Statistical analysis was performed using the gnomAD database as external controls. Through analysis of heterozygous premature truncating variants (PTV), we identified cancer-related genes and cancer-driver genes harboring multiple germline mutations. Association analysis detected 20 significant genes with false discovery rate, q < 0.05 including: PALD1, LRP1B, COL4A2, CYLC2, ZFYVE9, BRD3, AHDC1, etc. Functional annotation showed that the significant genes were enriched by gene clusters encoding for extracellular matrix and associated proteins. PTV genes were over-represented by functions related to transport of small molecules, innate immune system, ion channel transport, and stimuli-sensing channels. In conclusion, FDRs of FPC patients carry rare germline variants related to cancer pathogenesis that may contribute to increased susceptibility to PDAC. The identified variants may potentially be useful for risk prediction of high-risk individuals in predisposed families.

Ammonia exposure induces oxidative stress and inflammation by destroying the microtubule structures and the balance of solute carriers in the trachea of pigs

Ammonia (NH₃) is the most alkaline gaseous compound in the atmosphere and the primary gas pollutant in the piggery. It can cause irritation and damage to the airway after inhalation. However, the effects and toxicity mechanism of NH₃ on the trachea are still unclear. In order to evaluate the toxic effects of NH₃ inhalation on pig trachea, the changes of oxidative stress parameters (SOD, GSH, GSH-Px, and MDA), tissue structure and transcriptome in the trachea of pigs were examined after 30 days of exposure to NH₃. Our results showed SOD, GSH-Px and GSH in the trachea in the NH₃-treatment group were significantly decreased (P < 0.05) compared with the control group, on the contrary, MDA content was significantly higher (P < 0.05). The analysis of differentially expressed genes (DEGs) showed that 2542 DEGs (1109 up-regulated DEGs and 1433 down-regulated DEGs) were significantly changed under NH₃ exposure, including many DEGs associated with inflammation, oxidative stress, microtubule activity and SLC family, and the qRT-PCR verification results of these DEGs were consistent with the transcriptome results. The results indicated that NH₃ exposure could break down the mucosal barrier of the respiratory tract, induce oxidative stress and inflammation, reduce the activity of microtubules and disrupt the balance of SLC transporters. In this study, transcriptome analysis was used for the first time to explore the toxic mechanism of NH₃ on pig trachea, providing new insights for better assessing the toxicity mechanism of NH₃, as well as references for comparative medicine.

Immunological classification of gliomas based on immunogenomic profiling

Background

Gliomas are heterogeneous in the tumor immune microenvironment (TIM). However, a classification of gliomas based on immunogenomic profiling remains lacking.

Methods

We hierarchically clustered gliomas based on the enrichment levels of 28 immune cells in the TIM in five datasets and obtained three clusters: immunity-high, immunity-medium, and immunity-low.

Results

Glioblastomas were mainly distributed in immunity-high and immunity-medium, while lower-grade gliomas were distributed in all the three subtypes and predominated in immunity-low. Immunity-low displayed a better survival than other subtypes, indicating a negative correlation between immune infiltration and survival prognosis in gliomas. IDH mutations had a negative correlation with glioma immunity. Immunity-high had higher tumor stemness and epithelial-mesenchymal transition scores and included more high-grade tumors than immunity-low, suggesting that elevated immunity is associated with tumor progression in gliomas. Immunity-high had higher tumor mutation burden and more frequent somatic copy number alterations, suggesting a positive association between tumor immunity and genomic instability in gliomas.

Conclusions

The identification of immune-specific glioma subtypes has potential clinical implications for the immunotherapy of gliomas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12974-020-02030-w.

Discovery of increased epidermal DNAH10 expression after regeneration of dermis in a randomized with-in person trial - reflections on psoriatic inflammation

Because molecular memories of past inflammatory events can persist in epidermal cells, we evaluated the long-term epidermal protein expression landscapes after dermal regeneration and in psoriatic inflammation. We first characterized the effects of two dermal regeneration strategies on transplants of indicator split-thickness skin grafts (STSGs) in ten adult patients with deep burns covering more than 20% of their body surface area. After fascial excision, three adjacent areas within the wound were randomized to receive a permanent dermal matrix, a temporary granulation-tissue-inducing dressing or no dermal component as control. Control areas were covered with STSG immediately, and treated areas after two-weeks of dermis formation. Epidermis-dermis-targeted proteomics of one-year-follow-up samples were performed for protein expression profiling. Epidermal expression of axonemal dynein heavy chain 10 (DNAH10) was increased 20-fold in samples having had regenerating dermis vs control. Given the dermal inflammatory component found in our dermal regeneration samples as well as in early psoriatic lesions, we hypothesized that DNAH10 protein expression also would be affected in psoriatic skin samples. We discovered increased DNAH10 expression in inflammatory lesions when compared to unaffected skin. Our results associate DNAH10 expression with cell proliferation and inflammation as well as with the epidermal memory resulting from the previous regenerative signals of dermis. This study (ISRCTN14499986) was funded by the Finnish Ministry of Defense and by government subsidies for medical research.