Genetic changes in the FH gene cause vagal paraganglioma.

Vagal paraganglioma (VPGL) is a rare neuroendocrine tumor that originates from the paraganglion associated with the vagus nerve. VPGLs present challenges in terms of diagnostics and treatment. VPGL can occur as a hereditary tumor and, like other head and neck paragangliomas, is most frequently associated with mutations in the SDHx genes. However, data regarding the genetics of VPGL are limited. Herein, we report a rare case of a 41-year-old woman with VPGL carrying a germline variant in the FH gene. Using whole-exome sequencing, a variant, FH p.S249R, was identified; no variants were found in other PPGL susceptibility and candidate genes. Loss of heterozygosity analysis revealed the loss of the wild-type allele of the FH gene in the tumor. The pathogenic effect of the p.S249R variant on FH activity was confirmed by immunohistochemistry for S-(2-succino)cysteine (2SC). Potentially deleterious somatic variants were found in three genes, SLC7A7, ZNF225, and MED23. The latter two encode transcriptional regulators that can impact gene expression deregulation and are involved in tumor development and progression. Moreover, FH-mutated VPGL was characterized by a molecular phenotype different from SDHx-mutated PPGLs. In conclusion, the association of genetic changes in the FH gene with the development of VPGL was demonstrated. The germline variant FH: p.S249R and somatic deletion of the second allele can lead to biallelic gene damage that promotes tumor initiation. These results expand the clinical and mutation spectra of FH-related disorders and improve our understanding of the molecular genetic mechanisms underlying the pathogenesis of VPGL.

Matrisome Transcriptome Dynamics during Tissue Aging.

The extracellular matrix (ECM) is a complex three-dimensional network of macromolecules that provides structural support for the cells and plays a significant role in tissue homeostasis and repair. Growing evidence indicates that dysregulation of ECM remodeling contributes to various pathological conditions in the body, including age-associated diseases. In this work, gene expression data of normal human tissues obtained from the Genotype-Tissue Expression project, as well as data from MatrisomeDB 2.0, the ECM-protein knowledge database, are used to estimate the age-dependent matrisome transcriptome dynamics in the blood, heart, brain, liver, kidneys, lungs, and muscle. Differential gene expression (DE) analysis revealed dozens of matrisome genes encoding both structural elements of the ECM and ECM-associated proteins, which had a tissue-specific expression profile with age. Among common DE genes that changed expression with age in at least three tissues, COL18A1, MFAP1, IGFBP7, AEBP1, LTBP2, LTBP4, LG14, EFEMP1, PRELP, BGN, FAM20B, CTSC, CTSS, and CLEC2B were observed. The findings of the study also reveal that there are sex-specific alterations during aging in the matrisome gene expression. Taken together, the results obtained in this work may help in understanding the role of the ECM in tissue aging and might prove valuable for the future development of the field of ECM research in general.

Decelerated Epigenetic Aging in Long Livers.

Epigenetic aging is a hot topic in the field of aging research. The present study estimated epigenetic age in long-lived individuals, who are currently actively being studied worldwide as an example of successful aging due to their longevity. We used Bekaert's blood-based age prediction model to estimate the epigenetic age of 50 conditionally "healthy" and 45 frail long-livers over 90 years old. Frailty assessment in long-livers was conducted using the Frailty Index. The control group was composed of 32 healthy individuals aged 20-60 years. The DNA methylation status of 4 CpG sites (ASPA CpG1, PDE4C CpG1, ELOVL2 CpG6, and EDARADD CpG1) included in the epigenetic clock was assessed through pyrosequencing. According to the model calculations, the epigenetic age of long-livers was significantly lower than their chronological age (on average by 21 years) compared with data from the group of people aged 20 to 60 years. This suggests a slowing of epigenetic and potentially biological aging in long livers. At the same time, the obtained results showed no statistically significant differences in delta age (difference between the predicted and chronological age) between "healthy" long livers and long livers with frailty. We also failed to detect sex differences in epigenetic age either in the group of long livers or in the control group. It is possible that the predictive power of epigenetic clocks based on a small number of CpG sites is insufficient to detect such differences. Nevertheless, this study underscores the need for further research on the epigenetic status of centenarians to gain a deeper understanding of the factors contributing to delayed aging in this population.

