CellDrift: inferring perturbation responses in temporally sampled single-cell data.

Cells and tissues respond to perturbations in multiple ways that can be sensitively reflected in the alterations of gene expression. Current approaches to finding and quantifying the effects of perturbations on cell-level responses over time disregard the temporal consistency of identifiable gene programs. To leverage the occurrence of these patterns for perturbation analyses, we developed CellDrift (https://github.com/KANG-BIOINFO/CellDrift), a generalized linear model-based functional data analysis method that is capable of identifying covarying temporal patterns of various cell types in response to perturbations. As compared to several other approaches, CellDrift demonstrated superior performance in the identification of temporally varied perturbation patterns and the ability to impute missing time points. We applied CellDrift to multiple longitudinal datasets, including COVID-19 disease progression and gastrointestinal tract development, and demonstrated its ability to identify specific gene programs associated with sequential biological processes, trajectories and outcomes.

Predicting lung function decline in cystic fibrosis: the impact of initiating ivacaftor therapy.

BACKGROUND: Modulator therapies that seek to correct the underlying defect in cystic fibrosis (CF) have revolutionized the clinical landscape. Given the heterogeneous nature of lung disease progression in the post-modulator era, there is a need to develop prediction models that are robust to modulator uptake. METHODS: We conducted a retrospective longitudinal cohort study of the CF Foundation Patient Registry (N = 867 patients carrying the G551D mutation who were treated with ivacaftor from 2003 to 2018). The primary outcome was lung function (percent predicted forced expiratory volume in 1 s or FEV1pp). To characterize the association between ivacaftor initiation and lung function, we developed a dynamic prediction model through covariate selection of demographic and clinical characteristics. The ability of the selected model to predict a decline in lung function, clinically known as an FEV1-indicated exacerbation signal (FIES), was evaluated both at the population level and individual level. RESULTS: Based on the final model, the estimated improvement in FEV1pp after ivacaftor initiation was 4.89% predicted (95% confidence interval [CI]: 3.90 to 5.89). The rate of decline was reduced with ivacaftor initiation by 0.14% predicted/year (95% CI: 0.01 to 0.27). More frequent outpatient visits prior to study entry and being male corresponded to a higher overall FEV1pp. Pancreatic insufficiency, older age at study entry, a history of more frequent pulmonary exacerbations, lung infections, CF-related diabetes, and use of Medicaid insurance corresponded to lower FEV1pp. The model had excellent predictive accuracy for FIES events with an area under the receiver operating characteristic curve of 0.83 (95% CI: 0.83 to 0.84) for the independent testing cohort and 0.90 (95% CI: 0.89 to 0.90) for 6-month forecasting with the masked cohort. The root-mean-square errors of the FEV1pp predictions for these cohorts were 7.31% and 6.78% predicted, respectively, with standard deviations of 0.29 and 0.20. The predictive accuracy was robust across different covariate specifications. CONCLUSIONS: The methods and applications of dynamic prediction models developed using data prior to modulator uptake have the potential to inform post-modulator projections of lung function and enhance clinical surveillance in the new era of CF care.

Bayesian two-stage modeling of longitudinal and time-to-event data with an integrated fractional Brownian motion covariance structure.

It is difficult to characterize complex variations of biological processes, often longitudinally measured using biomarkers that yield noisy data. While joint modeling with a longitudinal submodel for the biomarker measurements and a survival submodel for assessing the hazard of events can alleviate measurement error issues, the continuous longitudinal submodel often uses random intercepts and slopes to estimate both between- and within-patient heterogeneity in biomarker trajectories. To overcome longitudinal submodel challenges, we replace random slopes with scaled integrated fractional Brownian motion (IFBM). As a more generalized version of integrated Brownian motion, IFBM reasonably depicts noisily measured biological processes. From this longitudinal IFBM model, we derive novel target functions to monitor the risk of rapid disease progression as real-time predictive probabilities. Predicted biomarker values from the IFBM submodel are used as inputs in a Cox submodel to estimate event hazard. This two-stage approach to fit the submodels is performed via Bayesian posterior computation and inference. We use the proposed approach to predict dynamic lung disease progression and mortality in women with a rare disease called lymphangioleiomyomatosis who were followed in a national patient registry. We compare our approach to those using integrated Ornstein-Uhlenbeck or conventional random intercepts-and-slopes terms for the longitudinal submodel. In the comparative analysis, the IFBM model consistently demonstrated superior predictive performance.

Automated, machine learning-based alerts increase epilepsy surgery referrals: A randomized controlled trial.

OBJECTIVE: To determine whether automated, electronic alerts increased referrals for epilepsy surgery. METHODS: We conducted a prospective, randomized controlled trial of a natural language processing-based clinical decision support system embedded in the electronic health record (EHR) at 14 pediatric neurology outpatient clinic sites. Children with epilepsy and at least two prior neurology visits were screened by the system prior to their scheduled visit. Patients classified as a potential surgical candidate were randomized 2:1 for their provider to receive an alert or standard of care (no alert). The primary outcome was referral for a neurosurgical evaluation. The likelihood of referral was estimated using a Cox proportional hazards regression model. RESULTS: Between April 2017 and April 2019, at total of 4858 children were screened by the system, and 284 (5.8%) were identified as potential surgical candidates. Two hundred four patients received an alert, and 96 patients received standard care. Median follow-up time was 24 months (range: 12-36 months). Compared to the control group, patients whose provider received an alert were more likely to be referred for a presurgical evaluation (3.1% vs 9.8%; adjusted hazard ratio [HR] = 3.21, 95% confidence interval [CI]: 0.95-10.8; one-sided p = .03). Nine patients (4.4%) in the alert group underwent epilepsy surgery, compared to none (0%) in the control group (one-sided p = .03). SIGNIFICANCE: Machine learning-based automated alerts may improve the utilization of referrals for epilepsy surgery evaluations.

Bayesian causal inference for observational studies with missingness in covariates and outcomes.

Missing data are a pervasive issue in observational studies using electronic health records or patient registries. It presents unique challenges for statistical inference, especially causal inference. Inappropriately handling missing data in causal inference could potentially bias causal estimation. Besides missing data problems, observational health data structures typically have mixed-type variables - continuous and categorical covariates - whose joint distribution is often too complex to be modeled by simple parametric models. The existence of missing values in covariates and outcomes makes the causal inference even more challenging, while most standard causal inference approaches assume fully observed data or start their works after imputing missing values in a separate preprocessing stage. To address these problems, we introduce a Bayesian nonparametric causal model to estimate causal effects with missing data. The proposed approach can simultaneously impute missing values, account for multiple outcomes, and estimate causal effects under the potential outcomes framework. We provide three simulation studies to show the performance of our proposed method under complicated data settings whose features are similar to our case studies. For example, Simulation Study 3 assumes the case where missing values exist in both outcomes and covariates. Two case studies were conducted applying our method to evaluate the comparative effectiveness of treatments for chronic disease management in juvenile idiopathic arthritis and cystic fibrosis.

Multilevel joint model of longitudinal continuous and binary outcomes for hierarchically structured data.

Joint modeling has been a useful strategy for incorporating latent associations between different types of outcomes simultaneously, often focusing on a longitudinal continuous outcome characterized by an LME submodel and a terminal event subject to a Cox proportional hazard or parametric survival submodel. Applications to hierarchical longitudinal studies have been less frequent, particularly with respect to a binary process, which is commonly specified by a GLMM. Furthermore, many of the joint model developments have not allowed for investigations of nested effects, such as those arising from multicenter studies. To fill this gap, we propose a multilevel joint model that encompasses the LME submodel and GLMM through a Bayesian approach. Motivated by the need for timely detection of pulmonary exacerbation and characterization of irregularly observed lung function measurements in people living with cystic fibrosis (CF) receiving care across multiple centers, we apply the model to the data arising from US CF Foundation Patient Registry. In parallel, we examine the extent of bias induced by a non-hierarchical model. Our simulation study and application results show that incorporating the center effect along with individual stochastic variation over time within the LME submodel improves model estimation and prediction. Given that the center effect is evident in lung function observed in the CF population, accounting for center-specific power parameters by incorporating the symmetric power exponential power (spep) link function in the GLMM can facilitate more accurate conclusions in clinical studies.

Lung function in children with cystic fibrosis in the USA and UK: a comparative longitudinal analysis of national registry data.

RATIONALE: A previous analysis found significantly higher lung function in the US paediatric cystic fibrosis (CF) population compared with the UK with this difference apparently decreasing in adolescence and adulthood. However, the cross-sectional nature of the study makes it hard to interpret these results. OBJECTIVES: To compare longitudinal trajectories of lung function in children with CF between the USA and UK and to explore reasons for any differences. METHODS: We used mixed effects regression analysis to model lung function trajectories in the study populations. Using descriptive statistics, we compared early growth and nutrition (height, weight, body mass index), infections (Pseudomonas aeruginosa, Staphylococcus aureus) and treatments (rhDnase, hypertonic saline, inhaled antibiotics). RESULTS: We included 9463 children from the USA and 3055 children from the UK with homozygous F508del genotype. Lung function was higher in the USA than in the UK when first measured at age six and remained higher throughout childhood. We did not find important differences in early growth and nutrition, or P.aeruginosa infection. Prescription of rhDNase and hypertonic saline was more common in the USA. Inhaled antibiotics were prescribed at similar levels in both countries, but Tobramycin was prescribed more in the USA and colistin in the UK. S. aureus infection was more common in the USA than the UK. CONCLUSIONS: Children with CF and homozygous F508del genotype in the USA had better lung function than UK children. These differences do not appear to be explained by early growth or nutrition, but differences in the use of early treatments need further investigation.

Projecting the impact of triple CFTR modulator therapy on intravenous antibiotic requirements in cystic fibrosis using patient registry data combined with treatment effects from randomised trials.

BACKGROUND: Cystic fibrosis (CF) is a life-threatening genetic disease, affecting around 10 500 people in the UK. Precision medicines have been developed to treat specific CF-gene mutations. The newest, elexacaftor/tezacaftor/ivacaftor (ELEX/TEZ/IVA), has been found to be highly effective in randomised controlled trials (RCTs) and became available to a large proportion of UK CF patients in 2020. Understanding the potential health economic impacts of ELEX/TEZ/IVA is vital to planning service provision. METHODS: We combined observational UK CF Registry data with RCT results to project the impact of ELEX/TEZ/IVA on total days of intravenous (IV) antibiotic treatment at a population level. Registry data from 2015 to 2017 were used to develop prediction models for IV days over a 1-year period using several predictors, and to estimate 1-year population total IV days based on standards of care pre-ELEX/TEZ/IVA. We considered two approaches to imposing the impact of ELEX/TEZ/IVA on projected outcomes using effect estimates from RCTs: approach 1 based on effect estimates on FEV1% and approach 2 based on effect estimates on exacerbation rate. RESULTS: ELEX/TEZ/IVA is expected to result in significant reductions in population-level requirements for IV antibiotics of 16.1% (~17 800 days) using approach 1 and 43.6% (~39 500 days) using approach 2. The two approaches require different assumptions. Increased understanding of the mechanisms through which ELEX/TEZ/IVA acts on these outcomes would enable further refinements to our projections. CONCLUSIONS: This work contributes to increased understanding of the changing healthcare needs of people with CF and illustrates how Registry data can be used in combination with RCT evidence to estimate population-level treatment impacts.

Bayesian regularization for a nonstationary Gaussian linear mixed effects model.

In omics experiments, estimation and variable selection can involve thousands of proteins/genes observed from a relatively small number of subjects. Many regression regularization procedures have been developed for estimation and variable selection in such high-dimensional problems. However, approaches have predominantly focused on linear regression models that ignore correlation arising from long sequences of repeated measurements on the outcome. Our work is motivated by the need to identify proteomic biomarkers that improve the prediction of rapid lung-function decline for individuals with cystic fibrosis (CF) lung disease. We extend four Bayesian penalized regression approaches for a Gaussian linear mixed effects model with nonstationary covariance structure to account for the complicated structure of longitudinal lung function data while simultaneously estimating unknown parameters and selecting important protein isoforms to improve predictive performance. Different types of shrinkage priors are evaluated to induce variable selection in a fully Bayesian framework. The approaches are studied with simulations. We apply the proposed method to real proteomics and lung-function outcome data from our motivating CF study, identifying a set of relevant clinical/demographic predictors and a proteomic biomarker for rapid decline of lung function. We also illustrate the methods on CD4 yeast cell-cycle genomic data, confirming that the proposed method identifies transcription factors that have been highlighted in the literature for their importance as cell cycle transcription factors.

Aligning Provider Prescribing With Guidelines for Soft Tissue Infections.

OBJECTIVE: Despite evidence-based guidelines, antibiotics prescribed for uncomplicated skin and soft tissue infections can involve inappropriate microbial coverage. Our aim was to evaluate the appropriateness of antibiotic prescribing practices for mild nonpurulent cellulitis in a pediatric tertiary academic medical center over a 1-year period. METHODS: Eligible patients treated in the emergency department or urgent care settings for mild nonpurulent cellulitis from January 2017 to December 2017 were identified by an International Classification of Diseases, Tenth Revision, code for cellulitis. The primary outcome was appropriateness of prescribed antibiotics as delineated by adherence with the Infectious Diseases Society of America guidelines. Secondary outcomes include reutilization rate as defined by revisit to the emergency department/urgent cares within 14 days of the initial encounter. RESULTS: A total of 967 encounters were evaluated with 60.0% overall having guideline-adherent care. Common reasons for nonadherence included inappropriate coverage of MRSA with clindamycin (n = 217, 56.1%) and single-agent coverage with sulfamethoxazole-trimethoprim (n = 129, 33.3%). There were 29 revisits within 14 days of initial patient encounters or a reutilization rate of 3.0%, which was not significantly associated with the Infectious Diseases Society of America adherence. CONCLUSIONS: Our data show antibiotic prescription for nonpurulent cellulitis as a potential area of standardization and optimization of care at our center.

Flexible link functions in a joint hierarchical Gaussian process model.

Many longitudinal studies often require jointly modeling a biomarker and an event outcome, in order to provide more accurate inference and dynamic prediction of disease progression. Cystic fibrosis (CF) studies have illustrated the benefits of these models, primarily examining the joint evolution of lung-function decline and survival. We propose a novel joint model within the shared-parameter framework that accommodates nonlinear lung-function trajectories, in order to provide more accurate inference on lung-function decline over time and to examine the association between evolution of lung function and risk of a pulmonary exacerbation (PE) event recurrence. Specifically, a two-level Gaussian process (GP) is used to estimate the nonlinear longitudinal trajectories and a flexible link function is introduced for a more accurate depiction of the binary process on the event outcome. Bayesian model assessment is used to evaluate each component of the joint model in simulation studies and an application to longitudinal data on patients receiving care from a CF center. A nonlinear structure is suggested by both longitudinal continuous and binary evaluations. Including a flexible link function improves model fit to these data. The proposed hierarchical GP model with a flexible power link function where Laplace distribution is the baseline (spep) has the best fit of all joint models considered, characterizing how accelerated lung-function decline corresponds to increased odds of experiencing another PE.

Flexible multivariate joint model of longitudinal intensity and binary process for medical monitoring of frequently collected data.

A frequent problem in longitudinal studies is that data may be assessed at subject-selected, irregularly spaced time-points, resulting in highly unbalanced outcome data, inducing bias, especially if availability of data is directly related to outcome. Our aim was to develop a multivariate joint model in a mixed outcomes framework to minimize irregular sampling bias. We demonstrate using blood glucose monitoring throughout pregnancy and risk of preterm birth among women with type 1 diabetes mellitus. Blood glucose measurements were unequally spaced and intensity of sampling varied between and within individuals over time. Multivariate linear mixed effects submodel for the longitudinal outcome (blood glucose), Poisson model for the intensity of glucose sampling, and logistic regression model for binary process (preterm birth) were specified. Association between models is captured through shared random effects. Markov chain Monte Carlo methods were used to fit the model. The multivariate joint model provided better prediction, compared with a joint model with a multivariate linear mixed effects submodel (ignoring intensity of glucose sampling) and a two-stage model. Most association parameters were significant in the preterm birth outcome model, signifying improvement of predictive ability of the binary endpoint by sharing random effects between glucose monitoring and preterm birth. A simulation study is presented to illustrate the effectiveness of the multivariate joint modeling approach.

Early identification of epilepsy surgery candidates: A multicenter, machine learning study.

OBJECTIVES: Epilepsy surgery is underutilized. Automating the identification of potential surgical candidates may facilitate earlier intervention. Our objective was to develop site-specific machine learning (ML) algorithms to identify candidates before they undergo surgery. MATERIALS & METHODS: In this multicenter, retrospective, longitudinal cohort study, ML algorithms were trained on n-grams extracted from free-text neurology notes, EEG and MRI reports, visit codes, medications, procedures, laboratories, and demographic information. Site-specific algorithms were developed at two epilepsy centers: one pediatric and one adult. Cases were defined as patients who underwent resective epilepsy surgery, and controls were patients with epilepsy with no history of surgery. The output of the ML algorithms was the estimated likelihood of candidacy for resective epilepsy surgery. Model performance was assessed using 10-fold cross-validation. RESULTS: There were 5880 children (n = 137 had surgery [2.3%]) and 7604 adults with epilepsy (n = 56 had surgery [0.7%]) included in the study. Pediatric surgical patients could be identified 2.0 years (range: 0-8.6 years) before beginning their presurgical evaluation with AUC =0.76 (95% CI: 0.70-0.82) and PR-AUC =0.13 (95% CI: 0.07-0.18). Adult surgical patients could be identified 1.0 year (range: 0-5.4 years) before beginning their presurgical evaluation with AUC =0.85 (95% CI: 0.78-0.93) and PR-AUC =0.31 (95% CI: 0.14-0.48). By the time patients began their presurgical evaluation, the ML algorithms identified pediatric and adult surgical patients with AUC =0.93 and 0.95, respectively. The mean squared error of the predicted probability of surgical candidacy (Brier scores) was 0.018 in pediatrics and 0.006 in adults. CONCLUSIONS: Site-specific machine learning algorithms can identify candidates for epilepsy surgery early in the disease course in diverse practice settings.

Prospective validation of a machine learning model that uses provider notes to identify candidates for resective epilepsy surgery.

OBJECTIVE: Delay to resective epilepsy surgery results in avoidable disease burden and increased risk of mortality. The objective was to prospectively validate a natural language processing (NLP) application that uses provider notes to assign epilepsy surgery candidacy scores. METHODS: The application was trained on notes from (1) patients with a diagnosis of epilepsy and a history of resective epilepsy surgery and (2) patients who were seizure-free without surgery. The testing set included all patients with unknown surgical candidacy status and an upcoming neurology visit. Training and testing sets were updated weekly for 1 year. One- to three-word phrases contained in patients' notes were used as features. Patients prospectively identified by the application as candidates for surgery were manually reviewed by two epileptologists. Performance metrics were defined by comparing NLP-derived surgical candidacy scores with surgical candidacy status from expert chart review. RESULTS: The training set was updated weekly and included notes from a mean of 519 ± 67 patients. The area under the receiver operating characteristic curve (AUC) from 10-fold cross-validation was 0.90 ± 0.04 (range = 0.83-0.96) and improved by 0.002 per week (P < .001) as new patients were added to the training set. Of the 6395 patients who visited the neurology clinic, 4211 (67%) were evaluated by the model. The prospective AUC on this test set was 0.79 (95% confidence interval [CI] = 0.62-0.96). Using the optimal surgical candidacy score threshold, sensitivity was 0.80 (95% CI = 0.29-0.99), specificity was 0.77 (95% CI = 0.64-0.88), positive predictive value was 0.25 (95% CI = 0.07-0.52), and negative predictive value was 0.98 (95% CI = 0.87-1.00). The number needed to screen was 5.6. SIGNIFICANCE: An electronic health record-integrated NLP application can accurately assign surgical candidacy scores to patients in a clinical setting.

Dynamic predictive probabilities to monitor rapid cystic fibrosis disease progression.

Cystic fibrosis (CF) is a progressive, genetic disease characterized by frequent, prolonged drops in lung function. Accurately predicting rapid underlying lung-function decline is essential for clinical decision support and timely intervention. Determining whether an individual is experiencing a period of rapid decline is complicated due to its heterogeneous timing and extent, and error component of the measured lung function. We construct individualized predictive probabilities for "nowcasting" rapid decline. We assume each patient's true longitudinal lung function, S(t), follows a nonlinear, nonstationary stochastic process, and accommodate between-patient heterogeneity through random effects. Corresponding lung-function decline at time t is defined as the rate of change, S'(t). We predict S'(t) conditional on observed covariate and measurement history by modeling a measured lung function as a noisy version of S(t). The method is applied to data on 30 879 US CF Registry patients. Results are contrasted with a currently employed decision rule using single-center data on 212 individuals. Rapid decline is identified earlier using predictive probabilities than the center's currently employed decision rule (mean difference: 0.65 years; 95% confidence interval (CI): 0.41, 0.89). We constructed a bootstrapping algorithm to obtain CIs for predictive probabilities. We illustrate real-time implementation with R Shiny. Predictive accuracy is investigated using empirical simulations, which suggest this approach more accurately detects peak decline, compared with a uniform threshold of rapid decline. Median area under the ROC curve estimates (Q1-Q3) were 0.817 (0.814-0.822) and 0.745 (0.741-0.747), respectively, implying reasonable accuracy for both. This article demonstrates how individualized rate of change estimates can be coupled with probabilistic predictive inference and implementation for a useful medical-monitoring approach.

Data driven decision making to characterize clinical personas of parents of children with cystic fibrosis: a mixed methods study.

BACKGROUND: Beginning at a young age, children with cystic fibrosis (CF) embark on demanding care regimens that pose challenges to parents. We examined the extent to which clinical, demographic and psychosocial features inform patterns of adherence to pulmonary therapies and how these patterns can be used to develop clinical personas, defined as aspects of adherence barriers that are presented by parents and/or perceived by clinicians, in order to enhance personalized CF care delivery. METHODS: We undertook an explanatory sequential mixed-methods study consisting of i) multivariate clustering to create clusters corresponding to parental adherence patterns (quantitative phase); ii) parental participant interviews to create clinical personas interpreted from clustering (qualitative phase). Clinical, demographic and psychosocial features were used in supervised clustering against clinical endpoints, which included adherence to airway clearance and aerosolized medications and self-efficacy score, which was used as a feature for modeling adherence. Clinical implications were developed for each persona by combing quantitative and qualitative data (integration phase). RESULTS: The quantitative phase showed that the 87 parent participants were segmented into three distinct patterns of adherence based on use of aerosolized medication and practice of airway clearance. Patterns were primarily influenced by self-efficacy, distance to CF care center and child BMI percentile. The two key patterns that emerged for the self-efficacy model were most heavily influenced by distance to CF care center and child BMI percentile. Eight clinical personas were developed in the qualitative phase from parent and clinician participant feedback of latent components from these models. Findings from the integration phase include recommendations to overcome specific challenges with maintaining treatment regimens and increasing support from social networks. CONCLUSIONS: Adherence patterns from multivariate models and resulting parent personas with their corresponding clinical implications have utility as clinical decision support tools and capabilities for tailoring intervention study designs that promote adherence.

