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.

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.

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.

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.

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.

Phenotypes of Rapid Cystic Fibrosis Lung Disease Progression during Adolescence and Young Adulthood.

RATIONALE: Individuals with cystic fibrosis are at risk for prolonged drops in lung function, clinically termed rapid decline, during discreet periods of the disease. OBJECTIVES: To identify phenotypes of rapid pulmonary decline and determine how these phenotypes are related to patient characteristics. METHODS: A longitudinal cohort study of patients with cystic fibrosis aged 6-21 years was conducted using the Cystic Fibrosis Foundation Patient Registry. A statistical approach for clustering longitudinal profiles, sparse functional principal components analysis, was used to classify patients into distinct phenotypes by evaluating trajectories of FEV1 decline. Phenotypes were compared with respect to baseline and mortality characteristics. MEASUREMENTS AND MAIN RESULTS: Three distinct phenotypes of rapid decline were identified, corresponding to early, middle, and late timing of maximal FEV1 loss, in the overall cohort (n = 18,387). The majority of variation (first functional principal component, 94%) among patient profiles was characterized by differences in mean longitudinal FEV1 trajectories. Average degree of rapid decline was similar among phenotypes (roughly -3% predicted/yr); however, average timing differed, with early, middle, and late phenotypes experiencing rapid decline at 12.9, 16.3, and 18.5 years of age, respectively. Individuals with the late phenotype had the highest initial FEV1 but experienced the greatest loss of lung function. The early phenotype was more likely to have respiratory infections and acute exacerbations at baseline or to develop them subsequently, compared with other phenotypes. CONCLUSIONS: By identifying phenotypes and associated risk factors, timing of interventions may be more precisely targeted for subgroups at highest risk of lung function loss.

Flexible semiparametric joint modeling: an application to estimate individual lung function decline and risk of pulmonary exacerbations in cystic fibrosis.

BACKGROUND: Epidemiologic surveillance of lung function is key to clinical care of individuals with cystic fibrosis, but lung function decline is nonlinear and often impacted by acute respiratory events known as pulmonary exacerbations. Statistical models are needed to simultaneously estimate lung function decline while providing risk estimates for the onset of pulmonary exacerbations, in order to identify relevant predictors of declining lung function and understand how these associations could be used to predict the onset of pulmonary exacerbations. METHODS: Using longitudinal lung function (FEV1) measurements and time-to-event data on pulmonary exacerbations from individuals in the United States Cystic Fibrosis Registry, we implemented a flexible semiparametric joint model consisting of a mixed-effects submodel with regression splines to fit repeated FEV1 measurements and a time-to-event submodel for possibly censored data on pulmonary exacerbations. We contrasted this approach with methods currently used in epidemiological studies and highlight clinical implications. RESULTS: The semiparametric joint model had the best fit of all models examined based on deviance information criterion. Higher starting FEV1 implied more rapid lung function decline in both separate and joint models; however, individualized risk estimates for pulmonary exacerbation differed depending upon model type. Based on shared parameter estimates from the joint model, which accounts for the nonlinear FEV1 trajectory, patients with more positive rates of change were less likely to experience a pulmonary exacerbation (HR per one standard deviation increase in FEV1 rate of change = 0.566, 95% CI 0.516-0.619), and having higher absolute FEV1 also corresponded to lower risk of having a pulmonary exacerbation (HR per one standard deviation increase in FEV1 = 0.856, 95% CI 0.781-0.937). At the population level, both submodels indicated significant effects of birth cohort, socioeconomic status and respiratory infections on FEV1 decline, as well as significant effects of gender, socioeconomic status and birth cohort on pulmonary exacerbation risk. CONCLUSIONS: Through a flexible joint-modeling approach, we provide a means to simultaneously estimate lung function trajectories and the risk of pulmonary exacerbations for individual patients; we demonstrate how this approach offers additional insights into the clinical course of cystic fibrosis that were not possible using conventional approaches.

Morphological and quantitative evaluation of emphysema in chronic obstructive pulmonary disease patients: A comparative study of MRI with CT.

PURPOSE: To further validate the ability of ultrashort echo-time (UTE) magnetic resonance imaging (MRI) in quantifying lung density in patients diagnosed with chronic obstructive pulmonary disease (COPD) and to develop an MRI-based emphysema index (EI). MATERIALS AND METHODS: Ten subjects clinically diagnosed with COPD (5M/5F, age 62.6 ± 8.5 years) and ten healthy subjects (2M/8F, age 48.9 ± 19.2 years) were imaged via UTE MRI at 3T (4 mm slices, 1.39 × 1.39 mm2 pixels). Chest computed tomography (CT) images (generally 5 mm slices, ≈0.55 × 0.55 mm2 pixels), acquired retrospectively, were compared to UTE MRI. CT lung densities, MR lung-signal density, and EI were quantified from both CT and UTE MR images via a quantitative automated analysis and compared to the percent predicted forced expiratory volume in 1 second (FEV1 % predicted). RESULTS: EI quantified in controls via CT and UTE MRI was 0.23 ± 0.78% and 2.40 ± 1.50%, respectively; in COPD subjects it was 13.3 ± 14.9% (P = 0.021) and 12.0 ± 9.8% (P = 0.013), respectively. Bland-Altman determined the mean differences and 95% limits of agreement for COPD subjects and healthy controls were 0.06 (12.50 to -12.38). Strong correlation (R2 = 0.79, P < 0.0001) existed between EIs quantified from both CT and UTE MRI. There was a slightly higher correlation between FEV1 % predicted and the UTE MRI EI (R2 = 0.65, P < 0.0001) compared to CT EI (R2 = 0.49, P < 0.0001). CONCLUSION: Our results demonstrate a significant positive correlation between lung density and EI assessed with CT and MRI. Furthermore, UTE MRI exhibits its potential as a diagnostic alternative to CT for assessing the extent and the severity of emphysema, particularly for longitudinal studies. J. Magn. Reson. Imaging 2016;44:1656-1663.

