Clonally Selected Lines After CRISPR-Cas Editing Are Not Isogenic.

The CRISPR-Cas9 system has enabled researchers to precisely modify/edit the sequence of a genome. A typical editing experiment consists of two steps: (1) editing cultured cells; (2) cell cloning and selection of clones with and without intended edit, presumed to be isogenic. The application of CRISPR-Cas9 system may result in off-target edits, whereas cloning will reveal culture-acquired mutations. We analyzed the extent of the former and the latter by whole genome sequencing in three experiments involving separate genomic loci and conducted by three independent laboratories. In all experiments we hardly found any off-target edits, whereas detecting hundreds to thousands of single nucleotide mutations unique to each clone after relatively short culture of 10-20 passages. Notably, clones also differed in copy number alterations (CNAs) that were several kb to several mb in size and represented the largest source of genomic divergence among clones. We suggest that screening of clones for mutations and CNAs acquired in culture is a necessary step to allow correct interpretation of DNA editing experiments. Furthermore, since culture associated mutations are inevitable, we propose that experiments involving derivation of clonal lines should compare a mix of multiple unedited lines and a mix of multiple edited lines.

A noncoding single-nucleotide polymorphism at 8q24 drives IDH1-mutant glioma formation.

Establishing causal links between inherited polymorphisms and cancer risk is challenging. Here, we focus on the single-nucleotide polymorphism rs55705857, which confers a sixfold greater risk of isocitrate dehydrogenase (IDH)-mutant low-grade glioma (LGG). We reveal that rs55705857 itself is the causal variant and is associated with molecular pathways that drive LGG. Mechanistically, we show that rs55705857 resides within a brain-specific enhancer, where the risk allele disrupts OCT2/4 binding, allowing increased interaction with the Myc promoter and increased Myc expression. Mutating the orthologous mouse rs55705857 locus accelerated tumor development in an Idh1R132H-driven LGG mouse model from 472 to 172 days and increased penetrance from 30% to 75%. Our work reveals mechanisms of the heritable predisposition to lethal glioma in ~40% of LGG patients.

Inherited genetics of adult diffuse glioma and polygenic risk scores-a review.

Knowledge about inherited and acquired genetics of adult diffuse glioma has expanded significantly over the past decade. Genomewide association studies (GWAS) stratified by histologic subtype identified six germline variants that were associated specifically with glioblastoma (GBM) and 12 that were associated with lower grade glioma. A GWAS performed using the 2016 WHO criteria, stratifying patients by IDH mutation and 1p/19q codeletion (as well as TERT promoter mutation), discovered that many of the known variants are associated with specific WHO glioma subtypes. In addition, the GWAS stratified by molecular group identified two additional novel regions: variants in D2HGDH that were associated with tumors that had an IDH mutation and a variant near FAM20C that was associated with tumors that had both IDH mutation and 1p/19q codeletion. The results of these germline associations have been used to calculate polygenic risk scores, from which to estimate relative and absolute risk of overall glioma and risk of specific glioma subtypes. We will review the concept of polygenic risk models and their potential clinical utility, as well as discuss the published adult diffuse glioma polygenic risk models. To date, these prior genetic studies have been done on European populations. Using the published glioma polygenic risk model, we show that the genetic associations published to date do not generalize across genetic ancestries, demonstrating that genetic studies need to be done on more diverse populations.

Functional analysis of low-grade glioma genetic variants predicts key target genes and transcription factors.

