A risk-reward examination of sample multiplexing reagents for single cell RNA-Seq

Single-cell RNA sequencing (scRNA-Seq) has emerged as a powerful tool for understanding cellular heterogeneity and function. However the choice of sample multiplexing reagents can impact data quality and experimental outcomes. In this study, we compared various multiplexing reagents, including MULTI-Seq, Hashtag antibody, and CellPlex, across diverse sample types such as human peripheral blood mononuclear cells (PBMCs), mouse embryonic brain and patient-derived xenografts (PDXs). We found that all multiplexing reagents worked well in cell types robust to ex vivo manipulation but suffered from signal-to-noise issues in more delicate sample types. We compared multiple demultiplexing algorithms which differed in performance depending on data quality. We find that minor improvements to laboratory workflows such as titration and rapid processing are critical to optimal performance. We also compared the performance of fixed scRNA-Seq kits and highlight the advantages of the Parse Biosciences kit for fragile samples. Highly multiplexed scRNA-Seq experiments require more sequencing resources, therefore we evaluated CRISPR-based destruction of non-informative genes to enhance sequencing value. Our comprehensive analysis provides insights into the selection of appropriate sample multiplexing reagents and protocols for scRNA-Seq experiments, facilitating more accurate and cost-effective studies.

The microtubule inhibitor eribulin demonstrates efficacy in platinum-resistant and refractory high-grade serous ovarian cancer patient-derived xenograft models

Background

Despite initial response to platinum-based chemotherapy and PARP inhibitor therapy (PARPi), nearly all recurrent high-grade serous ovarian cancer (HGSC) will acquire lethal drug resistance; indeed, ~15% of individuals have de novo platinum-refractory disease.

Objectives

To determine the potential of anti-microtubule agent (AMA) therapy (paclitaxel, vinorelbine and eribulin) in platinum-resistant or refractory (PRR) HGSC by assessing response in patient-derived xenograft (PDX) models of HGSC.

Design and methods

Of 13 PRR HGSC PDX, six were primary PRR, derived from chemotherapy-naïve samples (one was BRCA2 mutant) and seven were from samples obtained following chemotherapy treatment in the clinic (five were mutant for either BRCA1 or BRCA2 (BRCA1/2), four with prior PARPi exposure), recapitulating the population of individuals with aggressive treatment-resistant HGSC in the clinic. Molecular analyses and in vivo treatment studies were undertaken.

Results

Seven out of thirteen PRR PDX (54%) were sensitive to treatment with the AMA, eribulin (time to progressive disease (PD) ⩾100 days from the start of treatment) and 11 out of 13 PDX (85%) derived significant benefit from eribulin [time to harvest (TTH) for each PDX with p < 0.002]. In 5 out of 10 platinum-refractory HGSC PDX (50%) and one out of three platinum-resistant PDX (33%), eribulin was more efficacious than was cisplatin, with longer time to PD and significantly extended TTH (each PDX p < 0.02). Furthermore, four of these models were extremely sensitive to all three AMA tested, maintaining response until the end of the experiment (120d post-treatment start). Despite harbouring secondary BRCA2 mutations, two BRCA2-mutant PDX models derived from heavily pre-treated individuals were sensitive to AMA. PRR HGSC PDX models showing greater sensitivity to AMA had high proliferative indices and oncogene expression. Two PDX models, both with prior chemotherapy and/or PARPi exposure, were refractory to all AMA, one of which harboured the SLC25A40-ABCB1 fusion, known to upregulate drug efflux via MDR1.

Conclusion

The efficacy observed for eribulin in PRR HGSC PDX was similar to that observed for paclitaxel, which transformed ovarian cancer clinical practice. Eribulin is therefore worthy of further consideration in clinical trials, particularly in ovarian carcinoma with early failure of carboplatin/paclitaxel chemotherapy.

BRCA1 secondary splice-site mutations drive exon-skipping and PARP inhibitor resistance

BRCA1 splice isoforms Δ11 and Δ11q can contribute to PARP inhibitor (PARPi) resistance by splicing-out the mutation-containing exon, producing truncated, partially-functional proteins. However, the clinical impact and underlying drivers of BRCA1 exon skipping remain undetermined. We analyzed nine ovarian and breast cancer patient derived xenografts (PDX) with BRCA1 exon 11 frameshift mutations for exon skipping and therapy response, including a matched PDX pair derived from a patient pre- and post-chemotherapy/PARPi. BRCA1 exon 11 skipping was elevated in PARPi resistant PDX tumors. Two independent PDX models acquired secondary BRCA1 splice site mutations (SSMs), predicted in silico to drive exon skipping. Predictions were confirmed using qRT-PCR, RNA sequencing, western blots and BRCA1 minigene modelling. SSMs were also enriched in post-PARPi ovarian cancer patient cohorts from the ARIEL2 and ARIEL4 clinical trials. We demonstrate that SSMs drive BRCA1 exon 11 skipping and PARPi resistance, and should be clinically monitored, along with frame-restoring secondary mutations.

Interferon-ε is a tumour suppressor and restricts ovarian cancer

