Transgenic ferret models define pulmonary ionocyte diversity and function.

Speciation leads to adaptive changes in organ cellular physiology and creates challenges for studying rare cell-type functions that diverge between humans and mice. Rare cystic fibrosis transmembrane conductance regulator (CFTR)-rich pulmonary ionocytes exist throughout the cartilaginous airways of humans1,2, but limited presence and divergent biology in the proximal trachea of mice has prevented the use of traditional transgenic models to elucidate ionocyte functions in the airway. Here we describe the creation and use of conditional genetic ferret models to dissect pulmonary ionocyte biology and function by enabling ionocyte lineage tracing (FOXI1-CreERT2::ROSA-TG), ionocyte ablation (FOXI1-KO) and ionocyte-specific deletion of CFTR (FOXI1-CreERT2::CFTRL/L). By comparing these models with cystic fibrosis ferrets3,4, we demonstrate that ionocytes control airway surface liquid absorption, secretion, pH and mucus viscosity-leading to reduced airway surface liquid volume and impaired mucociliary clearance in cystic fibrosis, FOXI1-KO and FOXI1-CreERT2::CFTRL/L ferrets. These processes are regulated by CFTR-dependent ionocyte transport of Cl- and HCO3-. Single-cell transcriptomics and in vivo lineage tracing revealed three subtypes of pulmonary ionocytes and a FOXI1-lineage common rare cell progenitor for ionocytes, tuft cells and neuroendocrine cells during airway development. Thus, rare pulmonary ionocytes perform critical CFTR-dependent functions in the proximal airway that are hallmark features of cystic fibrosis airway disease. These studies provide a road map for using conditional genetics in the first non-rodent mammal to address gene function, cell biology and disease processes that have greater evolutionary conservation between humans and ferrets.

Revisiting the recombinant history of HIV-1 group M with dynamic network community detection.

The prevailing abundance of full-length HIV type 1 (HIV-1) genome sequences provides an opportunity to revisit the standard model of HIV-1 group M (HIV-1/M) diversity that clusters genomes into largely nonrecombinant subtypes, which is not consistent with recent evidence of deep recombinant histories for simian immunodeficiency virus (SIV) and other HIV-1 groups. Here we develop an unsupervised nonparametric clustering approach, which does not rely on predefined nonrecombinant genomes, by adapting a community detection method developed for dynamic social network analysis. We show that this method (dynamic stochastic block model [DSBM]) attains a significantly lower mean error rate in detecting recombinant breakpoints in simulated data (quasibinomial generalized linear model (GLM), P<8×10−8), compared to other reference-free recombination detection programs (genetic algorithm for recombination detection [GARD], recombination detection program 4 [RDP4], and RDP5). When this method was applied to a representative sample of n = 525 actual HIV-1 genomes, we determined k = 29 as the optimal number of DSBM clusters and used change-point detection to estimate that at least 95% of these genomes are recombinant. Further, we identified both known and undocumented recombination hotspots in the HIV-1 genome and evidence of intersubtype recombination in HIV-1 subtype reference genomes. We propose that clusters generated by DSBM can provide an informative framework for HIV-1 classification.

Multiple cohort study of hospitalized SARS-CoV-2 in-host infection dynamics: Parameter estimates, identifiability, sensitivity and the eclipse phase profile.