Dynamic Prediction of Survival in Cystic Fibrosis: A Landmarking Analysis Using UK Patient Registry Data.

BACKGROUND: Cystic fibrosis (CF) is an inherited, chronic, progressive condition affecting around 10,000 individuals in the United Kingdom and over 70,000 worldwide. Survival in CF has improved considerably over recent decades, and it is important to provide up-to-date information on patient prognosis. METHODS: The UK Cystic Fibrosis Registry is a secure centralized database, which collects annual data on almost all CF patients in the United Kingdom. Data from 43,592 annual records from 2005 to 2015 on 6181 individuals were used to develop a dynamic survival prediction model that provides personalized estimates of survival probabilities given a patient's current health status using 16 predictors. We developed the model using the landmarking approach, giving predicted survival curves up to 10 years from 18 to 50 years of age. We compared several models using cross-validation. RESULTS: The final model has good discrimination (C-indexes: 0.873, 0.843, and 0.804 for 2-, 5-, and 10-year survival prediction) and low prediction error (Brier scores: 0.036, 0.076, and 0.133). It identifies individuals at low and high risk of short- and long-term mortality based on their current status. For patients 20 years of age during 2013-2015, for example, over 80% had a greater than 95% probability of 2-year survival and 40% were predicted to survive 10 years or more. CONCLUSIONS: Dynamic personalized prediction models can guide treatment decisions and provide personalized information for patients. Our application illustrates the utility of the landmarking approach for making the best use of longitudinal and survival data and shows how models can be defined and compared in terms of predictive performance.

