Heat selection enables highly scalable methylome profiling in cell-free DNA for noninvasive monitoring of cancer patients.

Genome-wide analysis of cell-free DNA methylation profile is a promising approach for sensitive and specific detection of many cancers. However, scaling such assays for clinical translation is impractical because of the high cost of whole-genome bisulfite sequencing. We show that the small fraction of GC-rich genome is highly enriched in CpG sites and disproportionately harbors most of the cancer-specific methylation signature. Here, we report on the simple and effective heat enrichment of CpG-rich regions for bisulfite sequencing (Heatrich-BS) platform that allows for focused methylation profiling in these highly informative regions. Our novel method and bioinformatics algorithm enable accurate tumor burden estimation and quantitative tracking of colorectal cancer patient's response to treatment at much reduced sequencing cost suitable for frequent monitoring. We also show tumor epigenetic subtyping using Heatrich-BS, which could enable patient stratification. Heatrich-BS holds great potential for highly scalable screening and monitoring of cancer using liquid biopsy.

Identification of differential RNA modifications from nanopore direct RNA sequencing with xPore.

RNA modifications, such as N6-methyladenosine (m6A), modulate functions of cellular RNA species. However, quantifying differences in RNA modifications has been challenging. Here we develop a computational method, xPore, to identify differential RNA modifications from nanopore direct RNA sequencing (RNA-seq) data. We evaluate our method on transcriptome-wide m6A profiling data, demonstrating that xPore identifies positions of m6A sites at single-base resolution, estimates the fraction of modified RNA species in the cell and quantifies the differential modification rate across conditions. We apply xPore to direct RNA-seq data from six cell lines and multiple myeloma patient samples without a matched control sample and find that many m6A sites are preserved across cell types, whereas a subset exhibit significant differences in their modification rates. Our results show that RNA modifications can be identified from direct RNA-seq data with high accuracy, enabling analysis of differential modifications and expression from a single high-throughput experiment.

Spatial profiling of gastric cancer patient-matched primary and locoregional metastases reveals principles of tumour dissemination.

OBJECTIVE: Endoscopic mucosal biopsies of primary gastric cancers (GCs) are used to guide diagnosis, biomarker testing and treatment. Spatial intratumoural heterogeneity (ITH) may influence biopsy-derived information. We aimed to study ITH of primary GCs and matched lymph node metastasis (LNmet). DESIGN: GC resection samples were annotated to identify primary tumour superficial (PTsup), primary tumour deep (PTdeep) and LNmet subregions. For each subregion, we determined (1) transcriptomic profiles (NanoString 'PanCancer Progression Panel', 770 genes); (2) next-generation sequencing (NGS, 225 gastrointestinal cancer-related genes); (3) DNA copy number profiles by multiplex ligation-dependent probe amplification (MLPA, 16 genes); and (4) histomorphological phenotypes. RESULTS: NanoString profiling of 64 GCs revealed no differences between PTsup1 and PTsup2, while 43% of genes were differentially expressed between PTsup versus PTdeep and 38% in PTsup versus LNmet. Only 16% of genes were differently expressed between PTdeep and LNmet. Several genes with therapeutic potential (eg IGF1, PIK3CD and TGFB1) were overexpressed in LNmet and PTdeep compared with PTsup. NGS data revealed orthogonal support of NanoString results with 40% mutations present in PTdeep and/or LNmet, but not in PTsup. Conversely, only 6% of mutations were present in PTsup and were absent in PTdeep and LNmet. MLPA demonstrated significant ITH between subregions and progressive genomic changes from PTsup to PTdeep/LNmet. CONCLUSION: In GC, regional lymph node metastases are likely to originate from deeper subregions of the primary tumour. Future clinical trials of novel targeted therapies must consider assessment of deeper subregions of the primary tumour and/or metastases as several therapeutically relevant genes are only mutated, overexpressed or amplified in these regions.

