Search

Whereas traditional histology and light microscopy require multiple steps of formalin fixation, paraffin embedding, and sectioning to generate images for pathologic diagnosis, Microscopy using Ultraviolet Surface Excitation (MUSE) operates through UV excitation on the cut surface of tissue, generating images of high resolution without the need to fix or section tissue and allowing for potential use for downstream molecular tests. Here, we present the first study of the use and suitability of MUSE microscopy for neuropathological samples. MUSE images were generated from surgical biopsy samples of primary and metastatic brain tumor biopsy samples (n = 27), and blinded assessments of diagnoses, tumor grades, and cellular features were compared to corresponding hematoxylin and eosin (H&E) images. A set of MUSE-treated samples subsequently underwent exome and targeted sequencing, and quality metrics were compared to those from fresh frozen specimens. Diagnostic accuracy was relatively high, and DNA and RNA integrity appeared to be preserved for this cohort. This suggests that MUSE may be a reliable method of generating high-quality diagnostic-grade histologic images for neuropathology on a rapid and sample-sparing basis and for subsequent molecular analysis of DNA and RNA.

In BCR-ABL1 lymphoblastic leukemia, treatment heterogeneity to tyrosine kinase inhibitors (TKIs), especially in the absence of kinase domain mutations in BCR-ABL1, is poorly understood. Through deep molecular profiling, we uncovered three transcriptomic subtypes of BCR-ABL1 lymphoblastic leukemia, each representing a maturation arrest at a stage of B-cell progenitor differentiation. An earlier arrest was associated with lineage promiscuity, treatment refractoriness and poor patient outcomes. A later arrest was associated with lineage fidelity, durable leukemia remissions and improved patient outcomes. Each maturation arrest was marked by specific genomic events that control different transition points in B-cell development. Interestingly, these events were absent in BCR-ABL1+ preleukemic stem cells isolated from patients regardless of subtype, which supports that transcriptomic phenotypes are determined downstream of the leukemia-initialing event. Overall, our data indicate that treatment response and TKI efficacy are unexpected outcomes of the differentiation stage at which this leukemia transforms.

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.

BackgroundThere 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.MethodsUsing 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.ResultsAmong 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.ConclusionsAlthough 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.

Early diagnosis of melanoma is critical for improved survival. However, the biomarkers of early melanoma evolution and their origin within the tumor and its microenvironment, including the keratinocytes, are poorly defined. To address this, we used spatial transcript profiling that maintains the morphological tumor context to measure the expression of >1,000 RNAs in situ in patient-derived formalin-fixed, paraffin-embedded tissue sections in primary melanoma and melanocytic nevi. We profiled 134 regions of interest (each 200 μm in diameter) enriched in melanocytes, neighboring keratinocytes, or immune cells. This approach captured distinct expression patterns across cell types and tumor types during melanoma development. Unexpectedly, we discovered that S100A8 is expressed by keratinocytes within the tumor microenvironment during melanoma growth. Immunohistochemistry of 252 tumors showed prominent keratinocyte-derived S100A8 expression in melanoma but not in benign tumors and confirmed the same pattern for S100A8's binding partner S100A9, suggesting that injury to the epidermis may be an early and readily detectable indicator of melanoma development. Together, our results establish a framework for high-plex, spatial, and cell type‒specific resolution of gene expression in archival tissue applicable to the development of biomarkers and characterization of tumor microenvironment interactions in tumor evolution.

Aberrant expression of the E26 transformation-specific (ETS) transcription factors characterizes numerous human malignancies. Many of these proteins, including EWS:FLI1 and EWS:ERG fusions in Ewing sarcoma (EwS) and TMPRSS2:ERG in prostate cancer (PCa), drive oncogenic programs via binding to GGAA repeats. We report here that both EWS:FLI1 and ERG bind and transcriptionally activate GGAA-rich pericentromeric heterochromatin. The respective pathogen-like HSAT2 and HSAT3 RNAs, together with LINE, SINE, ERV, and other repeat transcripts, are expressed in EwS and PCa tumors, secreted in extracellular vesicles (EVs), and are highly elevated in plasma of patients with EwS with metastatic disease. High human satellite 2 and 3 (HSAT2,3) levels in EWS:FLI1- or ERG- expressing cells and tumors were associated with induction of G2/M checkpoint, mitotic spindle, and DNA damage programs. These programs were also activated in EwS EV-treated fibroblasts, coincident with accumulation of HSAT2,3 RNAs, proinflammatory responses, mitotic defects, and senescence. Mechanistically, HSAT2,3-enriched cancer EVs induced cGAS-TBK1 innate immune signaling and formation of cytosolic granules positive for double-strand RNAs, RNA-DNA, and cGAS. Hence, aberrantly expressed ETS proteins derepress pericentromeric heterochromatin, yielding pathogenic RNAs that transmit genotoxic stress and inflammation to local and distant sites. Monitoring HSAT2,3 plasma levels and preventing their dissemination may thus improve therapeutic strategies and blood-based diagnostics., Introduction Ewing sarcoma (EwS) is a poorly differentiated and highly aggressive childhood and adolescent malignancy of bone and soft tissues, thought to originate from primitive bone marrow-derived mesenchymal stem cells [...]

