Dr.Nod: computational framework for discovery of regulatory non-coding drivers in tissue-matched distal regulatory elements.

The discovery of cancer driver mutations is a fundamental goal in cancer research. While many cancer driver mutations have been discovered in the protein-coding genome, research into potential cancer drivers in the non-coding regions showed limited success so far. Here, we present a novel comprehensive framework Dr.Nod for detection of non-coding cis-regulatory candidate driver mutations that are associated with dysregulated gene expression using tissue-matched enhancer-gene annotations. Applying the framework to data from over 1500 tumours across eight tissues revealed a 4.4-fold enrichment of candidate driver mutations in regulatory regions of known cancer driver genes. An overarching conclusion that emerges is that the non-coding driver mutations contribute to cancer by significantly altering transcription factor binding sites, leading to upregulation of tissue-matched oncogenes and down-regulation of tumour-suppressor genes. Interestingly, more than half of the detected cancer-promoting non-coding regulatory driver mutations are over 20 kb distant from the cancer-associated genes they regulate. Our results show the importance of tissue-matched enhancer-gene maps, functional impact of mutations, and complex background mutagenesis model for the prediction of non-coding regulatory drivers. In conclusion, our study demonstrates that non-coding mutations in enhancers play a previously underappreciated role in cancer and dysregulation of clinically relevant target genes.

Sensitive tumour detection and classification using plasma cell-free DNA methylomes.

The use of liquid biopsies for cancer detection and management is rapidly gaining prominence1. Current methods for the detection of circulating tumour DNA involve sequencing somatic mutations using cell-free DNA, but the sensitivity of these methods may be low among patients with early-stage cancer given the limited number of recurrent mutations2-5. By contrast, large-scale epigenetic alterations-which are tissue- and cancer-type specific-are not similarly constrained6 and therefore potentially have greater ability to detect and classify cancers in patients with early-stage disease. Here we develop a sensitive, immunoprecipitation-based protocol to analyse the methylome of small quantities of circulating cell-free DNA, and demonstrate the ability to detect large-scale DNA methylation changes that are enriched for tumour-specific patterns. We also demonstrate robust performance in cancer detection and classification across an extensive collection of plasma samples from several tumour types. This work sets the stage to establish biomarkers for the minimally invasive detection, interception and classification of early-stage cancers based on plasma cell-free DNA methylation patterns.

Prediction of acute myeloid leukaemia risk in healthy individuals.

The incidence of acute myeloid leukaemia (AML) increases with age and mortality exceeds 90% when diagnosed after age 65. Most cases arise without any detectable early symptoms and patients usually present with the acute complications of bone marrow failure1. The onset of such de novo AML cases is typically preceded by the accumulation of somatic mutations in preleukaemic haematopoietic stem and progenitor cells (HSPCs) that undergo clonal expansion2,3. However, recurrent AML mutations also accumulate in HSPCs during ageing of healthy individuals who do not develop AML, a phenomenon referred to as age-related clonal haematopoiesis (ARCH)4-8. Here we use deep sequencing to analyse genes that are recurrently mutated in AML to distinguish between individuals who have a high risk of developing AML and those with benign ARCH. We analysed peripheral blood cells from 95 individuals that were obtained on average 6.3 years before AML diagnosis (pre-AML group), together with 414 unselected age- and gender-matched individuals (control group). Pre-AML cases were distinct from controls and had more mutations per sample, higher variant allele frequencies, indicating greater clonal expansion, and showed enrichment of mutations in specific genes. Genetic parameters were used to derive a model that accurately predicted AML-free survival; this model was validated in an independent cohort of 29 pre-AML cases and 262 controls. Because AML is rare, we also developed an AML predictive model using a large electronic health record database that identified individuals at greater risk. Collectively our findings provide proof-of-concept that it is possible to discriminate ARCH from pre-AML many years before malignant transformation. This could in future enable earlier detection and monitoring, and may help to inform intervention.

Tracing the origins of relapse in acute myeloid leukaemia to stem cells.