Singapore Undiagnosed Disease Program: Genomic Analysis aids Diagnosis and Clinical Management.

OBJECTIVE: Use next-generation sequencing (NGS) technology to improve our diagnostic yield in patients with suspected genetic disorders in the Asian setting. DESIGN: A diagnostic study conducted between 2014 and 2019 (and ongoing) under the Singapore Undiagnosed Disease Program. Date of last analysis was 1 July 2019. SETTING: Inpatient and outpatient genetics service at two large academic centres in Singapore. PATIENTS: Inclusion criteria: patients suspected of genetic disorders, based on abnormal antenatal ultrasound, multiple congenital anomalies and developmental delay. EXCLUSION CRITERIA: patients with known genetic disorders, either after clinical assessment or investigations (such as karyotype or chromosomal microarray). INTERVENTIONS: Use of NGS technology-whole exome sequencing (WES) or whole genome sequencing (WGS). MAIN OUTCOME MEASURES: (1) Diagnostic yield by sequencing type, (2) diagnostic yield by phenotypical categories, (3) reduction in time to diagnosis and (4) change in clinical outcomes and management. RESULTS: We demonstrate a 37.8% diagnostic yield for WES (n=172) and a 33.3% yield for WGS (n=24). The yield was higher when sequencing was conducted on trios (40.2%), as well as for certain phenotypes (neuromuscular, 54%, and skeletal dysplasia, 50%). In addition to aiding genetic counselling in 100% of the families, a positive result led to a change in treatment in 27% of patients. CONCLUSION: Genomic sequencing is an effective method for diagnosing rare disease or previous 'undiagnosed' disease. The clinical utility of WES/WGS is seen in the shortened time to diagnosis and the discovery of novel variants. Additionally, reaching a diagnosis significantly impacts families and leads to alteration in management of these patients.

Detection of genomic alterations in breast cancer with circulating tumour DNA sequencing.

Analysis of circulating cell-free DNA (cfDNA) has opened new opportunities for characterizing tumour mutational landscapes with many applications in genomic-driven oncology. We developed a customized targeted cfDNA sequencing approach for breast cancer (BC) using unique molecular identifiers (UMIs) for error correction. Our assay, spanning a 284.5 kb target region, is combined with a novel freely-licensed bioinformatics pipeline that provides detection of low-frequency variants, and reliable identification of copy number variations (CNVs) directly from plasma DNA. We first evaluated our pipeline on reference samples. Then in a cohort of 35 BC patients our approach detected actionable driver and clonal variants at low variant frequency levels in cfDNA that were concordant (77%) with sequencing of primary and/or metastatic solid tumour sites. We also detected ERRB2 gene CNVs used for HER2 subtype classification with 80% precision compared to immunohistochemistry. Further, we evaluated fragmentation profiles of cfDNA in BC and observed distinct differences compared to data from healthy individuals. Our results show that the developed assay addresses the majority of tumour associated aberrations directly from plasma DNA, and thus may be used to elucidate genomic alterations in liquid biopsy studies.

Breast cancer in women with neurofibromatosis type 1 (NF1): a comprehensive case series with molecular insights into its aggressive phenotype.

