Kathleen E. Houlahan

Combining tumor genome simulation with crowdsourcing to benchmark somatic single-nucleotide-variant detection

The detection of somatic mutations from cancer genome sequences is key to understanding the genetic basis of disease progression, patient survival and response to therapy. Benchmarking is needed for tool assessment and improvement but is complicated by a lack of gold standards, by extensive resource requirements and by difficulties in sharing personal genomic information. To resolve these issues, we launched the ICGC-TCGA DREAM Somatic Mutation Calling Challenge, a crowdsourced benchmark of somatic mutation detection algorithms. Here we report the BAMSurgeon tool for simulating cancer genomes and the results of 248 analyses of three in silico tumors created with it. Different algorithms exhibit characteristic error profiles, and, intriguingly, false positives show a trinucleotide profile very similar to one found in human tumors. Although the three simulated tumors differ in sequence contamination (deviation from normal cell sequence) and in subclonality, an ensemble of pipelines outperforms the best individual pipeline in all cases. BAMSurgeon is available at https://github.com/adamewing/bamsurgeon/.

Genomic hallmarks of localized, non-indolent prostate cancer

Prostate tumours are highly variable in their response to therapies, but clinically available prognostic factors can explain only a fraction of this heterogeneity. Here we analysed 200 whole-genome sequences and 277 additional whole-exome sequences from localized, non-indolent prostate tumours with similar clinical risk profiles, and carried out RNA and methylation analyses in a subset. These tumours had a paucity of clinically actionable single nucleotide variants, unlike those seen in metastatic disease. Rather, a significant proportion of tumours harboured recurrent non-coding aberrations, large-scale genomic rearrangements, and alterations in which an inversion repressed transcription within its boundaries. Local hypermutation events were frequent, and correlated with specific genomic profiles. Numerous molecular aberrations were prognostic for disease recurrence, including several DNA methylation events, and a signature comprised of these aberrations outperformed well-described prognostic biomarkers. We suggest that intensified treatment of genomically aggressive localized prostate cancer may improve cure rates.

Transcriptional profiling of rat white adipose tissue response to 2,3,7,8-tetrachlorodibenzo-ρ-dioxin

Polychlorinated dibenzodioxins are environmental contaminants commonly produced as a by-product of industrial processes. The most potent of these, 2,3,7,8-tetrachlorodibenzo-ρ-dioxin (TCDD), is highly lipophilic, leading to bioaccumulation. White adipose tissue (WAT) is a major site for energy storage, and is one of the organs in which TCDD accumulates. In laboratory animals, exposure to TCDD causes numerous metabolic abnormalities, including a wasting syndrome. We therefore investigated the molecular effects of TCDD exposure on WAT by profiling the transcriptomic response of WAT to 100μg/kg of TCDD at 1 or 4days in TCDD-sensitive Long-Evans (Turku/AB; L-E) rats. A comparative analysis was conducted simultaneously in identically treated TCDD-resistant Han/Wistar (Kuopio; H/W) rats one day after exposure to the same dose. We sought to identify transcriptomic changes coinciding with the onset of toxicity, while gaining additional insight into later responses. More transcriptional responses to TCDD were observed at 4days than at 1day post-exposure, suggesting WAT shows mostly secondary responses. Two classic AHR-regulated genes, Cyp1a1 and Nqo1, were significantly induced by TCDD in both strains, while several genes involved in the immune response, including Ms4a7 and F13a1 were altered in L-E rats alone. We compared genes affected by TCDD in rat WAT and human adipose cells, and observed little overlap. Interestingly, very few genes involved in lipid metabolism exhibited altered expression levels despite the pronounced lipid mobilization from peripheral fat pads by TCDD in L-E rats. Of these genes, the lipolysis-associated Lpin1 was induced slightly over 2-fold in L-E rat WAT on day 4.

BPG: Seamless, Automated and Interactive Visualization of Scientific Data

We introduce BPG, an easy-to-use framework for generating publication-quality, highly-customizable plots in the R statistical environment. This open-source package includes novel methods of displaying high-dimensional datasets and facilitates generation of complex multi-panel figures, making it ideal for complex datasets. A web-based interactive tool allows online figure customization, from which R code can be downloaded for seamless integration with computational pipelines. BPG is available at http://labs.oicr.on.ca/boutros-lab/software/bpg

Cross-species transcriptomic analysis elucidates constitutive aryl hydrocarbon receptor activity

BackgroundResearch on the aryl hydrocarbon receptor (AHR) has largely focused on variations in toxic outcomes resulting from its activation by halogenated aromatic hydrocarbons. But the AHR also plays key roles in regulating pathways critical for development, and after decades of research the mechanisms underlying physiological regulation by the AHR remain poorly characterized. Previous studies identified several core genes that respond to xenobiotic AHR ligands across a broad range of species and tissues. However, only limited inferences have been made regarding its role in regulating constitutive gene activity, i.e. in the absence of exogenous ligands. To address this, we profiled transcriptomic variations between AHR-active and AHR-less-active animals in the absence of an exogenous agonist across five tissues, three of which came from rats (hypothalamus, white adipose and liver) and two of which came from mice (kidney and liver). Because AHR status alone has been shown sufficient to alter transcriptomic responses, we reason that by contrasting profiles amongst AHR-variant animals, we may elucidate effects of the AHR on constitutive mRNA abundances.ResultsWe found significantly more overlap in constitutive mRNA abundances amongst tissues within the same species than from tissues between species and identified 13 genes (Agt, Car3, Creg1, Ctsc, E2f6, Enpp1, Gatm, Gstm4, Kcnj8, Me1, Pdk1, Slc35a3, and Sqrdl) that are affected by AHR-status in four of five tissues. One gene, Creg1, was significantly up-regulated in all AHR-less-active animals. We also find greater overlap between tissues at the pathway level than at the gene level, suggesting coherency to the AHR signalling response within these processes. Analysis of regulatory motifs suggests that the AHR mostly mediates transcriptional regulation via direct binding to response elements.ConclusionsThese findings, though preliminary, present a platform for further evaluating the role of the AHR in regulation of constitutive mRNA levels and physiologic function.