High-grade serous ovarian cancers have low survival rates because of their late presentation with extensive peritoneal metastases and frequent chemoresistance1, and require new treatments guided by novel insights into pathogenesis. Here we describe the intrinsic tumour-suppressive activities of interferon-ε (IFNε). IFNε is constitutively expressed in epithelial cells of the fallopian tube, the cell of origin of high-grade serous ovarian cancers, and is then lost during development of these tumours. We characterize its anti-tumour activity in several preclinical models: ovarian cancer patient-derived xenografts, orthotopic and disseminated syngeneic models, and tumour cell lines with or without mutations in Trp53 and Brca genes. We use manipulation of the IFNε receptor IFNAR1 in different cell compartments, differential exposure status to IFNε and global measures of IFN signalling to show that the mechanism of the anti-tumour activity of IFNε involves direct action on tumour cells and, crucially, activation of anti-tumour immunity. IFNε activated anti-tumour T and natural killer cells and prevented the accumulation and activation of myeloid-derived suppressor cells and regulatory T cells. Thus, we demonstrate that IFNε is an intrinsic tumour suppressor in the female reproductive tract whose activities in models of established and advanced ovarian cancer, distinct from other type I IFNs, are compelling indications of potential new therapeutic approaches for ovarian cancer.

Targeting homologous recombination deficiency in uterine leiomyosarcoma

BACKGROUND

Uterine leiomyosarcoma (uLMS) is a rare and aggressive gynaecological malignancy, with individuals with advanced uLMS having a five-year survival of < 10%. Mutations in the homologous recombination (HR) DNA repair pathway have been observed in ~ 10% of uLMS cases, with reports of some individuals benefiting from poly (ADP-ribose) polymerase (PARP) inhibitor (PARPi) therapy, which targets this DNA repair defect. In this report, we screened individuals with uLMS, accrued nationally, for mutations in the HR repair pathway and explored new approaches to therapeutic targeting.

METHODS

A cohort of 58 individuals with uLMS were screened for HR Deficiency (HRD) using whole genome sequencing (WGS), whole exome sequencing (WES) or NGS panel testing. Individuals identified to have HRD uLMS were offered PARPi therapy and clinical outcome details collected. Patient-derived xenografts (PDX) were generated for therapeutic targeting.

RESULTS

All 13 uLMS samples analysed by WGS had a dominant COSMIC mutational signature 3; 11 of these had high genome-wide loss of heterozygosity (LOH) (> 0.2) but only two samples had a CHORD score > 50%, one of which had a homozygous pathogenic alteration in an HR gene (deletion in BRCA2). A further three samples harboured homozygous HRD alterations (all deletions in BRCA2), detected by WES or panel sequencing, with 5/58 (9%) individuals having HRD uLMS. All five individuals gained access to PARPi therapy. Two of three individuals with mature clinical follow up achieved a complete response or durable partial response (PR) with the subsequent addition of platinum to PARPi upon minor progression during initial PR on PARPi. Corresponding PDX responses were most rapid, complete and sustained with the PARP1-specific PARPi, AZD5305, compared with either olaparib alone or olaparib plus cisplatin, even in a paired sample of a BRCA2-deleted PDX, derived following PARPi therapy in the patient, which had developed PARPi-resistance mutations in PRKDC, encoding DNA-PKcs.

CONCLUSIONS

Our work demonstrates the value of identifying HRD for therapeutic targeting by PARPi and platinum in individuals with the aggressive rare malignancy, uLMS and suggests that individuals with HRD uLMS should be included in trials of PARP1-specific PARPi.

Epithelial-to-Mesenchymal Transition Supports Ovarian Carcinosarcoma Tumorigenesis and Confers Sensitivity to Microtubule Targeting with Eribulin

Ovarian carcinosarcoma (OCS) is an aggressive and rare tumor type with limited treatment options. OCS is hypothesized to develop via the combination theory, with a single progenitor resulting in carcinomatous and sarcomatous components, or alternatively via the conversion theory, with the sarcomatous component developing from the carcinomatous component through epithelial-to-mesenchymal transition (EMT). In this study, we analyzed DNA variants from isolated carcinoma and sarcoma components to show that OCS from 18 women is monoclonal. RNA sequencing indicated that the carcinoma components were more mesenchymal when compared with pure epithelial ovarian carcinomas, supporting the conversion theory and suggesting that EMT is important in the formation of these tumors. Preclinical OCS models were used to test the efficacy of microtubule-targeting drugs, including eribulin, which has previously been shown to reverse EMT characteristics in breast cancers and induce differentiation in sarcomas. Vinorelbine and eribulin more effectively inhibited OCS growth than standard-of-care platinum-based chemotherapy, and treatment with eribulin reduced mesenchymal characteristics and N-MYC expression in OCS patient-derived xenografts. Eribulin treatment resulted in an accumulation of intracellular cholesterol in OCS cells, which triggered a downregulation of the mevalonate pathway and prevented further cholesterol biosynthesis. Finally, eribulin increased expression of genes related to immune activation and increased the intratumoral accumulation of CD8+ T cells, supporting exploration of immunotherapy combinations in the clinic. Together, these data indicate that EMT plays a key role in OCS tumorigenesis and support the conversion theory for OCS histogenesis. Targeting EMT using eribulin could help improve OCS patient outcomes.