Within-host SARS-CoV-2 modelling studies have been published throughout the COVID-19 pandemic. These studies contain highly variable numbers of individuals and capture varying timescales of pathogen dynamics; some studies capture the time of disease onset, the peak viral load and subsequent heterogeneity in clearance dynamics across individuals, while others capture late-time post-peak dynamics. In this study, we curate multiple previously published SARS-CoV-2 viral load data sets, fit these data with a consistent modelling approach, and estimate the variability of in-host parameters including the basic reproduction number, R0, as well as the best-fit eclipse phase profile. We find that fitted dynamics can be highly variable across data sets, and highly variable within data sets, particularly when key components of the dynamic trajectories (e.g. peak viral load) are not represented in the data. Further, we investigated the role of the eclipse phase time distribution in fitting SARS-CoV-2 viral load data. By varying the shape parameter of an Erlang distribution, we demonstrate that models with either no eclipse phase, or with an exponentially-distributed eclipse phase, offer significantly worse fits to these data, whereas models with less dispersion around the mean eclipse time (shape parameter two or more) offered the best fits to the available data across all data sets used in this work. This manuscript was submitted as part of a theme issue on "Modelling COVID-19 and Preparedness for Future Pandemics".

Current benchtop protocols are not appropriate for the evaluation of distraction-based growing rods: a literature review to justify a new protocol and its development.

PURPOSE: Although distraction-based growing rods (GR) are the gold standard for the treatment of early onset scoliosis, they suffer from high failure rates. We have (1) performed a literature search to understand the deficiencies of the current protocols, (2) in vitro evaluation of GRs using our proposed protocol and performed a finite element (FE) model validation, and (3) identified key features which should be considered in mechanical testing setups. METHODS: PubMed, Embase, and Web of Science databases were searched for articles published on (a) in vivo animal, in vitro cadaveric, and biomechanical studies analyzing the use of GRs as well as (b) failure mechanisms and risk factors for GRs. Both FE and benchtop models of a proposed TGR test construct were developed and evaluated for two cases, long tandem connectors (LT), and side-by-side connectors (SBS). The test construct consisted of five polymer blocks representing vertebral bodies, joined with springs to simulate spinal stiffness. The superior and inferior blocks accepted the pedicle screw anchors, while the three middle blocks were floating. After the pedicle screws, rods, and connectors were assembled onto this construct, distraction was performed, mimicking scoliosis surgery. The resulting distracted constructs were then subjected to static compression-bending loading. Yield load and stiffness were calculated and used to verify/validate the FE results. RESULTS: From the literature search, key features identified as significant were axial and transverse connectors, contoured rods, and distraction, distraction being the most challenging feature to incorporate in testing. The in silico analyses, once they are validated, can be used as a complementing technique to investigate other anatomical features which are not possible in the mechanical setup (like growth/scoliosis curvature). Based on our experiment, the LT constructs showed higher stiffness and yield load compared to SBS (78.85 N/mm vs. 59.68 N/mm and 838.84 N vs. 623.3 N). The FE predictions were in agreement with the experimental outcomes (within 10% difference). The maximum von Mises stresses were predicted adjacent to the distraction site, consistent with the location of observed failures in vivo. CONCLUSION: The two-way approach presented in this study can lead to a robust prediction of the contributing factors to the in vivo failure.

Goods, causes and intentions: problems with applying the doctrine of double effect to palliative sedation.

BACKGROUND: Palliative sedation and analgesia are employed in patients with refractory and intractable symptoms at the end of life to reduce their suffering by lowering their level of consciousness. The doctrine of double effect, a philosophical principle that justifies doing a "good action" with a potentially "bad effect," is frequently employed to provide an ethical justification for this practice. MAIN TEXT: We argue that palliative sedation and analgesia do not fulfill the conditions required to apply the doctrine of double effect, and therefore its use in this domain is inappropriate. Furthermore, we argue that the frequent application of the doctrine of double effect to palliative sedation and analgesia reflects physicians' discomfort with the complex moral, intentional, and causal aspects of end-of-life care. CONCLUSIONS: We are concerned that this misapplication of the doctrine of double effect can consequently impair physicians' ethical reasoning and relationships with patients at the end of life.

An investigation of listener envelopment utilizing a spherical microphone array and third-order ambisonics reproduction.

