Epidemiologic Factors in Relation to Colorectal Cancer Risk and Survival by Genotoxic Colibactin Mutational Signature

BACKGROUND

The genotoxin colibactin causes a tumor single-base substitution (SBS) mutational signature, SBS88. It is unknown whether epidemiologic factors' association with colorectal cancer risk and survival differs by SBS88.

METHODS

Within the Genetic Epidemiology of Colorectal Cancer Consortium and Colon Cancer Family Registry, we measured SBS88 in 4,308 microsatellite stable/microsatellite instability low tumors. Associations of epidemiologic factors with colorectal cancer risk by SBS88 were assessed using multinomial regression (N = 4,308 cases, 14,192 controls; cohort-only cases N = 1,911), and with colorectal cancer-specific survival using Cox proportional hazards regression (N = 3,465 cases).

RESULTS

392 (9%) tumors were SBS88 positive. Among all cases, the highest quartile of fruit intake was associated with lower risk of SBS88-positive colorectal cancer than SBS88-negative colorectal cancer [odds ratio (OR) = 0.53, 95% confidence interval (CI) 0.37-0.76; OR = 0.75, 95% CI 0.66-0.85, respectively, Pheterogeneity = 0.047]. Among cohort studies, associations of body mass index (BMI), alcohol, and fruit intake with colorectal cancer risk differed by SBS88. BMI ≥30 kg/m2 was associated with worse colorectal cancer-specific survival among those SBS88-positive [hazard ratio (HR) = 3.40, 95% CI 1.47-7.84], but not among those SBS88-negative (HR = 0.97, 95% CI 0.78-1.21, Pheterogeneity = 0.066).

CONCLUSIONS

Most epidemiologic factors did not differ by SBS88 for colorectal cancer risk or survival. Higher BMI may be associated with worse colorectal cancer-specific survival among those SBS88-positive; however, validation is needed in samples with whole-genome or whole-exome sequencing available.

IMPACT

This study highlights the importance of identification of tumor phenotypes related to colorectal cancer and understanding potential heterogeneity for risk and survival.

Genotoxic colibactin mutational signature in colorectal cancer is associated with clinicopathological features, specific genomic alterations and better survival

Background and Aims

The microbiome has long been suspected of a role in colorectal cancer (CRC) tumorigenesis. The mutational signature SBS88 mechanistically links CRC development with the strain of Escherichia coli harboring the pks island that produces the genotoxin colibactin, but the genomic, pathological and survival characteristics associated with SBS88-positive tumors are unknown.

Methods

SBS88-positive CRCs were identified from targeted sequencing data from 5,292 CRCs from 17 studies and tested for their association with clinico-pathological features, oncogenic pathways, genomic characteristics and survival.

Results

In total, 7.5% (398/5,292) of the CRCs were SBS88-positive, of which 98.7% (392/398) were microsatellite stable/microsatellite instability low (MSS/MSI-L), compared with 80% (3916/4894) of SBS88 negative tumors (p=1.5x10-28). Analysis of MSS/MSI-L CRCs demonstrated that SBS88 positive CRCs were associated with the distal colon (OR=1.84, 95% CI=1.40-2.42, p=1x10-5) and rectum (OR=1.90, 95% CI=1.44-2.51, p=6x10-6) tumor sites compared with the proximal colon. The top seven recurrent somatic mutations associated with SBS88-positive CRCs demonstrated mutational contexts associated with colibactin-induced DNA damage, the strongest of which was the APC:c.835-8A>G mutation (OR=65.5, 95%CI=39.0-110.0, p=3x10-80). Large copy number alterations (CNAs) including CNA loss on 14q and gains on 13q, 16q and 20p were significantly enriched in SBS88-positive CRCs. SBS88-positive CRCs were associated with better CRC-specific survival (p=0.007; hazard ratio of 0.69, 95% CI=0.52-0.90) when stratified by age, sex, study, and by stage.

Conclusion

SBS88-positivity, a biomarker of colibactin-induced DNA damage, can identify a novel subtype of CRC characterized by recurrent somatic mutations, copy number alterations and better survival. These findings provide new insights for treatment and prevention strategies for this subtype of CRC.

Body mass index and molecular subtypes of colorectal cancer

BACKGROUND

Obesity is an established risk factor for colorectal cancer (CRC), but the evidence for the association is inconsistent across molecular subtypes of the disease.

METHODS

We pooled data on body mass index (BMI), tumor microsatellite instability status, CpG island methylator phenotype status, BRAF and KRAS mutations, and Jass classification types for 11 872 CRC cases and 11 013 controls from 11 observational studies. We used multinomial logistic regression to estimate odds ratios (OR) and 95% confidence intervals (CI) adjusted for covariables.

RESULTS

Higher BMI was associated with increased CRC risk (OR per 5 kg/m2 = 1.18, 95% CI = 1.15 to 1.22). The positive association was stronger for men than women but similar across tumor subtypes defined by individual molecular markers. In analyses by Jass type, higher BMI was associated with elevated CRC risk for types 1-4 cases but not for type 5 CRC cases (considered familial-like/Lynch syndrome microsatellite instability-H, CpG island methylator phenotype-low or negative, BRAF-wild type, KRAS-wild type, OR = 1.04, 95% CI = 0.90 to 1.20). This pattern of associations for BMI and Jass types was consistent by sex and design of contributing studies (cohort or case-control).

CONCLUSIONS

In contrast to previous reports with fewer study participants, we found limited evidence of heterogeneity for the association between BMI and CRC risk according to molecular subtype, suggesting that obesity influences nearly all major pathways involved in colorectal carcinogenesis. The null association observed for the Jass type 5 suggests that BMI is not a risk factor for the development of CRC for individuals with Lynch syndrome.

A brain precursor atlas reveals the acquisition of developmental-like states in adult cerebral tumours