Cutaneous melanoma results from the malignant transformation of melanocytes in skin and accounts for 75% of deaths related to skin cancer. Between 5 to 12% of melanoma cases can be attributed to hereditary melanoma, melanoma caused by inherited germline mutations in melanoma predisposition genes. Although several inherited melanoma predisposition syndromes have been identified, a subset of melanoma families lack pathogenic mutations in known highly penetrant predisposition genes, including CDKN2A, CDK4, and BAP1. Here we report the case of a woman with a history of melanoma and multiple primary tumors with one pathogenic germline mutation in BUB1B

Early diagnosis of melanoma is critical for improved survival as melanoma is the deadliest of the common forms of skin cancer. The gold standard for the diagnosis of melanoma is a biopsy followed by histopathological analysis. Melanocytic nevi, which are very common benign neoplasms of the melanocytes, are often biopsied because they are mimics and precursor for melanoma. Annually, 4.5 million pigmented lesions are biopsied in the United States. A subset of melanoma is difficult to distinguish from melanocytic nevi, resulting in diagnostic errors and worsened patient outcomes. Therefore, improved diagnostic tests, including the utilization of novel biomarkers, are being developed to improve clinical and histological diagnostic accuracy of melanoma. Moreover, non-invasive tests prior to surgical biopsy have been introduced, including an adhesive patch-based assay that tests the expression of melanoma biomarkers PRAME and noncoding long RNA LINC00518. Our prior work identified that S100A8, a member of the calcium-binding S100 family, is differentially expressed in melanomas versus nevi4 . Specifically, S100A8 is expressed by the keratinocyte microenvironment of melanomas but not nevi. S100A8 is a melanoma biomarker of interest for an adhesive patch-based assay, because it is expressed in the epidermis, the most superficial layer of the skin.

