Long-term Healthcare Resource Utilization and Costs among Patients with Myasthenia Gravis: A Swedish Nationwide Population-based Study

INTRODUCTION

Health care costs and societal impact of myasthenia gravis (MG), a potentially life-threatening rare, chronic neuromuscular disease are sparsely studied. We assessed healthcare resource utilization (HCRU) and associated costs among patients with newly diagnosed (ND) and pre-existing (PE) MG in Sweden.

METHODS

This observational, retrospective cohort study used data from four linkable Swedish nationwide population-based registries. Adult MG patients receiving pharmacological treatment for MG and having ≥ 24-month follow-up during the period 1/1/2010 to 12/31/2017 were included.

RESULTS

A total of 1,275 patients were included in the analysis, of which 554 patients were categorized into the ND MG group and 721 into the PE MG group. Mean (±SD) age was 61.3 (±17.4) years and 52.3% were female. In first year post-diagnosis, ND patients had significantly higher utilization of acetylcholinesterase inhibitors (96.0% vs 83.9%), corticosteroids (59.6% vs 45.8%), thymectomy (12.1% vs 0.7%) and plasma exchange (3.8% vs 0.6%); had higher all-cause (70.9% vs 35.8%) and MG-related (62.5% vs 18.4%) hospitalization rates with 11 more hospitalization days (all p<0.01) and an increased risk of hospitalization (odds ratio [95% CI] = 4.4 [3.43, 5.64]) than PE MG. In year 1 post-diagnosis, ND MG patients incurred €7302 (p<0.01) higher total all-cause costs than PE MG, of which 84% were estimated to be MG-related and the majority (86%) were related to inpatient care. These results remained significant also after controlling for baseline demographics and comorbidities (p<0.01). In year 2 post-diagnosis, the all-cause medical costs decreased by ~55% for ND MG from year 1 and were comparable with PE MG.

CONCLUSION

In this population-based study, MG patients required significantly more healthcare resources in year 1 post-diagnosis than PE MG primarily due to more pharmacological treatments, thymectomies and associated hospitalizations. These findings highlight the need to better understand potential factors including disease characteristics associated with increased health resource use and costs and need for more efficacious treatments early in the disease course.

Tumour Growth Models of Breast Cancer for Evaluating Early Detection-A Summary and a Simulation Study

With the advent of nationwide mammography screening programmes, a number of natural history models of breast cancers have been developed and used to assess the effects of screening. The first half of this article provides an overview of a class of these models and describes how they can be used to study latent processes of tumour progression from observational data. The second half of the article describes a simulation study which applies a continuous growth model to illustrate how effects of extending the maximum age of the current Swedish screening programme from 74 to 80 can be evaluated. Compared to no screening, the current and extended programmes reduced breast cancer mortality by 18.5% and 21.7%, respectively. The proportion of screen-detected invasive cancers which were overdiagnosed was estimated to be 1.9% in the current programme and 2.9% in the extended programme. With the help of these breast cancer natural history models, we can better understand the latent processes, and better study the effects of breast cancer screening.

A unifying framework for continuous tumour growth modelling of breast cancer screening data

The aim of the current article is to present theory that can help unify continuous growth approaches for modelling breast cancer tumour growth based on human data. We present a framework that has three main features: a general likelihood function to account for patient specific screening regiments; stable disease assumptions describing tumour population dynamics; and mathematical models describing tumour growth, individual variation in tumour growth, a hazard for symptomatic detection, and screening test sensitivity. The framework is able to incorporate any random effects distributions for the tumour growth rate parameter, any hazard functions for symptomatic tumour detection, as well as any monotonously increasing function for the tumour growth model. Based on a sample of 1902 incident breast cancer cases with data on mammography screening, we show how the framework can be used to estimate tumour growth based on different growth functions.

Lymph node metastases in breast cancer: Investigating associations with tumor characteristics, molecular subtypes and polygenic risk score using a continuous growth model

We investigate the association between rate of breast cancer lymph node spread and grade, estrogen receptor (ER) status, progesteron receptor status, decision tree derived PAM50 molecular subtype and a polygenic risk score (PRS), using data on 10 950 women included from two different data sources. Lymph node spread was analyzed using a novel continuous tumor progression model that adjusts for tumor volume in a biologically motivated way and that incorporates covariates of interest. Grades 2 and 3 tumors, respectively, were associated with 1.63 and 2.17 times faster rates of lymph node spread than Grade 1 tumors (P < 10^-16 ). ER/PR negative breast cancer was associated with a 1.25/1.19 times faster spread than ER/PR positive breast cancer, respectively (P = .0011 and .0012). Among the molecular subtypes luminal A, luminal B, Her2-enriched and basal-like, Her2-enriched breast cancer was associated with 1.53 times faster spread than luminal A cancer (P = .00072). PRS was not associated with the rate of lymph node spread. Continuous growth models are useful for quantifying associations between lymph node spread and tumor characteristics. These may be useful for building realistic progression models for microsimulation studies used to design individualized screening programs.