In acute myeloid leukaemia, long-term survival is poor as most patients relapse despite achieving remission. Historically, the failure of therapy has been thought to be due to mutations that produce drug resistance, possibly arising as a consequence of the mutagenic properties of chemotherapy drugs. However, other lines of evidence have pointed to the pre-existence of drug-resistant cells. For example, deep sequencing of paired diagnosis and relapse acute myeloid leukaemia samples has provided direct evidence that relapse in some cases is generated from minor genetic subclones present at diagnosis that survive chemotherapy, suggesting that resistant cells are generated by evolutionary processes before treatment and are selected by therapy. Nevertheless, the mechanisms of therapy failure and capacity for leukaemic regeneration remain obscure, as sequence analysis alone does not provide insight into the cell types that are fated to drive relapse. Although leukaemia stem cells have been linked to relapse owing to their dormancy and self-renewal properties, and leukaemia stem cell gene expression signatures are highly predictive of therapy failure, experimental studies have been primarily correlative and a role for leukaemia stem cells in acute myeloid leukaemia relapse has not been directly proved. Here, through combined genetic and functional analysis of purified subpopulations and xenografts from paired diagnosis/relapse samples, we identify therapy-resistant cells already present at diagnosis and two major patterns of relapse. In some cases, relapse originated from rare leukaemia stem cells with a haematopoietic stem/progenitor cell phenotype, while in other instances relapse developed from larger subclones of immunophenotypically committed leukaemia cells that retained strong stemness transcriptional signatures. The identification of distinct patterns of relapse should lead to improved methods for disease management and monitoring in acute myeloid leukaemia. Moreover, the shared functional and transcriptional stemness properties that underlie both cellular origins of relapse emphasize the importance of developing new therapeutic approaches that target stemness to prevent relapse.

SCALP syndrome with a germline heterozygous DOCK6 mutation and somatic mosaic NRAS Q61R mutation.

We present a case of SCALP syndrome, which was diagnosed in a male infant with the characteristic findings of sebaceous nevi, central nervous system malformations, aplasia cutis congenita, limbal dermoid, and giant congenital melanocytic nevi, or pigmented nevi. We identified a germline compound heterozygous DOCK6 mutation and a somatic mosaic NRAS Q61R mutation in the giant congenital melanocytic nevus. This report will increase clinician awareness of SCALP syndrome and augment the literature in characterizing this rare syndrome, including its genetic background.

Environmental pesticide exposure and non-Hodgkin lymphoma survival: a population-based study.

BACKGROUND: There is evidence indicating that pesticide exposure is a risk factor for non-Hodgkin lymphoma (NHL) development. However, the association between pesticide exposure and NHL survival is not well-established. METHODS: Using the California Cancer Registry, we identified patients with a first primary diagnosis of NHL from 2010 to 2016 and linked these patients with CalEnviroScreen 3.0 to obtain production agriculture pesticide exposure to 70 chemicals from the state-mandated Pesticide Use Reporting (PUR) by census tract from 2012 to 2014. In addition, data from PUR was integrated into a geographic information system that employs land-use data to estimate cumulative exposure to specific pesticides previously associated with NHL (glyphosate, organophosphorus, carbamate, phenoxyherbicide, and 2,4-dimethylamine salt) between 10 years prior up to 1 year after NHL diagnosis. Multivariable Cox proportional hazards regression models were used to evaluate the association between total pesticide exposure from CalEnviroScreen 3.0 and individual pesticide exposure from geographic land use data and lymphoma-specific and overall survival. RESULTS: Among 35,808 NHL patients identified, 44.2% were exposed to pesticide in their census tract of residence. Glyphosate, organophosphorus, carbamate, phenoxyherbicide, and 2,4-dimethylamine salt exposure was observed in 34.1%, 26.0%, 10.6%, 14.0%, and 12.8% of NHL patients, respectively. Total pesticide exposure at the time of diagnosis was not associated with lymphoma-specific or overall survival. In addition, no association was consistently found between glyphosate, organophosphorus, carbamate, phenoxyherbicide, and 2,4 dimethylamine salt exposure and lymphoma-specific or overall survival. CONCLUSIONS: Although we found no consistent associations between agricultural pesticide exposure at the neighborhood level and worse survival, these results provide a platform for designing future studies to determine the association between pesticide and NHL.