In leukemia, a distinct subpopulation of cancer-initiating cells called leukemia stem cells (LSCs) is believed to drive population expansion and tumor growth. Failing to eliminate LSCs may result in disease relapse regardless of the amount of non-LSCs destroyed. The first step in targeting and eliminating LSCs is to identify and characterize them. Acute precursor B lymphoblastic leukemia (B-ALL) cells derived from patients were incubated with fluorescent glucose analog 2-(N-(7-Nitrobenz-2-oxa-1, 3-diazol-4-yl) Amino)-2-Deoxyglucose (NBDG) and sorted based on NBDG uptake. Cell subpopulations defined by glucose uptake were then serially transplanted into mice and evaluated for leukemia initiating capacity. Gene expression profiles of these cells were characterized using RNA-Sequencing (RNA-Seq). A distinct population of NBDG-low cells was identified in patient B-ALL samples. These cells are a small population (1.92% of the entire leukemia population), have lower HLA expression, and are smaller in size (4.0 to 7.0 μm) than the rest of the leukemia population. All mice transplanted with NBDG-low cells developed leukemia between 5 and 14 weeks, while those transplanted with NBDG-high cells did not develop leukemia (p ≤ 0.0001-0.002). Serial transplantation of the NBDG-low mouse model resulted in successful leukemia development. NBDG-medium (NBDG-med) populations also developed leukemia. Interestingly, comprehensive molecular characterization of NBDG-low and NBDG-med cells from patient-derived xenograft (PDX) models using RNA-Seq revealed a distinct profile of 2,162 differentially-expressed transcripts (DETs) (p

Abstract: Introduction: While previous studies propose pesticide exposure to be a risk factor for non-Hodgkin lymphoma (NHL) development, results are inconclusive. In addition, the impact of pesticide exposure on NHL survival is not well-established. Therefore, we identified NHL patients from the California Cancer Registry and linked these patients with the statewide pesticide use reporting database to determine the impact of pesticide exposure on NHL-related incidence and outcomes. Methods: Using the California Cancer Registry, we identified patients with a first primary diagnosis of NHL from 2010-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-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. SEER*Stat software was used to calculate NHL subtype incidence rates by census tract pesticide use level. Multivariable cox proportional hazards regression models were used to evaluate the impact of total pesticide exposure from CalEnviroScreen 3.0 and individual pesticide exposure from geographic land use data on lymphoma-specific and overall survival. Results: Among 35,808 NHL patients identified, 44.2% were exposed to pesticide in their census tract of residence. Pesticide exposure was higher in Hispanic/Latino (46.5%) and non-Hispanic white (45.6%) then Asian/Pacific Islander (37.2%) and African American (34.9%) patients with NHL. Glyphosate, organophosphorus, carbamate, phenoxyherbicide and 2,4-dimethylamine salt exposure was reported in 34.1%, 26.0%, 10.6%, 14.0% and 12.8% of NHL patients, respectively. Pesticide exposure was not associated with increased NHL incidence by NHL subtype or subgroups defined by sociodemographic factors. Total pesticide exposure at 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. Conclusion: In this large population-based study of neighborhood agricultural pesticide exposure, pesticide exposure was noted to be prevalent among patients diagnosed with NHL, with high pesticide exposure particularly observed in Hispanics/Latinos and non-Hispanic whites. However, pesticide exposure was not consistently associated with increased NHL incidence or worse NHL lymphoma-specific or overall survival. Disclosures: Poh: Acrotech: Honoraria; Incyte: Research Funding; Morphosys: Consultancy. Tuscano: BMS: Research Funding; Seattle Genetics: Research Funding; Takeda: Research Funding; Acrotech: Research Funding; Genentech: Research Funding; Pharmacyclics: Research Funding; Abbvie: Research Funding.

Colorectal cancer (CRC) is a biologically heterogeneous disease. To characterize its mutational profile, we conduct targeted sequencing of 205 genes for 2,105 CRC cases with survival data. Our data shows several findings in addition to enhancing the existing knowledge of CRC. We identify PRKCI, SPZ1, MUTYH, MAP2K4, FETUB, and TGFBR2 as additional genes significantly mutated in CRC. We find that among hypermutated tumors, an increased mutation burden is associated with improved CRC-specific survival (HR = 0.42, 95% CI: 0.21-0.82). Mutations in TP53 are associated with poorer CRC-specific survival, which is most pronounced in cases carrying TP53 mutations with predicted 0% transcriptional activity (HR = 1.53, 95% CI: 1.21-1.94). Furthermore, we observe differences in mutational frequency of several genes and pathways by tumor location, stage, and sex. Overall, this large study provides deep insights into somatic mutations in CRC, and their potential relationships with survival and tumor features.

Although DNA and RNA sequencing has a history spanning five decades, large-scale massively parallel sequencing, or next-generation sequencing (NGS), has only been commercially available for about 10 years. Nonetheless, the meteoric increase in sequencing throughput with NGS has dramatically changed our understanding of our genome and ourselves. Sequencing the first human genome as a haploid reference took nearly 10 years but now a full diploid human genome sequence can be accomplished in just a few days. NGS has also reduced the cost of generating sequence data and a plethora of sequence-based methods for probing a genome have emerged using NGS as the readout and have been applied to many species. NGS methods have also entered the medical realm and will see an increasing use in diagnosis and treatment. NGS has largely been driven by short-read generation (150 bp) but new platforms have emerged and are now capable of generating long multikilobase reads. These latter platforms enable reference-independent genome assemblies and long-range haplotype generation. Rapid DNA and RNA sequencing is now mainstream and will continue to have an increasing impact on biology and medicine.

Oncogenesis is driven by germline, environmental and stochastic factors. It is unknown how these interact to produce the molecular phenotypes of tumors. We therefore quantified the influence of germline polymorphisms on the somatic epigenome of 589 localized prostate tumors. Predisposition risk loci influence a tumor's epigenome, uncovering a mechanism for cancer susceptibility. We identified and validated 1,178 loci associated with altered methylation in tumoral but not nonmalignant tissue. These tumor methylation quantitative trait loci influence chromatin structure, as well as RNA and protein abundance. One prominent tumor methylation quantitative trait locus is associated with AKT1 expression and is predictive of relapse after definitive local therapy in both discovery and validation cohorts. These data reveal intricate crosstalk between the germ line and the epigenome of primary tumors, which may help identify germline biomarkers of aggressive disease to aid patient triage and optimize the use of more invasive or expensive diagnostic assays.

Since the first draft of the human genome was completed, next-generation DNA sequencing technology has dramatically reduced the cost of sequencing a genome. Computational analysis has not advanced as fast as the instruments that generate the data, and storing all the data remains a challenge. Nevertheless, personal genomics has arrived and is already being used in the clinic. Significant privacy issues remain, however, and these are not widely understood. The Genetic Information Non-Discrimination Act (GINA) needs to be extended and the probabilistic nature of genetic predisposition must be better explained to both the public and physicians. We must also be wary that this promising new technology and its applications do not amplify existing healthcare disparities.

Genomic rearrangements are a hallmark of cancer biology and progression, allowing cells to rapidly transform through alterations in regulatory structures, changes in expression patterns, reprogramming of signaling pathways, and creation of novel transcripts via gene fusion events. Though functional gene fusions encoding oncogenic proteins are the most dramatic outcomes of genomic rearrangements, we investigated the relationship between rearrangements evidenced by fusion transcripts and local expression changes in cancer using transcriptome data alone. 9,953 gene fusion predictions from 418 primary serious ovarian cancer tumors were analyzed, identifying depletions of gene fusion breakpoints within coding regions of fused genes as well as an N-terminal enrichment of breakpoints within fused genes. We identified 48 genes with significant fusion-associated upregulation and furthermore demonstrate that significant regional overexpression of intact genes in patient transcriptomes occurs within 1 megabase of 78 novel gene fusions that function as central markers of these regions. We reveal that cancer transcriptomes select for gene fusions that preserve protein and protein domain coding potential. The association of gene fusion transcripts with neighboring gene overexpression supports rearrangements as mechanism through which cancer cells remodel their transcriptomes and identifies a new way to utilize gene fusions as indicators of regional expression changes in diseased cells with only transcriptomic data.

IntroductionEGFR exon 20 insertions (EGFRex20ins) comprise an uncommon subset of EGFR-activating alterations relatively insensitive to first- and second-generation EGFR tyrosine kinase inhibitors (TKIs). However, recent early clinical data suggests these patients may benefit from newer-generation EGFR-TKIs. Comprehensive genomic profiling (CGP) identifies a broad spectrum of EGFRex20ins and associated co-occurring genomic alterations (GAs) present in NSCLC.MethodsHybrid capture-based CGP was performed prospectively on 14,483 clinically annotated consecutive NSCLC specimens to a mean coverage depth of greater than 650X for 236 or 315 cancer-related genes.ResultsOf 14,483 NSCLC cases, CGP identified 263 (1.8%) cases with EGFRex20ins, representing 12% (263 of 2251) of cases with EGFR mutations. Sixty-four unique EGFRex20ins were identified, most commonly D770_N771>ASVDN (21%) and N771_P772>SVDNP (20%). EGFR amplification occurred in 22% (57 of 263). The most common co-occurring GAs effected tumor protein p53 (TP53) (56%), cyclin dependent kinase inhibitor 2A (CDKN2A) (22%), cyclin dependent kinase inhibitor 2B (CDKN2B) (16%), NK2 homeobox 1 (NKX2-1) (14%) and RB transcriptional corepressor 1 (RB1) (11%); co-occurring GAs in other known lung cancer drivers were rare (5%). Average tumor mutational burden was low (mean 4.3, range 0 to 40.3 mutations/Mb). Clinical outcomes to first- and second-generation EGFR TKIs were obtained for five patients and none responded.ConclusionsIn the largest series of EGFRex20ins NSCLC, diverse EGFRex20ins were detected in 12% of EGFR-mutant NSCLC, a higher frequency than previously reported in smaller single-institution studies. Clinical outcomes showed lack of response to EGFR TKIs. Tumor mutational burden was low, consistent with non-smoking associated NSCLC. Comprehensive sequencing revealed increased proportion and wide variety of EGFRex20ins, representing a population of patients significant enough for focused efforts on effective interventions.

BackgroundA 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.ObjectiveTo identify histomorphologic features correlating with BRAF V600E status in nevi.MethodsWe retrospectively identified melanocytic nevi from our laboratory reporting system. We performed a histomorphologic analysis and analysis of BRAF V600E expression by immunohistochemistry.ResultsThirteen 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.LimitationsLimitations include the study's retrospective design and the small sample size of nevi negative for BRAF V600E.ConclusionsBRAF V600E is associated with distinct histomorphologic features in nevi. These features may contribute to improving the accuracy of classification and diagnosis of melanocytic neoplasms.

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.

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.

BackgroundA 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.ResultsIn 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).ConclusionsIn 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 .