Listener envelopment (LEV), the sense of being surrounded by the sound field, is a perception that has been found to be related to the overall impression of a concert hall. The purpose of this study was to investigate the relationship between the perception of LEV and the direction and arrival time of energy from spatial room impulse responses (IRs). IRs were obtained in a 2000-seat concert hall using a 32-channel spherical microphone array and analyzed using a third-order plane wave decomposition. Additionally, the IRs were convolved with anechoic music and processed for third-order Ambisonic reproductions and presented to subjects over a 30-loudspeaker array. Instances were found in which the energy in the late sound field did not correlate with LEV ratings as well as energy in a 70-100 ms time window. Follow-up listening tests were conducted with hybrid IRs containing portions of an enveloping IR and an unenveloping IR with crossover times ranging from 40 to 140 ms. Additional hybrid IRs were studied wherein portions of the spatial IRs were collapsed into all frontal energy with crossover times ranging from 40 to 120 ms. The tests confirmed that much of the important LEV information exists in the early portion of these IRs.

Mutations in the SPG7 gene cause chronic progressive external ophthalmoplegia through disordered mitochondrial DNA maintenance.

Despite being a canonical presenting feature of mitochondrial disease, the genetic basis of progressive external ophthalmoplegia remains unknown in a large proportion of patients. Here we show that mutations in SPG7 are a novel cause of progressive external ophthalmoplegia associated with multiple mitochondrial DNA deletions. After excluding known causes, whole exome sequencing, targeted Sanger sequencing and multiplex ligation-dependent probe amplification analysis were used to study 68 adult patients with progressive external ophthalmoplegia either with or without multiple mitochondrial DNA deletions in skeletal muscle. Nine patients (eight probands) were found to carry compound heterozygous SPG7 mutations, including three novel mutations: two missense mutations c.2221G>A; p.(Glu741Lys), c.2224G>A; p.(Asp742Asn), a truncating mutation c.861dupT; p.Asn288*, and seven previously reported mutations. We identified a further six patients with single heterozygous mutations in SPG7, including two further novel mutations: c.184-3C>T (predicted to remove a splice site before exon 2) and c.1067C>T; p.(Thr356Met). The clinical phenotype typically developed in mid-adult life with either progressive external ophthalmoplegia/ptosis and spastic ataxia, or a progressive ataxic disorder. Dysphagia and proximal myopathy were common, but urinary symptoms were rare, despite the spasticity. Functional studies included transcript analysis, proteomics, mitochondrial network analysis, single fibre mitochondrial DNA analysis and deep re-sequencing of mitochondrial DNA. SPG7 mutations caused increased mitochondrial biogenesis in patient muscle, and mitochondrial fusion in patient fibroblasts associated with the clonal expansion of mitochondrial DNA mutations. In conclusion, the SPG7 gene should be screened in patients in whom a disorder of mitochondrial DNA maintenance is suspected when spastic ataxia is prominent. The complex neurological phenotype is likely a result of the clonal expansion of secondary mitochondrial DNA mutations modulating the phenotype, driven by compensatory mitochondrial biogenesis.

SARS-CoV-2 booster vaccine dose significantly extends humoral immune response half-life beyond the primary series.

SARS-CoV-2 lipid nanoparticle mRNA vaccines continue to be administered as the predominant prophylactic measure to reduce COVID-19 disease pathogenesis. Quantifying the kinetics of the secondary immune response from subsequent doses beyond the primary series and understanding how dose-dependent immune waning kinetics vary as a function of age, sex, and various comorbidities remains an important question. We study anti-spike IgG waning kinetics in 152 individuals who received an mRNA-based primary series (first two doses) and a subset of 137 individuals who then received an mRNA-based booster dose. We find the booster dose elicits a 71-84% increase in the median Anti-S half life over that of the primary series. We find the Anti-S half life for both primary series and booster doses decreases with age. However, we stress that although chronological age continues to be a good proxy for vaccine-induced humoral waning, immunosenescence is likely not the mechanism, rather, more likely the mechanism is related to the presence of noncommunicable diseases, which also accumulate with age, that affect immune regulation. We are able to independently reproduce recent observations that those with pre-existing asthma exhibit a stronger primary series humoral response to vaccination than compared to those that do not, and further, we find this result is sustained for the booster dose. Finally, via a single-variate Kruskal-Wallis test we find no difference between male and female humoral decay kinetics, however, a multivariate approach utilizing  Least Absolute Shrinkage and Selection Operator (LASSO) regression for feature selection reveals a statistically significant (p < 1 × 10 - 3 ), albeit small, bias in favour of longer-lasting humoral immunity amongst males.