Adolescents' Spirituality and Cystic Fibrosis Airway Clearance Treatment Adherence: Examining Mediators.

OBJECTIVE: Adolescent cystic fibrosis (CF) treatment adherence is a significant multidimensional issue. Using the Theory of Reasoned Action (TRA), this study examined the role of spiritual factors in adherence. METHODS: Forty-five 11-19-year-olds diagnosed with CF completed questionnaires concerning psychosocial, spiritual, and adherence-related constructs and Daily Phone Diaries to calculate treatment adherence. Exploratory Factor Analysis identified two spiritual factors used in subsequent analyses. The mediating roles of attitude toward the treatment's value (utility), subjective behavioral norms (the product of perceived behavioral norms and one's motivation to comply with them), self-efficacy for completing the treatments and treatment intentions in the relationship between spiritual factors and treatment adherence were tested with path analysis. RESULTS: Lower 'spiritual struggle' and greater 'engaged spirituality' predicted treatment attitude (utility) and subjective behavioral norms, which, together with self-efficacy, predicted treatment intentions. Finally, treatment intentions predicted airway clearance adherence. CONCLUSIONS: Findings were consistent with the TRA. Engaged spirituality supports pro-adherence determinants and behavior. Spiritual struggle's negative associations with outcomes warrant screening and intervention.

Longitudinal Patterns of Glycemic Control and Blood Pressure in Pregnant Women with Type 1 Diabetes Mellitus: Phenotypes from Functional Data Analysis.

Objective To identify phenotypes of type 1 diabetes control and associations with maternal/neonatal characteristics based on blood pressure (BP), glucose, and insulin curves during gestation, using a novel functional data analysis approach that accounts for sparse longitudinal patterns of medical monitoring during pregnancy. Methods We performed a retrospective longitudinal cohort study of women with type 1 diabetes whose BP, glucose, and insulin requirements were monitored throughout gestation as part of a program-project grant. Scores from sparse functional principal component analysis (fPCA) were used to classify gestational profiles according to the degree of control for each monitored measure. Phenotypes created using fPCA were compared with respect to maternal and neonatal characteristics and outcome. Results Most of the gestational profile variation in the monitored measures was explained by the first principal component (82-94%). Profiles clustered into three subgroups of high, moderate, or low heterogeneity, relative to the overall mean response. Phenotypes were associated with baseline characteristics, longitudinal changes in glycohemoglobin A1 and weight, and to pregnancy-related outcomes. Conclusion Three distinct longitudinal patterns of glucose, insulin, and BP control were found. By identifying these phenotypes, interventions can be targeted for subgroups at highest risk for compromised outcome, to optimize diabetes management during pregnancy.

Semiparametric Mixed Models for Nested Repeated Measures Applied to Ambulatory Blood Pressure Monitoring Data.

Semiparametric mixed models are increasingly popular for statistical analysis of medical device studies in which long sequences of repeated measurements are recorded. Monitoring these sequences at different periods over time on the same individual, such as before and after an intervention, results in nested repeated measures (NRM). Covariance models to account for NRM and simultaneously address mean profile estimation with penalized splines via semiparametric regression are considered with application to a prospective study of 24-hour ambulatory blood pressure and the impact of surgical intervention on obstructive sleep apnea.

Regional ventilation changes in severe asthma after bronchial thermoplasty with (3)He MR imaging and CT.

PURPOSE: To quantify regional lung ventilation in healthy volunteers and patients with severe asthma (both before and after thermoplasty) by using a combination of helium 3 ((3)He) magnetic resonance (MR) imaging and computed tomography (CT), with the intention of developing more effective image-guided treatments for obstructive lung diseases. MATERIALS AND METHODS: With approval of the local institutional review board, informed consent, and an Investigational New Drug Exemption, six healthy volunteers and 10 patients with severe asthma were imaged in compliance with HIPAA regulations by using both multidetector CT and (3)He MR imaging. Individual bronchopulmonary segments were labeled voxel by voxel from the CT images and then registered to the (3)He MR images by using custom software. The (3)He signal intensity was then analyzed by evaluating the volume-weighted fraction of total-lung signal intensity present in each segment (segmental ventilation percentage [ SVP segmental ventilation percentage ]) and by identifying the whole-lung defect percentage and the segmental defect percentage. Of the 10 patients with asthma, seven received treatment with bronchial thermoplasty and were imaged with (3)He MR a second time. Changes in segmental defect percentages and whole-lung defect percentages are presented. RESULTS: Ventilation measures for healthy volunteers yielded smaller segment-to-segment variation (mean SVP segmental ventilation percentage , 100% ± 18 [standard deviation]) than did the measures for patients with severe asthma (mean SVP segmental ventilation percentage , 97% ± 23). Patients with asthma also demonstrated larger segmental defect percentages (median, 13.5%; interquartile range, 8.9%-17.8%) than healthy volunteers (median, 6%; interquartile range, 5.6%-6.3%). These quantitative results confirm what is visually observed on the (3)He images. A Spearman correlation of r = -0.82 was found between the change in whole-lung defect percentage and the number of days between final treatment and second (3)He imaging. CONCLUSION: Regional quantification of lung ventilation is indeed feasible and may be a useful technique for image-guided treatment of obstructive lung diseases, such as bronchial thermoplasty for severe asthma. In these patients, ventilation defects decreased as a function of time after treatment.