Human cerebral cancers are known to contain cell types resembling the varying stages of neural development. However, the basis of this association remains unclear. Here, we map the development of mouse cerebrum across the developmental time-course, from embryonic day 12.5 to postnatal day 365, performing single-cell transcriptomics on >100,000 cells. By comparing this reference atlas to single-cell data from >100 glial tumours of the adult and paediatric human cerebrum, we find that tumour cells have an expression signature that overlaps with temporally restricted, embryonic radial glial precursors (RGPs) and their immediate sublineages. Further, we demonstrate that prenatal transformation of RGPs in a genetic mouse model gives rise to adult cerebral tumours that show an embryonic/juvenile RGP identity. Together, these findings implicate the acquisition of embryonic-like states in the genesis of adult glioma, providing insight into the origins of human glioma, and identifying specific developmental cell types for therapeutic targeting.

Association between germline variants and somatic mutations in colorectal cancer

Colorectal cancer (CRC) is a heterogeneous disease with evidence of distinct tumor types that develop through different somatically altered pathways. To better understand the impact of the host genome on somatically mutated genes and pathways, we assessed associations of germline variations with somatic events via two complementary approaches. We first analyzed the association between individual germline genetic variants and the presence of non-silent somatic mutations in genes in 1375 CRC cases with genome-wide SNPs data and a tumor sequencing panel targeting 205 genes. In the second analysis, we tested if germline variants located within previously identified regions of somatic allelic imbalance were associated with overall CRC risk using summary statistics from a recent large scale GWAS (n≃125 k CRC cases and controls). The first analysis revealed that a variant (rs78963230) located within a CNA region associated with TLR3 was also associated with a non-silent mutation within gene FBXW7. In the secondary analysis, the variant rs2302274 located in CDX1/PDGFRB frequently gained/lost in colorectal tumors was associated with overall CRC risk (OR = 0.96, p = 7.50e-7). In summary, we demonstrate that an integrative analysis of somatic and germline variation can lead to new insights about CRC.

Identifying colorectal cancer caused by biallelic MUTYH pathogenic variants using tumor mutational signatures

Carriers of germline biallelic pathogenic variants in the MUTYH gene have a high risk of colorectal cancer. We test 5649 colorectal cancers to evaluate the discriminatory potential of a tumor mutational signature specific to MUTYH for identifying biallelic carriers and classifying variants of uncertain clinical significance (VUS). Using a tumor and matched germline targeted multi-gene panel approach, our classifier identifies all biallelic MUTYH carriers and all known non-carriers in an independent test set of 3019 colorectal cancers (accuracy = 100% (95% confidence interval 99.87-100%)). All monoallelic MUTYH carriers are classified with the non-MUTYH carriers. The classifier provides evidence for a pathogenic classification for two VUS and a benign classification for five VUS. Somatic hotspot mutations KRAS p.G12C and PIK3CA p.Q546K are associated with colorectal cancers from biallelic MUTYH carriers compared with non-carriers (p = 2 × 10-23 and p = 6 × 10-11, respectively). Here, we demonstrate the potential application of mutational signatures to tumor sequencing workflows to improve the identification of biallelic MUTYH carriers.

Evaluating the predictive accuracy of curated biological pathways in a public knowledgebase

Abstract

Reactome is a database of human biological pathways manually curated from the primary literature and peer-reviewed by experts. To evaluate the utility of Reactome pathways for predicting functional consequences of genetic perturbations, we compared predictions of perturbation effects based on Reactome pathways against published empirical observations. Ten cancer-relevant Reactome pathways, representing diverse biological processes such as signal transduction, cell division, DNA repair and transcriptional regulation, were selected for testing. For each pathway, root input nodes and key pathway outputs were defined. We then used pathway-diagram-derived logic graphs to predict, either by inspection by biocurators or using a novel algorithm MP-BioPath, the effects of bidirectional perturbations (upregulation/activation or downregulation/inhibition) of single root inputs on the status of key outputs. These predictions were then compared to published empirical tests. In total, 4968 test cases were analyzed across 10 pathways, of which 847 were supported by published empirical findings. Out of the 847 test cases, curators’ predictions agreed with the experimental evidence in 670 and disagreed in 177 cases, resulting in ∼81% overall accuracy. MP-BioPath predictions agreed with experimental evidence for 625 and disagreed for 222 test cases, resulting in ∼75% overall accuracy. The expected accuracy of random guessing was 33%. Per-pathway accuracy did not correlate with the number of pathway edges nor the number of pathway nodes but varied across pathways, ranging from 56% (curator)/44% (MP-BioPath) for ‘Mitotic G1 phase and G1/S transition’ to 100% (curator)/94% (MP-BioPath) for ‘RAF/MAP kinase cascade’. This study highlights the potential of pathway databases such as Reactome in modeling genetic perturbations, promoting standardization of experimental pathway activity readout and supporting hypothesis-driven research by revealing relationships between pathway inputs and outputs that have not yet been directly experimentally tested.

Database URL

www.reactome.org

Molecular and Pathology Features of Colorectal Tumors and Patient Outcomes Are Associated with Fusobacterium nucleatum and Its Subspecies animalis

BACKGROUND

Fusobacterium nucleatum (F. nucleatum) activates oncogenic signaling pathways and induces inflammation to promote colorectal carcinogenesis.

METHODS

We characterized F. nucleatum and its subspecies in colorectal tumors and examined associations with tumor characteristics and colorectal cancer-specific survival. We conducted deep sequencing of nusA, nusG, and bacterial 16s rRNA genes in tumors from 1,994 patients with colorectal cancer and assessed associations between F. nucleatum presence and clinical characteristics, colorectal cancer-specific mortality, and somatic mutations.

RESULTS

F. nucleatum, which was present in 10.3% of tumors, was detected in a higher proportion of right-sided and advanced-stage tumors, particularly subspecies animalis. Presence of F. nucleatum was associated with higher colorectal cancer-specific mortality (HR, 1.97; P = 0.0004). This association was restricted to nonhypermutated, microsatellite-stable tumors (HR, 2.13; P = 0.0002) and those who received chemotherapy [HR, 1.92; confidence interval (CI), 1.07-3.45; P = 0.029). Only F. nucleatum subspecies animalis, the main subspecies detected (65.8%), was associated with colorectal cancer-specific mortality (HR, 2.16; P = 0.0016), subspecies vincentii and nucleatum were not (HR, 1.07; P = 0.86). Additional adjustment for tumor stage suggests that the effect of F. nucleatum on mortality is partly driven by a stage shift. Presence of F. nucleatum was associated with microsatellite instable tumors, tumors with POLE exonuclease domain mutations, and ERBB3 mutations, and suggestively associated with TP53 mutations.