Transcriptome Profiling of Prostate Cancer, Considering Risk Groups and the TMPRSS2-ERG Molecular Subtype.

Molecular heterogeneity in prostate cancer (PCa) is one of the key reasons underlying the differing likelihoods of recurrence after surgical treatment in individual patients of the same clinical category. In this study, we performed RNA-Seq profiling of 58 localized PCa and 43 locally advanced PCa tissue samples obtained as a result of radical prostatectomy on a cohort of Russian patients. Based on bioinformatics analysis, we examined features of the transcriptome profiles within the high-risk group, including within the most commonly represented molecular subtype, TMPRSS2-ERG. The most significantly affected biological processes in the samples were also identified, so that they may be further studied in the search for new potential therapeutic targets for the categories of PCa under consideration. The highest predictive potential was found with the EEF1A1P5, RPLP0P6, ZNF483, CIBAR1, HECTD2, OGN, and CLIC4 genes. We also reviewed the main transcriptome changes in the groups at intermediate risk of PCa-Gleason Score 7 (groups 2 and 3 according to the ISUP classification)-on the basis of which the LPL, MYC, and TWIST1 genes were identified as promising additional prognostic markers, the statistical significance of which was confirmed using qPCR validation.

Lymphatic Dissemination in Prostate Cancer: Features of the Transcriptomic Profile and Prognostic Models.

Radical prostatectomy is the gold standard treatment for prostate cancer (PCa); however, it does not always completely cure PCa, and patients often experience a recurrence of the disease. In addition, the clinical and pathological parameters used to assess the prognosis and choose further tactics for treating a patient are insufficiently informative and need to be supplemented with new markers. In this study, we performed RNA-Seq of PCa tissue samples, aimed at identifying potential prognostic markers at the level of gene expression and miRNAs associated with one of the key signs of cancer aggressiveness-lymphatic dissemination. The relative expression of candidate markers was validated by quantitative PCR, including an independent sample of patients based on archival material. Statistically significant results, derived from an independent set of samples, were confirmed for miR-148a-3p and miR-615-3p, as well as for the CST2, OCLN, and PCAT4 genes. Considering the obtained validation data, we also analyzed the predictive value of models based on various combinations of identified markers using algorithms based on machine learning. The highest predictive potential was shown for the "CST2 + OCLN + pT" model (AUC = 0.863) based on the CatBoost Classifier algorithm.

ALDH3A2, ODF2, QSOX2, and MicroRNA-503-5p Expression to Forecast Recurrence in TMPRSS2-ERG-Positive Prostate Cancer.

Following radical surgery, patients may suffer a relapse. It is important to identify such patients so that therapy tactics can be modified appropriately. Existing stratification schemes do not display the probability of recurrence with enough precision since locally advanced prostate cancer (PCa) is classified as high-risk but is not ranked in greater detail. Between 40 and 50% of PCa cases belong to the TMPRSS2-ERG subtype that is a sufficiently homogeneous group for high-precision prognostic marker search to be possible. This study includes two independent cohorts and is based on high throughput sequencing and qPCR data. As a result, we have been able to suggest a perspective-trained model involving a deep neural network based on both qPCR data for mRNA and miRNA and clinicopathological criteria that can be used for recurrence risk forecasts in patients with TMPRSS2-ERG-positive, locally advanced PCa (the model uses ALDH3A2 + ODF2 + QSOX2 + hsa-miR-503-5p + ISUP + pT, with an AUC = 0.944). In addition to the prognostic model's use of identified differentially expressed genes and miRNAs, miRNA-target pairs were found that correlate with the prognosis and can be presented as an interactome network.

Dynamic Profiling of Exosomal microRNAs in Blood Plasma of Patients with Castration-Resistant Prostate Cancer.

Prostate cancer is one of the most common and socially significant cancers among men. The aim of this study was to identify significant changes in the expression of exosomal miRNAs associated with an increase in the level of prostate specific antigen in castration-resistant prostate cancer during therapy and to evaluate them as potential prognostic markers for this category of disease. High-throughput miRNA sequencing was performed on 49 blood plasma samples taken from 11 Russian patients with castration-resistant cancer during therapy. Bioinformatic analysis of the obtained miRNA-seq data was carried out. Additionally, miRNA-seq data from the PRJNA562276 project were analyzed to identify exosomal miRNAs associated with castration-resistant prostate cancer. We found 34 differentially expressed miRNAs associated with the progression of castration-resistant prostate cancer during therapy in Russian patients. It was also shown that hsa-miRNA-148a-3p expression can serve as a potential prognostic marker. We found the exosomal miRNA expression signature associated with castration-resistant prostate cancer progression, in particular on the Russian patient cohort. Many of these miRNAs are well-known players in either oncogenic transformation or tumor suppression. Further experimental studies with extended sampling are required to validate these results.