PURPOSE: The purpose of the study was to improve the understanding of NF1-associated breast cancer, given the increased risk of breast cancer in this tumour predisposition syndrome and the limited data. METHODS: We identified 18 women with NF1 and breast cancer at our institution. Clinical and pathologic characteristics of NF1-associated breast cancers were compared with 7132 breast cancers in patients without NF1 from our institutional database. Next generation sequencing was performed on DNA from blood and breast cancer specimens available. Blood specimens negative for NF1 mutation were subjected to multiplex ligation-dependent probe amplification (MLPA) to identify complete/partial deletions or duplications. Expression of neurofibromin in the NF1-associated breast cancers was evaluated using immunohistochemistry. RESULTS: There was a higher frequency of grade 3 (83.3% vs 45.4%, p = 0.005), oestrogen receptor (ER) negative (66.7% vs 26.3%, p < 0.001) and human epidermal growth factor receptor 2 (HER2)-positive (66.7% vs 23.4%, p < 0.001) tumours among NF1 patients compared to non-NF1 breast cancers. Overall survival was inferior in NF1 patients in multivariable analysis (hazard ratio 2.25, 95% CI 1.11-4.60; p = 0.025). Apart from germline NF1 mutations (11/16; 69%), somatic mutations in TP53 (8/10; 80%), second-hit NF1 (2/10; 20%), KMT2C (4/10; 40%), KMT2D (2/10; 20%), and PIK3CA (2/10; 20%) were observed. Immunohistochemical expression of neurofibromin was seen in the nuclei and/or cytoplasm of all specimens, but without any consistent pattern in the intensity or extent. CONCLUSIONS: This comprehensive series of NF1-associated breast cancers suggests that their aggressive features are related to germline NF1 mutations in cooperation with somatic mutations in TP53, KMT2C and other genes.

Individualised multiplexed circulating tumour DNA assays for monitoring of tumour presence in patients after colorectal cancer surgery.

Circulating tumour DNA (ctDNA) has the potential to be a specific biomarker for the monitoring of tumours in patients with colorectal cancer (CRC). Here, our aim was to develop a personalised surveillance strategy to monitor the clinical course of CRC after surgery. We developed patient-specific ctDNA assays based on multiplexed detection of somatic mutations identified from patient primary tumours, and applied them to detect ctDNA in 44 CRC patients, analysing a total of 260 plasma samples. We found that ctDNA detection correlated with clinical events - it is detectable in pre-operative but not post-operative plasma, and also in patients with recurrent CRC. We also detected ctDNA in 11 out of 15 cases at or before clinical or radiological recurrence of CRC, indicating the potential of our assay for early detection of metastasis. We further present data from a patient with multiple primary cancers to demonstrate the specificity of our assays to distinguish between CRC recurrence and a second primary cancer. Our approach can complement current methods for surveillance of CRC by adding an individualised biological component, allowing us not only to point to the presence of residual or recurrent disease, but also attribute it to the original cancer.

A non-active-site SET domain surface crucial for the interaction of MLL1 and the RbBP5/Ash2L heterodimer within MLL family core complexes.

The mixed lineage leukemia-1 (MLL1) enzyme is a histone H3 lysine 4 (H3K4) monomethyltransferase and has served as a paradigm for understanding the mechanism of action of the human SET1 family of enzymes that include MLL1-MLL4 and SETd1a,b. Dimethylation of H3K4 requires a sub-complex including WRAD (WDR5, RbBP5, Ash2L, and DPY-30), which binds to each SET1 family member forming a minimal core complex that is required for multiple lysine methylation. We recently demonstrated that WRAD is a novel histone methyltransferase that preferentially catalyzes H3K4 dimethylation in a manner that is dependent on an unknown non-active-site surface from the MLL1 SET domain. Recent genome sequencing studies have identified a number of human disease-associated missense mutations that localize to the SET domains of several MLL family members. In this investigation, we mapped many of these mutations onto the three-dimensional structure of the SET domain and noticed that a subset of MLL2 (KMT2D, ALR, MLL4)-associated Kabuki syndrome missense mutations map to a common solvent-exposed surface that is not expected to alter enzymatic activity. We introduced these mutations into the MLL1 SET domain and observed that all are defective for H3K4 dimethylation by the MLL1 core complex, which is associated with a loss of the ability of MLL1 to interact with WRAD or with the RbBP5/Ash2L heterodimer. Our results suggest that amino acids from this surface, which we term the Kabuki interaction surface or KIS, are required for formation of a second active site within SET1 family core complexes.

Exome sequencing identifies a spectrum of mutation frequencies in advanced and lethal prostate cancers.