CONCLUSIONS

F. nucleatum, and particularly subspecies animalis, was associated with a higher colorectal cancer-specific mortality and specific somatic mutated genes.

IMPACT

Our findings identify the F. nucleatum subspecies animalis as negatively impacting colorectal cancer mortality, which may occur through a stage shift and its effect on chemoresistance.

Abstract LB090: Associations of somatically mutated genes and pathways with colorectal cancer specific survival in 4,500 colorectal cancer patients

Colorectal cancer (CRC) is a heterogenous disease that develops through somatic mutations in driver genes, leading to activation of diverse neoplastic pathways. To systematically examine if somatically mutated genes and pathways impact survival, we sequenced tumor and normal DNA samples for 4,512 CRC cases using a targeted panel. We performed Cox regression to estimate hazard ratios (HR) and 95% confidence intervals (CIs) for associations of disease-specific (DS)-survival with somatically mutated genes and pathways, adjusting for age and sex, and stratifying baseline hazards by study population. We assessed statistical significance using Bonferroni p-value thresholds to account for multiple testing of 214 genes (2.34x10-4) and 6 key CRC pathways (8.3x10-3). We limited analyses to non-silent mutations.

We observed that DS-survival was significantly more favorable among individuals with hypermutated (HM) tumors, primarily consisting of microsatellite unstable and POLE-mutated tumors, as compared to those with non-hypermutated (NHM) tumors (HR=0.4, 95% CI: 0.3-0.5, p=1.2x10-17). BRAF V600E mutations were associated with poorer survival (HR=2.0, 95% CI: 1.6-2.5, p=1.6x10-10). This association was more pronounced among NHM tumors (HR=2.3, 95% CI: 1.8-2.9, p=4.3x10-12). We identified suggestive associations (p-values < 5.0x10-3 in overall or stratified analyses) between DS-survival and mutations in B2M, TP53, and SMAD4. Mutations in B2M may provide more favorable prognosis for survival (HR=0.5, 95% CI: 0.3-0.8, p=4.4x10-3), with similar effect sizes in HM and NHM tumors. Poorer survival may be associated with mutations in TP53 (HR=1.2, 95% CI:1.2-1.4, p=9.3x10-4) and SMAD4 (HR=1.3, 95% CI:1.1-1.6, p=3.2x10-3) in NHM tumors.

We further observed statistically significant associations between survival and three mutated pathways: TP53/ATM (HR=1.2, 95% CI:1.1-1.4, p=9.0x10-4), receptor tyrosine kinases (RTK) and RAS (HR=1.3, 95% CI:1.1-1.5, p=4.5x10-5), and TGF-beta (HR=1.3, 95% CI:1.1-1.5, p=2.9x10-4). Findings for TP53/ATM and RTK/RAS were primarily due to mutations in one gene within each pathway (TP53 and BRAF, respectively). However, the TGF-beta pathway finding, which was more pronounced in NHM tumors, included two genes with p-values less than 0.05 (SMAD4 p=3.2x10-3 and TGFBR2 p=0.01).

Despite our large sample size, relatively few somatically mutated genes were significantly associated with DS-survival. It may be that pathway-level testing affords more power for analyses and that larger sample sizes are needed. It is of interest that non-silent mutations in B2M, whose product is a component of the class I major histocompatibility complex, was associated with improved DS-survival, though our finding was not statistically significant after multiple testing correction so further exploration and replication is needed.

