Genomic-transcriptomic evolution in lung cancer and metastasis.

Intratumour heterogeneity (ITH) fuels lung cancer evolution, which leads to immune evasion and resistance to therapy1. Here, using paired whole-exome and RNA sequencing data, we investigate intratumour transcriptomic diversity in 354 non-small cell lung cancer tumours from 347 out of the first 421 patients prospectively recruited into the TRACERx study2,3. Analyses of 947 tumour regions, representing both primary and metastatic disease, alongside 96 tumour-adjacent normal tissue samples implicate the transcriptome as a major source of phenotypic variation. Gene expression levels and ITH relate to patterns of positive and negative selection during tumour evolution. We observe frequent copy number-independent allele-specific expression that is linked to epigenomic dysfunction. Allele-specific expression can also result in genomic-transcriptomic parallel evolution, which converges on cancer gene disruption. We extract signatures of RNA single-base substitutions and link their aetiology to the activity of the RNA-editing enzymes ADAR and APOBEC3A, thereby revealing otherwise undetected ongoing APOBEC activity in tumours. Characterizing the transcriptomes of primary-metastatic tumour pairs, we combine multiple machine-learning approaches that leverage genomic and transcriptomic variables to link metastasis-seeding potential to the evolutionary context of mutations and increased proliferation within primary tumour regions. These results highlight the interplay between the genome and transcriptome in influencing ITH, lung cancer evolution and metastasis.

The evolution of non-small cell lung cancer metastases in TRACERx.

Metastatic disease is responsible for the majority of cancer-related deaths1. We report the longitudinal evolutionary analysis of 126 non-small cell lung cancer (NSCLC) tumours from 421 prospectively recruited patients in TRACERx who developed metastatic disease, compared with a control cohort of 144 non-metastatic tumours. In 25% of cases, metastases diverged early, before the last clonal sweep in the primary tumour, and early divergence was enriched for patients who were smokers at the time of initial diagnosis. Simulations suggested that early metastatic divergence more frequently occurred at smaller tumour diameters (less than 8 mm). Single-region primary tumour sampling resulted in 83% of late divergence cases being misclassified as early, highlighting the importance of extensive primary tumour sampling. Polyclonal dissemination, which was associated with extrathoracic disease recurrence, was found in 32% of cases. Primary lymph node disease contributed to metastatic relapse in less than 20% of cases, representing a hallmark of metastatic potential rather than a route to subsequent recurrences/disease progression. Metastasis-seeding subclones exhibited subclonal expansions within primary tumours, probably reflecting positive selection. Our findings highlight the importance of selection in metastatic clone evolution within untreated primary tumours, the distinction between monoclonal versus polyclonal seeding in dictating site of recurrence, the limitations of current radiological screening approaches for early diverging tumours and the need to develop strategies to target metastasis-seeding subclones before relapse.

The evolution of lung cancer and impact of subclonal selection in TRACERx.

Lung cancer is the leading cause of cancer-associated mortality worldwide1. Here we analysed 1,644 tumour regions sampled at surgery or during follow-up from the first 421 patients with non-small cell lung cancer prospectively enrolled into the TRACERx study. This project aims to decipher lung cancer evolution and address the primary study endpoint: determining the relationship between intratumour heterogeneity and clinical outcome. In lung adenocarcinoma, mutations in 22 out of 40 common cancer genes were under significant subclonal selection, including classical tumour initiators such as TP53 and KRAS. We defined evolutionary dependencies between drivers, mutational processes and whole genome doubling (WGD) events. Despite patients having a history of smoking, 8% of lung adenocarcinomas lacked evidence of tobacco-induced mutagenesis. These tumours also had similar detection rates for EGFR mutations and for RET, ROS1, ALK and MET oncogenic isoforms compared with tumours in never-smokers, which suggests that they have a similar aetiology and pathogenesis. Large subclonal expansions were associated with positive subclonal selection. Patients with tumours harbouring recent subclonal expansions, on the terminus of a phylogenetic branch, had significantly shorter disease-free survival. Subclonal WGD was detected in 19% of tumours, and 10% of tumours harboured multiple subclonal WGDs in parallel. Subclonal, but not truncal, WGD was associated with shorter disease-free survival. Copy number heterogeneity was associated with extrathoracic relapse within 1 year after surgery. These data demonstrate the importance of clonal expansion, WGD and copy number instability in determining the timing and patterns of relapse in non-small cell lung cancer and provide a comprehensive clinical cancer evolutionary data resource.