Coming of age: ten years of next-generation sequencing technologies.

Since the completion of the human genome project in 2003, extraordinary progress has been made in genome sequencing technologies, which has led to a decreased cost per megabase and an increase in the number and diversity of sequenced genomes. An astonishing complexity of genome architecture has been revealed, bringing these sequencing technologies to even greater advancements. Some approaches maximize the number of bases sequenced in the least amount of time, generating a wealth of data that can be used to understand increasingly complex phenotypes. Alternatively, other approaches now aim to sequence longer contiguous pieces of DNA, which are essential for resolving structurally complex regions. These and other strategies are providing researchers and clinicians a variety of tools to probe genomes in greater depth, leading to an enhanced understanding of how genome sequence variants underlie phenotype and disease.

A renewed model of pancreatic cancer evolution based on genomic rearrangement patterns.

Pancreatic cancer, a highly aggressive tumour type with uniformly poor prognosis, exemplifies the classically held view of stepwise cancer development. The current model of tumorigenesis, based on analyses of precursor lesions, termed pancreatic intraepithelial neoplasm (PanINs) lesions, makes two predictions: first, that pancreatic cancer develops through a particular sequence of genetic alterations (KRAS, followed by CDKN2A, then TP53 and SMAD4); and second, that the evolutionary trajectory of pancreatic cancer progression is gradual because each alteration is acquired independently. A shortcoming of this model is that clonally expanded precursor lesions do not always belong to the tumour lineage, indicating that the evolutionary trajectory of the tumour lineage and precursor lesions can be divergent. This prevailing model of tumorigenesis has contributed to the clinical notion that pancreatic cancer evolves slowly and presents at a late stage. However, the propensity for this disease to rapidly metastasize and the inability to improve patient outcomes, despite efforts aimed at early detection, suggest that pancreatic cancer progression is not gradual. Here, using newly developed informatics tools, we tracked changes in DNA copy number and their associated rearrangements in tumour-enriched genomes and found that pancreatic cancer tumorigenesis is neither gradual nor follows the accepted mutation order. Two-thirds of tumours harbour complex rearrangement patterns associated with mitotic errors, consistent with punctuated equilibrium as the principal evolutionary trajectory. In a subset of cases, the consequence of such errors is the simultaneous, rather than sequential, knockout of canonical preneoplastic genetic drivers that are likely to set-off invasive cancer growth. These findings challenge the current progression model of pancreatic cancer and provide insights into the mutational processes that give rise to these aggressive tumours.

Complex rearrangements and oncogene amplifications revealed by long-read DNA and RNA sequencing of a breast cancer cell line.

The SK-BR-3 cell line is one of the most important models for HER2+ breast cancers, which affect one in five breast cancer patients. SK-BR-3 is known to be highly rearranged, although much of the variation is in complex and repetitive regions that may be underreported. Addressing this, we sequenced SK-BR-3 using long-read single molecule sequencing from Pacific Biosciences and develop one of the most detailed maps of structural variations (SVs) in a cancer genome available, with nearly 20,000 variants present, most of which were missed by short-read sequencing. Surrounding the important ERBB2 oncogene (also known as HER2), we discover a complex sequence of nested duplications and translocations, suggesting a punctuated progression. Full-length transcriptome sequencing further revealed several novel gene fusions within the nested genomic variants. Combining long-read genome and transcriptome sequencing enables an in-depth analysis of how SVs disrupt the genome and sheds new light on the complex mechanisms involved in cancer genome evolution.

Activation of hedgehog signaling associates with early disease progression in chronic lymphocytic leukemia.

Targeted sequencing of 103 leukemia-associated genes in leukemia cells from 841 treatment-naive patients with chronic lymphocytic leukemia (CLL) identified 89 (11%) patients as having CLL cells with mutations in genes encoding proteins that putatively are involved in hedgehog (Hh) signaling. Consistent with this finding, there was a significant association between the presence of these mutations and the expression of GLI1 (χ2 test, P < .0001), reflecting activation of the Hh pathway. However, we discovered that 38% of cases without identified mutations also were GLI1+ Patients with GLI1+ CLL cells had a shorter median treatment-free survival than patients with CLL cells lacking expression of GLI1 independent of IGHV mutation status. We found that GANT61, a small molecule that can inhibit GLI1, was highly cytotoxic for GLI1+ CLL cells relative to that of CLL cells without GLI1. Collectively, this study shows that a large proportion of patients have CLL cells with activated Hh signaling, which is associated with early disease progression and enhanced sensitivity to inhibition of GLI1.

