The MexTAg collaborative cross: host genetics affects asbestos related disease latency, but has little influence once tumours develop.

Objectives: This study combines two innovative mouse models in a major gene discovery project to assess the influence of host genetics on asbestos related disease (ARD). Conventional genetics studies provided evidence that some susceptibility to mesothelioma is genetic. However, the identification of host modifier genes, the roles they may play, and whether they contribute to disease susceptibility remain unknown. Here we report a study designed to rapidly identify genes associated with mesothelioma susceptibility by combining the Collaborative Cross (CC) resource with the well-characterised MexTAg mesothelioma mouse model. Methods: The CC is a powerful mouse resource that harnesses over 90% of common genetic variation in the mouse species, allowing rapid identification of genes mediating complex traits. MexTAg mice rapidly, uniformly, and predictably develop mesothelioma, but only after asbestos exposure. To assess the influence of host genetics on ARD, we crossed 72 genetically distinct CC mouse strains with MexTAg mice and exposed the resulting CC-MexTAg (CCMT) progeny to asbestos and monitored them for traits including overall survival, the time to ARD onset (latency), the time between ARD onset and euthanasia (disease progression) and ascites volume. We identified phenotype-specific modifier genes associated with these traits and we validated the role of human orthologues in asbestos-induced carcinogenesis using human mesothelioma datasets. Results: We generated 72 genetically distinct CCMT strains and exposed their progeny (2,562 in total) to asbestos. Reflecting the genetic diversity of the CC, there was considerable variation in overall survival and disease latency. Surprisingly, however, there was no variation in disease progression, demonstrating that host genetic factors do have a significant influence during disease latency but have a limited role once disease is established. Quantitative trait loci (QTL) affecting ARD survival/latency were identified on chromosomes 6, 12 and X. Of the 97-protein coding candidate modifier genes that spanned these QTL, eight genes (CPED1, ORS1, NDUFA1, HS1BP3, IL13RA1, LSM8, TES and TSPAN12) were found to significantly affect outcome in both CCMT and human mesothelioma datasets. Conclusion: Host genetic factors affect susceptibility to development of asbestos associated disease. However, following mesothelioma establishment, genetic variation in molecular or immunological mechanisms did not affect disease progression. Identification of multiple candidate modifier genes and their human homologues with known associations in other advanced stage or metastatic cancers highlights the complexity of ARD and may provide a pathway to identify novel therapeutic targets.

Time-course RNAseq data of murine AB1 mesothelioma and Renca renal cancer following immune checkpoint therapy.

Time-critical transcriptional events in the immune microenvironment are important for response to immune checkpoint blockade (ICB), yet these events are difficult to characterise and remain incompletely understood. Here, we present whole tumor RNA sequencing data in the context of treatment with ICB in murine models of AB1 mesothelioma and Renca renal cell cancer. We sequenced 144 bulk RNAseq samples from these two cancer types across 4 time points prior and after treatment with ICB. We also performed single-cell sequencing on 12 samples of AB1 and Renca tumors an hour before ICB administration. Our samples were equally distributed between responders and non-responders to treatment. Additionally, we sequenced AB1-HA mesothelioma tumors treated with two sample dissociation protocols to assess the impact of these protocols on the quality transcriptional information in our samples. These datasets provide time-course information to transcriptionally characterize the ICB response and provide detailed information at the single-cell level of the early tumor microenvironment prior to ICB therapy.

Immune checkpoint therapy responders display early clonal expansion of tumor infiltrating lymphocytes.

Immune checkpoint therapy (ICT) causes durable tumour responses in a subgroup of patients, but it is not well known how T cell receptor beta (TCRß) repertoire dynamics contribute to the therapeutic response. Using murine models that exclude variation in host genetics, environmental factors and tumour mutation burden, limiting variation between animals to naturally diverse TCRß repertoires, we applied TCRseq, single cell RNAseq and flow cytometry to study TCRß repertoire dynamics in ICT responders and non-responders. Increased oligoclonal expansion of TCRß clonotypes was observed in responding tumours. Machine learning identified TCRß CDR3 signatures unique to each tumour model, and signatures associated with ICT response at various timepoints before or during ICT. Clonally expanded CD8+ T cells in responding tumours post ICT displayed effector T cell gene signatures and phenotype. An early burst of clonal expansion during ICT is associated with response, and we report unique dynamics in TCRß signatures associated with ICT response.