Antibodies against endogenous retroviruses promote lung cancer immunotherapy.

B cells are frequently found in the margins of solid tumours as organized follicles in ectopic lymphoid organs called tertiary lymphoid structures (TLS)1,2. Although TLS have been found to correlate with improved patient survival and response to immune checkpoint blockade (ICB), the underlying mechanisms of this association remain elusive1,2. Here we investigate lung-resident B cell responses in patients from the TRACERx 421 (Tracking Non-Small-Cell Lung Cancer Evolution Through Therapy) and other lung cancer cohorts, and in a recently established immunogenic mouse model for lung adenocarcinoma3. We find that both human and mouse lung adenocarcinomas elicit local germinal centre responses and tumour-binding antibodies, and further identify endogenous retrovirus (ERV) envelope glycoproteins as a dominant anti-tumour antibody target. ERV-targeting B cell responses are amplified by ICB in both humans and mice, and by targeted inhibition of KRAS(G12C) in the mouse model. ERV-reactive antibodies exert anti-tumour activity that extends survival in the mouse model, and ERV expression predicts the outcome of ICB in human lung adenocarcinoma. Finally, we find that effective immunotherapy in the mouse model requires CXCL13-dependent TLS formation. Conversely, therapeutic CXCL13 treatment potentiates anti-tumour immunity and synergizes with ICB. Our findings provide a possible mechanistic basis for the association of TLS with immunotherapy response.

Lung adenocarcinoma promotion by air pollutants.

A complete understanding of how exposure to environmental substances promotes cancer formation is lacking. More than 70 years ago, tumorigenesis was proposed to occur in a two-step process: an initiating step that induces mutations in healthy cells, followed by a promoter step that triggers cancer development1. Here we propose that environmental particulate matter measuring ≤2.5 µm (PM2.5), known to be associated with lung cancer risk, promotes lung cancer by acting on cells that harbour pre-existing oncogenic mutations in healthy lung tissue. Focusing on EGFR-driven lung cancer, which is more common in never-smokers or light smokers, we found a significant association between PM2.5 levels and the incidence of lung cancer for 32,957 EGFR-driven lung cancer cases in four within-country cohorts. Functional mouse models revealed that air pollutants cause an influx of macrophages into the lung and release of interleukin-1ß. This process results in a progenitor-like cell state within EGFR mutant lung alveolar type II epithelial cells that fuels tumorigenesis. Ultradeep mutational profiling of histologically normal lung tissue from 295 individuals across 3 clinical cohorts revealed oncogenic EGFR and KRAS driver mutations in 18% and 53% of healthy tissue samples, respectively. These findings collectively support a tumour-promoting role for  PM2.5 air pollutants  and provide impetus for public health policy initiatives to address air pollution to reduce disease burden.

Molecular and evolutionary processes generating variation in gene expression.

Heritable variation in gene expression is common within and between species. This variation arises from mutations that alter the form or function of molecular gene regulatory networks that are then filtered by natural selection. High-throughput methods for introducing mutations and characterizing their cis- and trans-regulatory effects on gene expression (particularly, transcription) are revealing how different molecular mechanisms generate regulatory variation, and studies comparing these mutational effects with variation seen in the wild are teasing apart the role of neutral and non-neutral evolutionary processes. This integration of molecular and evolutionary biology allows us to understand how the variation in gene expression we see today came to be and to predict how it is most likely to evolve in the future.

Genetic and immune landscape evolution in MMR-deficient colorectal cancer.