Nuclear mutations are well known to drive tumor incidence, aggression and response to therapy. By contrast, the frequency and roles of mutations in the maternally inherited mitochondrial genome are poorly understood. Here we sequence the mitochondrial genomes of 384 localized prostate cancer patients, and identify a median of one mitochondrial single-nucleotide variant (mtSNV) per patient. Some of these mtSNVs occur in recurrent mutational hotspots and associate with aggressive disease. Younger patients have fewer mtSNVs than those who diagnosed at an older age. We demonstrate strong links between mitochondrial and nuclear mutational profiles, with co-occurrence between specific mutations. For example, certain control region mtSNVs co-occur with gain of the MYC oncogene, and these mutations are jointly associated with patient survival. These data demonstrate frequent mitochondrial mutation in prostate cancer, and suggest interplay between nuclear and mitochondrial mutational profiles in prostate cancer.In prostate cancer, the role of mutations in the maternally-inherited mitochondrial genome are not well known. Here, the authors demonstrate frequent, age-dependent mitochondrial mutation in prostate cancer. Strong links between mitochondrial and nuclear mutational profiles are associated with clinical aggressivity.

The contributions of coding mutations to tumorigenesis are relatively well known; however, little is known about somatic alterations in noncoding DNA. Here we describe GECCO (Genomic Enrichment Computational Clustering Operation) to analyze somatic noncoding alterations in 308 pancreatic ductal adenocarcinomas (PDAs) and identify commonly mutated regulatory regions. We find recurrent noncoding mutations to be enriched in PDA pathways, including axon guidance and cell adhesion, and newly identified processes, including transcription and homeobox genes. We identified mutations in protein binding sites correlating with differential expression of proximal genes and experimentally validated effects of mutations on expression. We developed an expression modulation score that quantifies the strength of gene regulation imposed by each class of regulatory elements, and found the strongest elements were most frequently mutated, suggesting a selective advantage. Our detailed single-cancer analysis of noncoding alterations identifies regulatory mutations as candidates for diagnostic and prognostic markers, and suggests new mechanisms for tumor evolution.