BACKGROUND: Large-scale genome-wide association studies (GWAS) have implicated thousands of germline genetic variants in modulating individuals' risk to various diseases, including cancer. At least 25 risk loci have been identified for low-grade gliomas (LGGs), but their molecular functions remain largely unknown. METHODS: We hypothesized that GWAS loci contain causal single nucleotide polymorphisms (SNPs) that reside in accessible open chromatin regions and modulate the expression of target genes by perturbing the binding affinity of transcription factors (TFs). We performed an integrative analysis of genomic and epigenomic data from The Cancer Genome Atlas and other public repositories to identify candidate causal SNPs within linkage disequilibrium blocks of LGG GWAS loci. We assessed their potential regulatory role via in silico TF binding sequence perturbations, convolutional neural network trained on TF binding data, and simulated annealing-based interpretation methods. RESULTS: We built an interactive website (http://education.knoweng.org/alg3/) summarizing the functional footprinting of 280 variants in 25 LGG GWAS regions, providing rich information for further computational and experimental scrutiny. We identified as case studies PHLDB1 and SLC25A26 as candidate target genes of rs12803321 and rs11706832, respectively, and predicted the GWAS variant rs648044 to be the causal SNP modulating ZBTB16, a known tumor suppressor in multiple cancers. We showed that rs648044 likely perturbed the binding affinity of the TF MAFF, as supported by RNA interference and in vitro MAFF binding experiments. CONCLUSIONS: The identified candidate (causal SNP, target gene, TF) triplets and the accompanying resource will help accelerate our understanding of the molecular mechanisms underlying genetic risk factors for gliomas.

Adult diffuse glioma GWAS by molecular subtype identifies variants in D2HGDH and FAM20C.

BACKGROUND: Twenty-five germline variants are associated with adult diffuse glioma, and some of these variants have been shown to be associated with particular subtypes of glioma. We hypothesized that additional germline variants could be identified if a genome-wide association study (GWAS) were performed by molecular subtype. METHODS: A total of 1320 glioma cases and 1889 controls were used in the discovery set and 799 glioma cases and 808 controls in the validation set. Glioma cases were classified into molecular subtypes based on combinations of isocitrate dehydrogenase (IDH) mutation, telomerase reverse transcriptase (TERT) promoter mutation, and 1p/19q codeletion. Logistic regression was applied to the discovery and validation sets to test for associations of variants with each of the subtypes. A meta-analysis was subsequently performed using a genome-wide P-value threshold of 5 × 10-8. RESULTS: Nine variants in or near D-2-hydroxyglutarate dehydrogenase (D2HGDH) on chromosome 2 were genome-wide significant in IDH-mutated glioma (most significant was rs5839764, meta P = 2.82 × 10-10). Further stratifying by 1p/19q codeletion status, one variant in D2HGDH was genome-wide significant in IDH-mutated non-codeleted glioma (rs1106639, meta P = 4.96 × 10-8). Further stratifying by TERT mutation, one variant near FAM20C (family with sequence similarity 20, member C) on chromosome 7 was genome-wide significant in gliomas that have IDH mutation, TERT mutation, and 1p/19q codeletion (rs111976262, meta P = 9.56 × 10-9). Thirty-six variants in or near GMEB2 on chromosome 20 near regulator of telomere elongation helicase 1 (RTEL1) were genome-wide significant in IDH wild-type glioma (most significant was rs4809313, meta P = 2.60 × 10-10). CONCLUSIONS: Performing a GWAS by molecular subtype identified 2 new regions and a candidate independent region near RTEL1, which were associated with specific glioma molecular subtypes.

Using germline variants to estimate glioma and subtype risks.

BACKGROUND: Twenty-five single nucleotide polymorphisms (SNPs) are associated with adult diffuse glioma risk. We hypothesized that the inclusion of these 25 SNPs with age at diagnosis and sex could estimate risk of glioma as well as identify glioma subtypes. METHODS: Case-control design and multinomial logistic regression were used to develop models to estimate the risk of glioma development while accounting for histologic and molecular subtypes. Case-case design and logistic regression were used to develop models to predict isocitrate dehydrogenase (IDH) mutation status. A total of 1273 glioma cases and 443 controls from Mayo Clinic were used in the discovery set, and 852 glioma cases and 231 controls from UCSF were used in the validation set. All samples were genotyped using a custom Illumina OncoArray. RESULTS: Patients in the highest 5% of the risk score had more than a 14-fold increase in relative risk of developing an IDH mutant glioma. Large differences in lifetime absolute risk were observed at the extremes of the risk score percentile. For both IDH mutant 1p/19q non-codeleted glioma and IDH mutant 1p/19q codeleted glioma, the lifetime risk increased from almost null to 2.3% and almost null to 1.7%, respectively. The SNP-based model that predicted IDH mutation status had a validation concordance index of 0.85. CONCLUSIONS: These results suggest that germline genotyping can provide new tools for the initial management of newly discovered brain lesions. Given the low lifetime risk of glioma, risk scores will not be useful for population screening; however, they may be useful in certain clinically defined high-risk groups.

