Jenkins, Robert B.

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.

Abstract Background The shared inherited genetic contribution to risk of different cancers is not fully known. In this study, we leverage results from 12 cancer genome-wide association studies (GWAS) to quantify pairwise genome-wide genetic correlations across cancers and identify novel cancer susceptibility loci. Methods We collected GWAS summary statistics for 12 solid cancers based on 376 759 participants with cancer and 532 864 participants without cancer of European ancestry. The included cancer types were breast, colorectal, endometrial, esophageal, glioma, head and neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancers. We conducted cross-cancer GWAS and transcriptome-wide association studies to discover novel cancer susceptibility loci. Finally, we assessed the extent of variant-specific pleiotropy among cancers at known and newly identified cancer susceptibility loci. Results We observed widespread but modest genome-wide genetic correlations across cancers. In cross-cancer GWAS and transcriptome-wide association studies, we identified 15 novel cancer susceptibility loci. Additionally, we identified multiple variants at 77 distinct loci with strong evidence of being associated with at least 2 cancer types by testing for pleiotropy at known cancer susceptibility loci. Conclusions Overall, these results suggest that some genetic risk variants are shared among cancers, though much of cancer heritability is cancer-specific and thus tissue-specific. The increase in statistical power associated with larger sample sizes in cross-disease analysis allows for the identification of novel susceptibility regions. Future studies incorporating data on multiple cancer types are likely to identify additional regions associated with the risk of multiple cancer types.