SIGNIFICANCE

Genomic analyses and preclinical models of ovarian carcinosarcoma support the conversion theory for disease development and indicate that microtubule inhibitors could be used to suppress EMT and stimulate antitumor immunity.

Improving PARP inhibitor efficacy in high-grade serous ovarian carcinoma: A focus on the immune system

High-grade serous ovarian carcinoma (HGSOC) is a genomically unstable malignancy responsible for over 70% of all deaths due to ovarian cancer. With roughly 50% of all HGSOC harboring defects in the homologous recombination (HR) DNA repair pathway (e.g., BRCA1/2 mutations), the introduction of poly ADP-ribose polymerase inhibitors (PARPi) has dramatically improved outcomes for women with HR defective HGSOC. By blocking the repair of single-stranded DNA damage in cancer cells already lacking high-fidelity HR pathways, PARPi causes the accumulation of double-stranded DNA breaks, leading to cell death. Thus, this synthetic lethality results in PARPi selectively targeting cancer cells, resulting in impressive efficacy. Despite this, resistance to PARPi commonly develops through diverse mechanisms, such as the acquisition of secondary BRCA1/2 mutations. Perhaps less well documented is that PARPi can impact both the tumour microenvironment and the immune response, through upregulation of the stimulator of interferon genes (STING) pathway, upregulation of immune checkpoints such as PD-L1, and by stimulating the production of pro-inflammatory cytokines. Whilst targeted immunotherapies have not yet found their place in the clinic for HGSOC, the evidence above, as well as ongoing studies exploring the synergistic effects of PARPi with immune agents, including immune checkpoint inhibitors, suggests potential for targeting the immune response in HGSOC. Additionally, combining PARPi with epigenetic-modulating drugs may improve PARPi efficacy, by inducing a BRCA-defective phenotype to sensitise resistant cancer cells to PARPi. Finally, invigorating an immune response during PARPi therapy may engage anti-cancer immune responses that potentiate efficacy and mitigate the development of PARPi resistance. Here, we will review the emerging PARPi literature with a focus on PARPi effects on the immune response in HGSOC, as well as the potential of epigenetic combination therapies. We highlight the potential of transforming HGSOC from a lethal to a chronic disease and increasing the likelihood of cure.

Acquired RAD51C promoter methylation loss causes PARP inhibitor resistance in high grade serous ovarian carcinoma

In high-grade serous ovarian carcinoma (HGSC), deleterious mutations in DNA repair gene RAD51C are established drivers of defective homologous recombination and are emerging biomarkers of PARP inhibitor (PARPi) sensitivity. RAD51C promoter methylation (meRAD51C) is detected at similar frequencies to mutations, yet its effects on PARPi responses remain unresolved. In this study, three HGSC patient-derived xenograft (PDX) models with methylation at most or all examined CpG sites in the RAD51C promoter show responses to PARPi. Both complete and heterogeneous methylation patterns were associated with RAD51C gene silencing and homologous recombination deficiency (HRD). PDX models lost meRAD51C following treatment with PARPi rucaparib or niraparib, where a single unmethylated copy of RAD51C was sufficient to drive PARPi resistance. Genomic copy number profiling of one of the PDX models using SNP arrays revealed that this resistance was acquired independently in two genetically distinct lineages. In a cohort of 11 patients with RAD51C-methylated HGSC, various patterns of meRAD51C were associated with genomic 'scarring', indicative of HRD history, but exhibited no clear correlations with clinical outcome. Differences in methylation stability under treatment pressure were also observed between patients, where one HGSC was found to maintain meRAD51C after 6 lines of therapy (4 platinum-based), whilst another HGSC sample was found to have heterozygous meRAD51C and elevated RAD51C gene expression (relative to homozygous meRAD51C controls) after only neo-adjuvant chemotherapy. As meRAD51C loss in a single gene copy was sufficient to cause PARPi resistance in PDX, methylation zygosity should be carefully assessed in previously treated patients when considering PARPi therapy.

Abstract 2057: BRCA1 D11q isoform expression impacts PARP inhibitor responses in high grade serous ovarian carcinoma patient derived xenografts

Background: Acquired PARP inhibitor (PARPi) resistance in High Grade Serous Ovarian Carcinoma (HGSOC) commonly occurs via restored homologous recombination DNA repair due to secondary mutations in BRCA1/2. However, overexpression of certain BRCA1 splice isoforms can also contribute to PARPi resistance in HGSOC. This includes the D11q isoform of BRCA1, where deleterious mutations in the 11q region of BRCA1 (exon 10) can be spliced out resulting in a truncated but partially functional BRCA1 protein. Here we describe four HGSOC Patient Derived Xenografts (PDX) with primary BRCA1 11q mutations, where D11q isoform expression has been shown to correlate with platinum and PARPi responses in vivo.

Methods: PDX were derived from chemo-naive and pre-treated patients (Table 1). PDX were treated with PARPi rucaparib (300mg/kg) and cisplatin (4mg/kg). D11q isoform expression was assessed using two-step RT-qPCR. DNA sequencing was performed using the BROCA Cancer Risk Panel.