Mismatch repair-deficient (MMRd) colorectal cancers (CRCs) have high mutation burdens, which make these tumours immunogenic and many respond to immune checkpoint inhibitors. The MMRd hypermutator phenotype may also promote intratumour heterogeneity (ITH) and cancer evolution. We applied multiregion sequencing and CD8 and programmed death ligand 1 (PD-L1) immunostaining to systematically investigate ITH and how genetic and immune landscapes coevolve. All cases had high truncal mutation burdens. Despite pervasive ITH, driver aberrations showed a clear hierarchy. Those in WNT/ß-catenin, mitogen-activated protein kinase, and TGF-ß receptor family genes were almost always truncal. Immune evasion (IE) drivers, such as inactivation of genes involved in antigen presentation or IFN-γ signalling, were predominantly subclonal and showed parallel evolution. These IE drivers have been implicated in immune checkpoint inhibitor resistance or sensitivity. Clonality assessments are therefore important for the development of predictive immunotherapy biomarkers in MMRd CRCs. Phylogenetic analysis identified three distinct patterns of IE driver evolution: pan-tumour evolution, subclonal evolution, and evolutionary stasis. These, but neither mutation burdens nor heterogeneity metrics, significantly correlated with T-cell densities, which were used as a surrogate marker of tumour immunogenicity. Furthermore, this revealed that genetic and T-cell infiltrates coevolve in MMRd CRCs. Low T-cell densities in the subgroup without any known IE drivers may indicate an, as yet unknown, IE mechanism. PD-L1 was expressed in the tumour microenvironment in most samples and correlated with T-cell densities. However, PD-L1 expression in cancer cells was independent of T-cell densities but strongly associated with loss of the intestinal homeobox transcription factor CDX2. This explains infrequent PD-L1 expression by cancer cells and may contribute to a higher recurrence risk of MMRd CRCs with impaired CDX2 expression. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Endothelial cell death after ionizing radiation does not impair vascular structure in mouse tumor models.

The effect of radiation therapy on tumor vasculature has long been a subject of debate. Increased oxygenation and perfusion have been documented during radiation therapy. Conversely, apoptosis of endothelial cells in irradiated tumors has been proposed as a major contributor to tumor control. To examine these contradictions, we use multiphoton microscopy in two murine tumor models: MC38, a highly vascularized, and B16F10, a moderately vascularized model, grown in transgenic mice with tdTomato-labeled endothelium before and after a single (15 Gy) or fractionated (5 × 3 Gy) dose of radiation. Unexpectedly, even these high doses lead to little structural change of the perfused vasculature. Conversely, non-perfused vessels and blind ends are substantially impaired after radiation accompanied by apoptosis and reduced proliferation of their endothelium. RNAseq analysis of tumor endothelial cells confirms the modification of gene expression in apoptotic and cell cycle regulation pathways after irradiation. Therefore, we conclude that apoptosis of tumor endothelial cells after radiation does not impair vascular structure.

Low doses of ionising radiation: definitions and contexts.

The term 'low dose' is applied to different levels of dose depending on the circumstances of exposure, with the potential for confusion unless the reasoning is clear. The United Nations Scientific Committee on the Effects of Ionising Radiation has defined low absorbed doses of ionising radiation as below about 100 mGy, and low dose rates as below 0.1 mGy min-1 (6 mGy h-1). These values relate to the interpretation of scientific evidence from epidemiological and biological studies. The International Commission on Radiological Protection has used similar values of 100 mSv and 5 mSv h-1 and applied this categorisation directly to the specific situation of patients undergoing diagnostic procedures: doses below 100 mSv were referred to as 'low' and doses below 10 mSv as 'very low'. Consideration of other exposure situations suggest that the same terms can be used for exposures received by emergency workers. However, for workers and members of the public in planned exposure situations, it is suggested that the term 'low dose' applies to doses below 10 mSv and 1 mSv, respectively - that is, below the dose limits. In each case, dose is being used as a surrogate for risk - risks at low doses are uncertain and estimates may change, but order of magnitude considerations are sufficient in most cases. Doses of < 100 mSv, < 10 mSv and < 1 mSv correspond to life-time cancer risk estimates of the order of < 10-2, < 10-3 and < 10-4, respectively. .

The Cellular Response to Complex DNA Damage Induced by Ionising Radiation.

Radiotherapy (ionising radiation; IR) is utilised in the treatment of ~50% of all human cancers, and where the therapeutic effect is largely achieved through DNA damage induction. In particular, complex DNA damage (CDD) containing two or more lesions within one to two helical turns of the DNA is a signature of IR and contributes significantly to the cell killing effects due to the difficult nature of its repair by the cellular DNA repair machinery. The levels and complexity of CDD increase with increasing ionisation density (linear energy transfer, LET) of the IR, such that photon (X-ray) radiotherapy is deemed low-LET whereas some particle ions (such as carbon ions) are high-LET radiotherapy. Despite this knowledge, there are challenges in the detection and quantitative measurement of IR-induced CDD in cells and tissues. Furthermore, there are biological uncertainties with the specific DNA repair proteins and pathways, including components of DNA single and double strand break mechanisms, that are engaged in CDD repair, which very much depends on the radiation type and associated LET. However, there are promising signs that advancements are being made in these areas and which will enhance our understanding of the cellular response to CDD induced by IR. There is also evidence that targeting CDD repair, particularly through inhibitors against selected DNA repair enzymes, can exacerbate the impact of higher LET, which could be explored further in a translational context.