Germline mutations in the BRCA2 tumour suppressor are associated with both an increased lifetime risk of developing prostate cancer (PCa) and increased risk of aggressive disease. To understand this aggression, here we profile the genomes and methylomes of localized PCa from 14 carriers of deleterious germline BRCA2 mutations (BRCA2-mutant PCa). We show that BRCA2-mutant PCa harbour increased genomic instability and a mutational profile that more closely resembles metastastic than localized disease. BRCA2-mutant PCa shows genomic and epigenomic dysregulation of the MED12L/MED12 axis, which is frequently dysregulated in metastatic castration-resistant prostate cancer (mCRPC). This dysregulation is enriched in BRCA2-mutant PCa harbouring intraductal carcinoma (IDC). Microdissection and sequencing of IDC and juxtaposed adjacent non-IDC invasive carcinoma in 10 patients demonstrates a common ancestor to both histopathologies. Overall we show that localized castration-sensitive BRCA2-mutant tumours are uniquely aggressive, due to de novo aberration in genes usually associated with metastatic disease, justifying aggressive initial treatment.

Histologic examination of tissues is central to the diagnosis and management of neoplasms and many other diseases, and is a foundational technique for preclinical and basic research. However, commonly used bright-field microscopy requires prior preparation of micrometre-thick tissue sections mounted on glass slides, a process that can require hours or days, that contributes to cost, and that delays access to critical information. Here, we introduce a simple, non-destructive slide-free technique that within minutes provides high-resolution diagnostic histological images resembling those obtained from conventional haematoxylin-and-eosin-histology. The approach, which we named microscopy with ultraviolet surface excitation (MUSE), can also generate shape and colour-contrast information. MUSE relies on ~280-nm ultraviolet light to restrict the excitation of conventional fluorescent stains to tissue surfaces, and it has no significant effects on downstream molecular assays (including fluorescence in situ hybridization and RNA-seq). MUSE promises to improve the speed and efficiency of patient care in both state-of-the-art and low-resource settings, and to provide opportunities for rapid histology in research.

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.

Chromosomal rearrangements are a hallmark of acute lymphoblastic leukemia (ALL) and are important ALL initiating events. We describe four different rearrangements of the erythropoietin receptor gene EPOR in Philadelphia chromosome-like (Ph-like) ALL. All of these rearrangements result in truncation of the cytoplasmic tail of EPOR at residues similar to those mutated in primary familial congenital polycythemia, with preservation of the proximal tyrosine essential for receptor activation and loss of distal regulatory residues. This resulted in deregulated EPOR expression, hypersensitivity to erythropoietin stimulation, and heightened JAK-STAT activation. Expression of truncated EPOR in mouse B cell progenitors induced ALL in vivo. Human leukemic cells with EPOR rearrangements were sensitive to JAK-STAT inhibition, suggesting a therapeutic option in high-risk ALL.