Mucociliary Clearance is Impaired in Small Airways of Cystic Fibrosis Pigs.

Rationale: Cystic fibrosis is a genetic disorder characterized by recurrent airway infections, inflammation, and progressive decline in lung function. Autopsy and spirometry data suggest that cystic fibrosis may start in the small airways which, due to the fractal nature of the airways, account for most of the airway tree surface area. However, they are not easily accessible for testing. Objectives: Here, we tested the hypothesis that mucociliary clearance is abnormal in the small airways of newborn cystic fibrosis pigs. Methods: Current mucociliary clearance assays are limited therefore we developed a dynamic positron emission tomography scan assay with high spatial and temporal resolution. Each study was accompanied by a high-resolution computed tomography scan that helped identify the thin outer region of the lung that contained small airways. Measurements and Main Results: Clearance of aerosolized [ 68 Ga]macro aggregated albumin from distal airways occurred within minutes after delivery and followed a two-phase process. In cystic fibrosis pigs, both early and late clearance rates were slower. Stimulation of the cystic fibrosis airways with the purinergic agonist UTP further impaired late clearance. Only 1 cystic fibrosis pig treated with UTP out of 6 cleared more than 20% of the delivered dose. Conclusions: These data indicate that mucociliary transport in the small airways is fast and can easily be missed if the acquisition is not fast enough. The data also indicate that mucociliary transport is impaired in small airways of cystic fibrosis pigs. This defect is exacerbated by stimulation of mucus secretions with purinergic agonists.

Titin mutation segregates with hereditary myopathy with early respiratory failure.

In 2001, we described an autosomal dominant myopathy characterized by neuromuscular ventilatory failure in ambulant patients. Here we describe the underlying genetic basis for the disorder, and we define the neuromuscular, respiratory and radiological phenotype in a study of 31 mutation carriers followed for up to 31 years. A combination of genome-wide linkage and whole exome sequencing revealed the likely causal genetic variant in the titin (TTN) gene (g.274375T>C; p.Cys30071Arg) within a shared haplotype of 2.93 Mbp on chromosome 2. This segregated with the phenotype in 21 individuals from the original family, nine subjects in a second family with the same highly selective pattern of muscle involvement on magnetic resonance imaging and a third familial case with a similar phenotype. Comparing the mutation carriers revealed novel features not apparent in our original report. The clinical presentation included predominant distal, proximal or respiratory muscle weakness. The age of onset was highly variable, from early adulthood, and including a mild phenotype in advanced age. Muscle weakness was earlier onset and more severe in the lower extremities in nearly all patients. Seven patients also had axial muscle weakness. Respiratory function studies demonstrated a gradual deterioration over time, reflecting the progressive nature of this condition. Cardiomyopathy was not present in any of our patients despite up to 31 years of follow-up. Magnetic resonance muscle imaging was performed in 21 affected patients and revealed characteristic abnormalities with semitendinosus involvement in 20 of 21 patients studied, including 3 patients who were presymptomatic. Diagnostic muscle histopathology most frequently revealed eosinophilic inclusions (inclusion bodies) and rimmed vacuoles, but was non-specific in a minority of patients. These findings have important clinical implications. This disease should be considered in patients with adult-onset proximal or distal myopathy and early respiratory failure, even in the presence of non-specific muscle pathology. Muscle magnetic resonance imaging findings are characteristic and should be considered as an initial investigation, and if positive should prompt screening for mutations in TTN. With 363 exons, screening TTN presented a major challenge until recently. However, whole exome sequencing provides a reliable cost-effective approach, providing the gene of interest is adequately captured.