Germline and somatic DNA repair gene alterations in prostate cancer.

BACKGROUND: Emerging evidence has suggested that DNA repair gene alterations may be important in prostate cancer pathogenesis. In the current study, the authors sought to characterize alterations in DNA repair pathway genes in both primary and metastatic prostate tumors with attention to tissue distribution as well as specific genomic alterations. METHODS: The authors studied the distribution and type of alterations in 24 genes that are considered important for DNA repair in 944 prostate cancers harvested from localized and metastatic tumors. Tumor DNA underwent hybrid capture for all coding exons of 287 or 395 cancer-related genes plus select introns from 19 or 31 genes frequently rearranged in cancer. Captured libraries were sequenced to a median exon coverage depth of >×500. Specific genomic alterations were characterized and the frequencies of mutations by tissue site (prostate vs metastases) were compared using logistic regression. RESULTS: A total of 152 patients from the cohort of 944 men (16%) harbored a germline or somatic mutation in ≥1 DNA repair genes. The most frequently mutated genes were BRCA2 (11.4%) and ATM (5.8%), followed by MSH6 (2.5%) and MSH2 (2.1%). Mutations were identified in approximately 20.1% of primary prostate tumors compared with 18.8% of bone metastases. When stratified by tissue site, the highest rates of DNA repair mutations were found in solid organ metastases, including brain and visceral metastases, compared with prostate. CONCLUSIONS: DNA repair gene mutations are more common in metastatic than localized prostate tumors. Visceral and other solid organ metastases appear enriched for these mutations compared with localized tumors or bone and lymph node metastases.

Identification of pre-leukaemic haematopoietic stem cells in acute leukaemia.

In acute myeloid leukaemia (AML), the cell of origin, nature and biological consequences of initiating lesions, and order of subsequent mutations remain poorly understood, as AML is typically diagnosed without observation of a pre-leukaemic phase. Here, highly purified haematopoietic stem cells (HSCs), progenitor and mature cell fractions from the blood of AML patients were found to contain recurrent DNMT3A mutations (DNMT3A(mut)) at high allele frequency, but without coincident NPM1 mutations (NPM1c) present in AML blasts. DNMT3A(mut)-bearing HSCs showed a multilineage repopulation advantage over non-mutated HSCs in xenografts, establishing their identity as pre-leukaemic HSCs. Pre-leukaemic HSCs were found in remission samples, indicating that they survive chemotherapy. Therefore DNMT3A(mut) arises early in AML evolution, probably in HSCs, leading to a clonally expanded pool of pre-leukaemic HSCs from which AML evolves. Our findings provide a paradigm for the detection and treatment of pre-leukaemic clones before the acquisition of additional genetic lesions engenders greater therapeutic resistance.

Comparative toxicoproteogenomics of mouse and rat liver identifies TCDD-resistance genes.