Impact TMPRSS2-ERG Molecular Subtype on Prostate Cancer Recurrence.

Currently, seven molecular subtypes of prostate cancer (PCa) are known, the most common of which being the subtype characterized by the presence of the TMPRSS2-ERG fusion transcript. While there is a considerable amount of work devoted to the influence of this transcript on the prognosis of the disease, data on its role in the progression and prognosis of PCa remain controversial. The present study is devoted to the analysis of the association between the TMPRSS2-ERG transcript and the biochemical recurrence of PCa. The study included two cohorts: the RNA-Seq sample of Russian patients with PCa (n = 72) and the TCGA-PRAD data (n = 203). The results of the analysis of the association between the TMPRSS2-ERG transcript and biochemical recurrence were contradictory. The differential expression analysis (biochemical recurrence cases versus biochemical recurrence-free) and the gene set enrichment analysis revealed a list of genes involved in major cellular pathways. The GNL3, QSOX2, SSPO, and SYS1 genes were selected as predictors of the potential prognostic model (AUC = 1.000 for a cohort of Russian patients with PCa and AUC = 0.779 for a TCGA-PRAD cohort).

NETO2 Is Deregulated in Breast, Prostate, and Colorectal Cancer and Participates in Cellular Signaling.

The NETO2 gene (neuropilin and tolloid-like 2) encodes a protein that acts as an accessory subunit of kainate receptors and is predominantly expressed in the brain. Upregulation of NETO2 has been observed in several tumors; however, its role in tumorigenesis remains unclear. In this study, we investigated NETO2 expression in breast, prostate, and colorectal cancer using quantitative PCR (qPCR), as well as the effect of shRNA-mediated NETO2 silencing on transcriptome changes in colorectal cancer cells. In the investigated tumors, we observed both increased and decreased NETO2 mRNA levels, presenting no correlation with the main clinicopathological characteristics. In HCT116 cells, NETO2 knockdown resulted in the differential expression of 17 genes and 2 long non-coding RNAs (lncRNAs), associated with the upregulation of circadian rhythm and downregulation of several cancer-associated pathways, including Wnt, transforming growth factor (TGF)-ß, Janus kinase (JAK)-signal transducer and activator of transcription (STAT), mitogen-activated protein kinase (MAPK), and phosphatidylinositol 3-kinase (PI3K)/protein kinase B (AKT) pathways. Furthermore, we demonstrated the possibility to utilize a novel model organism, short-lived fish Nothobranchius furzeri, for evaluating NETO2 functions. The ortholog neto2b in N. furzeri demonstrated a high similarity in nucleotide and amino acid sequences with human NETO2, as well as was characterized by stable expression in various fish tissues. Collectively, our findings demonstrate the deregulation of NETO2 in the breast, prostate, and colorectal cancer and its participation in the tumor development primarily through cellular signaling.

Mutation profiling in eight cases of vagal paragangliomas.

BACKGROUND: Vagal paragangliomas (VPGLs) belong to a group of rare head and neck neuroendocrine tumors. VPGLs arise from the vagus nerve and are less common than carotid paragangliomas. Both diagnostics and therapy of the tumors raise significant challenges. Besides, the genetic and molecular mechanisms behind VPGL pathogenesis are poorly understood. METHODS: The collection of VPGLs obtained from 8 patients of Russian population was used in the study. Exome library preparation and high-throughput sequencing of VPGLs were performed using an Illumina technology. RESULTS: Based on exome analysis, we identified pathogenic/likely pathogenic variants of the SDHx genes, frequently mutated in paragangliomas/pheochromocytomas. SDHB variants were found in three patients, whereas SDHD was mutated in two cases. Moreover, likely pathogenic missense variants were also detected in SDHAF3 and SDHAF4 genes encoding for assembly factors for the succinate dehydrogenase (SDH) complex. In a patient, we found a novel variant of the IDH2 gene that was predicted as pathogenic by a series of algorithms used (such as SIFT, PolyPhen2, FATHMM, MutationTaster, and LRT). Additionally, pathogenic/likely pathogenic variants were determined for several genes, including novel genes and some genes previously reported as associated with different types of tumors. CONCLUSIONS: Results indicate a high heterogeneity among VPGLs, however, it seems that driver events in most cases are associated with mutations in the SDHx genes and SDH assembly factor-coding genes that lead to disruptions in the SDH complex.