Compendium of TCDD-mediated transcriptomic response datasets in mammalian model systems

Background2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is the most potent congener of the dioxin class of environmental contaminants. Exposure to TCDD causes a wide range of toxic outcomes, ranging from chloracne to acute lethality. The severity of toxicity is highly dependent on the aryl hydrocarbon receptor (AHR). Binding of TCDD to the AHR leads to changes in transcription of numerous genes. Studies evaluating the transcriptional changes brought on by TCDD may provide valuable insight into the role of the AHR in human health and disease. We therefore compiled a collection of transcriptomic datasets that can be used to aid the scientific community in better understanding the transcriptional effects of ligand-activated AHR.ResultsSpecifically, we have created a datasets package – TCDD.Transcriptomics – for the R statistical environment, consisting of 63 unique experiments comprising 377 samples, including various combinations of 3 species (human derived cell lines, mouse and rat), 4 tissue types (liver, kidney, white adipose tissue and hypothalamus) and a wide range of TCDD exposure times and doses. These datasets have been fully standardized using consistent preprocessing and annotation packages (available as of September 14, 2015). To demonstrate the utility of this R package, a subset of “AHR-core” genes were evaluated across the included datasets. Ahrr, Nqo1 and members of the Cyp family were significantly induced following exposure to TCDD across the studies as expected while Aldh3a1 was induced specifically in rat liver. Inmt was altered only in liver tissue and primarily by rat-AHR.ConclusionsAnalysis of the “AHR-core” genes demonstrates a continued need for studies surrounding the impact of AHR-activity on the transcriptome; genes believed to be consistently regulated by ligand-activated AHR show surprisingly little overlap across species and tissues. Until now, a comprehensive assessment of the transcriptome across these studies was challenging due to differences in array platforms, processing methods and annotation versions. We believe that this package, which is freely available for download (http://labs.oicr.on.ca/boutros-lab/tcdd-transcriptomics) will prove to be a highly beneficial resource to the scientific community evaluating the effects of TCDD exposure as well as the variety of functions of the AHR.

Germline contamination and leakage in whole genome somatic single nucleotide variant detection

BackgroundThe clinical sequencing of cancer genomes to personalize therapy is becoming routine across the world. However, concerns over patient re-identification from these data lead to questions about how tightly access should be controlled. It is not thought to be possible to re-identify patients from somatic variant data. However, somatic variant detection pipelines can mistakenly identify germline variants as somatic ones, a process called “germline leakage”. The rate of germline leakage across different somatic variant detection pipelines is not well-understood, and it is uncertain whether or not somatic variant calls should be considered re-identifiable. To fill this gap, we quantified germline leakage across 259 sets of whole-genome somatic single nucleotide variant (SNVs) predictions made by 21 teams as part of the ICGC-TCGA DREAM Somatic Mutation Calling Challenge.ResultsThe median somatic SNV prediction set contained 4325 somatic SNVs and leaked one germline polymorphism. The level of germline leakage was inversely correlated with somatic SNV prediction accuracy and positively correlated with the amount of infiltrating normal cells. The specific germline variants leaked differed by tumour and algorithm. To aid in quantitation and correction of leakage, we created a tool, called GermlineFilter, for use in public-facing somatic SNV databases.ConclusionsThe potential for patient re-identification from leaked germline variants in somatic SNV predictions has led to divergent open data access policies, based on different assessments of the risks. Indeed, a single, well-publicized re-identification event could reshape public perceptions of the values of genomic data sharing. We find that modern somatic SNV prediction pipelines have low germline-leakage rates, which can be further reduced, especially for cloud-sharing, using pre-filtering software.

Accurate Reference-Free Somatic Variant-Calling by Integrating Genomic, Sequencing and Population Data

The detection of somatic single nucleotide variants (SNVs) is critical in both research and clinical applications. Studies of human cancer typically use matched normal (reference) samples from a distant tissue to increase SNV prediction accuracy. This process both doubles sequencing costs and poses challenges when reference samples are not readily available, such as for many cell-lines. To address these challenges, we created S22S: an approach for the prediction of somatic mutations without need for matched reference tissue. S22S takes underlying sequence data, augments them with genomic background context and population frequency information, and classifies SNVs as somatic or non-somatic. We validated S22S using primary tumor/normal pairs from four tumor types, spanning two different sequencing technologies. S22S robustly identifies somatic SNVs, with the area under the precision recall curve reaching 0.97 in kidney clear cell carcinoma, comparable to the best tumor/normal analysis pipelines. S22S is freely available at http://labs.oicr.on.ca/Boutros-lab/software/s22s.