To catalog protein-altering mutations that may drive the development of prostate cancers and their progression to metastatic disease systematically, we performed whole-exome sequencing of 23 prostate cancers derived from 16 different lethal metastatic tumors and three high-grade primary carcinomas. All tumors were propagated in mice as xenografts, designated the LuCaP series, to model phenotypic variation, such as responses to cancer-directed therapeutics. Although corresponding normal tissue was not available for most tumors, we were able to take advantage of increasingly deep catalogs of human genetic variation to remove most germline variants. On average, each tumor genome contained ~200 novel nonsynonymous variants, of which the vast majority was specific to individual carcinomas. A subset of genes was recurrently altered across tumors derived from different individuals, including TP53, DLK2, GPC6, and SDF4. Unexpectedly, three prostate cancer genomes exhibited substantially higher mutation frequencies, with 2,000-4,000 novel coding variants per exome. A comparison of castration-resistant and castration-sensitive pairs of tumor lines derived from the same prostate cancer highlights mutations in the Wnt pathway as potentially contributing to the development of castration resistance. Collectively, our results indicate that point mutations arising in coding regions of advanced prostate cancers are common but, with notable exceptions, very few genes are mutated in a substantial fraction of tumors. We also report a previously undescribed subtype of prostate cancers exhibiting "hypermutated" genomes, with potential implications for resistance to cancer therapeutics. Our results also suggest that increasingly deep catalogs of human germline variation may challenge the necessity of sequencing matched tumor-normal pairs.

Spectrum of MLL2 (ALR) mutations in 110 cases of Kabuki syndrome.

Kabuki syndrome is a rare, multiple malformation disorder characterized by a distinctive facial appearance, cardiac anomalies, skeletal abnormalities, and mild to moderate intellectual disability. Simplex cases make up the vast majority of the reported cases with Kabuki syndrome, but parent-to-child transmission in more than a half-dozen instances indicates that it is an autosomal dominant disorder. We recently reported that Kabuki syndrome is caused by mutations in MLL2, a gene that encodes a Trithorax-group histone methyltransferase, a protein important in the epigenetic control of active chromatin states. Here, we report on the screening of 110 families with Kabuki syndrome. MLL2 mutations were found in 81/110 (74%) of families. In simplex cases for which DNA was available from both parents, 25 mutations were confirmed to be de novo, while a transmitted MLL2 mutation was found in two of three familial cases. The majority of variants found to cause Kabuki syndrome were novel nonsense or frameshift mutations that are predicted to result in haploinsufficiency. The clinical characteristics of MLL2 mutation-positive cases did not differ significantly from MLL2 mutation-negative cases with the exception that renal anomalies were more common in MLL2 mutation-positive cases. These results are important for understanding the phenotypic consequences of MLL2 mutations for individuals and their families as well as for providing a basis for the identification of additional genes for Kabuki syndrome.

Exome sequencing identifies MLL2 mutations as a cause of Kabuki syndrome.

We demonstrate the successful application of exome sequencing to discover a gene for an autosomal dominant disorder, Kabuki syndrome (OMIM%147920). We subjected the exomes of ten unrelated probands to massively parallel sequencing. After filtering against existing SNP databases, there was no compelling candidate gene containing previously unknown variants in all affected individuals. Less stringent filtering criteria allowed for the presence of modest genetic heterogeneity or missing data but also identified multiple candidate genes. However, genotypic and phenotypic stratification highlighted MLL2, which encodes a Trithorax-group histone methyltransferase: seven probands had newly identified nonsense or frameshift mutations in this gene. Follow-up Sanger sequencing detected MLL2 mutations in two of the three remaining individuals with Kabuki syndrome (cases) and in 26 of 43 additional cases. In families where parental DNA was available, the mutation was confirmed to be de novo (n = 12) or transmitted (n = 2) in concordance with phenotype. Our results strongly suggest that mutations in MLL2 are a major cause of Kabuki syndrome.