miRNAs expression signature potentially associated with lymphatic dissemination in locally advanced prostate cancer.

BACKGROUND: Prostate cancer is one of the most common and socially significant cancers among men. The aim of our study was to reveal changes in miRNA expression profiles associated with lymphatic dissemination in prostate cancer and to identify the most prominent miRNAs as potential prognostic markers for future studies. METHODS: High-throughput miRNA sequencing was performed for 44 prostate cancer specimens taken from Russian patients, with and without lymphatic dissemination (N1 - 20 samples; N0 - 24 samples). RESULTS: We found at least 18 microRNAs with differential expression between N0 and N1 sample groups: miR-182-5p, miR-183-5p, miR-96-5p, miR-25-3p, miR-93-5p, miR-7-5p, miR-615-3p, miR-10b, miR-1248 (N1-miRs; elevated expression in N1 cohort; p < 0.05); miR-1271-5p, miR-184, miR-222-3p, miR-221-5p, miR-221-3p, miR-455-3p, miR-143-5p, miR-181c-3p and miR-455-5p (N0-miRs; elevated expression in N0; p < 0.05). The expression levels of N1-miRs were highly correlated between each other (the same is applied for N0-miRs) and the expression levels of N0-miRs and N1-miRs were anti-correlated. The tumor samples can be divided into two groups depending on the expression ratio between N0-miRs and N1-miRs. CONCLUSIONS: We found the miRNA expression signature associated with lymphatic dissemination, in particular on the Russian patient cohort. Many of these miRNAs are well-known players in either oncogenic transformation or tumor suppression. Further experimental studies with extended sampling are required to validate these results.

Gene Expression Changes and Associated Pathways Involved in the Progression of Prostate Cancer Advanced Stages.

Prostate cancer (PC) is one of the most common cancers among men worldwide, and advanced PCs, such as locally advanced PC (LAPC) and castration-resistant PC (CRPC), present the greatest challenges in clinical management. Current indicators have limited capacity to predict the disease course; therefore, better prognostic markers are greatly needed. In this study, we performed a bioinformatic analysis of The Cancer Genome Atlas (TCGA) datasets, including RNA-Seq data from the prostate adenocarcinoma (PRAD; n = 55) and West Coast Dream Team - metastatic CRPC (WCDT-MCRPC; n = 84) projects, to evaluate the transcriptome changes associated with progression-free survival (PFS) for LAPC and CRPC, respectively. We identified the genes whose expression was positively/negatively correlated with PFS. In LAPC, the genes with the most significant negative correlations were ZC2HC1A, SQLE, and KIF11, and the genes with the most significant positive correlations were SOD3, LRRC26, MIR22HG, MEG3, and MIR29B2CHG. In CRPC, the most significant positive correlations were found for BET1, CTAGE5, IFNGR1, and GIMAP6, and the most significant negative correlations were found for CLPB, PRPF19, ZNF610, MPST, and LINC02001. In addition, we performed a gene network interaction analysis using STRINGdb, which revealed a significant relationship between genes predominantly involved in the cell cycle and characterized by upregulated expression in early recurrence. Based on the results, we propose several genes that can be used as potential prognostic markers.

Differentially Expressed Genes Associated With Prognosis in Locally Advanced Lymph Node-Negative Prostate Cancer.

