Explore the vast range of titles available.

One of the most common analysis tasks in genomic research is to identify genes that are differentially expressed (DE) between experimental conditions. Empirical Bayes (EB) statistical tests using moderated genewise variances have been very effective for this purpose, especially when the number of biological replicate samples is small. The EB procedures can, however, be heavily influenced by a small number of genes with very large or very small variances. This article improves the differential expression tests by robustifying the hyperparameter estimation procedure. The robust procedure has the effect of decreasing the informativeness of the prior distribution for outlier genes while increasing its informativeness for other genes. This effect has the double benefit of reducing the chance that hypervariable genes will be spuriously identified as DE while increasing statistical power for the main body of genes. The robust EB algorithm is fast and numerically stable. The procedure allows exact small-sample null distributions for the test statistics and reduces exactly to the original EB procedure when no outlier genes are present. Simulations show that the robustified tests have similar performance to the original tests in the absence of outlier genes but have greater power and robustness when outliers are present. The article includes case studies for which the robust method correctly identifies and downweights genes associated with hidden covariates and detects more genes likely to be scientifically relevant to the experimental conditions. The new procedure is implemented in the limma software package freely available from the Bioconductor repository.

Analysing multiple cancer samples from an individual patient can provide insight into the way the disease evolves. Monitoring the expansion and contraction of distinct clones helps to reveal the mutations that initiate the disease and those that drive progression. Existing approaches for clonal tracking from sequencing data typically require the user to combine multiple tools that are not purpose-built for this task. Furthermore, most methods require a matched normal (non-tumour) sample, which limits the scope of application. We developed SuperFreq, a cancer exome sequencing analysis pipeline that integrates identification of somatic single nucleotide variants (SNVs) and copy number alterations (CNAs) and clonal tracking for both. SuperFreq does not require a matched normal and instead relies on unrelated controls. When analysing multiple samples from a single patient, SuperFreq cross checks variant calls to improve clonal tracking, which helps to separate somatic from germline variants, and to resolve overlapping CNA calls. To demonstrate our software we analysed 304 cancer-normal exome samples across 33 cancer types in The Cancer Genome Atlas (TCGA) and evaluated the quality of the SNV and CNA calls. We simulated clonal evolution through in silico mixing of cancer and normal samples in known proportion. We found that SuperFreq identified 93% of clones with a cellular fraction of at least 50% and mutations were assigned to the correct clone with high recall and precision. In addition, SuperFreq maintained a similar level of performance for most aspects of the analysis when run without a matched normal. SuperFreq is highly versatile and can be applied in many different experimental settings for the analysis of exomes and other capture libraries. We demonstrate an application of SuperFreq to leukaemia patients with diagnosis and relapse samples.

Venetoclax is a first-in-class cancer therapy that interacts with the cellular apoptotic machinery promoting apoptosis. Treatment of patients suffering chronic lymphocytic leukaemia with this BCL-2 antagonist has revealed emergence of a drug-selected BCL-2 mutation (G101V) in some patients failing therapy. To understand the molecular basis of this acquired resistance we describe the crystal structures of venetoclax bound to both BCL-2 and the G101V mutant. The pose of venetoclax in its binding site on BCL-2 reveals small but unexpected differences as compared to published structures of complexes with venetoclax analogues. The G101V mutant complex structure and mutant binding assays reveal that resistance is acquired by a knock-on effect of V101 on an adjacent residue, E152, with venetoclax binding restored by a E152A mutation. This provides a framework for considering analogues of venetoclax that might be effective in combating this mutation. The BCL-2 mutation G101V reduces venetoclax affinity and confers drug resistance in patients with chronic lymphocytic leukaemia. Here, the authors present crystal structures and biochemical analyses of venetoclax bound to BCL-2 and the G101V mutant, revealing the structural basis for venetoclax resistance.