The aryl hydrocarbon receptor (AHR) mediates many toxic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). However, the AHR alone does not explain the widely different outcomes among organisms. To identify the other factors involved, we evaluated three transgenic mouse lines, each expressing a different rat AHR isoform (rWT, DEL, and INS) providing widely different resistance to TCDD toxicity, as well as C57BL/6 and DBA/2 mice which exhibit a ~ tenfold divergence in TCDD sensitivity (exposures of 5-1000 µg/kg TCDD). We supplement these with whole-genome sequencing, together with transcriptomic and proteomic analyses of the corresponding rat models, Long-Evans (L-E) and Han/Wistar (H/W) rats (having a ~ 1000-fold difference in their TCDD sensitivities; 100 µg/kg TCDD), to identify genes associated with TCDD-response phenotypes. Overall, we identified up to 50% of genes with altered mRNA abundance following TCDD exposure are associated with a single AHR isoform (33.8%, 11.7%, 5.2% and 0.3% of 3076 genes altered unique to rWT, DEL, C57BL/6 and INS respectively following 1000 µg/kg TCDD). Hepatic Pxdc1 was significantly repressed in all three TCDD-sensitive animal models (C57BL/6 and rWT mice, and L-E rat) after TCDD exposure. Three genes, including Cxxc5, Sugp1 and Hgfac, demonstrated different AHRE-1 (full) motif occurrences within their promoter regions between rat strains, as well as different patterns of mRNA abundance. Several hepatic proteins showed parallel up- or downward alterations with their RNAs, with three genes (SNRK, IGTP and IMPA2) showing consistent, strain-dependent changes. These data show the value of integrating genomic, transcriptomic and proteomic evidence across multi-species models in toxicologic studies.

Improving classification of melanocytic nevi: Association of BRAF V600E expression with distinct histomorphologic features.

BACKGROUND: A subset of melanomas carrying a B-Raf proto-oncogene, serine/threonine kinase gene (BRAF) V600E mutation, which is the most common targetable mutation in melanoma, arise in association with a melanocytic nevus that is also harboring a BRAF V600E mutation. The detailed histomorphologic characteristics of nevi positive for BRAF V600E have not been systematically documented. OBJECTIVE: To identify histomorphologic features correlating with BRAF V600E status in nevi. METHODS: We retrospectively identified melanocytic nevi from our laboratory reporting system. We performed a histomorphologic analysis and analysis of BRAF V600E expression by immunohistochemistry. RESULTS: Thirteen nevi (14.8%) were negative and 76 (86.4%) were positive for BRAF V600E. The nevi positive for BRAF V600E were predominantly dermal (predominantly dermal growth in 55.3% of nevi positive for BRAF V600E and 15.4% of nevi negative for BRAF V600E [P = .01]) and showed a congenital growth pattern (congenital growth pattern in 51.3% of nevi positive for BRAF V600E and 15.4% of nevi negative for BRAF V600E [P = .02]). Compared with nevi negative for BRAF V600E, those that were positive for BRAF V600E often exhibited predominantly nested intraepidermal melanocytes, larger junctional nests, abrupt lateral circumscription, and larger cell size. Architectural disorder and inflammatory infiltrates were seen more often in nevi negative for BRAF V600E. BRAF sequencing of a subset of nevi confirmed the immunohistochemical results. LIMITATIONS: Limitations include the study's retrospective design and the small sample size of nevi negative for BRAF V600E. CONCLUSIONS: BRAF V600E is associated with distinct histomorphologic features in nevi. These features may contribute to improving the accuracy of classification and diagnosis of melanocytic neoplasms.

Molecular heterogeneity of non-small cell lung carcinoma patient-derived xenografts closely reflect their primary tumors.

Availability of lung cancer models that closely mimic human tumors remains a significant gap in cancer research, as tumor cell lines and mouse models may not recapitulate the spectrum of lung cancer heterogeneity seen in patients. We aimed to establish a patient-derived tumor xenograft (PDX) resource from surgically resected non-small cell lung cancer (NSCLC). Fresh tumor tissue from surgical resection was implanted and grown in the subcutaneous pocket of non-obese severe combined immune deficient (NOD SCID) gamma mice. Subsequent passages were in NOD SCID mice. A subset of matched patient and PDX tumors and non-neoplastic lung tissues were profiled by whole exome sequencing, single nucleotide polymorphism (SNP) and methylation arrays, and phosphotyrosine (pY)-proteome by mass spectrometry. The data were compared to published NSCLC datasets of NSCLC primary and cell lines. 127 stable PDXs were established from 441 lung carcinomas representing all major histological subtypes: 52 adenocarcinomas, 62 squamous cell carcinomas, one adeno-squamous carcinoma, five sarcomatoid carcinomas, five large cell neuroendocrine carcinomas, and two small cell lung cancers. Somatic mutations, gene copy number and expression profiles, and pY-proteome landscape of 36 PDXs showed greater similarity with patient tumors than with established cell lines. Novel somatic mutations on cancer associated genes were identified but only in PDXs, likely due to selective clonal growth in the PDXs that allows detection of these low allelic frequency mutations. The results provide the strongest evidence yet that PDXs established from lung cancers closely mimic the characteristics of patient primary tumors.