The Thrombin Receptor Restricts Subventricular Zone Neural Stem Cell Expansion and Differentiation.

Thrombin is frequently increased in the CNS after injury yet little is known regarding its effects on neural stem cells. Here we show that the subventricular zone (SVZ) of adult mice lacking the high affinity receptor for thrombin, proteinase activated receptor 1 (PAR1), show increased numbers of Sox2+ and Ki-67+ self-renewing neural stem cells (NSCs) and Olig2+ oligodendrocyte progenitors. SVZ NSCs derived from PAR1-knockout mice, or treated with a PAR1 small molecule inhibitor (SCH79797), exhibited enhanced capacity for self-renewal in vitro, including increases in neurosphere formation and BrdU incorporation. PAR1-knockout SVZ monolayer cultures contained more Nestin, NG2+ and Olig2+ cells indicative of enhancements in expansion and differentiation towards the oligodendrocyte lineage. Cultures of NSCs lacking PAR1 also expressed higher levels of myelin basic protein, proteolipid protein and glial fibrillary acidic protein upon differentiation. Complementing these findings, the corpus callosum and anterior commissure of adult PAR1-knockout mice contained greater numbers of Olig2+ progenitors and CC1+ mature oligodendrocytes. Together these findings highlight PAR1 inhibition as a means to expand adult SVZ NSCs and to promote an increased number of mature myelinating oligodendrocytes in vivo that may be of particular benefit in the context of neural injury where PAR1 agonists such as thrombin are deregulated.

Prognostic significance of multiple kallikreins in high-grade astrocytoma.

BACKGROUND: Kallikreins have clinical value as prognostic markers in a subset of malignancies examined to date, including kallikrein 3 (prostate specific antigen) in prostate cancer. We previously demonstrated that kallikrein 6 is expressed at higher levels in grade IV compared to grade III astrocytoma and is associated with reduced survival of GBM patients. METHODS: In this study we determined KLK1, KLK6, KLK7, KLK8, KLK9 and KLK10 protein expression in two independent tissue microarrays containing 60 grade IV and 8 grade III astrocytoma samples. Scores for staining intensity, percent of tumor stained and immunoreactivity scores (IR, product of intensity and percent) were determined and analyzed for correlation with patient survival. RESULTS: Grade IV glioma was associated with higher levels of kallikrein-immunostaining compared to grade III specimens. Univariable Cox proportional hazards regression analysis demonstrated that elevated KLK6- or KLK7-IR was associated with poor patient prognosis. In addition, an increased percent of tumor immunoreactive for KLK6 or KLK9 was associated with decreased survival in grade IV patients. Kaplan-Meier survival analysis indicated that patients with KLK6-IR < 10, KLK6 percent tumor core stained < 3, or KLK7-IR < 9 had a significantly improved survival. Multivariable analysis indicated that the significance of these parameters was maintained even after adjusting for gender and performance score. CONCLUSIONS: These data suggest that elevations in glioblastoma KLK6, KLK7 and KLK9 protein have utility as prognostic markers of patient survival.

The thrombin receptor is a critical extracellular switch controlling myelination.