In a classical view of hematopoiesis, the various blood cell lineages arise via a hierarchical scheme starting with multipotent stem cells that become increasingly restricted in their differentiation potential through oligopotent and then unipotent progenitors. We developed a cell-sorting scheme to resolve myeloid (My), erythroid (Er), and megakaryocytic (Mk) fates from single CD34(+) cells and then mapped the progenitor hierarchy across human development. Fetal liver contained large numbers of distinct oligopotent progenitors with intermingled My, Er, and Mk fates. However, few oligopotent progenitor intermediates were present in the adult bone marrow. Instead, only two progenitor classes predominate, multipotent and unipotent, with Er-Mk lineages emerging from multipotent cells. The developmental shift to an adult "two-tier" hierarchy challenges current dogma and provides a revised framework to understand normal and disease states of human hematopoiesis.

Chromosome 5p15.33 has been identified as a lung cancer susceptibility locus, however the underlying causal mechanisms were not fully elucidated. Previous fine-mapping studies of this locus have relied on imputation or investigated a small number of known, common variants. This study represents a significant advance over previous research by investigating a large number of novel, rare variants, as well as their underlying mechanisms through telomere length. Variants for this fine-mapping study were identified through a targeted deep sequencing (average depth of coverage greater than 4000×) of 576 individuals. Subsequently, 4652 SNPs, including 1108 novel SNPs, were genotyped in 5164 cases and 5716 controls of European ancestry. After adjusting for known risk loci, rs2736100 and rs401681, we identified a new, independent lung cancer susceptibility variant in LPCAT1: rs139852726 (OR = 0.46, P = 4.73×10(-9)), and three new adenocarcinoma risk variants in TERT: rs61748181 (OR = 0.53, P = 2.64×10(-6)), rs112290073 (OR = 1.85, P = 1.27×10(-5)), rs138895564 (OR = 2.16, P = 2.06×10(-5); among young cases, OR = 3.77, P = 8.41×10(-4)). In addition, we found that rs139852726 (P = 1.44×10(-3)) was associated with telomere length in a sample of 922 healthy individuals. The gene-based SKAT-O analysis implicated TERT as the most relevant gene in the 5p15.33 region for adenocarcinoma (P = 7.84×10(-7)) and lung cancer (P = 2.37×10(-5)) risk. In this largest fine-mapping study to investigate a large number of rare and novel variants within 5p15.33, we identified novel lung and adenocarcinoma susceptibility loci with large effects and provided support for the role of telomere length as the potential underlying mechanism.

BackgroundAccurate detection of somatic single nucleotide variants and small insertions and deletions from DNA sequencing experiments of tumour-normal pairs is a challenging task. Tumour samples are often contaminated with normal cells confounding the available evidence for the somatic variants. Furthermore, tumours are heterogeneous so sub-clonal variants are observed at reduced allele frequencies. We present here a cell-line titration series dataset that can be used to evaluate somatic variant calling pipelines with the goal of reliably calling true somatic mutations at low allele frequencies.ResultsCell-line DNA was mixed with matched normal DNA at 8 different ratios to generate samples with known tumour cellularities, and exome sequenced on Illumina HiSeq to depths of >300×. The data was processed with several different variant calling pipelines and verification experiments were performed to assay >1500 somatic variant candidates using Ion Torrent PGM as an orthogonal technology. By examining the variants called at varying cellularities and depths of coverage, we show that the best performing pipelines are able to maintain a high level of precision at any cellularity. In addition, we estimate the number of true somatic variants undetected as cellularity and coverage decrease.ConclusionsOur cell-line titration series dataset, along with the associated verification results, was effective for this evaluation and will serve as a valuable dataset for future somatic calling algorithm development. The data is available for further analysis at the European Genome-phenome Archive under accession number EGAS00001001016. Data access requires registration through the International Cancer Genome Consortium's Data Access Compliance Office (ICGC DACO).

As whole-genome sequencing for cancer genome analysis becomes a clinical tool, a full understanding of the variables affecting sequencing analysis output is required. Here using tumour-normal sample pairs from two different types of cancer, chronic lymphocytic leukaemia and medulloblastoma, we conduct a benchmarking exercise within the context of the International Cancer Genome Consortium. We compare sequencing methods, analysis pipelines and validation methods. We show that using PCR-free methods and increasing sequencing depth to ∼ 100 × shows benefits, as long as the tumour:control coverage ratio remains balanced. We observe widely varying mutation call rates and low concordance among analysis pipelines, reflecting the artefact-prone nature of the raw data and lack of standards for dealing with the artefacts. However, we show that, using the benchmark mutation set we have created, many issues are in fact easy to remedy and have an immediate positive impact on mutation detection accuracy.