Canadian Covid-19 pandemic public health mitigation measures at the province level.

The Covid-19 pandemic has prompted governments across the world to enforce a range of public health interventions. We introduce the Covid-19 Policy Response Canadian tracker (CPRCT) database that tracks and records implemented public health measures in every province and territory in Canada. The implementations are recorded on a four-level ordinal scale (0-3) for three domains, (Schools, Work, and Other), capturing differences in degree of response. The data-set allows the exploration of the effects of public health mitigation on the spread of Covid-19, as well as provides a near-real-time record in an accessible format that is useful for a diverse range of modeling and research questions.

A mathematical model of protein subunits COVID-19 vaccines.

We consider a general mathematical model for protein subunit vaccine with a focus on the MF59-adjuvanted spike glycoprotein-clamp vaccine for SARS-CoV-2, and use the model to study immunological outcomes in the humoral and cell-mediated arms of the immune response from vaccination. The mathematical model is fit to vaccine clinical trial data. We elucidate the role of Interferon-γ and Interleukin-4 in stimulating the immune response of the host. Model results, and results from a sensitivity analysis, show that a balance between the TH1 and TH2 arms of the immune response is struck, with the TH1 response being dominant. The model predicts that two-doses of the vaccine at 28 days apart will result in approximately 85% humoral immunity loss relative to peak immunity approximately 6 months post dose 1.

Acceptance and Commitment Therapy for people living with motor neuron disease: an uncontrolled feasibility study.

BACKGROUND: Motor neuron disease (MND) is a fatal, progressive neurodegenerative disease that causes progressive weakening and wasting of limb, bulbar, thoracic and abdominal muscles. Clear evidence-based guidance on how psychological distress should be managed in people living with MND (plwMND) is lacking. Acceptance and Commitment Therapy (ACT) is a form of psychological therapy that may be particularly suitable for this population. However, to the authors' knowledge, no study to date has evaluated ACT for plwMND. Consequently, the primary aim of this uncontrolled feasibility study was to examine the feasibility and acceptability of ACT for improving the psychological health of plwMND. METHODS: PlwMND aged ≥ 18 years were recruited from 10 UK MND Care Centres/Clinics. Participants received up to 8 one-to-one ACT sessions, developed specifically for plwMND, plus usual care. Co-primary feasibility and acceptability outcomes were uptake (≥ 80% of the target sample [N = 28] recruited) and initial engagement with the intervention (≥ 70% completing ≥ 2 sessions). Secondary outcomes included measures of quality of life, anxiety, depression, disease-related functioning, health status and psychological flexibility in plwMND and quality of life and burden in caregivers. Outcomes were assessed at baseline and 6 months. RESULTS: Both a priori indicators of success were met: 29 plwMND (104%) were recruited and 76% (22/29) attended ≥ 2 sessions. Attrition at 6-months was higher than anticipated (8/29, 28%), but only two dropouts were due to lack of acceptability of the intervention. Acceptability was further supported by good satisfaction with therapy and session attendance. Data were possibly suggestive of small improvements in anxiety and psychological quality of life from baseline to 6 months in plwMND, despite a small but expected deterioration in disease-related functioning and health status. CONCLUSIONS: There was good evidence of acceptability and feasibility. Limitations included the lack of a control group and small sample size, which complicate interpretation of findings. A fully powered RCT to evaluate the clinical and cost-effectiveness of ACT for plwMND is underway. TRIAL REGISTRATION: The study was pre-registered with the ISRCTN Registry (ISRCTN12655391).