Hemorrhagic white matter injuries in the perinatal period are a growing cause of cerebral palsy yet no neuroprotective strategies exist to prevent the devastating motor and cognitive deficits that ensue. We demonstrate that the thrombin receptor (protease-activated receptor 1, PAR1) exhibits peak expression levels in the spinal cord at term and is a critical regulator of the myelination continuum from initiation to the final levels achieved. Specifically, PAR1 gene deletion resulted in earlier onset of spinal cord myelination, including substantially more Olig2-positive oligodendrocytes, more myelinated axons, and higher proteolipid protein (PLP) levels at birth. In vitro, the highest levels of PAR1 were observed in oligodendrocyte progenitor cells (OPCs), being reduced with differentiation. In parallel, the expression of PLP and myelin basic protein (MBP), in addition to Olig2, were all significantly higher in cultures of PAR1-/- oligodendroglia. Moreover, application of a small molecule inhibitor of PAR1 (SCH79797) to OPCs in vitro increased PLP and MBP expression. Enhancements in myelination associated with PAR1 genetic deletion were also observed in adulthood as evidenced by higher amounts of MBP and thickened myelin sheaths across large, medium, and small diameter axons. Enriched spinal cord myelination in PAR1-/- mice was coupled to increases in extracellular-signal-regulated kinase 1/2 and AKT signaling developmentally. Nocturnal ambulation and rearing activity were also elevated in PAR1-/- mice. These studies identify the thrombin receptor as a powerful extracellular regulatory switch that could be readily targeted to improve myelin production in the face of white matter injury and disease.

Clinical significance and novel mechanism of action of kallikrein 6 in glioblastoma.

BACKGROUND: Kallikreins have prognostic value in specific malignancies, but few studies have addressed their clinical significance to glioblastoma multiforme (GBM). Kallikrein 6 (KLK6) is of potential high relevance to GBM, since it is upregulated at sites of CNS pathology and linked to reactive astrogliosis. Here we examine the clinical value of KLK6 as a prognostic indicator of GBM patient survival and its activity in promoting resistance to cytotoxic agents. METHODS: The association between patient survival and levels of KLK6 immunoreactivity were investigated in 60 grade IV astrocytoma tumor specimens. Levels of KLK6 RNA were also evaluated in a separate set of GBM patient tumors (n = 23). Recombinant KLK6 or enforced KLK6 overexpression in GBM cell lines was used to evaluate effects on astrocytoma cell survival. RESULTS: A range of KLK6 expression was observed across grade IV tumors, with higher levels a poor prognostic indicator of patient survival (P = .02) even after adjusting for gender and Eastern Cooperative Oncology Group performance scores (P = .01). KLK6 reduced the sensitivity of GBM cell lines to cytotoxic agents, including staurosporine and cisplatin, and to the current standard of patient care: radiotherapy or temozolomide alone or in combination. The ability of KLK6 to promote resistance to apoptosis was dependent on activation of the thrombin receptor, protease activated receptor 1. CONCLUSIONS: Taken together, these results indicate that elevated levels of KLK6 in GBM are likely to promote the resistance of tumor cells to cytotoxic agents and are an indicator of reduced patient postsurgical survival times.

Characterization and gene expression profiling in glioma cell lines with deletion of chromosome 19 before and after microcell-mediated restoration of normal human chromosome 19.

Nearly 10% of human gliomas are oligodendrogliomas. Deletion of chromosome arm 19q, often in conjunction with deletion of 1p, has been observed in 65-80% of these tumors. This has suggested the presence of a tumor suppressor gene located on the 19q arm. Chromosome 19 deletion is also of interest due to the better prognosis of patients with deletion, including longer survival and better response to chemotherapy, compared with patients without deletion. Two glioma cell lines with deletion of 19q were used for chromosome 19 microcell-mediated transfer, to assess the effect of replacing the deleted segment. Complementation with chromosome 19 significantly reduced the growth rate of the hybrid cells compared with the parental cell lines. Affymetrix U133 Plus 2.0 Gene Chip analysis was performed to measure and compare the expression of the chromosome 19 genes in the chromosome 19 hybrid cell lines to the parental cell line. Probes were considered significantly different when a P value <0.01 was seen in all of the cell line comparisons. Of 345 probes within the commonly deleted 19q region, seven genes (APOE, RCN3, FLJ10781, SAE1, STRN4, CCDC8, and BCL2L12) were identified as potential candidate genes. RT-PCR analysis of primary tumor specimens showed that several genes had significant differences when stratified by tumor morphology or deletion status. This suggests that one or more of these candidates may play a role in glioma formation or progression.