The oDGal Mouse: A Novel, Physiologically Relevant Rodent Model of Sporadic Alzheimer's Disease.

Sporadic Alzheimer's disease (sAD) represents a serious and growing worldwide economic and healthcare burden. Almost 95% of current AD patients are associated with sAD as opposed to patients presenting with well-characterized genetic mutations that lead to AD predisposition, i.e., familial AD (fAD). Presently, the use of transgenic (Tg) animals overexpressing human versions of these causative fAD genes represents the dominant research model for AD therapeutic development. As significant differences in etiology exist between sAD and fAD, it is perhaps more appropriate to develop novel, more sAD-reminiscent experimental models that would expedite the discovery of effective therapies for the majority of AD patients. Here we present the oDGal mouse model, a novel model of sAD that displays a range of AD-like pathologies as well as multiple cognitive deficits reminiscent of AD symptomology. Hippocampal cognitive impairment and pathology were delayed with N-acetyl-cysteine (NaC) treatment, which strongly suggests that reactive oxygen species (ROS) are the drivers of downstream pathologies such as elevated amyloid beta and hyperphosphorylated tau. These features demonstrate a desired pathophenotype that distinguishes our model from current transgenic rodent AD models. A preclinical model that presents a phenotype of non-genetic AD-like pathologies and cognitive deficits would benefit the sAD field, particularly when translating therapeutics from the preclinical to the clinical phase.

Differential Grainy head binding correlates with variation in chromatin structure and gene expression in Drosophila melanogaster.

Phenotypic evolution is often caused by variation in gene expression resulting from altered gene regulatory mechanisms. Genetic variation affecting chromatin remodeling has been identified as a potential source of variable gene expression; however, the roles of specific chromatin remodeling factors remain unclear. Here, we address this knowledge gap by examining the relationship between variation in gene expression, variation in chromatin structure, and variation in binding of the pioneer factor Grainy head between imaginal wing discs of two divergent strains of Drosophila melanogaster and their F1 hybrid. We find that (1) variation in Grainy head binding is mostly due to sequence changes that act in cis but are located outside of the canonical Grainy head binding motif, (2) variation in Grainy head binding correlates with changes in chromatin accessibility, and (3) this variation in chromatin accessibility, coupled with variation in Grainy head binding, correlates with variation in gene expression in some cases but not others. Interactions among these three molecular layers is complex, but these results suggest that genetic variation affecting the binding of pioneer factors contributes to variation in chromatin remodeling and the evolution of gene expression.

Genome-wide sexually antagonistic variants reveal long-standing constraints on sexual dimorphism in fruit flies.

The evolution of sexual dimorphism is constrained by a shared genome, leading to 'sexual antagonism', in which different alleles at given loci are favoured by selection in males and females. Despite its wide taxonomic incidence, we know little about the identity, genomic location, and evolutionary dynamics of antagonistic genetic variants. To address these deficits, we use sex-specific fitness data from 202 fully sequenced hemiclonal Drosophila melanogaster fly lines to perform a genome-wide association study (GWAS) of sexual antagonism. We identify approximately 230 chromosomal clusters of candidate antagonistic single nucleotide polymorphisms (SNPs). In contradiction to classic theory, we find no clear evidence that the X chromosome is a hot spot for sexually antagonistic variation. Characterising antagonistic SNPs functionally, we find a large excess of missense variants but little enrichment in terms of gene function. We also assess the evolutionary persistence of antagonistic variants by examining extant polymorphism in wild D. melanogaster populations and closely related species. Remarkably, antagonistic variants are associated with multiple signatures of balancing selection across the D. melanogaster distribution range and in their sister species D. simulans, indicating widespread and evolutionarily persistent (about 1 million years) genomic constraints on the evolution of sexual dimorphism. Based on our results, we propose that antagonistic variation accumulates because of constraints on the resolution of sexual conflict over protein coding sequences, thus contributing to the long-term maintenance of heritable fitness variation.