Older age is one of the main risk factors for cancer development. The incidence of prostate cancer, as a multifactorial disease, also depends upon demographic factors, race, and genetic predisposition. Prostate cancer most frequently occurs in men over 60 years of age, indicating a clear association between older age and disease onset. Carcinogenesis is followed by the deregulation of many genes, and some of these changes could serve as biomarkers for diagnosis, prognosis, prediction of drug therapy efficacy, as well as possible therapeutic targets. We have performed a bioinformatic analysis of a The Cancer Genome Atlas (TCGA) data and RNA-Seq profiling of a Russian patient cohort to reveal prognostic markers of locally advanced lymph node-negative prostate cancer (lymph node-negative LAPC). We also aimed to identify markers of the most common molecular subtype of prostate cancer carrying a fusion transcript TMPRSS2-ERG. We have found several genes that were differently expressed between the favorable and unfavorable prognosis groups and involved in the enriched KEGG pathways based on the TCGA (B4GALNT4, PTK6, and CHAT) and Russian patient cohort data (AKR1C1 and AKR1C3). Additionally, we revealed such genes for the TMPRSS2-ERG prostate cancer molecular subtype (B4GALNT4, ASRGL1, MYBPC1, RGS11, SLC6A14, GALNT13, and ST6GALNAC1). Obtained results contribute to a better understanding of the molecular mechanisms behind prostate cancer progression and could be used for further development of the LAPC prognosis marker panel.

The CIMP-high phenotype is associated with energy metabolism alterations in colon adenocarcinoma.

BACKGROUND: CpG island methylator phenotype (CIMP) is found in 15-20% of malignant colorectal tumors and is characterized by strong CpG hypermethylation over the genome. The molecular mechanisms of this phenomenon are not still fully understood. The development of CIMP is followed by global gene expression alterations and metabolic changes. In particular, CIMP-low colon adenocarcinoma (COAD), predominantly corresponded to consensus molecular subtype 3 (CMS3, "Metabolic") subgroup according to COAD molecular classification, is associated with elevated expression of genes participating in metabolic pathways. METHODS: We performed bioinformatics analysis of RNA-Seq data from The Cancer Genome Atlas (TCGA) project for CIMP-high and non-CIMP COAD samples with DESeq2, clusterProfiler, and topGO R packages. Obtained results were validated on a set of fourteen COAD samples with matched morphologically normal tissues using quantitative PCR (qPCR). RESULTS: Upregulation of multiple genes involved in glycolysis and related processes (ENO2, PFKP, HK3, PKM, ENO1, HK2, PGAM1, GAPDH, ALDOA, GPI, TPI1, and HK1) was revealed in CIMP-high tumors compared to non-CIMP ones. Most remarkably, the expression of the PKLR gene, encoding for pyruvate kinase participating in gluconeogenesis, was decreased approximately 20-fold. Up to 8-fold decrease in the expression of OGDHL gene involved in tricarboxylic acid (TCA) cycle was observed in CIMP-high tumors. Using qPCR, we confirmed the increase (4-fold) in the ENO2 expression and decrease (2-fold) in the OGDHL mRNA level on a set of COAD samples. CONCLUSIONS: We demonstrated the association between CIMP-high status and the energy metabolism changes at the transcriptomic level in colorectal adenocarcinoma against the background of immune pathway activation. Differential methylation of at least nine CpG sites in OGDHL promoter region as well as decreased OGDHL mRNA level can potentially serve as an additional biomarker of the CIMP-high status in COAD.

Mutational load in carotid body tumor.

BACKGROUND: Carotid body tumor (CBT) is a rare neoplasm arising from paraganglion located near the bifurcation of the carotid artery. There is great intra-tumor heterogeneity, and CBT development could be associated with both germline and somatic allelic variants. Studies on the molecular genetics of CBT are limited, and the molecular mechanisms of its pathogenesis are not fully understood. This work is focused on the estimation of mutational load (ML) in CBT. METHODS: Using the NextSeq 500 platform, we performed exome sequencing of tumors with matched lymph node tissues and peripheral blood obtained from six patients with CBT. To obtain reliable results in tumors with low ML, we developed and successfully applied a complex approach for the analysis of sequencing data. ML was evaluated as the number of somatic variants per megabase (Mb) of the target regions covered by the Illumina TruSeq Exome Library Prep Kit. RESULTS: The ML in CBT varied in the range of 0.09-0.28/Mb. Additionally, we identified several pathogenic/likely pathogenic somatic and germline allelic variants across six patients studied (including TP53 variants). CONCLUSIONS: Using the developed approach, we estimated the ML in CBT, which is much lower than in common malignant tumors. Identified variants in known paraganglioma/pheochromocytoma-causative genes and novel genes could be associated with the pathogenesis of CBT. The obtained results expand our knowledge of the mutation process in CBT as well as the biology of tumor development.