Investigation of bias in an epilepsy machine learning algorithm trained on physician notes.

Racial disparities in the utilization of epilepsy surgery are well documented, but it is unknown whether a natural language processing (NLP) algorithm trained on physician notes would produce biased recommendations for epilepsy presurgical evaluations. To assess this, an NLP algorithm was trained to identify potential surgical candidates using 1097 notes from 175 epilepsy patients with a history of resective epilepsy surgery and 268 patients who achieved seizure freedom without surgery (total N = 443 patients). The model was tested on 8340 notes from 3776 patients with epilepsy whose surgical candidacy status was unknown (2029 male, 1747 female, median age = 9 years; age range = 0-60 years). Multiple linear regression using demographic variables as covariates was used to test for correlations between patient race and surgical candidacy scores. After accounting for other demographic and socioeconomic variables, patient race, gender, and primary language did not influence surgical candidacy scores (P > .35 for all). Higher scores were given to patients >18 years old who traveled farther to receive care, and those who had a higher family income and public insurance (P < .001, .001, .001, and .01, respectively). Demographic effects on surgical candidacy scores appeared to reflect patterns in patient referrals.

Feasibility and Acceptability of a Telephone-Based Chaplaincy Intervention to Decrease Parental Spiritual Struggle.

Spiritual struggles (SSs) are distressing spiritual thoughts associated with poorer health outcomes. This study's purpose was to test feasibility, acceptability, and fidelity of an intervention to decrease SS of parents of children with CF. Parents screening positive for SS were enrolled and were randomized to intervention or attention-control condition. Intervention focused on intra-, inter-, and divine SS. Mixed linear modeling examined between-group differences. We present analyses of N = 23, and participants all showed decreased levels of SS. Acceptability was high; feasibility was higher in the intervention arm. GuideSS_CF is acceptable and feasible and warrants development as a potentially efficacious intervention.