Olaparib increases the therapeutic index of hemithoracic irradiation compared with hemithoracic irradiation alone in a mouse lung cancer model.

BACKGROUND: The radiosensitising effect of the poly(ADP-ribose) polymerase inhibitor olaparib on tumours has been reported. However, its effect on normal tissues in combination with radiation has not been well studied. Herein, we investigated the therapeutic index of olaparib combined with hemithoracic radiation in a urethane-induced mouse lung cancer model. METHODS: To assess tolerability, A/J mice were treated with olaparib plus whole thorax radiation (13 Gy), body weight changes were monitored and normal tissue effects were assessed by histology. In anti-tumour (intervention) studies, A/J mice were injected with urethane to induce lung tumours, and were then treated with olaparib alone, left thorax radiation alone or the combination of olaparib plus left thorax radiation at 8 weeks (early intervention) or 18 weeks (late intervention) after urethane injection. Anti-tumour efficacy and normal tissue effects were assessed by visual inspection, magnetic resonance imaging and histology. RESULTS: Enhanced body weight loss and oesophageal toxicity were observed when olaparib was combined with whole thorax but not hemithorax radiation. In both the early and late intervention studies, olaparib increased the anti-tumour effects of hemithoracic irradiation without increasing lung toxicity. CONCLUSIONS: The addition of olaparib increased the therapeutic index of hemithoracic radiation in a mouse model of lung cancer.

Cost-effectiveness of baloxavir marboxil compared to laninamivir for the treatment of influenza in Japan.

BACKGROUND: Baloxavir marboxil (baloxavir) is a new oral antiviral for influenza types A and B. OBJECTIVES: To determine the cost-effectiveness of baloxavir versus laninamivir in otherwise healthy (OwH) adults in Japan. METHODS: A decision tree was utilized to describe the course of influenza and predict associated costs and quality-adjusted life-years (QALYs) over one year by antiviral. Costs were valued from the public healthcare payer perspective, including influenza test, antiviral acquisition, other medications, physician visits, other outpatient costs associated with influenza or drug-related adverse events (DRAEs), and hospitalizations. Resource utilization and unit costs were obtained from the analysis of the JammNet claims database. Health state utilities were obtained from a clinical trial of baloxavir and previous models, and were driven by influenza symptoms, DRAEs, and complications caused by influenza. Sensitivity analyses were also performed. RESULTS: The total payer expenditure per patient for baloxavir versus laninamivir was ¥9383 versus ¥9132. The additional acquisition costs of baloxavir were partly offset by the DRAE costs avoided. Baloxavir showed a small gain in QALYs versus laninamivir and the incremental cost per QALY gained (¥2,231,260) was lower than the considered willingness-to-pay threshold (¥5,000,000/QALY). Key model drivers were the probability of DRAEs and the duration of symptoms. The probability of baloxavir being cost-effective was 64%. CONCLUSION: This cost-effectiveness study on baloxavir suggests that it would be cost-effective compared to laninamivir in OwH adults in Japan. Further studies are needed in different settings such as high-risk population and with different comparators.

Early metal use and crematory practices in the American Southeast.

Long-distance exchange of copper objects during the Archaic Period (ca. 8000-3000 cal B.P.) is a bellwether of emergent social complexity in the Eastern Woodlands. Originating from the Great Lakes, the Canadian Maritimes, and the Appalachian Mountains, Archaic-age copper is found in significant amounts as far south as Tennessee and in isolated pockets at major trade centers in Louisiana but is absent from most of the southeastern United States. Here we report the discovery of a copper band found with the cremated remains of at least seven individuals buried in the direct center of a Late Archaic shell ring located in coastal Georgia. Late Archaic shell rings are massive circular middens thought to be constructed, in part, during large-scale ritual gatherings and feasting events. The exotic copper and cremated remains are unique in coastal South Carolina and Georgia where Archaic-age cremations are conspicuously absent and no other Archaic copper objects have been reported. Elemental data produced through laser ablation inductively coupled plasma mass spectrometry shows the copper originated from the Great Lakes, effectively extending Archaic copper exchange almost 1,000 km beyond its traditional boundaries. Similarities in mortuary practices and the presence of copper originating from the Great Lakes reveal the presence of long-distance exchange relations spanning vast portions of the eastern United States and suggest an unexpected level of societal complexity at shell ring localities. These findings are consistent with the hypothesis that elite actors solidified their positions through ritual gatherings and the long-distance exchange of exotic objects during the Archaic.