Citation Format: Tabitha A. Harrison, Syed H. Zaidi, Conghui Qu, Amanda I. Phipps, Robert S. Steinfelder, Quang M. Trinh, Sonja I. Berndt, Daniel D. Buchanan, Peter T. Campbell, Andrew T. Chan, Kim F. Doheny, David A. Drew, Jane C. Figueiredo, Steven J. Gallinger, Marios Giannakis, Andrea Gsur, Marc J. Gunter, Michael Hoffmeister, Wen-Yi Huang, Paul J. Limburg, Victor Moreno, Polly A. Newcomb, Shuji Ogino, Ross L. Prentice, Tameka Shelford, Wei Sun, Stephen N. Thibodeau, Li Hsu, Ulrike Peters. Associations of somatically mutated genes and pathways with colorectal cancer specific survival in 4,500 colorectal cancer patients [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr LB090.

Abstract 3131: Tumour spatial heterogeneity in breast cancer and the impact on clinical management

We recognize that many cancers are highly complex mixtures of sub-populations of cells, which are also influenced by their microenvironment. This heterogeneity explains in part, why the majority of current therapeutic approaches for cancer work best when multiple agents are combined. Therefore, in an era of targeted therapeutics it becomes critical to understand the complexity of tumors both at diagnosis and over the course of therapy, including measures of heterogeneity. Here we present genomic, transcriptomic and in situ proteomic profiling of a cohort of breast cancer (BCa) lumpectomies to reveal the extent of heterogeneity in pathologically defined unifocal and multifocal cancers. Integration of molecular profiling, phylogenetic analyses and radiomics has the potential to significantly improve BCa clinical management and stratification to targeted therapies that are already available in the clinic. In this study, lumpectomies with imaging data were processed as whole mounts with serial blocks reviewed. Tissue cores were taken from at least three different regions through the lumpectomy for nucleic acid extraction and tissue microarray (TMA) construction, focusing on morphologic/histological differences in addition to the spatial orientation of the sampled region within the lumpectomy. Targeted sequencing using the Oncomine Comprehensive Assay v3 (OCAv3); transcriptional profiling using the NanoString Breast Cancer 360 Panel; in situ profiling by multiplex fluorescence immunohistochemistry (MxIF) and Digital Spatial Profiling (see abstract Spears et al) were performed. From this cohort of 60 patients, we present a subset of patients demonstrating integration of the molecular profiling to reveal the phylogenetic relationship between the multiple samplings and the impact on clinical decision making in BCa. Briefly, we identified differences in the molecular subtypes between the different sample regions from the same unifocal cancer as well as differences in the predicted responses to anti-PDL1 therapy by transcriptional profiling; while targeted sequencing of driver mutations suggested the likelihood of an ancestral tumor cell giving rise to the lesions in pathologically defined multifocal cancers. However it was evident that genes and pathways found to be aberrant in these different lesions from the same cancer could impact the response to standard BCa treatment, or the selection of targeted therapies. In situ proteomics demonstrated differences in the expression of standard BCa markers ER, PgR, HER2 and Ki67 in addition to immune markers in the tumor and tumor microenvironment. While there are clinically validated transcriptional risk test available for BCa, we have demonstrated that transcriptomics or genomics alone is insufficient for a rational stratification of patients to currently available targeted therapies; therefore supporting the need for an integrative approach.

Citation Format: Jane Bayani, Quang M. Trinh, Mary Anne Quintayo, Cheryl Crozier, Megan Hopkins, Jingping Qiao, Alison Cheung, James Mainprize, Quaid Morris, Melanie Spears, Martin Yaffe, Lincoln Stein, John M. Bartlett. Tumour spatial heterogeneity in breast cancer and the impact on clinical management [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 3131.

Abstract PO-002: Revealing tumour spatial heterogeneity in breast cancer and the impact on clinical management

Many cancers are highly complex mixtures of many sub-populations of cells, also influenced by their microenvironment. This heterogeneity explains, in part, why the majority of current therapeutic approaches for cancer work best when multiple agents are combined. Therefore, in an era of targeted therapeutics it becomes critical to understand the complexity of tumours both at diagnosis and over the course of therapy, including measures of heterogeneity. Here we present genomic, transcriptomic and in situ proteomic profiling of a cohort of breast cancer (BCa) lumpectomies with associated imaging data to reveal the extent of heterogeneity in pathologically defined unifocal and multifocal cancers. Integration of molecular profiling, phylogenetic analyses and radiomics has the potential to significantly improve BCa clinical management and stratification to targeted therapies that are already available in the clinic. In this study, lumpectomies were processed as whole mounts with serial blocks reviewed. Tissue cores were taken from at least three different regions throughout the lumpectomy for nucleic acid extraction and tissue microarray (TMA) construction, focusing on morphologic/histological differences in addition to the spatial orientation of the sampled region within the lumpectomy. Targeted sequencing using the Oncomine Comprehensive Assay v3 (OCAv3); transcriptional profiling using the NanoString Breast Cancer 360 Panel; and in situ profiling by multiplex fluorescence immunohistochemistry (MxIF) performed (see abstract Cheung et al). From this cohort of 60 patients, we present a subset of patients demonstrating integration of the molecular profiling to reveal the phylogenetic relationship between the multiple samplings and the potential impact on clinical decision making in BCa. We identified differences in the molecular subtypes between the different sample regions from the same unifocal cancer as well as differences in the predicted responses to anti-PDL1 therapy by transcriptional profiling. Targeted sequencing of driver mutations suggested the likelihood of an ancestral tumour cell giving rise to the lesions in pathologically defined multifocal cancers; however it was evident that genes and pathways found to be aberrant in these different lesions from the same cancer could impact the response to standard BCa treatment, or the selection of targeted therapies. In situ proteomics demonstrated differences in the expression of standard BCa markers ER, PgR, HER2 and Ki67, in addition to immune markers in the tumour and its microenvironment. While there are clinically validated transcriptional risk tests available for BCa, we have demonstrated that transcriptomics or genomics alone is insufficient for a rational stratification of patients to currently available targeted therapies, therefore, supporting the need for an integrative approach.

Citation Format: Jane Bayani, Quang M. Trinh, Mary Anne Quintayo, Cheryl Crozier, Megan Hopkins, Jingping Qiao, Alison Cheung, James G Mainprize, Quaid Morris, Melanie Spears, Martin J. Yaffe, Lincoln D. Stein, John M.S. Bartlett. Revealing tumour spatial heterogeneity in breast cancer and the impact on clinical management [abstract]. In: Proceedings of the AACR Virtual Special Conference on Tumor Heterogeneity: From Single Cells to Clinical Impact; 2020 Sep 17-18. Philadelphia (PA): AACR; Cancer Res 2020;80(21 Suppl):Abstract nr PO-002.

Abstract P3-08-22: The mutational landscape of cancer driver genes in matched primary ductal carcinoma in situ and recurrent ductal carcinoma in situ or recurrent invasive cancers

Due to breast screening, ductal carcinoma in situ (DCIS) accounts for approximately 25% of all newly diagnosed breast neoplasms. Believed to be a precursor to invasive carcinoma, a significant number of patients diagnosed with DCIS are effectively managed by surgery alone or in conjunction with radiotherapy and endocrine therapy. In general, there is an 8-11% relative risk for a subsequent invasive carcinoma over a period of 10 year, with 98% breast cancer-specific survival after 10 years of follow up. Although mastectomy, breast conserving surgery, and radiotherapy can reduce the risk of recurrence, there are ongoing lifetime consequences of treatment. Thus, there is a clinical need to identify those patients who are at risk of an invasive recurrence from those who might not recur or experience a subsequent DCIS recurrence. In this study, 60 patients with pure primary DCIS, treated with only with breast conserving therapy, across three different clinical outcome groups were examined: patients who did not experience a recurrence within 5 years (n=20); patients who experienced a recurrent DCIS within 5 years (n=20); and patients that recurred with an invasive cancer within 5 years (n=20). Pure primary DCIS lesions, as well as the matched DCIS recurrence or invasive recurrence, were macrodissected from formalin fixed paraffin embedded tissues and subjected to nucleic acid extraction. All samples were profiled using Thermo Fisher Scientific’s validated targeted sequencing panel, the Oncomine Comprehensive Assay v3.0 (OCAv3.0). This assay is comprised of common cancer driver genes shown to be prognostic and predictive to targeted therapies in use or in late-phase clinical trials, and is currently being used in the NCI-MATCH trial (NCT02465060). While the OCA panel and accompanying Oncomine Knowledgebase Reporter provides information regarding the targeted treatments linked to known actionable mutations, this study utilized all somatic mutations and copy number changes to reveal the genomic landscape of DCIS and their matched recurrences across these pan-cancer driver genes Amongst all primary DCIS samples across the three different clinical outcome groups, PIK3CA, was frequently found to be affected by SNV and Indels (32.8%) in addition to TP53 (26.2%), NF1 (21.3%), CREBBP (16.4%), ATM (14.8%), PALB2 (14.8%). DNA repair genes, including CHEK1, RAD50, RAD51B, MRE11A, BRCA1 and BRCA2, were found to be frequently subject to mutation in these primary DCIS samples ranging from 5%-15%. Similarly, copy number gains were frequently detected in HER2 (26.2%), CDK12 (18%), FGFR1 (4.9%), GNAS (4.9%), MYC, CCNE1, AR, RAD51C and RNF43 (1.6% each), and loses at H3F3A and KNSTRN (each 1.6%). Matched primary DCIS and their recurrent DCIS or invasive lesions exhibited similar changes with invasive cancers suggesting that in some cases, the primary DCIS gives rise to the invasive cancer. We will present the preliminary findings mapping the mutational and copy-number landscape of primary DCIS and matched recurrences to identify putative genomic changes defining these clinical outcome groups and to investigate the genomic progression of DCIS to invasive carcinoma.

Citation Format: Jane Bayani, Quang M Trinh, Mary Anne Quintayo, Cheryl Crozier, Ilinca Lungu, Dan Dion, Joema Felipe-Lima, Giancarlo Pruneri, Jonas Bergh, Fredrik Warnberg, Giuseppe Viale, Lincoln D Stein, John MS Bartlett. The mutational landscape of cancer driver genes in matched primary ductal carcinoma in situ and recurrent ductal carcinoma in situ or recurrent invasive cancers [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr P3-08-22.

Abstract P4-10-03: The genomic landscape of male breast cancers using the oncomine comprehensive assay for actionable mutations

Introduction: Male breast cancer (BCa) is a rare disease accounting for less than 1% of all breast cancers (BC) and 1% of all cancers in males. The clinical management is largely extrapolated from female BCa. Few studies have examined the genomic landscape of male BCas, with six male BCas included in The Cancer Genome Atlas (TCGA). Familial studies of male BCas have shown genomic changes similar to female BCa, while a larger targeted sequencing study of 59 male breast cancers identified recurrent mutations affecting PIK3CA and GATA3. To date, there is still limited information regarding the genomic landscape of male BCas; particularly in the context of identifying targeted treatments. To reveal genomic changes that characterize male BCas in the context of known cancer driver genes linked to prognosis and targeted agents, we performed a targeted sequencing study on 248 male BCas from the International Male Breast Cancer Program. Methods: 248 primary M0, ER+ve, HER2-ve male BCas enrolled in the Part 1 (retrospective joint analysis) International Male Breast Cancer Program of 1483 patients diagnosed between 1990-2010 (Cardoso et al. Annals of Oncology, 2018) were processed for nucleic acid extraction from formalin-fixed paraffin embedded (FFPE) tissues. Using the Thermo Fisher Scientific Oncomine Comprehensive Assay v3 (OCAv3), a validated targeted sequencing panel currently used in the NCI-MATCH trial (NCT02465060), we evaluated mutational and copy number variations (CNVs) of genes that are prognostic or predictive to targeted therapies currently in use in the clinic or late-stage clinical trials. The OCAv3 DNA pan-cancer panel assays 115 genes for determining mutational status (48 full coding and 67 hotspot) as well as copy-number assessment in 43 genes. The OCAv3 uses Ampliseq-based technology linked to the Oncomine NGS workflow to identify actionable mutations and CNVs Results: Of the 248 samples assayed, 216 passed strict quality control parameters (87.1%). Using the Oncomine NGS workflow, actionable mutations at ≥5% variant allele frequency (VAF) were most frequently identified in PIK3CA (29.2%), BRCA2 (11.1%), NF1 (11.6%), indels in TP53 (10.6%), ATR (5.6%), ATRX (5.1%), indel BRCA2 (5.1%), TP53 point mutations (4.6%), MET (4.6%), ATM (4.6%), NOTCH2 (4.6%), CHEK1 (4.2%), FANCI (4.2%), PTEN (3.2%) with a number of additional genes identified at lower frequencies. Gene amplifications were most frequently detected in MYC (24.5%), FGFR1 (14.8%), CCND1 (12%), FGF3 (9.7%), FGF19 (9.7%), MDM2 (6.5%), CDK4 (1.4%), FGFR3, MDM4, ERBB2 (0.9% each), and FLT3, AR, MYCL, CDK6, IGF1R, FGFR4, KRAS, AKT3 and ESR1 (0.5% each). Although the results here describe the mutations and copy-number changes deemed to be actionable, further analysis of all non-actionable somatic mutations and CNVs will be presented and compared to female BCas previously assayed using the same panel. Conclusion: In this targeted sequencing study of the largest series of male BCas to our knowledge, we have revealed that PIK3CA continues to be a frequently altered gene in both male and female BCas. However, there is an enrichment of mutations in genes related to DNA repair in male BCs. Interestingly, while MYC is commonly amplified in female BCa, a higher frequency of amplified cases were seen in male BCas, in contrast to female BCas. Together with our previously generated transcriptional profiling data in this data set, we believe that both common and unique biological processes comprising male and female BCas will ultimately improve their clinical management and move towards the goal of precision medicine.

This work has been funded by the Breast Cancer Research Foundation (BCRF).

Citation Format: Jane Bayani, Coralie Poncet, Cheryl Crozier, Quang M Trinh, Megan Hopkins, Aime Lambert Uwimana, Tammy Piper, Carrie Cunningham, Monika Sobol, Stefan Aebi, Kim Benstead, Oliver Bogler, Lissandra Dal Lago, Florentine Hilbers, Ingrid Hedenfalk, Larissa Korde, Barbro Linderholm, John Martens, Lavinia Middleton, Melissa Murray, Catherine Kelly, Cecilia Nilsson, Monika Nowaczyk, Stephanie Peeters, Aleksandra Peric, Peggy Porter, Carolien Schröder, Isabel T Rubio, Kathryn J Ruddy, Christi van Asperen, Danielle Van Den Weyngaert, Carolien HM van Deurzen, Elise van Leeuwen-Stok, Joanna Vermeij, Eric Winer, Lincoln D Stein, Sharon H Giordano, Fatima Cardoso, John MS Bartlett. The genomic landscape of male breast cancers using the oncomine comprehensive assay for actionable mutations [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr P4-10-03.

Identification of Distinct Prognostic Groups: Implications for Patient Selection to Targeted Therapies Among Anti-Endocrine Therapy–Resistant Early Breast Cancers

PURPOSE

Hormone receptor–positive breast cancer remains an ongoing therapeutic challenge, despite optimal anti-endocrine therapies. In this study, we assessed the prognostic ability of genomic signatures to identify patients at risk for recurrence after endocrine therapy. Analysis was performed on the basis of an a priori hypothesis related to molecular pathways, which might predict response to existing targeted therapies.

PATIENTS AND METHODS

A subset of patients from the Tamoxifen Versus Exemestane Adjuvant Multinational trial ( ClinicalTrials.gov identifiers: NCT00279448 and NCT00032136, and NCT00036270) pathology cohort were analyzed to determine the prognostic ability of mutational and copy number aberration biomarkers that represent the cyclin D/cyclin-dependent kinase (CCND/CDK), fibroblast growth factor receptor/fibroblast growth factor (FGFR/FGF), and phosphatidylinositol 3-kinase/protein kinase B (PI3K/ATK) pathways to inform the potential choice of additional therapies to standard endocrine treatment. Copy number analysis and targeted sequencing was performed. Pathways were identified as aberrant if there were copy number aberrations and/or mutations in any of the predetermined pathway genes: CCND1/CCND2/CCND3/CDK4/CDK6, FGFR1/FGFR2/FGFR2/FGFR4, and AKT1/AKT2/PIK3CA/PTEN.

RESULTS

The 390 of 420 samples that passed quality control were analyzed for distant metastasis–free survival between groups. Patients with no changes in the CCND/CDK pathway experienced a better distant metastasis–free survival (hazard ratio, 1.94; 95% CI, 1.45 to 2.61; P < .001) than those who possessed aberrations. In the FGFR/FGF and PI3K/AKT pathways, a similar outcome was observed (hazard ratio, 1.43 [95% CI, 1.07 to 1.92; P = .017] and 1.34 [95% CI, 1.00 to 1.81; P = .053], respectively).

CONCLUSION

We show that aberrations of genes in these pathways are independently linked to a higher risk of relapse after endocrine treatment. Improvement of the clinical management of early breast cancers could be made by identifying those for whom current endocrine therapies are sufficient, thus reducing unnecessary treatment, and secondly, by identifying those who are at high risk for recurrence and linking molecular features that drive these cancers to treatment with targeted therapies.

Abstract P2-10-04: Targeted sequencing in early breast cancer: Identification of novel candidate mutations predictive of anthracycline benefit

Background The use of chemotherapies such as anthracyclines and taxanes have improved overall and disease free survival in breast cancer. For all patients, anthracyclines can have significant toxicities including cardiotoxicity and leukemia. It is therefore essential to select the subset of patients who will receive the optimal overall benefit from anthracycline therapy and to identify molecular pathways driving resistance. To fully understand the impact of mutations in the context of current breast cancer therapy, requires a comprehensive mapping of key molecular events in the context of treatment. We sequenced 101 genes, that were prioritized based on not only gene frequency, but also taking into account the importance of amino acid substitution, type of mutation and network connectivity, in 692 primary tumours to both identify driver genes and pathway cassettes and to understand their clinical significance in response to anthracycline treatment.

Methods We performed targeted sequencing in patients from the BR9601 (n=374) and CCTG MA.5 (n=703) clinical trials. The BR9601 and MA.5 clinical trials examined the effectiveness of combination chemotherapy consisting of CMF (cyclophosphamide, methotrexate and 5-fluorouracil) with or without epirubicin. DNA was extracted, samples were sequenced using AmpliSeq Technology adapted to Illumina and somatic mutations were called using a novel mutation calling pipeline (ISOWN). A priori analyses were performed using distant recurrence free survival (DRFS) as the primary endpoint.

Results: In 692 successfully analysed samples 509 (73.6%) samples exhibited at least one single nucleotide mutation (range 0-54). 94/101 genes were mutated in at least one patient. Only variants in PIK3CA, TP53, CDH1, TLE6, MLL3 and USH2A were detected in 5% or more of samples. TSC22D1, RB1 and ZNF565 were associated with increased risk of distant relapse in multivariate analyses corrected for clinic-pathological variables. No single genes were predictive of anthracycline treatment compared to CMF in multivariate analyses corrected for clinic-pathological variables. Signaling cassettes/modules were designed based on the pathway database, Reactome. Within the signaling cassettes one module was predictive of anthracycline failure. Patients with one or more mutations in this module had an increased risk of distant relapse (HR 0.52, 95% CI 0.29-0.95, p=0.034) when treated with an anthracycline containing chemotherapy regimen compared to CMF (HR 1.34 95% CI 1.05-1.72, p=0.019).

Conclusions: We successfully performed a signaling pathway-based targeted sequencing analysis within predefined signaling modules. We identified a single signaling cassette linked to anthracycline resistance in early breast cancer. However, further work to validate this study in a separate clinical trial is warranted.

Citation Format: Spears M, Kalatskaya I, Trinh QM, Liao L, Chong TM, Crozier C, Dion D, Heisler L, Timms L, Stein LD, Pritchard KI, Levine MN, Shepherd L, Twelves CJ, Bartlett JMS. Targeted sequencing in early breast cancer: Identification of novel candidate mutations predictive of anthracycline benefit [abstract]. In: Proceedings of the 2017 San Antonio Breast Cancer Symposium; 2017 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2018;78(4 Suppl):Abstract nr P2-10-04.

ISOWN: accurate somatic mutation identification in the absence of normal tissue controls

Background

A key step in cancer genome analysis is the identification of somatic mutations in the tumor. This is typically done by comparing the genome of the tumor to the reference genome sequence derived from a normal tissue taken from the same donor. However, there are a variety of common scenarios in which matched normal tissue is not available for comparison.

Results

In this work, we describe an algorithm to distinguish somatic single nucleotide variants (SNVs) in next-generation sequencing data from germline polymorphisms in the absence of normal samples using a machine learning approach. Our algorithm was evaluated using a family of supervised learning classifications across six different cancer types and ~1600 samples, including cell lines, fresh frozen tissues, and formalin-fixed paraffin-embedded tissues; we tested our algorithm with both deep targeted and whole-exome sequencing data. Our algorithm correctly classified between 95 and 98% of somatic mutations with F1-measure ranges from 75.9 to 98.6% depending on the tumor type. We have released the algorithm as a software package called ISOWN (Identification of SOmatic mutations Without matching Normal tissues).

Conclusions

In this work, we describe the development, implementation, and validation of ISOWN, an accurate algorithm for predicting somatic mutations in cancer tissues in the absence of matching normal tissues. ISOWN is available as Open Source under Apache License 2.0 from https://github.com/ikalatskaya/ISOWN.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-017-0446-9) contains supplementary material, which is available to authorized users.

Targeted sequencing in a phase III trial of luminal breast cancer: Identification of novel targets

505

Background: The International Cancer Genome Consortium and The Cancer Genome Atlas have had a global transformative impact on our understanding of cancer. These programs have mapped the genomic landscape of common and rare tumors setting the scene for a comprehensive change in the approach to cancer diagnosis and treatment. However, the task remains incomplete until these mutational events are linked to clinical outcomes in the context of current therapeutic intervention to drive future stratified medicine approaches. Methods: We performed targeted sequencing in patients from the Tamoxifen Exemestane Adjuvant Multicentre trial. DNA was extracted and a 101 gene panel analysed using a novel mutation calling pipeline. Both a priori and machine learning analyses were performed using distant recurrence free survival as the primary endpoint. Results: In 1,491 successfully analyzed samples 1,070 (71.76%) samples exhibited at least one single nucleotide mutation (range 0-94, 1.828+/-0.133, mean+/-s.e.). 98/101 genes were mutated in at least one patient. Only variants in PIK3CA, TP53, MLL3, CDH1 were detected in 5% or more of samples. Twenty genes were associated with increased risk of recurrence in multivariate analyses corrected for clinic-pathological variables, 50% of these genes were involved in transcriptional regulation or RNA/protein processing. In a multivariate analysis, two combined signalling modules were independently prognostic for residual risk following hormone therapy (HRvalidation 3.10 95%CI 1.78-5.40 and HRvalidation 2.70 95%CI 1.57-4.64). Conclusions: We successfully performed a signalling pathway-based targeted sequencing analysis within predefined signalling modules. In supervised and unsupervised analyses we identified multiple signalling cassettes linked to poor outcome in patients with ER+ve breast cancers treated with modern endocrine therapy in the context of a phase III clinical trial. These results identify novel candidates as targets to treat endocrine refractory breast cancers.

Copy-number and targeted sequencing analyses to identify distinct prognostic groups: Implications for patient selection to targeted therapies amongst anti-endocrine therapy resistant early breast cancers

524

Background: Hormone receptor positive breast cancer is a therapeutic challenge. Despite optimal anti-endocrine therapies, most breast cancer deaths follow a diagnosis of early luminal cancer. To understand the impact of multiple aberrations in the context of current therapy, we assessed the prognostic ability of genomic signatures as a putative stratification tool to targeted therapies. Methods: This a priori study is based on molecular pathways which might predict response to targeted therapies. DNA from 420 patients from the phase III TEAM pathology cohort were used. Patients with a distant recurrence within 5 years were matched by clinical variables to those disease-free at follow up. Copy number analysis was performed using the Affymetrix Oncoscan Assay. Targeted sequencing was performed in a subset of samples for genes based on signaling cassettes mined from the ICGC. Pathways were identified as aberrant if there were copy number variations (CNVs) and/or mutations in any of the pre-determined pathway genes: 1) CCND1/CCND2/CCND3/CDK4/CDK6; 2) FGFR1/FGFR2/FGFR2/FGFR4; and 3) AKT1/AKT2/PIK3CA/PTEN. Kaplan-Meier and log-rank analyses were used for DFS between groups. Hazard ratios were calculated using the Cox proportional hazard models adjusted for age, tumour size, grade, lymph node and HER2 status. Results: 390/420 samples passed informatics QC filters. For the CCND/CDK pathway, patients with no CNV changes experienced a better DFS (HR = 1.7, 95% CI 1.3-2.3, p < 0.001). For the FGFR/FGF pathway, a similar outcome is seen among patients without CNVs (HR = 1.5, 95% CI 1.1-2.0; p = 0.005). For AKT/PIK3CA, a decrease in DFS was seen in those with aberrations (HR = 1.4, 95% CI 1.0-1.8, p = 0.03). Conclusions: We demonstrated that CNVs of genes within CDK4/CCND, PIK3CA/AKT and FGFR pathways are independently linked to high risk of relapse following endocrine treatment, with implications for identifying those patients who are at high-risk for recurrence despite optimal anti-endocrine therapy and linking molecular features driving these cancers to targeted therapies.

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.

Tuberculosis and HIV co-infection in Vietnam

Tuberculosis (TB) and human immunodeficiency virus (HIV) infection are leading causes of disease and death in Vietnam, but TB/HIV disease trends and the profile of co-infected patients are poorly described.

METHODS

We examined national TB and HIV notification data to provide a geographic overview and describe relevant disease trends within Vietnam. We also compared the demographic and clinical profiles of TB patients with and without HIV infection.

RESULTS

During the past 10 years (2005-2014) cumulative HIV case numbers and deaths increased to 298,151 and 71,332 respectively, but access to antiretroviral therapy (ART) improved and new infections and deaths declined. From 2011-2014 routine HIV testing of TB patients increased from 58.9% to 72.5% and of all TB patients diagnosed with HIV in 2014, 2,803 (72.4%) received ART. The number of multidrug resistant (MDR)-TB cases enrolled for treatment increased almost 3-fold (578 to 1,532) from 2011-2014. The rate of HIV co-infection in MDR and non-MDR TB cases (51/1,532; 3.3% vs 3,774/100,555; 3.8%; OR 0.77, 95% CI 0.7-1.2) was similar in 2014.

CONCLUSIONS

The care of TB/HIV co-infected patients have shown sustained improvement in Vietnam. Rising numbers of MDR-TB cases is a concern, but this is not "driven" by HIV co-infection.

Tuberculosis and HIV co-infection-focus on the Asia-Pacific region

Tuberculosis (TB) is the leading opportunistic disease and cause of death in patients with HIV infection. In 2013 there were 1.1 million new TB/HIV co-infected cases globally, accounting for 12% of incident TB cases and 360,000 deaths. The Asia-Pacific region, which contributes more than a half of all TB cases worldwide, traditionally reports low TB/HIV co-infection rates. However, routine testing of TB patients for HIV infection is not universally implemented and the estimated prevalence of HIV in new TB cases increased to 6.3% in 2013. Although HIV infection rates have not seen the rapid rise observed in Sub-Saharan Africa, indications are that rates are increasing among specific high-risk groups. This paper reviews the risks of TB exposure and progression to disease, including the risk of TB recurrence, in this vulnerable population. There is urgency to scale up interventions such as intensified TB case-finding, isoniazid preventive therapy, and TB infection control, as well as HIV testing and improved access to antiretroviral treatment. Increased awareness and concerted action is required to reduce TB/HIV co-infection rates in the Asia-Pacific region and to improve the outcomes of people living with HIV.

Identification of pre-leukaemic haematopoietic stem cells in acute leukaemia

In acute myeloid leukemia (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-leukemic phase. Here, highly purified hematopoietic stem cells (HSC), 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 HSC exhibited multilineage repopulation advantage over non-mutated HSC in xenografts, establishing their identity as pre-leukemic-HSC (preL-HSC). preL-HSC were found in remission samples indicating that they survive chemotherapy. Thus DNMT3a^mut arises early in AML evolution, likely in HSC, leading to a clonally expanded pool of preL-HSC from which AML evolves. Our findings provide a paradigm for the detection and treatment of pre-leukemic clones before the acquisition of additional genetic lesions engenders greater therapeutic resistance.

Identification of genes expressed by immune cells of the colon that are regulated by colorectal cancer-associated variants

A locus on human chromosome 11q23 tagged by marker rs3802842 was associated with colorectal cancer (CRC) in a genome-wide association study; this finding has been replicated in case–control studies worldwide. In order to identify biologic factors at this locus that are related to the etiopathology of CRC, we used microarray-based target selection methods, coupled to next-generation sequencing, to study 103 kb at the 11q23 locus. We genotyped 369 putative variants from 1,030 patients with CRC (cases) and 1,061 individuals without CRC (controls) from the Ontario Familial Colorectal Cancer Registry. Two previously uncharacterized genes, COLCA1 and COLCA2, were found to be co-regulated genes that are transcribed from opposite strands. Expression levels of COLCA1 and COLCA2 transcripts correlate with rs3802842 genotypes. In colon tissues, COLCA1 co-localizes with crystalloid granules of eosinophils and granular organelles of mast cells, neutrophils, macrophages, dendritic cells and differentiated myeloid-derived cell lines. COLCA2 is present in the cytoplasm of normal epithelial, immune and other cell lineages, as well as tumor cells. Tissue microarray analysis demonstrates the association of rs3802842 with lymphocyte density in the lamina propria (p = 0.014) and levels of COLCA1 in the lamina propria (p = 0.00016) and COLCA2 (tumor cells, p = 0.0041 and lamina propria, p = 6 × 10^–5). In conclusion, genetic, expression and immunohistochemical data implicate COLCA1 and COLCA2 in the pathogenesis of colon cancer. Histologic analyses indicate the involvement of immune pathways.

Cloud-based uniform ChIP-Seq processing tools for modENCODE and ENCODE

Background

Funded by the National Institutes of Health (NIH), the aim of the Model Organism ENCyclopedia of DNA Elements (modENCODE) project is to provide the biological research community with a comprehensive encyclopedia of functional genomic elements for both model organisms C. elegans (worm) and D. melanogaster (fly). With a total size of just under 10 terabytes of data collected and released to the public, one of the challenges faced by researchers is to extract biologically meaningful knowledge from this large data set. While the basic quality control, pre-processing, and analysis of the data has already been performed by members of the modENCODE consortium, many researchers will wish to reinterpret the data set using modifications and enhancements of the original protocols, or combine modENCODE data with other data sets. Unfortunately this can be a time consuming and logistically challenging proposition.

Results

In recognition of this challenge, the modENCODE DCC has released uniform computing resources for analyzing modENCODE data on Galaxy (https://github.com/modENCODE-DCC/Galaxy), on the public Amazon Cloud (http://aws.amazon.com), and on the private Bionimbus Cloud for genomic research (http://www.bionimbus.org). In particular, we have released Galaxy workflows for interpreting ChIP-seq data which use the same quality control (QC) and peak calling standards adopted by the modENCODE and ENCODE communities. For convenience of use, we have created Amazon and Bionimbus Cloud machine images containing Galaxy along with all the modENCODE data, software and other dependencies.

Conclusions

Using these resources provides a framework for running consistent and reproducible analyses on modENCODE data, ultimately allowing researchers to use more of their time using modENCODE data, and less time moving it around.

Exome sequencing identifies nonsegregating nonsense ATM and PALB2 variants in familial pancreatic cancer

We sequenced 11 germline exomes from five families with familial pancreatic cancer (FPC). One proband had a germline nonsense variant in ATM with somatic loss of the variant allele. Another proband had a nonsense variant in PALB2 with somatic loss of the variant allele. Both variants were absent in a relative with FPC. These findings question the causal mechanisms of ATM and PALB2 in these families and highlight challenges in identifying the causes of familial cancer syndromes using exome sequencing.

Spatiotemporal integration of molecular and anatomical data in virtual reality using semantic mapping

We have developed a computational framework for spatiotemporal integration of molecular and anatomical datasets in a virtual reality environment. Using two case studies involving gene expression data and pharmacokinetic data, respectively, we demonstrate how existing knowledge bases for molecular data can be semantically mapped onto a standardized anatomical context of human body. Our data mapping methodology uses ontological representations of heterogeneous biomedical datasets and an ontology reasoner to create complex semantic descriptions of biomedical processes. This framework provides a means to systematically combine an increasing amount of biomedical imaging and numerical data into spatiotemporally coherent graphical representations. Our work enables medical researchers with different expertise to simulate complex phenomena visually and to develop insights through the use of shared data, thus paving the way for pathological inference, developmental pattern discovery and biomedical hypothesis testing.