Documenting an epidemic of suffering: low health-related quality of life among transgender youth.

PURPOSE: To quantify HRQOL of TGN patients using the PedsQL 4.0 generic core scales, and to compare reported HRQOL of TGN adolescents with published data from comparison populations. METHODS: Transgender children and adolescents (N = 142; 68% natal females) ages 6-23 years (M = 15.9, SD = 3.7) attending an outpatient clinic for TGN care at an academic pediatric hospital and caregivers of children and adolescents (N = 95) completed the PedsQL 4.0 generic core scales. Scores were compared with published scores for healthy adolescents and adolescents with 10 chronic diseases. RESULTS: TGN youth reported significantly lower overall HRQOL (more than twice the clinically meaningful difference) compared to youth without chronic disease. Total self-reported TGN HRQOL (M(SD), 65.72(17.40)) was lower than all chronic disease comparison groups except for rheumatology and cerebral palsy. TGN youth reported physical functioning (M(SD), 75.33(22.87)) lower than or similar to chronically ill comparisons, but higher than rheumatology and cerebral palsy groups. Psychosocial functioning (M(SD), 59.87(17.83)) was lower than all comparison samples and similar to youth with cerebral palsy. Results were similar for parent proxy-reports of TGN youth HRQOL (LS means: 68.75; 95% CI 65.87-71.61 vs 66.16; 95% CI 62.87-69.45; p = 0.12). CONCLUSIONS: TGN youth reported low HRQOL across all domains; most were significantly lower than healthy peers or peers with chronic diseases. Clinicians should understand the magnitude of TGN youth's low HRQOL and offer them and their caregivers resources to maximize their ability to achieve their full potential for healthy and productive lives.