Results: PARPi resistant PDX #1032 and #1049 had high D11q expression relative to D11q-high cell line UWB1.289 and the PARPi-responsive models (Table 1). Interestingly, PDX #1049 harboured a second BRCA1 exon 11 donor splice site mutation, shown previously to enhance splicing and increase abundance of the D11q isoform. Indeed, other proximal splice site mutations have been found following relapse in patients with BRCA1 11q mutations post-PARPi.

Conclusions: Our work in PDX provides additional functional evidence that BRCA1 D11q expression impacts on PARPi and platinum responses, and is a clinically relevant mechanism of PARPi resistance in a subset of patients with BRCA1 mutations. Whilst high D11q expression can be associated with secondary BRCA1 splice site mutations (PDX #1049), it can also occur in their absence (PDX #1032). Thus, BRCA1 D11q isoform expression should be considered a potential mechanism of PARP-resistance in patients with BRCA1 11q mutations.

Table 1.Summary of results.PDX modelPrimary BRCA1 mutationBRCA1 D11q isoform expressionSecondary BRCA1 mutationsOther mutations (BROCA)PDX platinum responsePDX PARPi responsePatient treatments prior to PDX generation#1049Germline BRCA1: c.2475delCVery highBRCA1: c.4096+3A&gt;C (33%) - splice site mutation known to cause increased D11q expressionTP53: c.428_429insCAStable diseaseNon-responsive - progressive disease3 lines of platinum chemotherapy and PARPi#1032Somatic BRCA1: c.1251delTHighNone detectedTP53: c.97-1_97-6del; FANCA: c.3788_3790delResistant - progressive diseaseNon-responsive - progressive disease3 lines of platinum chemotherapy and PARPi#56Germline BRCA1: c.894_895delTGModerateNone detectedNone detectedSensitive - complete responseMixed responsesNo prior treatment#206Germline BRCA1: c.3817C&gt;TLowNone in PDX, one found in post-PARPi patient biopsy: BRCA1:c.3746_4081delTP53:c.451C&gt;T; BLM rearrangementSensitive - complete responseResponsive - complete response1 line of platinum chemotherapy

Citation Format: Ksenija Nesic, Cassandra J. Vandenberg, Olga Kondrashova, John J. Krais, Neil Johnson, Elizabeth M. Swisher, Matthew J. Wakefield, Clare L. Scott. BRCA1 D11q isoform expression impacts PARP inhibitor responses in high grade serous ovarian carcinoma patient derived xenografts [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 2057.

CHK2 Inhibition Provides a Strategy to Suppress Hematologic Toxicity from PARP Inhibitors

Patients with cancer treated with PARP inhibitors (PARPi) experience various side effects, with hematologic toxicity being most common. Short-term treatment of mice with olaparib resulted in depletion of reticulocytes, B-cell progenitors, and immature thymocytes, whereas longer treatment induced broader myelosuppression. We performed a CRISPR/Cas9 screen that targeted DNA repair genes in Eμ-Myc pre-B lymphoma cell lines as a way to identify strategies to suppress hematologic toxicity from PARPi. The screen revealed that single-guide RNAs targeting the serine/threonine kinase checkpoint kinase 2 (CHK2) were enriched following olaparib treatment. Genetic or pharmacologic inhibition of CHK2-blunted PARPi response in lymphoid and myeloid cell lines, and in primary murine pre-B/pro-B cells. Using a Cas9 base editor, we found that blocking CHK2-mediated phosphorylation of p53 also impaired olaparib response. Our results identify the p53 pathway as a major determinant of the acute response to PARPi in normal blood cells and demonstrate that targeting CHK2 can short circuit this response. Cotreatment with a CHK2 inhibitor did not antagonize olaparib response in ovarian cancer cell lines. Selective inhibition of CHK2 may spare blood cells from the toxic influence of PARPi and broaden the utility of these drugs. IMPLICATIONS: We reveal that genetic or pharmacologic inhibition of CHK2 may offer a way to alleviate the toxic influence of PARPi in the hematologic system.

Consequences of Zmat3 loss in c-MYC- and mutant KRAS-driven tumorigenesis

TP53 is a critical tumor suppressor that is mutated in approximately 50% of human cancers. Unveiling the downstream target genes of TP53 that fulfill its tumor suppressor function is an area of intense investigation. Zmat3 (also known as Wig-1 or PAG608) is one such downstream target of p53, whose loss in hemopoietic stem cells lacking the apoptosis and cell cycle regulators, Puma and p21, respectively, promotes the development of leukemia. The function of Zmat3 in tumorigenesis however remains unclear. Here, to investigate which oncogenic drivers co-operate with Zmat3 loss to promote neoplastic transformation, we utilized Zmat3 knockout mice in models of c-MYC-driven lymphomagenesis and Kras^G12D-driven lung adenocarcinoma development. Interestingly, unlike loss of p53, Zmat3 germline loss had little impact on the rate of tumor development or severity of malignant disease upon either the c-MYC or Kras^G12D oncogenic activation. Furthermore, loss of Zmat3 failed to rescue Kras^G12D primary lung tumor cells from oncogene-induced senescence. Taken together, we conclude that in the context of c-MYC-driven lymphomagenesis or mutant Kras^G12D-driven lung adenocarcinoma development, additional co-occurring mutations are required to resolve Zmat3 tumor suppressive activity.

Deep profiling of apoptotic pathways with mass cytometry identifies a synergistic drug combination for killing myeloma cells

Multiple myeloma is an incurable and fatal cancer of immunoglobulin-secreting plasma cells. Most conventional therapies aim to induce apoptosis in myeloma cells but resistance to these drugs often arises and drives relapse. In this study, we sought to identify the best adjunct targets to kill myeloma cells resistant to conventional therapies using deep profiling by mass cytometry (CyTOF). We validated probes to simultaneously detect 26 regulators of cell death, mitosis, cell signaling, and cancer-related pathways at the single-cell level following treatment of myeloma cells with dexamethasone or bortezomib. Time-resolved visualization algorithms and machine learning random forest models (RFMs) delineated putative cell death trajectories and a hierarchy of parameters that specified myeloma cell survival versus apoptosis following treatment. Among these parameters, increased amounts of phosphorylated cAMP response element-binding protein (CREB) and the pro-survival protein, MCL-1, were defining features of cells surviving drug treatment. Importantly, the RFM prediction that the combination of an MCL-1 inhibitor with dexamethasone would elicit potent, synergistic killing of myeloma cells was validated in other cell lines, in vivo preclinical models and primary myeloma samples from patients. Furthermore, CyTOF analysis of patient bone marrow cells clearly identified myeloma cells and their key cell survival features. This study demonstrates the utility of CyTOF profiling at the single-cell level to identify clinically relevant drug combinations and tracking of patient responses for future clinical trials.

Anti-apoptotic A1 is not essential for lymphoma development in Eµ-Myc mice but helps sustain transplanted Eµ-Myc tumour cells

The transcription factor c-MYC regulates a multiplicity of genes involved in cellular growth, proliferation, metabolism and DNA damage response and its overexpression is a hallmark of many tumours. Since MYC promotes apoptosis under conditions of stress, such as limited availability of nutrients or cytokines, MYC-driven cells are very much dependent on signals that inhibit cell death. Stress signals trigger apoptosis via the pathway regulated by opposing fractions of the BCL-2 protein family and previous genetic studies have shown that the development of B lymphoid tumours in Eµ-Myc mice is critically dependent on expression of pro-survival BCL-2 relatives MCL-1, BCL-W and, to a lesser extent, BCL-X_L, but not BCL-2 itself, and that sustained growth of these lymphomas is dependent on MCL-1. Using recently developed mice that lack expression of all three functional pro-survival A1 genes, we show here that the kinetics of lymphoma development in Eµ-Myc mice and the competitive repopulation capacity of Eµ-Myc haemopoietic stem and progenitor cells is unaffected by the absence of A1. However, conditional loss of a single remaining functional A1 gene from transplanted A1-a^−/−A1-b^fl/flA1-c^−/− Eµ-Myc lymphomas slowed their expansion, significantly extending the life of the transplant recipients. Thus, A1 contributes to the survival of malignant Eµ-Myc-driven B lymphoid cells. These results strengthen the case for BFL-1, the human homologue of A1, being a valid target for drug development for MYC-driven tumours.

Mnt modulates Myc-driven lymphomagenesis

The transcriptional represser Mnt is a functional antagonist of the proto-oncoprotein Myc. Both Mnt and Myc utilise Max as an obligate partner for DNA binding, and Myc/Max and Mnt/Max complexes compete for occupancy at E-box DNA sequences in promoter regions. We have previously shown in transgenic mouse models that the phenotype and kinetics of onset of haemopoietic tumours varies with the level of Myc expression. We reasoned that a decrease in the level of Mnt would increase the functional level of Myc and accelerate Myc-driven tumorigenesis. We tested the impact of reduced Mnt in three models of myc transgenic mice and in p53+/- mice. To our surprise, mnt heterozygosity actually slowed Myc-driven tumorigenesis in vavP-MYC10 and Eμ-myc mice, suggesting that Mnt facilitates Myc-driven oncogenesis. To explore the underlying cause of the delay in tumour development, we enumerated Myc-driven cell populations in healthy young vavP-MYC10 and Eμ-myc mice, expecting that the reduced rate of leukaemogenesis in mnt heterozygous mice would be reflected in a reduced number of preleukaemic cells, due to increased apoptosis or reduced proliferation or both. However, no differences were apparent. Furthermore, when mnt+/+ and mnt+/- pre-B cells from healthy young Eμ-myc mice were compared in vitro, no differences were seen in their sensitivity to apoptosis or in cell size or cell cycling. Moreover, the frequencies of apoptotic, senescent and proliferating cells were comparable in vivo in mnt+/- and mnt+/+ Eμ-myc lymphomas. Thus, although mnt heterozygosity clearly slowed lymphomagenesis in vavP-MYC10 and Eμ-myc mice, the change(s) in cellular properties responsible for this effect remain to be identified.

Overexpression of Mcl-1 exacerbates lymphocyte accumulation and autoimmune kidney disease in lpr mice

Cell death by apoptosis has a critical role during embryonic development and in maintaining tissue homeostasis. In mammals, there are two converging apoptosis pathways: the 'extrinsic' pathway, which is triggered by engagement of cell surface 'death receptors' such as Fas/APO-1; and the 'intrinsic' pathway, which is triggered by diverse cellular stresses, and is regulated by pro-survival and pro-apoptotic members of the Bcl-2 family of proteins. Pro-survival Mcl-1, which can block activation of the pro-apoptotic proteins, Bax and Bak, appears critical for the survival and maintenance of multiple haemopoietic cell types. To investigate the impact on haemopoiesis of simultaneously inhibiting both apoptosis pathways, we introduced the vavP-Mcl-1 transgene, which causes overexpression of Mcl-1 protein in all haemopoietic lineages, into Fas^lpr/lpr mice, which lack functional Fas and are prone to autoimmunity. The combined mutations had a modest impact on myelopoiesis, primarily an increase in the macrophage/monocyte population in Mcl-1tg/lpr mice compared with lpr or Mcl-1tg mice. The impact on lymphopoiesis was striking, with a marked elevation in all major lymphoid subsets, including the non-conventional double-negative (DN) T cells (TCRβ⁺CD4^-CD8^-B220⁺) characteristic of Fas^lpr/lpr mice. Of note, the onset of autoimmunity was markedly accelerated in Mcl-1tg/lpr mice compared with lpr mice, and this was preceded by an increase in immunoglobulin (Ig)-producing cells and circulating autoantibodies. This degree of impact was surprising, given the relatively mild phenotype conferred by the vavP-Mcl-1 transgene by itself: a two- to threefold elevation of peripheral B and T cells, no significant increase in the non-conventional DN T-cell population and no autoimmune disease. Comparison of the phenotype with that of other susceptible mice suggests that the development of autoimmune disease in Mcl-1tg/lpr mice may be influenced not only by Ig-producing cells but also other haemopoietic cell types.

RAG-induced DNA lesions activate proapoptotic BIM to suppress lymphomagenesis in p53-deficient mice

Delbridge, Strasser, and collaborators show that potentially oncogenic RAG1/2-dependent DNA lesions trigger apoptosis through the induction of BIM, which functions as an efficient tumor suppressor.

Neoplastic transformation is driven by oncogenic lesions that facilitate unrestrained cell expansion and resistance to antiproliferative signals. These oncogenic DNA lesions, acquired through errors in DNA replication, gene recombination, or extrinsically imposed damage, are thought to activate multiple tumor suppressive pathways, particularly apoptotic cell death. DNA damage induces apoptosis through well-described p53-mediated induction of PUMA and NOXA. However, loss of both these mediators (even together with defects in p53-mediated induction of cell cycle arrest and cell senescence) does not recapitulate the tumor susceptibility observed in p53^−/− mice. Thus, potentially oncogenic DNA lesions are likely to also trigger apoptosis through additional, p53-independent processes. We found that loss of the BH3-only protein BIM accelerated lymphoma development in p53-deficient mice. This process was negated by concomitant loss of RAG1/2-mediated antigen receptor gene rearrangement. This demonstrates that BIM is critical for the induction of apoptosis caused by potentially oncogenic DNA lesions elicited by RAG1/2-induced gene rearrangement. Furthermore, this highlights the role of a BIM-mediated tumor suppressor pathway that acts in parallel to the p53 pathway and remains active even in the absence of wild-type p53 function, suggesting this may be exploited in the treatment of p53-deficient cancers.

Loss of a Single Mcl-1 Allele Inhibits MYC-Driven Lymphomagenesis by Sensitizing Pro-B Cells to Apoptosis

MCL-1 is critical for progenitor cell survival during emergency hematopoiesis, but its role in sustaining cells undergoing transformation and in lymphomagenesis is only poorly understood. We investigated the importance of MCL-1 in the survival of B lymphoid progenitors undergoing MYC-driven transformation and its functional interactions with pro-apoptotic BIM and PUMA and the tumor suppressor p53 in lymphoma development. Loss of one Mcl-1 allele almost abrogated MYC-driven-lymphoma development owing to a reduction in lymphoma initiating pre-B cells. Although loss of the p53 target PUMA had minor impact, loss of one p53 allele substantially accelerated lymphoma development when MCL-1 was limiting, most likely because p53 loss also causes defects in non-apoptotic tumor suppressive processes. Remarkably, loss of BIM restored the survival of lymphoma initiating cells and rate of tumor development. Thus, MCL-1 has a major role in lymphoma initiating pro-B cells to oppose BIM, which is upregulated in response to oncogenic stress.

BCL-2 is dispensable for thrombopoiesis and platelet survival

Navitoclax (ABT-263), an inhibitor of the pro-survival BCL-2 family proteins BCL-2, BCL-XL and BCL-W, has shown clinical efficacy in certain BCL-2-dependent haematological cancers, but causes dose-limiting thrombocytopaenia. The latter effect is caused by Navitoclax directly inducing the apoptotic death of platelets, which are dependent on BCL-XL for survival. Recently, ABT-199, a selective BCL-2 antagonist, was developed. It has shown promising anti-leukaemia activity in patients whilst sparing platelets, suggesting that the megakaryocyte lineage does not require BCL-2. In order to elucidate the role of BCL-2 in megakaryocyte and platelet survival, we generated mice with a lineage-specific deletion of Bcl2, alone or in combination with loss of Mcl1 or Bclx. Platelet production and platelet survival were analysed. Additionally, we made use of BH3 mimetics that selectively inhibit BCL-2 or BCL-XL. We show that the deletion of BCL-2, on its own or in concert with MCL-1, does not affect platelet production or platelet lifespan. Thrombocytopaenia in Bclx-deficient mice was not affected by additional genetic loss or pharmacological inhibition of BCL-2. Thus, BCL-2 is dispensable for thrombopoiesis and platelet survival in mice.

DNA damage-induced primordial follicle oocyte apoptosis and loss of fertility require TAp63-mediated induction of Puma and Noxa

Trp63, a transcription factor related to the tumor suppressor p53, is activated by diverse stimuli and can initiate a range of cellular responses. TAp63 is the predominant Trp53 family member in primordial follicle oocyte nuclei and is essential for their apoptosis triggered by DNA damage in vivo. After γ-irradiation, induction of the proapoptotic BH3-only members Puma and Noxa was observed in primordial follicle oocytes from WT and Trp53(-/-) mice but not in those from TAp63-deficient mice. Primordial follicle oocytes from mice lacking Puma or both Puma and Noxa were protected from γ-irradiation-induced apoptosis and, remarkably, could produce healthy offspring. Hence, PUMA and NOXA are critical for DNA damage-induced, TAp63-mediated primordial follicle oocyte apoptosis. Thus, blockade of PUMA may protect fertility during cancer therapy and prevent premature menopause, improving women's health.

Apoptosis-promoted tumorigenesis: gamma-irradiation-induced thymic lymphomagenesis requires Puma-driven leukocyte death

Although tumor development requires impaired apoptosis, we describe a novel paradigm of apoptosis-dependent tumorigenesis. Because DNA damage triggers apoptosis through p53-mediated induction of BH3-only proteins Puma and Noxa, we explored their roles in gamma-radiation-induced thymic lymphomagenesis. Surprisingly, whereas Noxa loss accelerated it, Puma loss ablated tumorigenesis. Tumor suppression by Puma deficiency reflected its protection of leukocytes from gamma-irradiation-induced death, because their glucocorticoid-mediated decimation in Puma-deficient mice activated cycling of stem/progenitor cells and restored thymic lymphomagenesis. Our demonstration that cycles of cell attrition and repopulation by stem/progenitor cells can drive tumorigenesis has parallels in human cancers, such as therapy-induced malignancies.

The BH3 mimetic ABT-737 targets selective Bcl-2 proteins and efficiently induces apoptosis via Bak/Bax if Mcl-1 is neutralized

Since apoptosis is impaired in malignant cells overexpressing prosurvival Bcl-2 proteins, drugs mimicking their natural antagonists, BH3-only proteins, might overcome chemoresistance. Of seven putative BH3 mimetics tested, only ABT-737 triggered Bax/Bak-mediated apoptosis. Despite its high affinity for Bcl-2, Bcl-x(L), and Bcl-w, many cell types proved refractory to ABT-737. We show that this resistance reflects ABT-737's inability to target another prosurvival relative, Mcl-1. Downregulation of Mcl-1 by several strategies conferred sensitivity to ABT-737. Furthermore, enforced Mcl-1 expression in a mouse lymphoma model conferred resistance. In contrast, cells overexpressing Bcl-2 remained highly sensitive to ABT-737. Hence, ABT-737 should prove efficacious in tumors with low Mcl-1 levels, or when combined with agents that inactivate Mcl-1, even to treat those tumors that overexpress Bcl-2.

The BRIP1 helicase functions independently of BRCA1 in the Fanconi anemia pathway for DNA crosslink repair

BRIP1 (also called BACH1) is a DEAH helicase that interacts with the BRCT domain of BRCA1 (refs. 1-6) and has an important role in BRCA1-dependent DNA repair and checkpoint functions. We cloned the chicken ortholog of BRIP1 and established a homozygous knockout in the avian B-cell line DT40. The phenotype of these brip1 mutant cells in response to DNA damage differs from that of brca1 mutant cells and more closely resembles that of fancc mutant cells, with a profound sensitivity to the DNA-crosslinking agent cisplatin and acute cell-cycle arrest in late S-G2 phase. These defects are corrected by expression of human BRIP1 lacking the BRCT-interaction domain. Moreover, in human cells exposed to mitomycin C, short interfering RNA-mediated knock-down of BRIP1 leads to a substantial increase in chromosome aberrations, a characteristic phenotype of cells derived from individuals with Fanconi anemia. Because brip1 mutant cells are proficient for ubiquitination of FANCD2 protein, our data indicate that BRIP1 has a function in the Fanconi anemia pathway that is independent of BRCA1 and downstream of FANCD2 activation.

BRCA1-independent ubiquitination of FANCD2

Monoubiquitination of the FANCD2 protein is a key step in the Fanconi anemia (FA) tumor suppressor pathway, coinciding with this molecule's accumulation at sites of genome damage. Strong circumstantial evidence points to a requirement for the BRCA1 gene product in this step. Here, we show that the purified BRCA1/BARD1 complex, together with E1 and UbcH5a, is sufficient to reconstitute the monoubiquitination of FANCD2 in vitro. Although siRNA-mediated knockdown of BRCA1 in human cells results in defective targeting of FANCD2 to sites of DNA damage, it does not lead to a defect in FANCD2 ubiquitination. Furthermore, ablation of the RING finger domains of either BRCA1 or BARD1 in the chicken B cell line DT40 also leaves FANCD2 modification intact. Consequently, while BRCA1 affects the accumulation of FANCD2 at sites of DNA damage, BRCA1/BARD1 E3 ligase activity is not essential for the monoubiquitination of FANCD2.

Activation of the E3 ligase function of the BRCA1/BARD1 complex by polyubiquitin chains

Loss of the tumour suppressor BRCA1 results in profound chromosomal instability. The fundamental defect underlying this catastrophic phenotype is not yet known. In vivo, BRCA1 forms a heterodimeric complex with BARD1. Both proteins contain an N-terminal zinc RING-finger domain which confers E3 ubiquitin ligase activity. We have isolated full-length human BRCA1/BARD1 complex and have shown that it has a dual E3 ubiquitin ligase activity. First, it mediates the monoubiquitylation of nucleosome core histones in vitro, including the variant histone H2AX that co-localizes with BRCA1 at sites of DNA damage. Secondly, BRCA1/BARD1 catalyses the formation of multiple polyubiquitin chains on itself. Remarkably, this auto-polyubiquitylation potentiates the E3 ubiquitin ligase activity of the BRCA1/BARD1 complex >20-fold. Even though BRCA1 has been reported to associate with a C-terminal ubiquitin hydrolase, BAP1, this enzyme does not appear to function in the deubiquitylation of the BRCA1/BARD1 complex.

Variant subunit specificity in the quaternary structure of Artemia hemoglobin

The brine shrimp Artemia has three extracellular hemoglobins (Hbs) that are developmentally expressed and exhibit distinct oxygen-binding characteristics (Heip, Moens, and Kondo 1978; Heip et al. 1978 ). These Hbs are composed of two polymers, each of which comprises nine covalently linked globin domains. Although the cDNA sequences of two nine-domain globins from Artemia have been published, there is evidence for the existence of further expressed globin genes (Manning, Trotman, and Tate 1990 ). In the present study extensive analysis at the cDNA and genomic levels was performed in order to determine the globin gene copy number in Artemia. Sequence and Southern analysis suggest that four Hb polymers (T1, T2, C1, and C2) are expressed in Artemia. In addition, there is also at least one globin pseudogene. Protein sequencing of the native Hbs revealed that there are limitations on which two polymers can associate. The composition of the Hbs has been determined to be: Hb I, C1C2; Hb II, C1T2; and Hb III, T1T2. These pairings allow the levels of the three Artemia Hbs to be regulated independently by polymer expression alone, therefore explaining the previously inconsistent developmental and hypoxia-induced expression patterns.

Variable substitution rates of the 18 domain sequences in Artemia hemoglobin

The Artemia hemoglobin is a dimer comprising two nine-domain covalent polymers in quaternary association. Each polymer is encoded by a gene representing nine successive globin domains which have different sequences and are presumed to have been copied originally from a single-domain gene. Two different polymers exist as the result of a complete duplication of the nine-domain gene, allowing the formation of either homodimers or the heterodimer. The total population size of 18 domains comprising nine corresponding pairs, coupled with the probability that they reflect several hundred million years of evolution in the same lineage, provides a unique model in which the process of gene multiplication can be analyzed. The outcome has important implications for the reliability of local molecular clocks. The two polymers differ from each other at 11.7% of amino acid sites; however when corresponding individual domains are compared between polymers, amino acid substitution fluctuates by a factor of 2.7-fold from lowest to highest. This variation is not obvious at the DNA level: Domain pair identity values fluctuate by 1. 3-fold. Identity values are, however, uncorrected for multiple substitutions, and both silent and nonsilent changes are pooled. Therefore, to determine the variability in relative substitution rates at the DNA level, we have used the method of Li (1993, J Mol Evol 36:96-99) to determine estimates of nonsynonymous (KA) and synonymous (KS) substitutions per site for the nine pairs of domains. As expected, the overall level of silent substitutions (KS of 56. 9%) far exceeded nonsilent substitutions (KA of 6.7%); however, for corresponding domain pairs, KA fluctuates by 2.3-fold and KS by 1. 7-fold. The large discrepancies reflected in the expressed protein have accrued within a single lineage and the implication is that divergence dates of different genera based on amino acid sequences, even with well-studied proteins of reasonable size, can be wrong by a factor well in excess of 2.