Search

Abstract Noonan Syndrome (NS) is associated with an increased risk of low-grade central nervous system tumours in children but only very rarely associated with high-grade gliomas. Here we describe the first reported case of a spinal high-grade astrocytoma with piloid features (HGAP) in a child with NS. This case was a diagnostic and treatment dilemma, prior to whole-genome germline and tumour sequencing, tumour transcriptome sequencing and DNA methylation analysis. The methylation profile matched strongly with HGAP and sequencing identified somatic FGFR1 and NF1 variants and a PTPN11 germline pathogenic variant. Therapeutic targets were identified but also alterations novel to HGAP such as differential expression of VEGFA and PD-L1. The germline PTPN11 finding has not been previously described in individuals with HGAP. This case underscores the power of precision medicine from a diagnostic, therapeutic and clinical management perspective, and describes an association between HGAP and NS which has not previously been reported.

Abstract Autosomal dominant polycystic kidney disease (ADPKD) is the most common monogenic cause of kidney failure and is primarily associated with PKD1 or PKD2. Approximately 10% of patients remain undiagnosed after standard genetic testing. We aimed to utilise short and long-read genome sequencing and RNA studies to investigate undiagnosed families. Patients with typical ADPKD phenotype and undiagnosed after genetic diagnostics were recruited. Probands underwent short-read genome sequencing, PKD1 and PKD2 coding and non-coding analyses and then genome-wide analysis. Targeted RNA studies investigated variants suspected to impact splicing. Those undiagnosed then underwent Oxford Nanopore Technologies long-read genome sequencing. From over 172 probands, 9 met inclusion criteria and consented. A genetic diagnosis was made in 8 of 9 (89%) families undiagnosed on prior genetic testing. Six had variants impacting splicing, five in non-coding regions of PKD1. Short-read genome sequencing identified novel branchpoint, AG-exclusion zone and missense variants generating cryptic splice sites and a deletion causing critical intron shortening. Long-read sequencing confirmed the diagnosis in one family. Most undiagnosed families with typical ADPKD have splice-impacting variants in PKD1. We describe a pragmatic method for diagnostic laboratories to assess PKD1 and PKD2 non-coding regions and validate suspected splicing variants through targeted RNA studies.

Abstract Predicting the impact of coding and noncoding variants on splicing is challenging, particularly in non-canonical splice sites, leading to missed diagnoses in patients. Existing splice prediction tools are complementary but knowing which to use for each splicing context remains difficult. Here, we describe Introme, which uses machine learning to integrate predictions from several splice detection tools, additional splicing rules, and gene architecture features to comprehensively evaluate the likelihood of a variant impacting splicing. Through extensive benchmarking across 21,000 splice-altering variants, Introme outperformed all tools (auPRC: 0.98) for the detection of clinically significant splice variants. Introme is available at https://github.com/CCICB/introme .

Abstract Background Molecular profiling of the tumour immune microenvironment (TIME) has enabled the rational choice of immunotherapies in some adult cancers. In contrast, the TIME of paediatric cancers is relatively unexplored. We speculated that a more refined appreciation of the TIME in childhood cancers, rather than a reliance on commonly used biomarkers such as tumour mutation burden (TMB), neoantigen load and PD-L1 expression, is an essential prerequisite for improved immunotherapies in childhood solid cancers. Methods We combined immunohistochemistry (IHC) with RNA sequencing and whole-genome sequencing across a diverse spectrum of high-risk paediatric cancers to develop an alternative, expression-based signature associated with CD8+ T-cell infiltration of the TIME. Furthermore, we explored transcriptional features of immune archetypes and T-cell receptor sequencing diversity, assessed the relationship between CD8+ and CD4+ abundance by IHC and deconvolution predictions and assessed the common adult biomarkers such as neoantigen load and TMB. Results A novel 15-gene immune signature, Immune Paediatric Signature Score (IPASS), was identified. Using this signature, we estimate up to 31% of high-risk cancers harbour infiltrating T-cells. In addition, we showed that PD-L1 protein expression is poorly correlated with PD-L1 RNA expression and TMB and neoantigen load are not predictive of T-cell infiltration in paediatrics. Furthermore, deconvolution algorithms are only weakly correlated with IHC measurements of T-cells. Conclusions Our data provides new insights into the variable immune-suppressive mechanisms dampening responses in paediatric solid cancers. Effective immune-based interventions in high-risk paediatric cancer will require individualised analysis of the TIME.