Prospective comparison of combined 18F-FDG and 18F-NaF PET/CT vs. 18F-FDG PET/CT imaging for detection of malignancy.

PURPOSE: Typically, (18)F-FDG PET/CT and (18)F-NaF PET/CT scans are done as two separate studies on different days to allow sufficient time for the radiopharmaceutical from the first study to decay. This is inconvenient for the patients and exposes them to two doses of radiation from the CT component of the examinations. In the current study, we compared the clinical usefulness of a combined (18)F-FDG/(18)F-NaF PET/CT scan with that of a separate (18)F-FDG-only PET/CT scan. METHODS: There were 62 patients enrolled in this prospective trial. All had both an (18)F-FDG-alone PET/CT scan and a combined (18)F-FDG/(18)F-NaF PET/CT scan. Of the 62 patients, 53 (85%) received simultaneous tracer injections, while 9 (15%) received (18)F-NaF subsequent to the initial (18)F-FDG dose (average delay 2.2 h). Images were independently reviewed for PET findings by two Board-Certified nuclear medicine physicians, with discrepancies resolved by a third reader. Interpreters were instructed to only report findings that were concerning for malignancy. Reading the (18)F-FDG-only scan first for half of the patients controlled for order bias. RESULTS: In 15 of the 62 patients (24%) neither the (18)F-FDG-only PET/CT scan nor the combined (18)F-FDG/(18)F-NaF PET/CT scan identified malignancy. In the remaining 47 patients who had PET findings of malignancy, a greater number of lesions were detected in 16 of 47 patients (34%) using the combined (18)F-FDG/(18)F-NaF PET/CT scan compared to the (18)F-FDG-only PET/CT scan. In 2 of these 47 patients (4%), the (18)F-FDG-only scan demonstrated soft tissue lesions that were not prospectively identified on the combined study. In 29 of these 47 patients (62%), the combined scan detected an equal number of lesions compared to the (18)F-FDG-only scan. Overall, 60 of all the 62 patients (97%) showed an equal or greater number of lesions on the combined scan than on the (18)F-FDG-only scan. CONCLUSION: The current study demonstrated that (18)F-FDG and (18)F-NaF can be combined in a single PET/CT scan by administering the two radiopharmaceuticals simultaneously or in sequence on the same day. In addition to patient convenience and reduced radiation exposure from the CT component, the combined (18)F-FDG/(18)F-NaF PET/CT scan appeared to increase the sensitivity for detection of osseous lesions compared to the (18)F-FDG-only PET/CT scan in the studied population.

Modeling waning and boosting of COVID-19 in Canada with vaccination.

SARS-CoV-2, the causative agent of COVID-19, has caused devastating health and economic impacts around the globe since its appearance in late 2019. The advent of effective vaccines leads to open questions on how best to vaccinate the population. To address such questions, we developed a model of COVID-19 infection by age that includes the waning and boosting of immunity against SARS-CoV-2 in the context of infection and vaccination. The model also accounts for changes to infectivity of the virus, such as public health mitigation protocols over time, increases in the transmissibility of variants of concern, changes in compliance to mask wearing and social distancing, and changes in testing rates. The model is employed to study public health mitigation and vaccination of the COVID-19 epidemic in Canada, including different vaccination programs (rollout by age), and delays between doses in a two-dose vaccine. We find that the decision to delay the second dose of vaccine is appropriate in the Canadian context. We also find that the benefits of a COVID-19 vaccination program in terms of reductions in infections is increased if vaccination of 15-19 year olds are included in the vaccine rollout.

Long-term durability of immune responses to the BNT162b2 and mRNA-1273 vaccines based on dosage, age and sex.