BackgroundDormant leukemia stem cells (LSC) promote therapeutic resistance and leukemic progression as a result of unbridled activation of stem cell gene expression programs. Thus, we hypothesized that 1) deregulation of the hedgehog (Hh) stem cell self-renewal and cell cycle regulatory pathway would promote dormant human LSC generation and 2) that PF-04449913, a clinical antagonist of the GLI2 transcriptional activator, smoothened (SMO), would enhance dormant human LSC eradication.MethodsTo test these postulates, whole transcriptome RNA sequencing (RNA-seq), microarray, qRT-PCR, stromal co-culture, confocal fluorescence microscopic, nanoproteomic, serial transplantation and cell cycle analyses were performed on FACS purified normal, chronic phase (CP) chronic myeloid leukemia (CML), blast crisis (BC) phase CML progenitors with or without PF-04449913 treatment.ResultsNotably, RNA-seq analyses revealed that Hh pathway and cell cycle regulatory gene overexpression correlated with leukemic progression. While lentivirally enforced GLI2 expression enhanced leukemic progenitor dormancy in stromal co-cultures, this was not observed with a mutant GLI2 lacking a transactivation domain, suggesting that GLI2 expression prevented cell cycle transit. Selective SMO inhibition with PF-04449913 in humanized stromal co-cultures and LSC xenografts reduced downstream GLI2 protein and cell cycle regulatory gene expression. Moreover, SMO inhibition enhanced cell cycle transit and sensitized BC LSC to tyrosine kinase inhibition in vivo at doses that spare normal HSC.ConclusionIn summary, while GLI2, forms part of a core HH pathway transcriptional regulatory network that promotes human myeloid leukemic progression and dormant LSC generation, selective inhibition with PF-04449913 reduces the dormant LSC burden thereby providing a strong rationale for clinical trials predicated on SMO inhibition in combination with TKIs or chemotherapeutic agents with the ultimate aim of obviating leukemic therapeutic resistance, persistence and progression.

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.

BackgroundChromothripsis, a newly discovered type of complex genomic rearrangement, has been implicated in the evolution of several types of cancers. To date, it has been described in bone cancer, SHH-medulloblastoma and acute myeloid leukemia, amongst others, however there are still no formal or automated methods for detecting or annotating it in high throughput sequencing data. As such, findings of chromothripsis are difficult to compare and many cases likely escape detection altogether.ResultsWe introduce ShatterProof, a software tool for detecting and quantifying chromothriptic events. ShatterProof takes structural variation calls (translocations, copy-number variations, short insertions and loss of heterozygosity) produced by any algorithm and using an operational definition of chromothripsis performs robust statistical tests to accurately predict the presence and location of chromothriptic events. Validation of our tool was conducted using clinical data sets including matched normal, prostate cancer samples in addition to the colorectal cancer and SCLC data sets used in the original description of chromothripsis.ConclusionsShatterProof is computationally efficient, having low memory requirements and near linear computation time. This allows it to become a standard component of sequencing analysis pipelines, enabling researchers to routinely and accurately assess samples for chromothripsis. Source code and documentation can be found at http://search.cpan.org/~sgovind/Shatterproof.

Polymorphous low-grade adenocarcinoma (PLGA) is the second most frequent type of malignant tumor of the minor salivary glands. We identified PRKD1 hotspot mutations encoding p.Glu710Asp in 72.9% of PLGAs but not in other salivary gland tumors. Functional studies demonstrated that this kinase-activating alteration likely constitutes a driver of PLGA.

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.

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.

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.

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.

Somatic cells can be reprogrammed to induced pluripotent stem cells (iPSCs) using oncogenic transcription factors. However, this method leads to genetic aberrations in iPSCs via unknown mechanisms, which may limit their clinical use. Here, we demonstrate that the supplementation of growth media with antioxidants reduces the genome instability of cells transduced with the reprogramming factors. Antioxidant supplementation did not affect transgene expression level or silencing kinetics. Importantly, iPSCs made with antioxidants had significantly fewer de novo copy number variations, but not fewer coding point mutations, than iPSCs made without antioxidants. Our results suggest that the quality and safety of human iPSCs might be enhanced by using antioxidants in the growth media during the generation and maintenance of iPSCs.