Tretinoin improves the anti-cancer response to cyclophosphamide, in a model-selective manner.

BACKGROUND: Chemotherapy is included in treatment regimens for many solid cancers, but when administered as a single agent it is rarely curative. The addition of immune checkpoint therapy to standard chemotherapy regimens has improved response rates and increased survival in some cancers. However, most patients do not respond to treatment and immune checkpoint therapy can cause severe side effects. Therefore, there is a need for alternative immunomodulatory drugs that enhance chemotherapy. METHODS: We used gene expression data from cyclophosphamide (CY) responders and non-responders to identify existing clinically approved drugs that could phenocopy a chemosensitive tumor microenvironment (TME), and tested combination treatments in multiple murine cancer models. RESULTS: The vitamin A derivative tretinoin was the top predicted upstream regulator of response to CY. Tretinoin pre-treatment induced an inflammatory, interferon-associated TME, with increased infiltration of CD8 + T cells, sensitizing the tumor to subsequent chemotherapy. However, while combination treatment significantly improved survival and cure rate in a CD4+ and CD8+ T cell dependent manner in AB1-HA murine mesothelioma, this effect was model-selective, and could not be replicated using other cell lines. CONCLUSIONS: Despite the promising data in one model, the inability to validate the efficacy of combination treatment in multiple cancer models deprioritizes tretinoin/cyclophosphamide combination therapy for clinical translation.

Geldanamycin treatment does not result in anti-cancer activity in a preclinical model of orthotopic mesothelioma.

Mesothelioma is characterised by its aggressive invasive behaviour, affecting the surrounding tissues of the pleura or peritoneum. We compared an invasive pleural model with a non-invasive subcutaneous model of mesothelioma and performed transcriptomic analyses on the tumour samples. Invasive pleural tumours were characterised by a transcriptomic signature enriched for genes associated with MEF2C and MYOCD signaling, muscle differentiation and myogenesis. Further analysis using the CMap and LINCS databases identified geldanamycin as a potential antagonist of this signature, so we evaluated its potential in vitro and in vivo. Nanomolar concentrations of geldanamycin significantly reduced cell growth, invasion, and migration in vitro. However, administration of geldanamycin in vivo did not result in significant anti-cancer activity. Our findings show that myogenesis and muscle differentiation pathways are upregulated in pleural mesothelioma which may be related to the invasive behaviour. However, geldanamycin as a single agent does not appear to be a viable treatment for mesothelioma.

CD4+ T cells drive an inflammatory, TNF-α/IFN-rich tumor microenvironment responsive to chemotherapy.

While chemotherapy remains the first-line treatment for many cancers, it is still unclear what distinguishes responders from non-responders. Here, we characterize the chemotherapy-responsive tumor microenvironment in mice, using RNA sequencing on tumors before and after cyclophosphamide, and compare the gene expression profiles of responders with progressors. Responsive tumors have an inflammatory and highly immune infiltrated pre-treatment tumor microenvironment characterized by the enrichment of pathways associated with CD4+ T cells, interferons (IFNs), and tumor necrosis factor alpha (TNF-α). The same gene expression profile is associated with response to cyclophosphamide-based chemotherapy in patients with breast cancer. Finally, we demonstrate that tumors can be sensitized to cyclophosphamide and 5-FU chemotherapy by pre-treatment with recombinant TNF-α, IFNγ, and poly(I:C). Thus, a CD4+ T cell-inflamed pre-treatment tumor microenvironment is necessary for response to chemotherapy, and this state can be therapeutically attained by targeted immunotherapy.

Cancer chemotherapy: insights into cellular and tumor microenvironmental mechanisms of action.

Chemotherapy has historically been the mainstay of cancer treatment, but our understanding of what drives a successful therapeutic response remains limited. The diverse response of cancer patients to chemotherapy has been attributed principally to differences in the proliferation rate of the tumor cells, but there is actually very little experimental data supporting this hypothesis. Instead, other mechanisms at the cellular level and the composition of the tumor microenvironment appear to drive chemotherapy sensitivity. In particular, the immune system is a critical determinant of chemotherapy response with the depletion or knock-out of key immune cell populations or immunological mediators completely abrogating the benefits of chemotherapy in pre-clinical models. In this perspective, we review the literature regarding the known mechanisms of action of cytotoxic chemotherapy agents and the determinants of response to chemotherapy from the level of individual cells to the composition of the tumor microenvironment. We then summarize current work toward the development of dynamic biomarkers for response and propose a model for a chemotherapy sensitive tumor microenvironment.