Genomic instability is a hallmark of cancer and, as such, structural alterations and fusion genes are common events in the cancer landscape. RNA sequencing (RNA-Seq) is a powerful method for profiling cancers, but current methods for identifying fusion genes are optimised for short reads. JAFFA ( https://github.com/Oshlack/JAFFA/wiki ) is a sensitive fusion detection method that outperforms other methods with reads of 100 bp or greater. JAFFA compares a cancer transcriptome to the reference transcriptome, rather than the genome, where the cancer transcriptome is inferred using long reads directly or by de novo assembling short reads.

Calling fusion genes from RNA-seq data is well established, but other transcriptional variants are difficult to detect using existing approaches. To identify all types of variants in transcriptomes we developed MINTIE, an integrated pipeline for RNA-seq data. We take a reference-free approach, combining de novo assembly of transcripts with differential expression analysis to identify up-regulated novel variants in a case sample. We compare MINTIE with eight other approaches, detecting > 85% of variants while no other method is able to achieve this. We posit that MINTIE will be able to identify new disease variants across a range of disease types.

The gradual shift from cytotoxic drugs to highly selective, targeted therapeutic agents for cancer requires a parallel effort to characterize cancers at the molecular level to guide the choice of therapy for the individual patient. Here we review the genomic technologies that can be used to develop these drug response indicators, or biomarkers. We also discuss hurdles in their development and the implementation of biomarkers in clinical practice.

Background RNA sequencing allows the study of both gene expression changes and transcribed mutations, providing a highly effective way to gain insight into cancer biology. When planning the sequencing of a large cohort of samples, library size is a fundamental factor affecting both the overall cost and the quality of the results. Here we specifically address how overall library size influences the detection of somatic mutations in RNA-seq data in two acute myeloid leukaemia datasets. Results  We simulated shallower sequencing depths by downsampling 45 acute myeloid leukaemia samples (100 bp PE) that are part of the Leucegene project, which were originally sequenced at high depth. We compared the sensitivity of six methods of recovering validated mutations on the same samples. The methods compared are a combination of three popular callers (MuTect, VarScan, and VarDict) and two filtering strategies. We observed an incremental loss in sensitivity when simulating libraries of 80M, 50M, 40M, 30M and 20M fragments, with the largest loss detected with less than 30M fragments (below 90%, average loss of 7%). The sensitivity in recovering insertions and deletions varied markedly between callers, with VarDict showing the highest sensitivity (60%). Single nucleotide variant sensitivity is relatively consistent across methods, apart from MuTect, whose default filters need adjustment when using RNA-Seq. We also analysed 136 RNA-Seq samples from the TCGA-LAML cohort (50 bp PE) and assessed the change in sensitivity between the initial libraries (average 59M fragments) and after downsampling to 40M fragments. When considering single nucleotide variants in recurrently mutated myeloid genes we found a comparable performance, with a 6% average loss in sensitivity using 40M fragments. Conclusions Between 30M and 40M 100 bp PE reads are needed to recover 90–95% of the initial variants on recurrently mutated myeloid genes. To extend this result to another cancer type, an exploration of the characteristics of its mutations and gene expression patterns is suggested.

Background Patients with BRCA1 -like tumors correlate with improved response to DNA double-strand break-inducing therapy. A gene expression-based classifier was developed to distinguish between BRCA1 -like and non- BRCA1 -like tumors. We hypothesized that these tumors may also be more sensitive to PARP inhibitors than standard treatments. Methods A diagnostic gene expression signature ( BRCA1 ness) was developed using a centroid model with 128 triple-negative breast cancer samples from the EU FP7 RATHER project. This BRCA1 ness signature was then tested in HER2-negative patients ( n  = 116) from the I-SPY 2 TRIAL who received an oral PARP inhibitor veliparib in combination with carboplatin (V-C), or standard chemotherapy alone. We assessed the association between BRCA1 ness and pathologic complete response in the V-C and control arms alone using Fisher’s exact test, and the relative performance between arms (biomarker × treatment interaction, likelihood ratio p  < 0.05) using a logistic model and adjusting for hormone receptor status (HR). Results We developed a gene expression signature to identify BRCA1 -like status. In the I-SPY 2 neoadjuvant setting the BRCA1 ness signature associated significantly with response to V-C ( p  = 0.03), but not in the control arm ( p  = 0.45). We identified a significant interaction between BRCA1 ness and V-C ( p  = 0.023) after correcting for HR. Conclusions A genomic-based BRCA1 -like signature was successfully translated to an expression-based signature ( BRC1A ness). In the I-SPY 2 neoadjuvant setting, we determined that the BRCA1 ness signature is capable of predicting benefit of V-C added to standard chemotherapy compared to standard chemotherapy alone. Trial registration I-SPY 2 TRIAL beginning December 31, 2009: Neoadjuvant and Personalized Adaptive Novel Agents to Treat Breast Cancer (I-SPY 2), NCT01042379 .

Background Ganciclovir (GCV) is widely used in solid organ and haematopoietic stem cell transplant patients for prophylaxis and treatment of cytomegalovirus. It has long been considered a mutagen and carcinogen. However, the contribution of GCV to cancer incidence and other factors that influence its mutagenicity remains unknown. Methods This retrospective cohort study analysed genomics data for 121,771 patients who had undergone targeted sequencing compiled by the Genomics Evidence Neoplasia Information Exchange (GENIE) or Foundation Medicine (FM). A statistical approach was developed to identify patients with GCV-associated mutational signature (GCV sig ) from targeted sequenced data of tumour samples. Cell line exposure models were further used to quantify mutation burden and DNA damage caused by GCV and other antiviral and immunosuppressive drugs. Results Mutational profiles from 22 of 121,771 patient samples in the GENIE and FM cohorts showed evidence of GCV sig . A diverse range of cancers was represented. All patients with detailed clinical history available had previously undergone solid organ transplantation and received GCV and mycophenolate treatment. RAS hotspot mutations associated with GCV sig were present in 9 of the 22 samples, with all samples harbouring multiple GCV-associated protein-altering mutations in cancer driver genes. In vitro testing in cell lines showed that elevated DNA damage response and GCV sig are uniquely associated with GCV but not acyclovir, a structurally similar antiviral. Combination treatment of GCV with the immunosuppressant, mycophenolate mofetil (MMF), increased the misincorporation of GCV in genomic DNA and mutations attributed to GCV sig in cell lines and organoids. Conclusions In summary, GCV can cause a diverse range of cancers. Its mutagenicity may be potentiated by other therapies, such as mycophenolate, commonly co-prescribed with GCV for post-transplant patients. Further investigation of the optimal use of these drugs could help reduce GCV-associated mutagenesis in post-transplant patients.

Polycomb group proteins are transcriptional repressors that play a central role in the establishment and maintenance of gene expression patterns during development. Using mice with an N-ethyl-N-nitrosourea (ENU)-induced mutation in Suppressor of Zeste 12 (Suz12), a core component of Polycomb Repressive Complex 2 (PRC2), we show here that loss of Suz12 function enhances hematopoietic stem cell (HSC) activity. In addition to these effects on a wild-type genetic background, mutations in Suz12 are sufficient to ameliorate the stem cell defect and thrombocytopenia present in mice that lack the thrombopoietin receptor (c-Mpl). To investigate the molecular targets of the PRC2 complex in the HSC compartment, we examined changes in global patterns of gene expression in cells deficient in Suz12. We identified a distinct set of genes that are regulated by Suz12 in hematopoietic cells, including eight genes that appear to be highly responsive to PRC2 function within this compartment. These data suggest that PRC2 is required to maintain a specific gene expression pattern in hematopoiesis that is indispensable to normal stem cell function.