We describe a method for pooling and sequencing DNA from a large number of individual samples while preserving information regarding sample identity. DNA from 576 individuals was arranged into four 12 row by 12 column matrices and then pooled by row and by column resulting in 96 total pools with 12 individuals in each pool. Pooling of DNA was carried out in a two-dimensional fashion, such that DNA from each individual is present in exactly one row pool and exactly one column pool. By considering the variants observed in the rows and columns of a matrix we are able to trace rare variants back to the specific individuals that carry them. The pooled DNA samples were enriched over a 250 kb region previously identified by GWAS to significantly predispose individuals to lung cancer. All 96 pools (12 row and 12 column pools from 4 matrices) were barcoded and sequenced on an Illumina HiSeq 2000 instrument with an average depth of coverage greater than 4,000×. Verification based on Ion PGM sequencing confirmed the presence of 91.4% of confidently classified SNVs assayed. In this way, each individual sample is sequenced in multiple pools providing more accurate variant calling than a single pool or a multiplexed approach. This provides a powerful method for rare variant detection in regions of interest at a reduced cost to the researcher.

Short tandem repeat (STR) analysis, such as the AmpFlSTR® Identifiler® Plus kit, is a standard, PCR-based human genotyping method used in the field of forensics. Misidentification of cell line and tissue DNA can be costly if not detected early; therefore it is necessary to have quality control measures such as STR profiling in place. A major issue in large-scale research studies involving archival formalin-fixed paraffin embedded (FFPE) tissues is that varying levels of DNA degradation can result in failure to correctly identify samples using STR genotyping. PCR amplification of STRs of several hundred base pairs is not always possible when DNA is degraded. The Sample ID Plus® panel from Sequenom allows for human DNA identification and authentication using SNP genotyping. In comparison to lengthy STR amplicons, this multiplexing PCR assay requires amplification of only 76-139 base pairs, and utilizes 47 SNPs to discriminate between individual samples. In this study, we evaluated both STR and SNP genotyping methods of sample identification, with a focus on paired FFPE tumor/normal DNA samples intended for next-generation sequencing (NGS). The ability to successfully validate the identity of FFPE samples can enable cost savings by reducing rework.

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.

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.

The International Cancer Genome Consortium (ICGC) was launched to coordinate large-scale cancer genome studies in tumours from 50 different cancer types and/or subtypes that are of clinical and societal importance across the globe. Systematic studies of more than 25,000 cancer genomes at the genomic, epigenomic and transcriptomic levels will reveal the repertoire of oncogenic mutations, uncover traces of the mutagenic influences, define clinically relevant subtypes for prognosis and therapeutic management, and enable the development of new cancer therapies.

The International Cancer Genome Consortium (ICGC) was launched to coordinate large-scale cancer genome studies in tumours from 50 different cancer types and/or subtypes that are of clinical and societal importance across the globe. Systematic studies of more than 25,000 cancer genomes at the genomic, epigenomic and transcriptomic levels will reveal the repertoire of oncogenic mutations, uncover traces of the mutagenic influences, define clinically relevant subtypes for prognosis and therapeutic management, and enable the development of new cancer therapies.

Metagenomic studies characterize both the composition and diversity of uncultured viral and microbial communities. BLAST-based comparisons have typically been used for such analyses; however, sampling biases, high percentages of unknown sequences, and the use of arbitrary thresholds to find significant similarities can decrease the accuracy and validity of estimates. Here, we present Genome relative Abundance and Average Size (GAAS), a complete software package that provides improved estimates of community composition and average genome length for metagenomes in both textual and graphical formats. GAAS implements a novel methodology to control for sampling bias via length normalization, to adjust for multiple BLAST similarities by similarity weighting, and to select significant similarities using relative alignment lengths. In benchmark tests, the GAAS method was robust to both high percentages of unknown sequences and to variations in metagenomic sequence read lengths. Re-analysis of the Sargasso Sea virome using GAAS indicated that standard methodologies for metagenomic analysis may dramatically underestimate the abundance and importance of organisms with small genomes in environmental systems. Using GAAS, we conducted a meta-analysis of microbial and viral average genome lengths in over 150 metagenomes from four biomes to determine whether genome lengths vary consistently between and within biomes, and between microbial and viral communities from the same environment. Significant differences between biomes and within aquatic sub-biomes (oceans, hypersaline systems, freshwater, and microbialites) suggested that average genome length is a fundamental property of environments driven by factors at the sub-biome level. The behavior of paired viral and microbial metagenomes from the same environment indicated that microbial and viral average genome sizes are independent of each other, but indicative of community responses to stressors and environmental conditions.