Temporally restricted activation of IFNß signaling underlies response to immune checkpoint therapy in mice.

The biological determinants of the response to immune checkpoint blockade (ICB) in cancer remain incompletely understood. Little is known about dynamic biological events that underpin therapeutic efficacy due to the inability to frequently sample tumours in patients. Here, we map the transcriptional profiles of 144 responding and non-responding tumours within two mouse models at four time points during ICB. We find that responding tumours display on/fast-off kinetics of type-I-interferon (IFN) signaling. Phenocopying of this kinetics using time-dependent sequential dosing of recombinant IFNs and neutralizing antibodies markedly improves ICB efficacy, but only when IFNß is targeted, not IFNα. We identify Ly6C+/CD11b+ inflammatory monocytes as the primary source of IFNß and find that active type-I-IFN signaling in tumour-infiltrating inflammatory monocytes is associated with T cell expansion in patients treated with ICB. Together, our results suggest that on/fast-off modulation of IFNß signaling is critical to the therapeutic response to ICB, which can be exploited to drive clinical outcomes towards response.

Comprehensive Testing of Chemotherapy and Immune Checkpoint Blockade in Preclinical Cancer Models Identifies Additive Combinations.

Antibodies that target immune checkpoints such as cytotoxic T lymphocyte antigen 4 (CTLA-4) and the programmed cell death protein 1/ligand 1 (PD-1/PD-L1) are now a treatment option for multiple cancer types. However, as a monotherapy, objective responses only occur in a minority of patients. Chemotherapy is widely used in combination with immune checkpoint blockade (ICB). Although a variety of isolated immunostimulatory effects have been reported for several classes of chemotherapeutics, it is unclear which chemotherapeutics provide the most benefit when combined with ICB. We investigated 10 chemotherapies from the main canonical classes dosed at the clinically relevant maximum tolerated dose in combination with anti-CTLA-4/anti-PD-L1 ICB. We screened these chemo-immunotherapy combinations in two murine mesothelioma models from two different genetic backgrounds, and identified chemotherapies that produced additive, neutral or antagonistic effects when combined with ICB. Using flow cytometry and bulk RNAseq, we characterized the tumor immune milieu in additive chemo-immunotherapy combinations. 5-fluorouracil (5-FU) or cisplatin were additive when combined with ICB while vinorelbine and etoposide provided no additional benefit when combined with ICB. The combination of 5-FU with ICB augmented an inflammatory tumor microenvironment with markedly increased CD8+ T cell activation and upregulation of IFNγ, TNFα and IL-1ß signaling. The effective anti-tumor immune response of 5-FU chemo-immunotherapy was dependent on CD8+ T cells but was unaffected when TNFα or IL-1ß cytokine signaling pathways were blocked. Our study identified additive and non-additive chemotherapy/ICB combinations and suggests a possible role for increased inflammation in the tumor microenvironment as a basis for effective combination therapy.

Retinoic Acid Induces an IFN-Driven Inflammatory Tumour Microenvironment, Sensitizing to Immune Checkpoint Therapy.

With immune checkpoint therapy (ICT) having reshaped the treatment of many cancers, the next frontier is to identify and develop novel combination therapies to improve efficacy. Previously, we and others identified beneficial immunological effects of the vitamin A derivative tretinoin on anti-tumour immunity. Although it is known that tretinoin preferentially depletes myeloid derived suppressor cells in blood, little is known about the effects of tretinoin on the tumour microenvironment, hampering the rational design of clinical trials using tretinoin in combination with ICT. Here, we aimed to identify how tretinoin changed the tumour microenvironment in mouse tumour models, using flow cytometry and RNAseq, and we sought to use that information to establish optimal dosing and scheduling of tretinoin in combination with several ICT antibodies in multiple cancer models. We found that tretinoin rapidly induced an interferon dominated inflammatory tumour microenvironment, characterised by increased CD8+ T cell infiltration. This phenotype completely overlapped with the phenotype that was induced by ICT itself, and we confirmed that the combination further amplified this inflammatory milieu. The addition of tretinoin significantly improved the efficacy of anti-CTLA4/anti-PD-L1 combination therapy, and staggered scheduling was more efficacious than concomitant scheduling, in a dose-dependent manner. The positive effects of tretinoin could be extended to ICT antibodies targeting OX40, GITR and CTLA4 monotherapy in multiple cancer models. These data show that tretinoin induces an interferon driven, CD8+ T cell tumour microenvironment that is responsive to ICT.

