PANACEA: network-based methods for pharmacotherapy prioritization in personalized oncology.

MOTIVATION: Identifying appropriate pharmacotherapy options from genomics results is a significant challenge in personalized oncology. However, computational methods for prioritizing drugs are underdeveloped. With the hypothesis that network-based approaches can improve the performance by extending the use of potential drug targets beyond direct interactions, we devised two network-based methods for personalized pharmacotherapy prioritization in cancer. RESULTS: We developed novel personalized drug prioritization approaches, PANACEA: PersonAlized Network-based Anti-Cancer therapy EvaluAtion. In PANACEA, initially, the protein interaction network is extended with drugs, and a driverness score is assigned to each altered gene. For scoring drugs, either (i) the 'distance-based' method, incorporating the shortest distance between drugs and altered genes, and driverness scores, or (ii) the 'propagation' method involving the propagation of driverness scores via a random walk with restart framework is performed. We evaluated PANACEA using multiple datasets, and demonstrated that (i) the top-ranking drugs are relevant for cancer pharmacotherapy using TCGA data; (ii) drugs that cancer cell lines are sensitive to are identified using GDSC data; and (iii) PANACEA can perform adequately in the clinical setting using cases with known drug responses. We also illustrate that the proposed methods outperform iCAGES and PanDrugs, two previous personalized drug prioritization approaches. AVAILABILITY AND IMPLEMENTATION: The corresponding R package is available on GitHub. (https://github.com/egeulgen/PANACEA.git). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

driveR: a novel method for prioritizing cancer driver genes using somatic genomics data.

BACKGROUND: Cancer develops due to "driver" alterations. Numerous approaches exist for predicting cancer drivers from cohort-scale genomics data. However, methods for personalized analysis of driver genes are underdeveloped. In this study, we developed a novel personalized/batch analysis approach for driver gene prioritization utilizing somatic genomics data, called driveR. RESULTS: Combining genomics information and prior biological knowledge, driveR accurately prioritizes cancer driver genes via a multi-task learning model. Testing on 28 different datasets, this study demonstrates that driveR performs adequately, achieving a median AUC of 0.684 (range 0.651-0.861) on the 28 batch analysis test datasets, and a median AUC of 0.773 (range 0-1) on the 5157 personalized analysis test samples. Moreover, it outperforms existing approaches, achieving a significantly higher median AUC than all of MutSigCV (Wilcoxon rank-sum test p < 0.001), DriverNet (p < 0.001), OncodriveFML (p < 0.001) and MutPanning (p < 0.001) on batch analysis test datasets, and a significantly higher median AUC than DawnRank (p < 0.001) and PRODIGY (p < 0.001) on personalized analysis datasets. CONCLUSIONS: This study demonstrates that the proposed method is an accurate and easy-to-utilize approach for prioritizing driver genes in cancer genomes in personalized or batch analyses. driveR is available on CRAN: https://cran.r-project.org/package=driveR .

Endometrial gene expression profiling of recurrent implantation failure after in vitro fertilization.

Recurrent implantation failure (RIF) is diagnosed when good-quality embryos repeatedly fail to implant after transfer in several in vitro fertilization (IVF) treatment cycles. Different expression profiles in maternal mRNAs could be referring to many diseases including RIF. This study aimed to reveal significantly dysregulated selected genes expression between healthy fertile women and RIF patients in the implantation window days of the natural menstrual cycle. MME, WWC1, TNC, and FOXP3 genes were chosen as target genes regarding their possible relations with the implantation process. Pathways with these genes were identified and the relationship between these pathways and RIF was investigated. In this study, the endometrial biopsy samples were collected in the secretory phase (cycle day 20-24) of the menstrual cycle from RIF patients (n = 34) and healthy fertile controls (n = 34). After "Pathway and network-oriented GWAS analysis" (PANOGA) and "Kyoto Encyclopedia of Genes and Genomes" (KEGG) pathway analysis; "Membrane Metalloendopeptidase" (MME), "WW and C2 Domain Containing 1" (WWC1), "Tenascin C" (TNC) and "Forkhead Box P3" (FOXP3) genes were chosen as target genes by regarding their possible relation with implantation process. Detection of differences in mRNA expressions between the control group and RIF patients has been performed with the droplet digital PCR (ddPCR) method. Results of the study showed that MME and WWC1 genes expression levels are significantly (p < 0,05) up-regulated 4.9 and 5.2 times respectively and TNC gene expression level is significantly (p < 0,05) down-regulated 9 times in the RIF samples compared to the control group. However, no statistically significant difference was observed between the patient group and the control group in the expression of the FOXP3 gene (p < 0.05). Changes are observed in the expression of the renin-angiotensin system pathway in which the MME gene is involved in the implantation process. The increase in MME gene expression can be speculated to cause implantation failure by restricting the invasion of trophoblast cells. Increasing WWC1 gene expression in the Hippo signaling pathway inhibits "Yes-associated protein 1" (YAP) expression, which is a transcriptional cofactor. Inhibition of YAP protein expression may impair the implantation process by causing the failure of endometrial decidualization. The TNC gene is located in the focal adhesion pathway and this pathway reduces cell adhesion on the endometrial surface to facilitate the attachment of the embryo to the endometrium. The reason for implantation failure might be that the intercellular connections are not suitable for implantation as a result of decreased expression of the focal adhesion pathway in which the TNC gene is effective. Considering the relations between the pathways of the target genes and the implantation process, changes in the expression of target genes might be a cause of RIF.

Clinical, Electrodiagnostic, and Genetic Features of Tangier Disease in an Adolescent Girl with Presentation of Peripheral Neuropathy.

Tangier disease (TD) is a rare, autosomal recessive inherited disorder caused by a mutation in the adenosine triphosphate-binding cassette transporter 1 (ABCA1) gene, which results in a decrease in plasma high-density lipoprotein (HDL) levels. Peripheral neuropathy can be seen in approximately 50% of patients with TD, which usually occurs after the age of 15 years, and is characterized by relapsing-remitting mono- or polyneuropathy or syringomyelia-like neuropathy. Herein, we report a 16-year-old female patient who was initially diagnosed with peripheral neuropathy at the age of 13 years. Whole exome sequencing was performed, and a nonsense mutation (p.Arg1817X) in ABCA1 was identified. The patient was investigated for systemic findings of TD after the genetic diagnosis was made, and low (< 5 mg/dL) levels of HDL cholesterol were detected by lipid electrophoresis. Other family members were reexamined after the diagnosis of the proband, and asymptomatic sister of the proband was diagnosed with TD. We would like to emphasize that TD should be considered in the differential diagnosis of pediatric patients presenting with peripheral neuropathy; furthermore detection of HDL levels by lipid electrophoresis is a simple but indicative diagnostic test.

A microdeletion event at 19q13.43 in IDH-mutant astrocytomas is strongly correlated with MYC overexpression.

MYC dysregulation is pivotal in the onset and progression of IDH-mutant gliomas, mostly driven by copy-number alterations, regulatory element alterations, or epigenetic changes. Our pilot analysis uncovered instances of relative MYC overexpression without alterations in the proximal MYC network (PMN), prompting a deeper investigation into potential novel oncogenic mechanisms. Analysing comprehensive genomics profiles of 236 "IDH-mutant 1p/19q non-co-deleted" lower-grade gliomas from The Cancer Genome Atlas, we identified somatic genomic alterations within the PMN. In tumours without PMN-alterations but with MYC-overexpression, genes correlated with MYC-overexpression were identified. Our analyses yielded that 86/236 of astrocytomas exhibited no PMN-alterations, a subset of 21/86 displaying relative MYC overexpression. Within this subset, we discovered 42 genes inversely correlated with relative MYC expression, all on 19q. Further analysis pinpointed a minimal common region at 19q13.43, encompassing 15 genes. The inverse correlations of these 15 genes with relative MYC overexpression were re-confirmed using independent scRNAseq data. Further, the micro-deleted astrocytoma subset displayed significantly higher genomic instability compared to WT cases, but lower instability compared to PMN-hit cases. This newly identified 19q micro-deletion represents a potential novel mechanism underlying MYC dysregulation in astrocytomas. Given the prominence of 19q loss in IDH-mutant gliomas, our findings bear significant implications for understanding gliomagenesis.

Integrated bioinformatics analysis of validated and circulating miRNAs in ovarian cancer.

Recent studies have focused on the early detection of ovarian cancer (OC) using tumor materials by liquid biopsy. The mechanisms of microRNAs (miRNAs) to impact OC and signaling pathways are still unknown. This study aims to reliably perform functional analysis of previously validated circulating miRNAs' target genes by using pathfindR. Also, overall survival and pathological stage analyses were evaluated with miRNAs' target genes which are common in the The Cancer Genome Atlas and GTEx datasets. Our previous studies have validated three downregulated miRNAs (hsa-miR-885-5p, hsa-miR-1909-5p, and hsalet7d-3p) having a diagnostic value in OC patients' sera, with high-throughput techniques. The predicted target genes of these miRNAs were retrieved from the miRDB database (v6.0). Active-subnetwork-oriented Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was conducted by pathfindR using the target genes. Enrichment of KEGG pathways assessed by the analysis of pathfindR indicated that 24 pathways were related to the target genes. Ubiquitin-mediated proteolysis, spliceosome and Notch signaling pathway were the top three pathways with the lowest p-values (p < 0.001). Ninety-three common genes were found to be differentially expressed (p < 0.05) in the datasets. No significant genes were found to be significant in the analysis of overall survival analyses, but 24 genes were found to be significant with pathological stages analysis (p < 0.05). The findings of our study provide in-silico evidence that validated circulating miRNAs' target genes and enriched pathways are related to OC and have potential roles in theranostics applications. Further experimental investigations are required to validate our results which will ultimately provide a new perspective for translational applications in OC management.

Gene pathway analysis of the endometrium at the start of the window of implantation in women with unexplained infertility and unexplained recurrent pregnancy loss: is unexplained recurrent pregnancy loss a subset of unexplained infertility?

This study aims to understand differences/similarities in the genetic profile of the endometrium at the start of window of implantation (WOI) in women with unexplained infertility (UI) and unexplained recurrent pregnancy loss (uRPL). Differentially expressed genes (DEGs) from the endometrium were evaluated using gene expression array and pathway enrichment analysis was performed to analyse gene expression pathways involved in both conditions. We found 2,171 genes arranged in 117 pathways and 730 genes arranged in 33 pathways differentially expressed in endometrium of patients in UI and uRPL, respectively. Complement-coagulation cascades, morphine addiction pathway, and PI3K-Akt signalling pathway were predominantly differentially expressed in UI. Cancer pathways, NF-κB signalling pathway, and actin cytoskeleton regulation pathway showed significant changes in uRPL. Forty-eight percent of DEGs and 84% of differentially expressed pathways in uRPL were found in the endometrium of UI patients. Unexpected close association in gene expression pathways between UI and uRPL is observed supporting the hypothesis 'uRPL is a clinical subset of UI'. Yet 100% DEGs overlap wasn't found suggesting the endometrium has still some different gene expression patterns at start of WOI in UI and uRPL. Lastly, diagnostic tools may be developed for uRPL because more specific genes-pathways are involved compared with UI, which shows broader genetic expression profile.

Correlation of anatomical involvement patterns of insular gliomas with subnetworks of the limbic system.

OBJECTIVE: Gliomas frequently involve the insula both primarily and secondarily by invasion. Despite the high connectivity of the human insula, gliomas do not spread randomly to or from the insula but follow stereotypical anatomical involvement patterns. In the majority of cases, these patterns correspond to the intrinsic connectivity of the limbic system, except for tumors with aggressive biology. On the basis of these observations, the authors hypothesized that these different involvement patterns may be correlated with distinct outcomes and analyzed these correlations in an institutional cohort. METHODS: Fifty-nine patients who had undergone surgery for insular diffuse gliomas and had complete demographic, pre- and postoperative imaging, pathology, molecular genetics, and clinical follow-up data were included in the analysis (median age 37 years, range 21-71 years, M/F ratio 1.68). Patients with gliomatosis and those with only minor involvement of the insula were excluded. The presence of T2-hyperintense tumor infiltration was evaluated in 12 anatomical structures. Hierarchical biclustering was used to identify co-involved structures, and the findings were correlated with established functional anatomy knowledge. Overall survival was evaluated using Kaplan-Meier and Cox proportional hazards regression analysis (17 parameters). RESULTS: The tumors involved the anterior insula (98.3%), posterior insula (67.8%), temporal operculum (47.5%), amygdala (42.4%), frontal operculum (40.7%), temporal pole (39%), parolfactory area (35.6%), hypothalamus (23.7%), hippocampus (16.9%), thalamus (6.8%), striatum (5.1%), and cingulate gyrus (3.4%). A mean 4.2 ± 2.6 structures were involved. On the basis of hierarchical biclustering, 7 involvement patterns were identified and correlated with cortical functional anatomy (pure insular [11.9%], olfactocentric [15.3%], olfactoopercular [33.9%], operculoinsular [15.3%], striatoinsular [3.4%], translimbic [11.9%], and multifocal [8.5%] patterns). Cox regression identified hippocampal involvement (p = 0.006) and postoperative tumor volume (p = 0.027) as significant negative independent prognosticators of overall survival and extent of resection (p = 0.015) as a significant positive independent prognosticator. CONCLUSIONS: The study findings indicate that insular gliomas primarily involve the olfactocentric limbic girdle and that involvement in the hippocampocentric limbic girdle is associated with a worse prognosis.

CogNet: classification of gene expression data based on ranked active-subnetwork-oriented KEGG pathway enrichment analysis.

Most of the traditional gene selection approaches are borrowed from other fields such as statistics and computer science, However, they do not prioritize biologically relevant genes since the ultimate goal is to determine features that optimize model performance metrics not to build a biologically meaningful model. Therefore, there is an imminent need for new computational tools that integrate the biological knowledge about the data in the process of gene selection and machine learning. Integrative gene selection enables incorporation of biological domain knowledge from external biological resources. In this study, we propose a new computational approach named CogNet that is an integrative gene selection tool that exploits biological knowledge for grouping the genes for the computational modeling tasks of ranking and classification. In CogNet, the pathfindR serves as the biological grouping tool to allow the main algorithm to rank active-subnetwork-oriented KEGG pathway enrichment analysis results to build a biologically relevant model. CogNet provides a list of significant KEGG pathways that can classify the data with a very high accuracy. The list also provides the genes belonging to these pathways that are differentially expressed that are used as features in the classification problem. The list facilitates deep analysis and better interpretability of the role of KEGG pathways in classification of the data thus better establishing the biological relevance of these differentially expressed genes. Even though the main aim of our study is not to improve the accuracy of any existing tool, the performance of the CogNet outperforms a similar approach called maTE while obtaining similar performance compared to other similar tools including SVM-RCE. CogNet was tested on 13 gene expression datasets concerning a variety of diseases.

Sequential filtering for clinically relevant variants as a method for clinical interpretation of whole exome sequencing findings in glioma.

BACKGROUND: In the clinical setting, workflows for analyzing individual genomics data should be both comprehensive and convenient for clinical interpretation. In an effort for comprehensiveness and practicality, we attempted to create a clinical individual whole exome sequencing (WES) analysis workflow, allowing identification of genomic alterations and presentation of neurooncologically-relevant findings. METHODS: The analysis workflow detects germline and somatic variants and presents: (1) germline variants, (2) somatic short variants, (3) tumor mutational burden (TMB), (4) microsatellite instability (MSI), (5) somatic copy number alterations (SCNA), (6) SCNA burden, (7) loss of heterozygosity, (8) genes with double-hit, (9) mutational signatures, and (10) pathway enrichment analyses. Using the workflow, 58 WES analyses from matched blood and tumor samples of 52 patients were analyzed: 47 primary and 11 recurrent diffuse gliomas. RESULTS: The median mean read depths were 199.88 for tumor and 110.955 for normal samples. For germline variants, a median of 22 (14-33) variants per patient was reported. There was a median of 6 (0-590) reported somatic short variants per tumor. A median of 19 (0-94) broad SCNAs and a median of 6 (0-12) gene-level SCNAs were reported per tumor. The gene with the most frequent somatic short variants was TP53 (41.38%). The most frequent chromosome-/arm-level SCNA events were chr7 amplification, chr22q loss, and chr10 loss. TMB in primary gliomas were significantly lower than in recurrent tumors (p = 0.002). MSI incidence was low (6.9%). CONCLUSIONS: We demonstrate that WES can be practically and efficiently utilized for clinical analysis of individual brain tumors. The results display that NOTATES produces clinically relevant results in a concise but exhaustive manner.

Comparison of endometrial prostanoid profiles in three infertile subgroups: the missing part of receptivity?

OBJECTIVE: To study the prostanoid profile of the endometria of patients with recurrent implantation failure (RIF), unexplained infertility (UIF), and recurrent miscarriages (RM), and to compare them with the endometria of healthy fertile controls. DESIGN: Prospective cohort study. SETTING: University hospital. PATIENT(S): Fifteen patients with RIF, 18 patients with UIF, 16 patients with RM, and 23 fertile controls were recruited. INTERVENTION(S): Endometrial samples were taken during the window of implantation. After tissue homogenization and extraction, analysis with ultra-performance liquid chromatography diode array detector electrospray ionisation tandem mass spectrometry was performed. MAIN OUTCOME MEASURES: Concentrations of prostaglandin (PG) D1, PGE1, PGF1α, 6-ketoPGF1α, PGD2, PGE2, PGF2α, 15-deoxy-Δ12,14-PGJ2, PGD3, PGE3, PGF3α, thromboxane B2, 13,14-dihydro-PGE1, 13,14-dihydro-PGF1α, 13,14-dihydro-PGF2α, 13,14-dihydro-15-keto-PGE1, 13,14-dihydro-15-keto-PGE2, and 13,14-dihydro-15-keto-PGF2α were assessed. RESULT(S): Comparison of the endometria of patients with UIF and the controls showed no statistically significant differences. When the endometria of patients with RIF were compared with the controls, thromboxane B2 (TXB2) was found significantly higher (843.1 pg/mg vs. 133.5 pg/mg). When the endometria of patients with RM were compared with controls, 13,14-dihydro-15-keto PGF2α and TXB2 were found significantly higher (3907.30 pg/mg vs. 17.80 pg/mg and 858.7 pg/mg vs. 133.5 pg/mg respectively). CONCLUSION(S): We identified increased endometrial presence of TXB2 in patients with RM and RIF, and 13,14-dihydro-15-keto PGF2α in patients with RM. Although common ground is observed for RM and RIF, prostanoids, on the other hand, might make their own contribution to endometrial receptivity as important as genes and proteins. Attempts to normalize the prostaglandin profile of the endometrium via enzymatic activity can open new therapeutic options.

Mutations and Copy Number Alterations in IDH Wild-Type Glioblastomas Are Shaped by Different Oncogenic Mechanisms.

Little is known about the mutational processes that shape the genetic landscape of gliomas. Numerous mutational processes leave marks on the genome in the form of mutations, copy number alterations, rearrangements or their combinations. To explore gliomagenesis, we hypothesized that gliomas with different underlying oncogenic mechanisms would have differences in the burden of various forms of these genomic alterations. This was an analysis on adult diffuse gliomas, but IDH-mutant gliomas as well as diffuse midline gliomas H3-K27M were excluded to search for the possible presence of new entities among the very heterogenous group of IDH-WT glioblastomas. The cohort was divided into two molecular subsets: (1) Molecularly-defined GBM (mGBM) as those that carried molecular features of glioblastomas (including TERT promoter mutations, 7/10 pattern, or EGFR-amplification), and (2) those who did not (others). Whole exome sequencing was performed for 37 primary tumors and matched blood samples as well as 8 recurrences. Single nucleotide variations (SNV), short insertion or deletions (indels) and copy number alterations (CNA) were quantified using 5 quantitative metrics (SNV burden, indel burden, copy number alteration frequency-wGII, chromosomal arm event ratio-CAER, copy number amplitude) as well as 4 parameters that explored underlying oncogenic mechanisms (chromothripsis, double minutes, microsatellite instability and mutational signatures). Findings were validated in the TCGA pan-glioma cohort. mGBM and "Others" differed significantly in their SNV (only in the TCGA cohort) and CNA metrics but not indel burden. SNV burden increased with increasing age at diagnosis and at recurrences and was driven by mismatch repair deficiency. On the contrary, indel and CNA metrics remained stable over increasing age at diagnosis and with recurrences. Copy number alteration frequency (wGII) correlated significantly with chromothripsis while CAER and CN amplitude correlated significantly with the presence of double minutes, suggesting separate underlying mechanisms for different forms of CNA.

pathfindR: An R Package for Comprehensive Identification of Enriched Pathways in Omics Data Through Active Subnetworks.

Pathway analysis is often the first choice for studying the mechanisms underlying a phenotype. However, conventional methods for pathway analysis do not take into account complex protein-protein interaction information, resulting in incomplete conclusions. Previously, numerous approaches that utilize protein-protein interaction information to enhance pathway analysis yielded superior results compared to conventional methods. Hereby, we present pathfindR, another approach exploiting protein-protein interaction information and the first R package for active-subnetwork-oriented pathway enrichment analyses for class comparison omics experiments. Using the list of genes obtained from an omics experiment comparing two groups of samples, pathfindR identifies active subnetworks in a protein-protein interaction network. It then performs pathway enrichment analyses on these identified subnetworks. To further reduce the complexity, it provides functionality for clustering the resulting pathways. Moreover, through a scoring function, the overall activity of each pathway in each sample can be estimated. We illustrate the capabilities of our pathway analysis method on three gene expression datasets and compare our results with those obtained from three popular pathway analysis tools. The results demonstrate that literature-supported disease-related pathways ranked higher in our approach compared to the others. Moreover, pathfindR identified additional pathways relevant to the conditions that were not identified by other tools, including pathways named after the conditions.

Potential Marker Pathways in the Endometrium That May Cause Recurrent Implantation Failure.

The aim of this prospective cohort study was to identify altered biologic processes in the endometrium that may be potential markers of receptive endometrium in patients with repeated implantation failure (RIF) as compared with fertile controls. The study was conducted in a university-affiliated in vitro fertilization (IVF) gynecology clinic and molecular biology and genetics laboratory. Healthy fertile controls (n = 24) and patients with RIF (n = 24) were recruited. Window of implantation gene profiling associated with RIF was performed. Six hundred forty-one differentially expressed genes were identified, and 44 pathways were found enriched. Upon clustering of the enriched pathways, 9 representative pathways were established. The important pathways that were identified included circadian rhythm, pathways in cancer, proteasome, complement and coagulation cascades, citrate cycle, adherens junction, immune system and inflammation, cell cycle, and renin-angiotensin system. The involvement of the circadian rhythm pathway and other related pathways may alter the endometrium's functioning to ultimately cause RIF. Furthermore, we found that the pathogenesis of RIF was multifaceted and that numerous processes were involved. We believe that a better understanding of the underlying mechanisms of RIF will ultimately give rise to better treatment opportunities and to better outcomes in IVF.

Meningiomas Display a Specific Immunoexpression Pattern in a Rostrocaudal Gradient: An Analysis of 366 Patients.

BACKGROUND: Meningiomas are heterogeneous, with differences in anatomical, histopathological, and clinical characteristics. Such spatial variability in meningioma biology is thought to result from differences in the expression of critical developmental regulators. We hypothesized that the variability in meningioma biology would follow gradients such as in embryology and tested a cohort of 366 meningiomas for histopathological and immunohistochemical gradients. METHODS: The medical records from 366 patients treated for meningiomas from 2003 to 2016 were retrospectively analyzed for age, gender, anatomical localization, recurrence-free survival, overall survival, histopathological diagnosis, and immunohistochemistry findings for 6 markers: epithelial membrane antigen (EMA), progesterone receptor (PR), CD34, S100, p53, and Ki-67 labeling index. RESULTS: EMA, PR, S100, p53, and CD34 were expressed in 94%, 73%, 49%, 26%, and 23% of the tumors, respectively. p53 expression correlated positively with Ki-67 and World Health Organization (WHO) grade (rτ = 0.31 and rτ = 0.4, respectively). PR positivity correlated inversely with S100, p53, Ki-67, and WHO grade (rτ = -0.19, rτ = -0.14, rτ = -0.15, and rτ = -0.16, respectively). All secretory meningiomas were positive for EMA and PR and negative for S100, and this pattern exhibited a rostrocaudal gradient. The overall proportion of EMA+PR+S100- cases was significantly lower in the cranial vault (30.3%) than in the skull base (45.89%; P = 0.021). The proportion of WHO grade II-III tumors was greater in cranial vault than in skull base meningiomas. CONCLUSIONS: Unsupervised methods detected an association between the anatomical location and tumor biology in meningiomas. Unlike the categorical associations that former studies had indicated, the present study revealed a rostrocaudal gradient in both the cranial vault and the skull base, correlating with human developmental biology.

Whole exome sequencing-based analysis to identify DNA damage repair deficiency as a major contributor to gliomagenesis in adult diffuse gliomas.

OBJECTIVE: Processes that cause or contribute to cancer, such as aging, exposure to carcinogens, or DNA damage repair deficiency (DDRd), create predictable and traceable nucleotide alterations in one's genetic code (termed "mutational signatures"). Large studies have previously identified various such mutational signatures across cancers that can be attributed to the specific causative processes. To gain further insight into the processes in glioma development, the authors analyzed mutational signatures in adult diffuse gliomas (DGs). METHODS: Twenty-five DGs and paired blood samples were whole exome sequenced. Somatic mutational signatures were identified using 2 different methods. Associations of the signatures with age at diagnosis, molecular subset, and mutational load were investigated. As DDRd-related signatures were frequently observed, germline and somatic DDR gene mutations as well as microsatellite instability (MSI) status were determined for all samples. For validation of signature prevalence, publicly available data from The Cancer Genome Atlas (TCGA) were used. RESULTS: Each tumor had a unique combination of signatures. The most common signatures were signature 1 (88%, aging related), signature 3 (52%, homologous recombination related), and signature 15 (56%, mismatch repair related). Eighty-four percent of the tumors contained at least 1 DDRd signature. The findings were validated using public TCGA data. The weight of signature 1 positively correlated with age (r = 0.43) while cumulative weight of DDRd signatures negatively correlated with age (r = -0.16). Each subject had at least 1 germline/somatic alteration in a DDR gene, the most common being the risk single nucleotide polymorphism rs1800734 in MLH1. The rs1800734-AA genotype had a higher cumulative DDRd weight as well as higher mutational load; TP53 was the most common somatically altered DDR gene. MSI was observed in 24% of the tumors. No significant associations of MSI status with mutational load, rs1800734, or the cumulative weight of DDRd signatures were identified. CONCLUSIONS: Current findings suggest that DDRd may act as a fundamental mechanism in gliomagenesis rather than being a random, secondary event.

Use of telomerase promoter mutations to mark specific molecular subsets with reciprocal clinical behavior in IDH mutant and IDH wild-type diffuse gliomas.

OBJECTIVE Recent studies have established that hemispheric diffuse gliomas may be grouped into subsets on the basis of molecular markers; these subsets are loosely correlated with the histopathological diagnosis but are strong predictors of clinical tumor behavior. Based on an analysis of molecular and clinical parameters, the authors hypothesized that mutations of the telomerase promoter (TERTp-mut) mark separate oncogenic programs among isocitrate dehydrogenase 1 and/or 2 (IDH) mutant (IDH-mut) and IDH wild-type (IDH-wt) diffuse gliomas independent of histopathology or WHO grade. METHODS Four molecular subsets of the combined statuses of IDH and TERT-promoter mutations (double mutant, IDH only, TERT only, and double negative) were defined. Differences in age, anatomical location, molecular genetics, and survival rates in a surgical cohort of 299 patients with a total of 356 hemispheric diffuse gliomas (WHO Grade II, III, or IV) were analyzed. RESULTS TERTp-mut were present in 38.8% of IDH-mut and 70.2% of IDH-wt gliomas. The mutational status was stable in each patient at 57 recurrence events over a 2645-month cumulative follow-up period. Among patients with IDH-mut gliomas, those in the double-mutant subset had better survival and a lower incidence of malignant degeneration than those in the IDH-only subset. Of patients in the double-mutant subset, 96.3% were also positive for 1p/19q codeletions. All patients with 1p/19q codeletions had TERTp-mut. In patients with IDH-mut glioma, epidermal growth factor receptor or phosphatase and tensin homolog mutations were not observed, and copy-number variations were uncommon. Among IDH-wt gliomas, the TERT-only subset was associated with significantly higher age, higher Ki-67 labeling index, primary glioblastoma-specific oncogenic changes, and poor survival. The double-negative subset was genetically and biologically heterogeneous. Survival analyses (Kaplan-Meier, multivariate, and regression-tree analyses) confirmed that patients in the 4 molecular subsets had distinct prognoses. CONCLUSIONS Molecular subsets result in different tumor biology and clinical behaviors in hemispheric diffuse gliomas.

ALPK3 gene mutation in a patient with congenital cardiomyopathy and dysmorphic features.

Primary cardiomyopathy is one of the most common inherited cardiac diseases and harbors significant phenotypic and genetic heterogeneity. Because of this, genetic testing has become standard in treatment of this disease group. Indeed, in recent years, next-generation DNA sequencing has found broad applications in medicine, both as a routine diagnostic tool for genetic disorders and as a high-throughput discovery tool for identifying novel disease-causing genes. We describe a male infant with primary dilated cardiomyopathy who was diagnosed using intrauterine echocardiography and found to progress to hypertrophic cardiomyopathy after birth. This proband was born to a nonconsanguineous family with a past history of a male fetus that died because of cardiac abnormalities at 30 wk of gestation. Using whole-exome sequencing, a novel homozygous frameshift mutation (c.2018delC; p.Gln675SerfsX30) in ALPK3 was identified and confirmed with Sanger sequencing. Heterozygous family members were normal with echocardiographic examination. To date, only two studies have reported homozygous pathogenic variants of ALPK3, with a total of seven affected individuals with cardiomyopathy from four unrelated consanguineous families. We include a discussion of the patient's phenotypic features and a review of relevant literature findings.

IDH-mutant glioma specific association of rs55705857 located at 8q24.21 involves MYC deregulation.

The single nucleotide polymorphism rs55705857, located in a non-coding but evolutionarily conserved region at 8q24.21, is strongly associated with IDH-mutant glioma development and was suggested to be a causal variant. However, the molecular mechanism underlying this association has remained unknown. With a case control study in 285 gliomas, 316 healthy controls, 380 systemic cancers, 31 other CNS-tumors, and 120 IDH-mutant cartilaginous tumors, we identified that the association was specific to IDH-mutant gliomas. Odds-ratios were 9.25 (5.17-16.52; 95% CI) for IDH-mutated gliomas and 12.85 (5.94-27.83; 95% CI) for IDH-mutated, 1p/19q co-deleted gliomas. Decreasing strength with increasing anaplasia implied a modulatory effect. No somatic mutations were noted at this locus in 114 blood-tumor pairs, nor was there a copy number difference between risk-allele and only-ancestral allele carriers. CCDC26 RNA-expression was rare and not different between the two groups. There were only minor subtype-specific differences in common glioma driver genes. RNA sequencing and LC-MS/MS comparisons pointed to significantly altered MYC-signaling. Baseline enhancer activity of the conserved region specifically on the MYC promoter and its further positive modulation by the SNP risk-allele was shown in vitro. Our findings implicate MYC deregulation as the underlying cause of the observed association.

A rare case of congenital fibrosis of extraocular muscle type 1A due to KIF21A mutation with Marcus Gunn jaw-winking phenomenon.

BACKGROUND: Congenital fibrosis of the extraocular muscles (CFEOM1) is classically a congenital, non-progressive, restrictive strabismus syndrome characterized by bilateral ptosis and ophthalmoplegia with an infraducted position of the globes. This autosomal dominant syndrome is caused by mutations in the KIF21A gene. METHODS AND RESULTS: In this report we describe a 5-year-old boy, and his mother, both of whom have a mutation in the KIF21A gene, who possesses typical features of CFEOM1 syndrome. Besides displaying typical features of CFEOM1, he demonstrated Marcus Gunn jaw-winking phenomenon. The patient additionally had a positive family history of such features. CONCLUSION: This is first report of the coexistence of CFEOM and Marcus Gunn jaw-winking phenomenon in a patient with a KIF21A mutation from Turkey. We explain the phenotypic findings associated with mutations in KIF21A including CFEOM1A and Marcus Gunn jaw-winking phenomenon.