Abstract Diffuse midline gliomas (DMG) harbouring H3K27M mutation are paediatric tumours with a dismal outcome. Recently, a new subtype of midline gliomas has been described with similar features to DMG, including loss of H3K27 trimethylation, but lacking the canonical H3K27M mutation (H3-WT). Here, we report a cohort of five H3-WT tumours profiled by whole-genome sequencing, RNA sequencing and DNA methylation profiling and combine their analysis with previously published cases. We show that these tumours have recurrent and mutually exclusive mutations in either ACVR1 or EGFR and are characterised by high expression of EZHIP associated to its promoter hypomethylation. Affected patients share a similar poor prognosis as patients with H3K27M DMG. Global molecular analysis of H3-WT and H3K27M DMG reveal distinct transcriptome and methylome profiles including differential methylation of homeobox genes involved in development and cellular differentiation. Patients have distinct clinical features, with a trend demonstrating ACVR1 mutations occurring in H3-WT tumours at an older age. This in-depth exploration of H3-WT tumours further characterises this novel DMG, H3K27-altered sub-group, characterised by a specific immunohistochemistry profile with H3K27me3 loss, wild-type H3K27M and positive EZHIP. It also gives new insights into the possible mechanism and pathway regulation in these tumours, potentially opening new therapeutic avenues for these tumours which have no known effective treatment. This study has been retrospectively registered on clinicaltrial.gov on 8 November 2017 under the registration number NCT03336931 ( https://clinicaltrials.gov/ct2/show/NCT03336931 ).

IntroductionEpendymomas (EPN) are the third most common malignant brain cancer in children. Treatment strategies for pediatric EPN have remained unchanged over recent decades, with 10-year survival rates stagnating at just 67% for children aged 0-14 years. Moreover, a proportion of patients who survive treatment often suffer long-term neurological side effects as a result of therapy. It is evident that there is a need for safer, more effective treatments for pediatric EPN patients. There are ten distinct subgroups of EPN, each with their own molecular and prognostic features. To identify and facilitate the testing of new treatments for EPN, in vivo laboratory models representative of the diverse molecular subtypes are required. Here, we describe the establishment of a patient-derived orthotopic xenograft (PDOX) model of posterior fossa A (PFA) EPN, derived from a metastatic cranial lesion.MethodsPatient and PDOX tumors were analyzed using immunohistochemistry, DNA methylation profiling, whole genome sequencing (WGS) and RNA sequencing.ResultsBoth patient and PDOX tumors classified as PFA EPN by methylation profiling, and shared similar histological features consistent with this molecular subgroup. RNA sequencing revealed that gene expression patterns were maintained across the primary and metastatic tumors, as well as the PDOX. Copy number profiling revealed gains of chromosomes 7, 8 and 19, and loss of chromosomes 2q and 6q in the PDOX and matched patient tumor. No clinically significant single nucleotide variants were identified, consistent with the low mutation rates observed in PFA EPN. Overexpression of EZHIP RNA and protein, a common feature of PFA EPN, was also observed. Despite the aggressive nature of the tumor in the patient, this PDOX was unable to be maintained past two passages in vivo.DiscussionOthers who have successfully developed PDOX models report some of the lowest success rates for EPN compared to other pediatric brain cancer types attempted, with loss of tumorigenicity not uncommon, highlighting the challenges of propagating these tumors in the laboratory. Here, we discuss our collective experiences with PFA EPN PDOX model generation and propose potential approaches to improve future success in establishing preclinical EPN models.

Abstract Whole genome sequencing (WGS) has the potential to outperform clinical microarrays for the detection of structural variants (SV) including copy number variants (CNVs), but has been challenged by high false positive rates. Here we present ClinSV, a WGS based SV integration, annotation, prioritization, and visualization framework, which identified 99.8% of simulated pathogenic ClinVar CNVs > 10 kb and 11/11 pathogenic variants from matched microarrays. The false positive rate was low (1.5–4.5%) and reproducibility high (95–99%). In clinical practice, ClinSV identified reportable variants in 22 of 485 patients (4.7%) of which 35–63% were not detectable by current clinical microarray designs. ClinSV is available at https://github.com/KCCG/ClinSV .

Abstract The prognosis of recurrent malignant peripheral nerve sheath tumors (MPNST) is dismal, with surgical resection being the only definitive salvage therapy. Treatment with chemoradiation approaches has not significantly improved patient outcomes. Similarly, trials of therapies targeting MPNST genomic drivers have thus far been unsuccessful. Improved understanding of the molecular pathogenesis of MPNST indicates frequent activation of the mitogen-activated protein kinase (MAPK) cell signaling pathway. MEK inhibitors have shown activity in preclinical studies; however, their clinical efficacy has not been reported to date. We describe here a case of sustained complete response to MEK inhibition in an adolescent patient with a recurrent metastatic MPNST with multiple alterations in the MAPK pathway, guided by a precision oncology approach.

Summary: Complex somatic genomic rearrangements and copy number alterations are hallmarks of nearly all cancers. We have developed an algorithm, LINX, to aid interpretation of structural variant and copy number data derived from short-read, whole-genome sequencing. LINX classifies raw structural variant calls into distinct events and predicts their effect on the local structure of the derivative chromosome and the functional impact on affected genes. Visualizations facilitate further investigation of complex rearrangements. LINX allows insights into a diverse range of structural variation events and can reliably detect pathogenic rearrangements, including gene fusions, immunoglobulin enhancer rearrangements, intragenic deletions, and duplications. Uniquely, LINX also predicts chained fusions that we demonstrate account for 13% of clinically relevant oncogenic fusions. LINX also reports a class of inactivation events that we term homozygous disruptions that may be a driver mutation in up to 9% of tumors and may frequently affect PTEN, TP53, and RB1.

Abstract Background Basal-like breast cancer (BLBC) is a poorly characterised, heterogeneous disease. Patients are diagnosed with aggressive, high-grade tumours and often relapse with chemotherapy resistance. Detailed understanding of the molecular underpinnings of this disease is essential to the development of personalised therapeutic strategies. Inhibitor of differentiation 4 (ID4) is a helix-loop-helix transcriptional regulator required for mammary gland development. ID4 is overexpressed in a subset of BLBC patients, associating with a stem-like poor prognosis phenotype, and is necessary for the growth of cell line models of BLBC through unknown mechanisms. Methods Here, we have defined unique molecular insights into the function of ID4 in BLBC and the related disease high-grade serous ovarian cancer (HGSOC), by combining RIME proteomic analysis, ChIP-seq mapping of genomic binding sites and RNA-seq. Results These studies reveal novel interactions with DNA damage response proteins, in particular, mediator of DNA damage checkpoint protein 1 (MDC1). Through MDC1, ID4 interacts with other DNA repair proteins (γH2AX and BRCA1) at fragile chromatin sites. ID4 does not affect transcription at these sites, instead binding to chromatin following DNA damage. Analysis of clinical samples demonstrates that ID4 is amplified and overexpressed at a higher frequency in BRCA1-mutant BLBC compared with sporadic BLBC, providing genetic evidence for an interaction between ID4 and DNA damage repair deficiency. Conclusions These data link the interactions of ID4 with MDC1 to DNA damage repair in the aetiology of BLBC and HGSOC.

Healthspan and healthy aging are areas of research with potential socioeconomic impact. Here, the authors present the Medical Genome Reference Bank (MGRB) which consist of over 4,000 individuals aged 70 years and older without a history of the major age-related diseases and report on results from whole-genome sequencing and association analyses.

Pediatric high grade gliomas (HGG) are primary brain malignancies that result in significant morbidity and mortality. One of the challenges in their treatment is inter- and intra-tumoral heterogeneity. Precision medicine approaches have the potential to enhance diagnostic, prognostic and/or therapeutic information. In this case study we describe the molecular characterization of a pediatric HGG and the use of an integrated approach based on genomic, in vitro and in vivo testing to identify actionable targets and treatment options. Molecular analysis based on WGS performed on initial and recurrent tumor biopsies revealed mutations in TP53, TSC1 and CIC genes, focal amplification of MYCN, and copy number gains in SMO and c-MET. Transcriptomic analysis identified increased expression of MYCN, and genes involved in sonic hedgehog signaling proteins (SHH, SMO, GLI1, GLI2) and receptor tyrosine kinase pathways (PLK, AURKA, c-MET). HTS revealed no cytotoxic efficacy of SHH pathway inhibitors while sensitivity was observed to the mTOR inhibitor temsirolimus, the ALK inhibitor ceritinib, and the PLK1 inhibitor BI2536. Based on the integrated approach, temsirolimus, ceritinib, BI2536 and standard therapy temozolomide were selected for further in vivo evaluation. Using the PDX animal model (median survival 28 days) we showed significant in vivo activity for mTOR inhibition by temsirolimus and BI2536 (median survival 109 and 115.5 days respectively) while ceritinib and temozolomide had only a moderate effect (43 and 75.5 days median survival respectively). This case study demonstrates that an integrated approach based on genomic, in vitro and in vivo drug efficacy testing in a PDX model may be useful to guide the management of high risk pediatric brain tumor in a clinically meaningful timeframe.

Zellweger syndrome spectrum disorders are caused by mutations in any of at least 12 different PEX genes. This includes PEX16, an important regulator of peroxisome biogenesis. Using whole genome sequencing, we detected previously unreported, biallelic variants in PEX16 [NM_004813.2:c.658G>A, p.(Ala220Thr) and NM_004813.2:c.830G>A, p.(Arg277Gln)] in an individual with leukodystrophy, spastic paraplegia, cerebellar ataxia, and craniocervical dystonia with normal plasma very long chain fatty acids. Using olfactory-neurosphere derived cells, a population of neural stem cells, we showed patient cells had reduced peroxisome density and increased peroxisome size, replicating previously reported findings in PEX16 cell lines. Along with alterations in peroxisome morphology, patient cells also had impaired peroxisome function with reduced catalase activity. Furthermore, patient cells had reduced oxidative stress levels after exposure to hydrogen-peroxide (H2O2), which may be a result of compensation by H2O2 metabolising enzymes other than catalase to preserve peroxisome-related cell functions. Our findings of impaired catalase activity and altered oxidative stress response are novel. Our study expands the phenotype of PEX16 mutations by including dystonia and provides further insights into the pathological mechanisms underlying PEX16-associated disorders. Additional studies of the full spectrum of peroxisomal dysfunction could improve our understanding of the mechanism underlying PEX16-associated disorders. Keywords: Whole genome sequencing, PEX16, Peroxisomes, Leukodystrophy, Dystonia

Summary: Early infantile epileptic encephalopathies (EOEE) are a debilitating spectrum of disorders associated with cognitive impairments. We present a clinical report of a KCNT2 mutation in an EOEE patient. The de novo heterozygous variant Phe240Leu SLICK was identified by exome sequencing and confirmed by Sanger sequencing. Phe240Leu rSlick and hSLICK channels were electrophysiologically, heterologously characterized to reveal three significant alterations to channel function. First, [Cl−]i sensitivity was reversed in Phe240Leu channels. Second, predominantly K+-selective WT channels were made to favor Na+ over K+ by Phe240Leu. Third, and consequent to altered ion selectivity, Phe240Leu channels had larger inward conductance. Further, rSlick channels induced membrane hyperexcitability when expressed in primary neurons, resembling the cellular seizure phenotype. Taken together, our results confirm that Phe240Leu is a “change-of-function” KCNT2 mutation, demonstrating unusual altered selectivity in KNa channels. These findings establish pathogenicity of the Phe240Leu KCNT2 mutation in the reported EOEE patient. : Gururaj et al. report a KCNT2 mutation in a patient with epileptic encephalopathy and employ electrophysiological analyses to establish channel properties that could underlie epileptogenesis: namely, inhibition by high [Cl−]i and loss of exclusive selectivity to K+. Furthermore, primary neurons expressing Ph240Leu display a hyperexcitable phenotype. Keywords: seizures, epileptic encephalopathy, ion channels, potassium channels, Slick, KCNT2, selectivity, Xenopus oocytes, Slack, KCNT1

Abstract Cancer is a multifactorial disease driven by a combination of genetic and environmental factors. Many cancer driver mutations have been characterised in protein-coding regions of the genome. However, mutations in noncoding regions associated with cancer have been less investigated. G-quadruplex (G4) nucleic acids are four-stranded secondary structures formed in guanine-rich sequences and prevalent in the regulatory regions. In this study, we used published whole cancer genome sequence data to find mutations in cancer patients that overlap potential RNA G4-forming sequences in 5′ UTRs. Using RNAfold, we assessed the effect of these mutations on the thermodynamic stability of predicted RNA G4s in the context of full-length 5′ UTRs. Of the 217 identified mutations, we found that 33 are predicted to destabilise and 21 predicted to stabilise potential RNA G4s. We experimentally validated the effect of destabilising mutations in the 5′ UTRs of BCL2 and CXCL14 and one stabilising mutation in the 5′ UTR of TAOK2. These mutations resulted in an increase or a decrease in translation of these mRNAs, respectively. These findings suggest that mutations that modulate the G4 stability in the noncoding regions could act as cancer driver mutations, which present an opportunity for early cancer diagnosis using individual sequencing information.

Chronic myelomonocytic leukaemia is treated with agents that modify DNA methylation but whether they have direct cytotoxic effects is unclear. Here, the authors show that cells from treated patients show marked methylation changes without altered somatic mutation burden, suggesting that cytotoxicity is not a major factor in therapeutic efficacy.

We examined whether hypoxic exposure prior to the event of transplantation would have a positive or negative effect upon later islet graft function. Mouse islets exposed to hypoxic culture were transplanted into syngeneic recipients. Islet graft function, β-cell physiology, as well as molecular changes were examined. Expression of hypoxia-response genes in human islets pre- and posttransplant was examined by microarray. Hypoxia-preexposed murine islet grafts provided poor glycemic control in their syngeneic recipients, marked by persistent hyperglycemia and pronounced glucose intolerance with failed first- and second-phase glucose-stimulated insulin secretion in vivo. Mechanistically, hypoxic preexposure stabilized HIF-1α with a concomitant increase in hypoxic-response genes including LDHA, and a molecular gene set, which would favor glycolysis and lactate production and impair glucose sensing. Indeed, static incubation studies showed that hypoxia-exposed islets exhibited dysregulated glucose responsiveness with elevated basal insulin secretion. Isolated human islets, prior to transplantation, express a characteristic hypoxia-response gene expression signature, including high levels of LDHA, which is maintained posttransplant. Hypoxic preexposure of an islet graft drives a HIF-dependent switch to glycolysis with subsequent poor glycemic control and loss of glucose-stimulated insulin secretion (GSIS). Early intervention to reverse or prevent these hypoxia-induced metabolic gene changes may improve clinical islet transplantation.

In the context of islet transplantation, experimental models show that induction of islet intrinsic NF-κB-dependent proinflammatory genes can contribute to islet graft rejection. Isolation of human islets triggers activation of the NF-κB and mitogen-activated kinase (MAPK) stress response pathways. However, the downstream NF-κB target genes induced in human islets during the isolation process are poorly described. Therefore, in this study, using microarray, bioinformatic, and RTqPCR approaches, we determined the pattern of genes expressed by a set of 14 human islet preparations. We found that isolated human islets express a panel of genes reminiscent of cells undergoing a marked NF-κB-dependent proinflammatory response. Expressed genes included matrix metallopeptidase 1 (MMP1) and fibronectin 1 (FN1), factors involved in tissue remodeling, adhesion, and cell migration; inflammatory cytokines IL-1β and IL-8; genes regulating cell survival including A20 and ATF3; and notably high expression of a set of chemokines that would favor neutrophil and monocyte recruitment including CXCL2, CCL2, CXCL12, CXCL1, CXCL6, and CCL28. Of note, the inflammatory profile of isolated human islets was maintained after transplantation into RAG -/- recipients. Thus, human islets can provide a reservoir of NF-κB-dependent inflammatory factors that have the potential to contribute to the anti-islet-graft immune response. To test this hypothesis, we extracted rodent islets under optimal conditions, forced activation of NF-κB, and transplanted them into allogenic recipients. These NF-κB activated islets not only expressed the same chemokine profile observed in human islets but also struggled to maintain normoglycemia posttransplantation. Further, NF-κB-activated islets were rejected with a faster tempo as compared to non-NF-κB-activated rodent islets. Thus, isolated human islets can make cell autonomous contributions to the ensuing allograft response by elaborating inflammatory factors that contribute to their own demise. These data highlight the potential importance of islet intrinsic proinflammatory responses as targets for therapeutic intervention.

Human islets are subjected to a number of stresses before and during their isolation that may influence their survival and engraftment after transplantation. Apoptosis is likely to be activated in response to these stresses. Apoptosis due to intrinsic stresses is regulated by pro- and antiapoptotic members of the Bcl-2 family. While the role of the Bcl-2 family in apoptosis of rodent islets is becoming increasingly understood, little is known about which of these molecules are expressed or required for apoptosis of human islets. This study investigated the expression of the Bcl-2 family of molecules in isolated human islets. RNA and protein lysates were extracted from human islets immediately postisolation. At the same time, standard quality control assays including viability staining and β-cell content were performed on each islet preparation. Microarrays, RT-PCR, and Western blotting were performed on islet RNA and protein. The prosurvival molecules Bcl-xl and Mcl-1, but not Bcl-2, were highly expressed. The multidomain proapoptotic effector molecule Bax was expressed at higher levels than Bak. Proapoptotic BH3-only molecules were expressed at low levels, with Bid being the most abundant. The proapoptotic molecules BNIP3, BNIP3L, and Beclin-1 were all highly expressed, indicating exposure of islets to oxygen and nutrient deprivation during isolation. Our data provide a comprehensive analysis of expression levels of pro- and antiapoptotic Bcl-2 family members in isolated human islets. Knowledge of which molecules are expressed will guide future research to understand the apoptotic pathways activated during isolation or after transplantation. This is crucial for the design of methods to achieve improved transplantation outcomes.