PPARα and PPARγ activation is associated with pleural mesothelioma invasion but therapeutic inhibition is ineffective.

Mesothelioma is a cancer that typically originates in the pleura of the lungs. It rapidly invades the surrounding tissues, causing pain and shortness of breath. We compared cell lines injected either subcutaneously or intrapleurally and found that only the latter resulted in invasive and rapid growth. Pleural tumors displayed a transcriptional signature consistent with increased activity of nuclear receptors PPARα and PPARγ and with an increased abundance of endogenous PPAR-activating ligands. We found that chemical probe GW6471 is a potent, dual PPARα/γ antagonist with anti-invasive and anti-proliferative activity in vitro. However, administration of GW6471 at doses that provided sustained plasma exposure levels sufficient for inhibition of PPARα/γ transcriptional activity did not result in significant anti-mesothelioma activity in mice. Lastly, we demonstrate that the in vitro anti-tumor effect of GW6471 is off-target. We conclude that dual PPARα/γ antagonism alone is not a viable treatment modality for mesothelioma.

How to Better Understand the Influence of Host Genetics on Developing an Effective Immune Response to Thoracic Cancers.

Thoracic cancers pose a significant global health burden. Immune checkpoint blockade therapies have improved treatment outcomes, but durable responses remain limited. Understanding how the host immune system interacts with a developing tumor is essential for the rational development of improved treatments for thoracic malignancies. Recent technical advances have improved our understanding of the mutational burden of cancer cells and changes in cancer-specific gene expression, providing a detailed understanding of the complex biology underpinning tumor-host interactions. While there has been much focus on the genetic alterations associated with cancer cells and how they may impact treatment outcomes, how host genetics affects cancer development is also critical and will greatly determine treatment response. Genome-wide association studies (GWAS) have identified genetic variants associated with cancer predisposition. This approach has successfully identified host genetic risk factors associated with common thoracic cancers like lung cancer, but is less effective for rare cancers like malignant mesothelioma. To assess how host genetics impacts rare thoracic cancers, we used the Collaborative Cross (CC); a powerful murine genetic resource designed to maximize genetic diversity and rapidly identify genes associated with any biological trait. We are using the CC in conjunction with our asbestos-induced MexTAg mouse model, to identify host genes associated with mesothelioma development. Once genes that moderate tumor development and progression are known, human homologues can be identified and human datasets interrogated to validate their association with disease outcome. Furthermore, our CC-MexTAg animal model enables in-depth study of the tumor microenvironment, allowing the correlation of immune cell infiltration and gene expression signatures with disease development. This strategy provides a detailed picture of the underlying biological pathways associated with mesothelioma susceptibility and progression; knowledge that is crucial for the rational development of new diagnostic and therapeutic strategies. Here we discuss the influence of host genetics on developing an effective immune response to thoracic cancers. We highlight current knowledge gaps, and with a focus on mesothelioma, describe the development and application of the CC-MexTAg to overcome limitations and illustrate how the knowledge gained from this unique study will inform the rational design of future treatments of mesothelioma.

Malignant Pleural Effusions-A Window Into Local Anti-Tumor T Cell Immunity?

The success of immunotherapy that targets inhibitory T cell receptors for the treatment of multiple cancers has seen the anti-tumor immune response re-emerge as a promising biomarker of response to therapy. Longitudinal characterization of T cells in the tumor microenvironment (TME) helps us understand how to promote effective anti-tumor immunity. However, serial analyses at the tumor site are rarely feasible in clinical practice. Malignant pleural effusions (MPE) associated with thoracic cancers are an abnormal accumulation of fluid in the pleural space that is routinely drained for patient symptom control. This fluid contains tumor cells and immune cells, including lymphocytes, macrophages and dendritic cells, providing a window into the local tumor microenvironment. Recurrent MPE is common, and provides an opportunity for longitudinal analysis of the tumor site in a clinical setting. Here, we review the phenotype of MPE-derived T cells, comparing them to tumor and blood T cells. We discuss the benefits and limitations of their use as potential dynamic biomarkers of response to therapy.