The lipid nanoparticle (LNP)-formulated mRNA vaccines BNT162b2 and mRNA-1273 are a widely adopted multi vaccination public health strategy to manage the COVID-19 pandemic. Clinical trial data has described the immunogenicity of the vaccine, albeit within a limited study time frame. Here, we use a within-host mathematical model for LNP-formulated mRNA vaccines, informed by available clinical trial data from 2020 to September 2021, to project a longer term understanding of immunity as a function of vaccine type, dosage amount, age, and sex. We estimate that two standard doses of either mRNA-1273 or BNT162b2, with dosage times separated by the company-mandated intervals, results in individuals losing more than 99% humoral immunity relative to peak immunity by 8 months following the second dose. We predict that within an 8 month period following dose two (corresponding to the original CDC time-frame for administration of a third dose), there exists a period of time longer than 1 month where an individual has lost more than 99% humoral immunity relative to peak immunity, regardless of which vaccine was administered. We further find that age has a strong influence in maintaining humoral immunity; by 8 months following dose two we predict that individuals aged 18-55 have a four-fold humoral advantage compared to aged 56-70 and 70+ individuals. We find that sex has little effect on the immune response and long-term IgG counts. Finally, we find that humoral immunity generated from two low doses of mRNA-1273 decays at a substantially slower rate relative to peak immunity gained compared to two standard doses of either mRNA-1273 or BNT162b2. Our predictions highlight the importance of the recommended third booster dose in order to maintain elevated levels of antibodies.

COVID-19 Seroprevalence in Canada Modelling Waning and Boosting COVID-19 Immunity in Canada a Canadian Immunization Research Network Study.

COVID-19 seroprevalence changes over time, with infection, vaccination, and waning immunity. Seroprevalence estimates are needed to determine when increased COVID-19 vaccination coverage is needed, and when booster doses should be considered, to reduce the spread and disease severity of COVID-19 infection. We use an age-structured model including infection, vaccination and waning immunity to estimate the distribution of immunity to COVID-19 in the Canadian population. This is the first mathematical model to do so. We estimate that 60-80% of the Canadian population has some immunity to COVID-19 by late Summer 2021, depending on specific characteristics of the vaccine and the waning rate of immunity. Models results indicate that increased vaccination uptake in age groups 12-29, and booster doses in age group 50+ are needed to reduce the severity COVID-19 Fall 2021 resurgence.

Clinically relevant finite element technique based protocol to evaluate growing rods for early onset scoliosis correction.

OBJECTIVE: The emergence of distraction-based growing rods has provided the means to reduce the progression of spinal deformity in early onset scoliosis (EOS). The current protocols for evaluating spinal implants (ie, ASTM-F1717 and ISO-12189) were developed for fusion/dynamic devices. These protocols do not feature long unsupported rod lengths subjected to distraction. Due to the unsuitability of the existing guidelines for the evaluation of growing rods, a new distraction-based finite element protocol is presented herein for the first time. METHOD: A vertebrectomy (VO) model from current protocols was modified to accommodate multi-spinal segments (ie, MS model) in which springs with appropriate stiffness were simulated in between the plastic blocks. To assess the efficacy of the protocol, two different computational studies were conducted: (a) compression-bending (MS-CB) with no distraction, and (b) distraction followed by compression-bending (MS-D + CB). In each study, the model with no axial connector (rods-only) was modified to include a) 80-mm long tandem (LT) connectors, and b) side-by-side (SBS) connectors. Stiffness and yield loads were calculated as per ASTM-F1717 guidelines and compared with the corresponding VO models with no distraction. In the MS-D + CB models, distraction was applied at the top block, stretching the spring-block construct in a simulation of scoliosis surgery prior to locking the construct at the connector-rods' interface. RESULTS: MS-CB models predicted higher stiffness and yield loads, compared to the VO models. The locking mechanism produced pre-existing stresses on the rod-connector interface, which caused a shift in the location of high-stress regions to the distraction site. Distraction led to a decrease in the construct's stiffness and yield load. DISCUSSION: The proposed protocol enables the simulation of clinical parameters that are not feasible in the F1717 models and predicted stress patterns in the hardware consistent with observed clinical failures.