Random effects models of lymph node metastases in breast cancer: quantifying the roles of covariates and screening using a continuous growth model

We recently described a joint model of breast cancer tumor size and number of affected lymph nodes, which conditions on screening history, mammographic density, and mode of detection, and can be used to infer growth rates, time to symptomatic detection, screening sensitivity, and rates of lymph node spread. The model of lymph node spread can be estimated in isolation from measurements of tumor volume and number of affected lymph nodes, giving inference identical to the joint model. Here, we extend our model to include covariate effects. We also derive theoretical results in order to study the role of screening on lymph node metastases at diagnosis. We analyze the association between hormone replacement therapy (HRT) and breast cancer lymph node spread, using data from a case-control study designed specifically to study the effects of HRT on breast cancer. Using our method, we estimate that women using HRT at time of diagnosis have a 36% lower rate of lymph node spread than nonusers (95% confidence interval [CI] =(8%,58%)). This can be contrasted with the effect of HRT on the tumor growth rate, estimated here to be 15% slower in HRT users (95% CI = (-34%,+7%)). For screen-detected cancers, we illustrate how lead time can relate to lymph node spread; and using symptomatic cancers, we illustrate the potential consequences of false negative screens in terms of lymph node spread.

Continuous tumour growth models, lead time estimation and length bias in breast cancer screening studies

Comparisons of survival times between screen-detected and symptomatically detected breast cancer cases are subject to lead time and length biases. Whilst the existence of these biases is well known, correction procedures for these are not always clear, as are not the interpretation of these biases. In this paper we derive, based on a recently developed continuous tumour growth model, conditional lead time distributions, using information on each individual's tumour size, screening history and percent mammographic density. We show how these distributions can be used to obtain an individual-based (conditional) procedure for correcting survival comparisons. In stratified analyses, our correction procedure works markedly better than a previously used unconditional lead time correction, based on multi-state Markov modelling. In a study of postmenopausal invasive breast cancer patients, we estimate that, in large (>12 mm) tumours, the multi-state Markov model correction over-corrects five-year survival by 2–3 percentage points. The traditional view of length bias is that tumours being present in a woman's breast for a long time, due to being slow-growing, have a greater chance of being screen-detected. This gives a survival advantage for screening cases which is not due to the earlier detection by screening. We use simulated data to share the new insight that, not only the tumour growth rate but also the symptomatic tumour size will affect the sampling procedure, and thus be a part of the length bias through any link between tumour size and survival. We explain how this has a bearing on how observable breast cancer-specific survival curves should be interpreted. We also propose an approach for correcting survival comparisons for the length bias.

Joint models of tumour size and lymph node spread for incident breast cancer cases in the presence of screening

Continuous growth models show great potential for analysing cancer screening data. We recently described such a model for studying breast cancer tumour growth based on modelling tumour size at diagnosis, as a function of screening history, detection mode, and relevant patient characteristics. In this article, we describe how the approach can be extended to jointly model tumour size and number of lymph node metastases at diagnosis. We propose a new class of lymph node spread models which are biologically motivated and describe how they can be extended to incorporate random effects to allow for heterogeneity in underlying rates of spread. Our final model provides a dramatically better fit to empirical data on 1860 incident breast cancer cases than models in current use. We validate our lymph node spread model on an independent data set consisting of 3961 women diagnosed with invasive breast cancer.

Modelling breast cancer tumour growth for a stable disease population

Statistical models of breast cancer tumour progression have been used to further our knowledge of the natural history of breast cancer, to evaluate mammography screening in terms of mortality, to estimate overdiagnosis, and to estimate the impact of lead-time bias when comparing survival times between screen detected cancers and cancers found outside of screening programs. Multi-state Markov models have been widely used, but several research groups have proposed other modelling frameworks based on specifying an underlying biological continuous tumour growth process. These continuous models offer some advantages over multi-state models and have been used, for example, to quantify screening sensitivity in terms of mammographic density, and to quantify the effect of body size covariates on tumour growth and time to symptomatic detection. As of yet, however, the continuous tumour growth models are not sufficiently developed and require extensive computing to obtain parameter estimates. In this article, we provide a detailed description of the underlying assumptions of the continuous tumour growth model, derive new theoretical results for the model, and show how these results may help the development of this modelling framework. In illustrating the approach, we develop a model for mammography screening sensitivity, using a sample of 1901 post-menopausal women diagnosed with invasive breast cancer.