Tumour draining lymph node-generated CD8 T cells play a role in controlling lung metastases after a primary tumour is removed but not when adjuvant immunotherapy is used.

Surgical resection of cancer remains the frontline therapy for millions of patients annually, but post-operative recurrence is common, with a relapse rate of around 45% for non-small cell lung cancer. The tumour draining lymph nodes (dLN) are resected at the time of surgery for staging purposes, and this cannot be a null event for patient survival and future response to immune checkpoint blockade treatment. This project investigates cancer surgery, lymphadenectomy, onset of metastatic disease, and response to immunotherapy in a novel model that closely reflects the clinical setting. In a murine metastatic lung cancer model, primary subcutaneous tumours were resected with associated dLNs remaining intact, completely resected or partially resected. Median survival after surgery was significantly shorter with complete dLN resection at the time of surgery (49 days (95%CI)) compared to when lymph nodes remained intact (> 88 days; p < 0.05). Survival was partially restored with incomplete lymph node resection and CD8 T cell dependent. Treatment with aCTLA4 whilst effective against the primary tumour was ineffective for metastatic lung disease. Conversely, aPD-1/aCD40 treatment was effective in both the primary and metastatic disease settings and restored the detrimental effects of complete dLN resection on survival. In this pre-clinical lung metastatic disease model that closely reflects the clinical setting, we observe decreased frequency of survival after complete lymphadenectomy, which was ameliorated with partial lymph node removal or with early administration of aPD-1/aCD40 therapy. These findings have direct relevance to surgical lymph node resection and adjuvant immunotherapy in lung cancer, and perhaps other cancer, patients.

A phase 1b clinical trial optimizing regulatory T cell depletion in combination with platinum-based chemotherapy in thoracic cancers.

Background: Single-agent cyclophosphamide can deplete regulatory T-cells (Treg). We aimed to determine optimal dosing and scheduling of oral cyclophosphamide, alongside pemetrexed-based chemotherapy, to deplete Treg in mesothelioma or non-small-cell lung cancer patients.Methods: 31 Patients received pemetrexed ± cisplatin or carboplatin on day 1 of a 21-day cycle (maximum 6 cycles). From cycle two, patients received cyclophosphamide, 50 mg/day, with intrapatient escalation to maximum 100/150 mg/day alternately. Immunological changes were examined by flow cytometry. Primary endpoint was Treg proportion of CD4+ T-cells, with doses tailored to target Treg nadir <4%.Results: Reduction in Treg proportion was observed on day 8 of all cycles, and was not augmented by cyclophosphamide. Few patients achieved the <4% Treg target. Treg proliferation reached nadir one week after chemotherapy, and peaked on day 1 of the subsequent cycle. Efficacy parameters were similar to chemotherapy alone. Seventeen percent of patients ceased cyclophosphamide due to toxicity.Conclusions: Specific Treg depletion to the degree seen with single-agent cyclophosphamide was not observed during pemetrexed-based chemotherapy. This study highlights the poor evidence basis for use of cyclophosphamide as an immunotherapeutic in combination with chemotherapy, and the importance of detailed flow cytometry studies.Trial registration: Clinical trial registration: www.anzctr.org.au identifier is ACTRN12609000260224.

Tumor Infiltrating Effector Memory Antigen-Specific CD8+ T Cells Predict Response to Immune Checkpoint Therapy.

Immune checkpoint therapy (ICT) results in durable responses in individuals with some cancers, but not all patients respond to treatment. ICT improves CD8+ cytotoxic T lymphocyte (CTL) function, but changes in tumor antigen-specific CTLs post-ICT that correlate with successful responses have not been well characterized. Here, we studied murine tumor models with dichotomous responses to ICT. We tracked tumor antigen-specific CTL frequencies and phenotype before and after ICT in responding and non-responding animals. Tumor antigen-specific CTLs increased within tumor and draining lymph nodes after ICT, and exhibited an effector memory-like phenotype, expressing IL-7R (CD127), KLRG1, T-bet, and granzyme B. Responding tumors exhibited higher infiltration of effector memory tumor antigen-specific CTLs, but lower frequencies of regulatory T cells compared to non-responders. Tumor antigen-specific CTLs persisted in responding animals and formed memory responses against tumor antigens. Our results suggest that increased effector memory tumor antigen-specific CTLs, in the presence of reduced immunosuppression within tumors is part of a successful ICT response. Temporal and nuanced analysis of T cell subsets provides a potential new source of immune based biomarkers for response to ICT.

