Mixed lineage kinase domain-like (MLKL) is a component of the “necrosome,” the multiprotein complex that triggers tumor necrosis factor (TNF)-induced cell death by necroptosis. To define the specific ...role and molecular mechanism of MLKL action, we generated MLKL-deficient mice and solved the crystal structure of MLKL. Although MLKL-deficient mice were viable and displayed no hematopoietic anomalies or other obvious pathology, cells derived from these animals were resistant to TNF-induced necroptosis unless MLKL expression was restored. Structurally, MLKL comprises a four-helical bundle tethered to the pseudokinase domain, which contains an unusual pseudoactive site. Although the pseudokinase domain binds ATP, it is catalytically inactive and its essential nonenzymatic role in necroptotic signaling is induced by receptor-interacting serine-threonine kinase 3 (RIPK3)-mediated phosphorylation. Structure-guided mutation of the MLKL pseudoactive site resulted in constitutive, RIPK3-independent necroptosis, demonstrating that modification of MLKL is essential for propagation of the necroptosis pathway downstream of RIPK3.
•The crystal structure of MLKL reveals an unusual pseudoactive site•The MLKL pseudokinase domain binds ATP but is not catalytically active•MLKL-deficient cells are resistant to TNF-induced necroptosis•MLKL pseudoactive site mutations cause caspase- and RIPK3-independent necroptosis
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.
Abstract
Motivation
Calling copy number alterations (CNAs) from RNA sequencing (RNA-Seq) is challenging, because of the marked variability in coverage across genes and paucity of single nucleotide ...polymorphisms (SNPs). We have adapted SuperFreq to call absolute and allele sensitive CNAs from RNA-Seq. SuperFreq uses an error-propagation framework to combine and maximize information from read counts and B-allele frequencies.
Results
We used datasets from The Cancer Genome Atlas (TCGA) to assess the validity of CNA calls from RNA-Seq. When ploidy estimates were consistent, we found agreement with DNA SNP-arrays for over 98% of the genome for acute myeloid leukaemia (TCGA-AML, n = 116) and 87% for colorectal cancer (TCGA-CRC, n = 377). The sensitivity of CNA calling from RNA-Seq was dependent on gene density. Using RNA-Seq, SuperFreq detected 78% of CNA calls covering 100 or more genes with a precision of 94%. Recall dropped for focal events, but this also depended on signal intensity. For example, in the CRC cohort SuperFreq identified all cases (7/7) with high-level amplification of ERBB2, where the copy number was typically >20, but identified only 6% of cases (1/17) with moderate amplification of IGF2, which occurs over a smaller interval. SuperFreq offers an integrated platform for identification of CNAs and point mutations. As evidence of how SuperFreq can be applied, we used it to reproduce the established relationship between somatic mutation load and CNA profile in CRC using RNA-Seq alone.
Availability and implementation
SuperFreq is implemented in R and the code is available through GitHub: https://github.com/ChristofferFlensburg/SuperFreq/. Data and code to reproduce the figures are available at: https://gitlab.wehi.edu.au/flensburg.c/SuperFreq_RNA_paper. Data from TCGA (phs000178) was accessed from GDC following completion of a data access request through the database of Genotypes and Phenotypes (dbGaP). Data from the Leucegene consortium was downloaded from GEO (AML samples: GSE67040; normal CD34+ cells: GSE48846).
Supplementary information
Supplementary data are available at Bioinformatics online.
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.
Abstract
Molnupiravir (EIDD-2801) is an antiviral that received approval for the treatment of severe acute respiratory syndrome coronavirus 2 (SARS-CoV2) infection. Treatment of bacteria or cell ...lines with the active form of molnupiravir, β-d-N4-hydroxycytidine (NHC, or EIDD-1931), induces mutations in DNA. Yet these results contrast in vivo genotoxicity studies conducted during registration of the drug. Using a CRISPR screen, we found that inactivating the pyrimidine salvage pathway component uridine–cytidine kinase 2 (Uck2) renders cells more tolerant of NHC. Short-term exposure to NHC increased the mutation rate in a mouse myeloid cell line, with most mutations being T:A to C:G transitions. Inactivating Uck2 impaired the mutagenic activity of NHC, whereas over-expression of Uck2 enhanced mutagenesis. UCK2 is upregulated in many cancers and cell lines. Our results suggest differences in ribonucleoside metabolism contribute to the variable mutagenicity of NHC observed in cancer cell lines and primary tissues.
Graphical Abstract
Graphical Abstract
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.
The B-cell lymphoma 2 (BCL-2) inhibitor venetoclax has an emerging role in acute myeloid leukemia (AML), with promising response rates in combination with hypomethylating agents or low-dose ...cytarabine in older patients. The tolerability and efficacy of venetoclax in combination with intensive chemotherapy in AML is unknown.
Patients with AML who were ≥ 65 years (≥ 60 years if monosomal karyotype) and fit for intensive chemotherapy were allocated to venetoclax dose-escalation cohorts (range, 50-600 mg). Venetoclax was administered orally for 14 days each cycle. During induction, a 7-day prephase/dose ramp-up (days -6 to 0) was followed by an additional 7 days of venetoclax combined with infusional cytarabine 100 mg/m
on days 1-5 and idarubicin 12 mg/m
intravenously on days 2-3 (ie, 5 + 2). Consolidation (4 cycles) included 14 days of venetoclax (days -6 to 7) combined with cytarabine (days 1-2) and idarubicin (day 1). Maintenance venetoclax was permitted (7 cycles). The primary objective was to assess the optimal dose schedule of venetoclax with 5 + 2.
Fifty-one patients with a median age of 72 years (range, 63-80 years) were included. The maximum tolerated dose was not reached with venetoclax 600 mg/day. The main grade ≥ 3 nonhematologic toxicities during induction were febrile neutropenia (55%) and sepsis (35%). In contrast to induction, platelet recovery was notably delayed during consolidation cycles. The overall response rate (complete remission CR/CR with incomplete count recovery) was 72%; it was 97% in de novo AML and was 43% in secondary AML. During the venetoclax prephase, marrow blast reductions (≥ 50%) were noted in
-,
-, and
-mutant AML.
Venetoclax combined with 5 + 2 induction chemotherapy was safe and tolerable in fit older patients with AML. Although the optimal postremission therapy remains to be determined, the high remission rate in de novo AML warrants additional investigation (ANZ Clinical Trial Registry No. ACTRN12616000445471).