SeqControl: process control for DNA sequencing.

As high-throughput sequencing continues to increase in speed and throughput, routine clinical and industrial application draws closer. These 'production' settings will require enhanced quality monitoring and quality control to optimize output and reduce costs. We developed SeqControl, a framework for predicting sequencing quality and coverage using a set of 15 metrics describing overall coverage, coverage distribution, basewise coverage and basewise quality. Using whole-genome sequences of 27 prostate cancers and 26 normal references, we derived multivariate models that predict sequencing quality and depth. SeqControl robustly predicted how much sequencing was required to reach a given coverage depth (area under the curve (AUC) = 0.993), accurately classified clinically relevant formalin-fixed, paraffin-embedded samples, and made predictions from as little as one-eighth of a sequencing lane (AUC = 0.967). These techniques can be immediately incorporated into existing sequencing pipelines to monitor data quality in real time. SeqControl is available at http://labs.oicr.on.ca/Boutros-lab/software/SeqControl/.

Optimization of miRNA-seq data preprocessing.

The past two decades of microRNA (miRNA) research has solidified the role of these small non-coding RNAs as key regulators of many biological processes and promising biomarkers for disease. The concurrent development in high-throughput profiling technology has further advanced our understanding of the impact of their dysregulation on a global scale. Currently, next-generation sequencing is the platform of choice for the discovery and quantification of miRNAs. Despite this, there is no clear consensus on how the data should be preprocessed before conducting downstream analyses. Often overlooked, data preprocessing is an essential step in data analysis: the presence of unreliable features and noise can affect the conclusions drawn from downstream analyses. Using a spike-in dilution study, we evaluated the effects of several general-purpose aligners (BWA, Bowtie, Bowtie 2 and Novoalign), and normalization methods (counts-per-million, total count scaling, upper quartile scaling, Trimmed Mean of M, DESeq, linear regression, cyclic loess and quantile) with respect to the final miRNA count data distribution, variance, bias and accuracy of differential expression analysis. We make practical recommendations on the optimal preprocessing methods for the extraction and interpretation of miRNA count data from small RNA-sequencing experiments.

Next-generation sequencing identifies rare variants associated with Noonan syndrome.

Noonan syndrome (NS) is a relatively common genetic disorder, characterized by typical facies, short stature, developmental delay, and cardiac abnormalities. Known causative genes account for 70-80% of clinically diagnosed NS patients, but the genetic basis for the remaining 20-30% of cases is unknown. We performed next-generation sequencing on germ-line DNA from 27 NS patients lacking a mutation in the known NS genes. We identified gain-of-function alleles in Ras-like without CAAX 1 (RIT1) and mitogen-activated protein kinase kinase 1 (MAP2K1) and previously unseen loss-of-function variants in RAS p21 protein activator 2 (RASA2) that are likely to cause NS in these patients. Expression of the mutant RASA2, MAP2K1, or RIT1 alleles in heterologous cells increased RAS-ERK pathway activation, supporting a causative role in NS pathogenesis. Two patients had more than one disease-associated variant. Moreover, the diagnosis of an individual initially thought to have NS was revised to neurofibromatosis type 1 based on an NF1 nonsense mutation detected in this patient. Another patient harbored a missense mutation in NF1 that resulted in decreased protein stability and impaired ability to suppress RAS-ERK activation; however, this patient continues to exhibit a NS-like phenotype. In addition, a nonsense mutation in RPS6KA3 was found in one patient initially diagnosed with NS whose diagnosis was later revised to Coffin-Lowry syndrome. Finally, we identified other potential candidates for new NS genes, as well as potential carrier alleles for unrelated syndromes. Taken together, our data suggest that next-generation sequencing can provide a useful adjunct to RASopathy diagnosis and emphasize that the standard clinical categories for RASopathies might not be adequate to describe all patients.