Abstract Background Mendelian randomization (MR) studies are susceptible to metadata errors (e.g. incorrect specification of the effect allele column) and other analytical issues that can introduce substantial bias into analyses. We developed a quality control (QC) pipeline for the Fatty Acids in Cancer Mendelian Randomization Collaboration (FAMRC) that can be used to identify and correct for such errors. Methods We collated summary association statistics from fatty acid and cancer genome-wide association studies (GWAS) and subjected the collated data to a comprehensive QC pipeline. We identified metadata errors through comparison of study-specific statistics to external reference data sets (the National Human Genome Research Institute-European Bioinformatics Institute GWAS catalogue and 1000 genome super populations) and other analytical issues through comparison of reported to expected genetic effect sizes. Comparisons were based on three sets of genetic variants: (i) GWAS hits for fatty acids, (ii) GWAS hits for cancer and (iii) a 1000 genomes reference set. Results We collated summary data from 6 fatty acid and 54 cancer GWAS. Metadata errors and analytical issues with the potential to introduce substantial bias were identified in seven studies (11.6%). After resolving metadata errors and analytical issues, we created a data set of 219 842 genetic associations with 90 cancer types, generated in analyses of 566 665 cancer cases and 1 622 374 controls. Conclusions In this large MR collaboration, 11.6% of included studies were affected by a substantial metadata error or analytical issue. By increasing the integrity of collated summary data prior to their analysis, our protocol can be used to increase the reliability of downstream MR analyses. Our pipeline is available to other researchers via the CheckSumStats package (https://github.com/MRCIEU/CheckSumStats).

AbstractBackgroundMeningiomas are the most common primary central nervous system (CNS) tumor in adults and CNS World Health Organization grade 2 (atypical) meningiomas show an intermediate risk of recurrence/progression. Molecular parameters are needed to better inform management following gross total resection (GTR).MethodsWe performed comprehensive genomic analysis of tumor tissue from 63 patients who underwent radiologically confirmed GTR of a primary grade 2 meningioma, including a CLIA-certified target next-generation sequencing panel (n = 61), chromosomal microarray (n = 63), genome-wide methylation profiling (n = 62), H3K27me3 immunohistochemistry (n = 62), and RNA-sequencing (n = 19). Genomic features were correlated with long-term clinical outcomes (median follow-up: 10 years) using Cox proportional hazards regression modeling and published molecular prognostic signatures were evaluated.ResultsThe presence of specific copy number variants (CNVs), including -1p, -10q, -7p, and -4p, was the strongest predictor of decreased recurrence-free survival (RFS) within our cohort (P < .05). NF2 mutations were frequent (51%) but did not show a significant association with RFS. DNA methylation-based classification assigned tumors to DKFZ Heidelberg benign (52%) or intermediate (47%) meningioma subclasses and was not associated with RFS. H3K27 trimethylation (H3K27me3) was unequivocally lost in 4 tumors, insufficient for RFS analysis. Application of published integrated histologic/molecular grading systems did not improve prediction of recurrence risk over the presence of -1p or -10q alone.ConclusionsCNVs are strong predictors of RFS in grade 2 meningiomas following GTR. Our study supports incorporation of CNV profiling into clinical evaluation to better guide postoperative patient management, which can be readily implemented using existing, clinically validated technologies.

Abstract BACKGROUND The topologically associated domain (TAD) on 8q24 surrounding MYC is important for many types of cancer. The region is also important in the formation of IDH-mutant glioma; the risk allele rs55705857 being associated with the development of these tumors. The rs55705857 risk allele significantly increased the formation of tumors in an IDH-mutant mouse model. Cerebral organoids have been found to recapitulate the early development of the brain. We hypothesized that the rs55705857 risk allele may alter the phenotype of cerebral organoids. METHODS Isogenic induced pluripotent stem cells (iPSCs) were developed with (n=3) and without (n=2) the rs55705857 risk allele from the parental cell line PGP-1 (GM23338). These cells were cultured and differentiated into unpatterned cerebral organoids and observed throughout differentiation. Organoids were fixed, sectioned, and stained by H&E and immunofluorescence. RESULTS Organoids containing the risk allele were smaller throughout development, starting at the embryoid body stage and continuing up to over 75 days in culture, compared to two isogenic non-risk lines. At 4 weeks the risk allele containing organoids also contained a larger number of small rosette-like neuroepithelia, instead of the more continuous neuroepithelia observed in the non-risk allele organoids. At 75 days, the organoids with the risk allele maintained a larger proportion of Sox2 positive cells, which are enriched in the tight clusters of small cells with round nuclei. CONCLUSIONS The risk allele rs55705857 in the MYC TAD at 8q24 is altering the development of unpatterned cerebral organoids. This difference is visible as early as embryoid body development and persists through neuroepithelial development. Further studies of these differences may help our understanding of how rs55705857 accelerates IDH-mutant glioma development.

Abstract Missing MRI sequences are common in radiology due to differences in protocols, time constraints, and image artifacts. Missing sequences are a nuisance for subsequent algorithms, such as automated segmentation algorithms requiring multiple MRI sequences (e.g., HD-GLIO, which requires T2, T1, post-contrast T1 (T1c), and FLAIR). Accurate tumor segmentation is a fundamental step in extracting radiomics analysis to use, for example, to predict the IDH status of glioma or their overall survival. We validated a conditional (pix2pix) and super-resolution GAN (SOUP-GAN) for imputing MRI sequences on 285 adult glioma with complete MRI sequences available. We simulated missing T1 sequences and subsequently imputed them from T1c using pix2pix and from sagittal T1 using SOUP-GAN. Using HD-GLIO, we obtained the true segmentations from the complete MRI sequences. HD-GLIO was also used after replacing the original T1 sequence with the one imputed with pix2pix and SOUP-GAN. Thus, there were three sets of segmentations for each subject. We compared T2 and T1c radiomics features extracted with Pyradiomics using segmentations obtained with complete data with features extracted from segmentation obtained with the imputed T1 sequences using intraclass correlation coefficient (ICC) and Bland-Altman (BA) analysis. Median (interquartile range) ICC for shape features extracted from T2 and T1c were 0.82(0.76-0.92) using pix2pix and 0.89(0.86-0.94) using SOUP-GAN. Median ICC for first-order-statistics features were 0.90(0.83-0.93) for T2 and 0.90(0.88-0.92) for T1c using pix2pix, and 0.61(0.22-0.82) for T2 and 0.42(0.34-0.57) for T1c using SOUP-GAN. Median ICC for texture features were 0.86(0.82-0.87) for T2 and 0.87(0.79-0.90) for T1c using pix2pix, and 0.69(0.47-0.81) for T2 and 0.76(0.65-0.83) for T1c using SOUP-GAN. We show that T2 and T1c shape features (e.g., volumetric analysis) can be adequately extracted from segmentations obtained with an algorithm requiring multiple MRI sequences using T1 sequences imputed with GANs; however, the effect of imputation differed for first-order-statistics and texture features.

Abstract Diffuse hemispheric glioma, H3 G34-mutant (DHG_H3G34) is a novel tumor type in the 2021 WHO CNS classification. We describe a comprehensive sequencing and high-resolution genome-wide copy number analysis of a series of cases clinically tested by a single laboratory (2018-2022). Cases included tumors from 47 unique patients (1 reportedly recurrent) that had an H3-3A G34 mutation detected using an 187-gene mutation and fusion targeted neuro-oncology NGS panel (n=47). A subset (n=18) of cases was also tested Oncoscan chromosomal microarray. Median age at testing was 20 years (range, 12-50). H3-3A G34 mutations included G34R (n=44; 94%), G34V (n=2) and a novel G34E variant considered to be likely clinically relevant (n=1). All DHG_H3G34 were hemispheric tumors, and one tumor was multifocal with midline involvement. Concurrent mutations recurrently involved ATRX (n=38), TERT promoter (n=3; mutually exclusive with ATRX mutations), TP53 (n=43), PDGFRA (n=26), PTEN (n=5), NF1 (n=4, all patients over 20 years), PIK3CA (n=4) and CDKN2A (n=3). The single recurrent tumor also had an MSH6 mutation, and in addition to 3 other cases, had increased number of sequence variants suggestive of a high tumor mutational burden. A single case also showed an FGFR3::FAM184B (exon 17::exon 2) fusion. All cases with available chromosomal microarray data had unbalanced genomes with multiple chromosomal gains/losses, and regions of copy-neutral loss of heterozygosity. The most frequent recurrent (at least 50% of cases) abnormalities were losses involving 3q, 4q, 10q, 13q and 18q, and 17p copy-neutral loss of heterozygosity encompassing TP53 (all but one case had concomitant TP53 mutation), a pattern reminiscent of IDH-mutant astrocytomas. Oncogene amplifications (PDGFRA and PIK3CA) were detected in 3 (of 18) cases. CDKN2A/B deletion was observed in 8 cases (2 homozygous). In conclusion, we describe additional sequence and copy-number abnormalities in DHG_H3G34, expanding the spectrum of genetic changes associated with this novel tumor type.

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 Idh1 R132H -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.