Characteristics of TCR Repertoire Associated With Successful Immune Checkpoint Therapy Responses.

Immunotherapies have revolutionized cancer treatment. In particular, immune checkpoint therapy (ICT) leads to durable responses in some patients with some cancers. However, the majority of treated patients do not respond. Understanding immune mechanisms that underlie responsiveness to ICT will help identify predictive biomarkers of response and develop treatments to convert non-responding patients to responding ones. ICT primarily acts at the level of adaptive immunity. The specificity of adaptive immune cells, such as T and B cells, is determined by antigen-specific receptors. T cell repertoires can be comprehensively profiled by high-throughput sequencing at the bulk and single-cell level. T cell receptor (TCR) sequencing allows for sensitive tracking of dynamic changes in antigen-specific T cells at the clonal level, giving unprecedented insight into the mechanisms by which ICT alters T cell responses. Here, we review how the repertoire influences response to ICT and conversely how ICT affects repertoire diversity. We will also explore how changes to the repertoire in different anatomical locations can better correlate and perhaps predict treatment outcome. We discuss the advantages and limitations of current metrics used to characterize and represent TCR repertoire diversity. Discovery of predictive biomarkers could lie in novel analysis approaches, such as network analysis of amino acids similarities between TCR sequences. Single-cell sequencing is a breakthrough technology that can link phenotype with specificity, identifying T cell clones that are crucial for successful ICT. The field of immuno-sequencing is rapidly developing and cross-disciplinary efforts are required to maximize the analysis, application, and validation of sequencing data. Unravelling the dynamic behavior of the TCR repertoire during ICT will be highly valuable for tracking and understanding anti-tumor immunity, biomarker discovery, and ultimately for the development of novel strategies to improve patient outcomes.

Voluntary exercise in mesothelioma: effects on tumour growth and treatment response in a murine model.

OBJECTIVE: There is substantial evidence that exercise can safely reduce the risk of cancer and improve survival in different human cancer populations. Long latency periods associated with carcinogen-induced cancers like asbestos induced mesothelioma provide an opportunity to implement exercise as an intervention to delay or prevent disease development. However, there are limited studies investigating the ability of exercise to prevent or delay cancer, and exercise as a preventive strategy has never been assessed in models with a known carcinogen. We investigated the potential of voluntary exercise (VE) to delay development of asbestos related disease (ARD) in our well-characterised, asbestos induced MexTAg model of mesothelioma. RESULTS: Asbestos exposed MexTAg mice were given continuous or delayed access to VE and ARD assessed over time. We found that the addition of VE did not affect ARD development in asbestos exposed MexTAg mice. However, non-asbestos exposed, aged matched control mice participated in significantly more VE behaviours, suggesting subclinical development of ARD after asbestos exposure had a greater impact on VE participation than age alone. These data highlight the importance of model choice and the potential limitation that some pre-clinical studies may not accurately represent the clinical paradigm, particularly in the context of prevention studies.

Bilateral murine tumor models for characterizing the response to immune checkpoint blockade.

The therapeutic response to immune checkpoint blockade (ICB) is highly variable, not only between different cancers but also between patients with the same cancer type. The biological mechanisms underlying these differences in response are incompletely understood. Identifying correlates in patient tumor samples is challenging because of genetic and environmental variability. Murine studies usually compare different tumor models or treatments, introducing potential confounding variables. This protocol describes bilateral murine tumor models, derived from syngeneic cancer cell lines, that display a symmetrical yet dichotomous response to ICB. These models enable detailed analysis of whole tumors in a highly homogeneous background, combined with knowledge of the therapeutic outcome within a few weeks, and could potentially be used for mechanistic studies using other (immuno-)therapies. We discuss key considerations and describe how to use two cell lines as fully optimized models. We discuss experimental details, including proper inoculation technique to achieve symmetry and one-sided surgical tumor removal, which takes only 5 min per mouse. Furthermore, we outline the preparation of bulk tissue or single-cell suspensions for downstream analyses such as bulk RNA-seq, immunohistochemistry, single-cell RNA-seq and flow cytometry.