High-resolution µCT of a mouse embryo using a compact laser-driven X-ray betatron source.

In the field of X-ray microcomputed tomography (µCT) there is a growing need to reduce acquisition times at high spatial resolution (approximate micrometers) to facilitate in vivo and high-throughput operations. The state of the art represented by synchrotron light sources is not practical for certain applications, and therefore the development of high-brightness laboratory-scale sources is crucial. We present here imaging of a fixed embryonic mouse sample using a compact laser-plasma-based X-ray light source and compare the results to images obtained using a commercial X-ray µCT scanner. The radiation is generated by the betatron motion of electrons inside a dilute and transient plasma, which circumvents the flux limitations imposed by the solid or liquid anodes used in conventional electron-impact X-ray tubes. This X-ray source is pulsed (duration <30 fs), bright (>1010 photons per pulse), small (diameter <1 µm), and has a critical energy >15 keV. Stable X-ray performance enabled tomographic imaging of equivalent quality to that of the µCT scanner, an important confirmation of the suitability of the laser-driven source for applications. The X-ray flux achievable with this approach scales with the laser repetition rate without compromising the source size, which will allow the recording of high-resolution µCT scans in minutes.

Evolving Plastic Responses to External and Genetic Environments.

Phenotypic plasticity can mitigate adaptive trade-offs in fluctuating environments but how plasticity arises is little known. New research documents this process in a bacterial system. We highlight remarkable parallels to the evolution of sexual dimorphism and argue that their approach can aid our understanding of adaptive conflicts between the sexes.

ReN VM spheroids in matrix: A neural progenitor three-dimensional in vitro model reveals DYRK1A inhibitors as potential regulators of radio-sensitivity.

INTRODUCTION: Pre-clinical testing of small molecules for therapeutic development across many pathologies relies on the use of in-vitro and in-vivo models. When designed and implemented well, these models serve to predict the clinical outcome as well as the toxicity of the evaluated therapies. The two-dimensional (2D) reductionist approach where cells are incubated in a mono-layer on hard plastic microtiter plates is relatively inexpensive but not physiologically relevant. In contrast, well developed and applied three dimensional (3D) in vitro models could be employed to bridge the gap between 2D in vitro primary screening and expensive in vivo rodent models by incorporating key features of the tissue microenvironment to explore differentiation, cortical development, cancers and various neuronal dysfunctions. These features include an extracellular matrix, co-culture, tension and perfusion and could replace several hundred rodents in the drug screening validation cascade. METHODS: Human neural progenitor cells from middle brain (ReN VM, Merck Millipore, UK) were expanded as instructed by the supplier (Merck Millipore, UK), and then seeded in 96-well low-attachment plates (Corning, UK) to form multicellular spheroids followed by adding a Matrigel layer to mimic extracellular matrix around neural stem cell niche. ReN VM cells were then differentiated via EGF and bFGF deprivation for 7 days and were imaged at day 7. Radiotherapy was mimicked via gamma-radiation at 2Gy in the absence and presence of selected DYRK1A inhibitors Harmine, INDY and Leucettine 41 (L41). Cell viability was measured by AlamarBlue assay. Immunofluorescence staining was used to assess cell pluripotency marker SOX2 and differentiation marker GFAP. RESULTS: After 7 days of differentiation, neuron early differentiation marker (GFAP, red) started to be expressed among the cells expressing neural stem cell marker SOX2 (green). Radiation treatment caused significant morphology change including the reduced viability of the spheroids. These spheroids also revealed sensitizing potential of DYRK1A inhibitors tested in this study, including Harmine, INDY and L41. DISCUSSION & CONCLUSIONS: Combined with the benefit of greatly reducing the issues associated with in vivo rodent models, including reducing numbers of animals used in a drug screening cascade, cost, ethics, and potential animal welfare burden, we feel the well-developed and applied 3D neural spheroid model presented in this study will provide a crucial tool to evaluate combinatorial therapies, optimal drug concentrations and treatment dosages.