Assessing Clinician Utilization of Next-Generation Antibiotics Against Resistant Gram-Negative Infections in U.S. Hospitals : A Retrospective Cohort Study.

BACKGROUND: The U.S. antibiotic market failure has threatened future innovation and supply. Understanding when and why clinicians underutilize recently approved gram-negative antibiotics might help prioritize the patient in future antibiotic development and potential market entry rewards. OBJECTIVE: To determine use patterns of recently U.S. Food and Drug Administration (FDA)-approved gram-negative antibiotics (ceftazidime-avibactam, ceftolozane-tazobactam, meropenem-vaborbactam, plazomicin, eravacycline, imipenem-relebactam-cilastatin, and cefiderocol) and identify factors associated with their preferential use (over traditional generic agents) in patients with gram-negative infections due to pathogens displaying difficult-to-treat resistance (DTR; that is, resistance to all first-line antibiotics). DESIGN: Retrospective cohort. SETTING: 619 U.S. hospitals. PARTICIPANTS: Adult inpatients. MEASUREMENTS: Quarterly percentage change in antibiotic use was calculated using weighted linear regression. Machine learning selected candidate variables, and mixed models identified factors associated with new (vs. traditional) antibiotic use in DTR infections. RESULTS: Between quarter 1 of 2016 and quarter 2 of 2021, ceftolozane-tazobactam (approved 2014) and ceftazidime-avibactam (2015) predominated new antibiotic usage whereas subsequently approved gram-negative antibiotics saw relatively sluggish uptake. Among gram-negative infection hospitalizations, 0.7% (2551 [2631 episodes] of 362 142) displayed DTR pathogens. Patients were treated exclusively using traditional agents in 1091 of 2631 DTR episodes (41.5%), including "reserve" antibiotics such as polymyxins, aminoglycosides, and tigecycline in 865 of 1091 episodes (79.3%). Patients with bacteremia and chronic diseases had greater adjusted probabilities and those with do-not-resuscitate status, acute liver failure, and Acinetobacter baumannii complex and other nonpseudomonal nonfermenter pathogens had lower adjusted probabilities of receiving newer (vs. traditional) antibiotics for DTR infections, respectively. Availability of susceptibility testing for new antibiotics increased probability of usage. LIMITATION: Residual confounding. CONCLUSION: Despite FDA approval of 7 next-generation gram-negative antibiotics between 2014 and 2019, clinicians still frequently treat resistant gram-negative infections with older, generic antibiotics with suboptimal safety-efficacy profiles. Future antibiotics with innovative mechanisms targeting untapped pathogen niches, widely available susceptibility testing, and evidence demonstrating improved outcomes in resistant infections might enhance utilization. PRIMARY FUNDING SOURCE: U.S. Food and Drug Administration; NIH Intramural Research Program.

Safety signal detection with control of latent factors.

Postmarket drug safety database like vaccine adverse event reporting system (VAERS) collect thousands of spontaneous reports annually, with each report recording occurrences of any adverse events (AEs) and use of vaccines. We hope to identify signal vaccine-AE pairs, for which certain vaccines are statistically associated with certain adverse events (AE), using such data. Thus, the outcomes of interest are multiple AEs, which are binary outcomes and could be correlated because they might share certain latent factors; and the primary covariates are vaccines. Appropriately accounting for the complex correlation among AEs could improve the sensitivity and specificity of identifying signal vaccine-AE pairs. We propose a two-step approach in which we first estimate the shared latent factors among AEs using a working multivariate logistic regression model, and then use univariate logistic regression model to examine the vaccine-AE associations after controlling for the latent factors. Our simulation studies show that this approach outperforms current approaches in terms of sensitivity and specificity. We apply our approach in analyzing VAERS data and report our findings.

Reliability of Admission Procalcitonin Testing for Capturing Bacteremia Across the Sepsis Spectrum: Real-World Utilization and Performance Characteristics, 65 U.S. Hospitals, 2008-2017.

OBJECTIVES: Serum procalcitonin is often ordered at admission for patients with suspected sepsis and bloodstream infections (BSIs), although its performance characteristics in this setting remain contested. This study aimed to evaluate use patterns and performance characteristics of procalcitonin-on-admission in patients with suspected BSI, with or without sepsis. DESIGN: Retrospective cohort study. SETTING: Cerner HealthFacts Database (2008-2017). PATIENTS: Adult inpatients (≥ 18 yr) who had blood cultures and procalcitonin drawn within 24 hours of admission. INTERVENTIONS: None. MEASUREMENTS AND MAIN RESULTS: Testing frequency of procalcitonin was determined. Sensitivity of procalcitonin-on-admission for detecting BSI due to different pathogens was calculated. Area under the receiver operating characteristic curve (AUC) was calculated to assess discrimination by procalcitonin-on-admission for BSI in patients with and without fever/hypothermia, ICU admission and sepsis defined by Centers for Disease Control and Prevention Adult Sepsis Event criteria. AUCs were compared using Wald test and p values were adjusted for multiple comparisons. At 65 procalcitonin-reporting hospitals, 74,958 of 739,130 patients (10.1%) who had admission blood cultures also had admission procalcitonin testing. Most patients (83%) who had admission day procalcitonin testing did not have a repeat procalcitonin test. Median procalcitonin varied considerably by pathogen, BSI source, and acute illness severity. At a greater than or equal to 0.5 ng/mL cutoff, sensitivity for BSI detection was 68.2% overall, ranging between 58.0% for enterococcal BSI without sepsis and 96.4% for pneumococcal sepsis. Procalcitonin-on-admission displayed moderate discrimination at best for overall BSI (AUC, 0.73; 95% CI, 0.72-0.73) and showed no additional utility in key subgroups. Empiric antibiotic use proportions were not different between blood culture sampled patients with a positive procalcitonin (39.7%) and negative procalcitonin (38.4%) at admission. CONCLUSIONS: At 65 study hospitals, procalcitonin-on-admission demonstrated poor sensitivity in ruling out BSI, moderate-to-poor discrimination for both bacteremic sepsis and occult BSI and did not appear to meaningfully alter empiric antibiotic usage. Diagnostic stewardship of procalcitonin-on-admission and risk assessment of admission procalcitonin-guided clinical decisions is warranted.

Joint semiparametric models for case-cohort designs.

Two-phase studies such as case-cohort and nested case-control studies are widely used cost-effective sampling strategies. In the first phase, the observed failure/censoring time and inexpensive exposures are collected. In the second phase, a subgroup of subjects is selected for measurements of expensive exposures based on the information from the first phase. One challenging issue is how to utilize all the available information to conduct efficient regression analyses of the two-phase study data. This paper proposes a joint semiparametric modeling of the survival outcome and the expensive exposures. Specifically, we assume a class of semiparametric transformation models and a semiparametric density ratio model for the survival outcome and the expensive exposures, respectively. The class of semiparametric transformation models includes the proportional hazards model and the proportional odds model as special cases. The density ratio model is flexible in modeling multivariate mixed-type data. We develop efficient likelihood-based estimation and inference procedures and establish the large sample properties of the nonparametric maximum likelihood estimators. Extensive numerical studies reveal that the proposed methods perform well under practical settings. The proposed methods also appear to be reasonably robust under various model mis-specifications. An application to the National Wilms Tumor Study is provided.

The prevalence and clinical significance of HER2 expression in prostate adenocarcinoma.

AIMS: Abnormalities in HER2 are well-established oncogenic drivers and are approved therapeutic targets in various malignancies. Prior studies evaluating HER2 expression in prostate cancer (PCa) have yielded variable results. Most of these studies used immunohistochemistry scoring systems based on breast cancer data. The goal of this study was to determine the prevalence and clinical significance of HER2 expression using a scoring system that better reflects the HER2 staining pattern observed in PCa. METHODS: We randomly selected similar numbers of localized low risk (AJCC stage I), locally advanced (AJCC stages II & III), and metastatic (AJCC stage IV) PCa patients treated at the DC VA Medical Center between 2000 and 2022. Among them, we included patients who had sufficient PCa tissue samples and clinical information required for this analysis. Two experienced pathologists independently scored HER2 expression (Ventana Pathway anti-HER2) according to a modified gastric cancer HER2 scoring system. RESULTS: Out of the 231 patients included, 85 % self-identified as Black. 58.9 % of patients expressed HER2 (1+: 35.5 %; 2+: 18.2 %; 3+: 5.2 %). Validity of the results was confirmed using the HercepTest antibody. Higher HER2 expression was associated with a higher Gleason Score/Grade Group and advanced disease. CONCLUSIONS: Our findings support the adverse prognostic impact on HER2 in PCa. We propose the use of a modified scoring system to evaluate HER2 expression in PCa. The observed high prevalence of HER2 expression supports the consideration of novel HER2-targeted therapies addressing even low levels of HER2 expression in future PCa trials.

Synthesizing studies for comparing different treatment sequences in clinical trials.

With advances in cancer treatments and improved patient survival, more patients may go through multiple lines of treatment. It is of clinical importance to choose a sequence of effective treatments (eg, lines of treatment) for individual patients with the goal of optimizing their long-term clinical outcome (eg, survival). Several important issues arise in cancer studies. First, cancer clinical trials are usually conducted by each line of treatment. For a treatment sequence, we may have first line and second line treatment data from two different studies. Second, there is typically a treatment initiation period varying from patient to patient between progression of disease and the start of the second line treatment due to administrative reasons. Additionally, the choice of the second line treatment for patients with progression of disease may depend on their characteristics. We address all these issues and develop semiparametric methods under the potential outcome framework for the estimation of the overall survival probability for a treatment sequence and for comparing different treatment sequences. We establish the large sample properties of the proposed inferential procedures. Simulation studies and an application to a colorectal clinical trial are provided.

Mixture survival trees for cancer risk classification.

In oncology studies, it is important to understand and characterize disease heterogeneity among patients so that patients can be classified into different risk groups and one can identify high-risk patients at the right time. This information can then be used to identify a more homogeneous patient population for developing precision medicine. In this paper, we propose a mixture survival tree approach for direct risk classification. We assume that the patients can be classified into a pre-specified number of risk groups, where each group has distinct survival profile. Our proposed tree-based methods are devised to estimate latent group membership using an EM algorithm. The observed data log-likelihood function is used as the splitting criterion in recursive partitioning. The finite sample performance is evaluated by extensive simulation studies and the proposed method is illustrated by a case study in breast cancer.

Inferring latent heterogeneity using many feature variables supervised by survival outcome.

In cancer studies, it is important to understand disease heterogeneity among patients so that precision medicine can particularly target high-risk patients at the right time. Many feature variables such as demographic variables and biomarkers, combined with a patient's survival outcome, can be used to infer such latent heterogeneity. In this work, we propose a mixture model to model each patient's latent survival pattern, where the mixing probabilities for latent groups are modeled through a multinomial distribution. The Bayesian information criterion is used for selecting the number of latent groups. Furthermore, we incorporate variable selection with the adaptive lasso into inference so that only a few feature variables will be selected to characterize the latent heterogeneity. We show that our adaptive lasso estimator has oracle properties when the number of parameters diverges with the sample size. The finite sample performance is evaluated by the simulation study, and the proposed method is illustrated by two datasets.

Statistically efficient association analysis of quantitative traits with haplotypes and untyped SNPs in family studies.

BACKGROUND: Associations between haplotypes and quantitative traits provide valuable information about the genetic basis of complex human diseases. Haplotypes also provide an effective way to deal with untyped SNPs. Two major challenges arise in haplotype-based association analysis of family data. First, haplotypes may not be inferred with certainty from genotype data. Second, the trait values within a family tend to be correlated because of common genetic and environmental factors. RESULTS: To address these challenges, we present an efficient likelihood-based approach to analyzing associations of quantitative traits with haplotypes or untyped SNPs. This approach properly accounts for within-family trait correlations and can handle general pedigrees with arbitrary patterns of missing genotypes. We characterize the genetic effects on the quantitative trait by a linear regression model with random effects and develop efficient likelihood-based inference procedures. Extensive simulation studies are conducted to examine the performance of the proposed methods. An application to family data from the Childhood Asthma Management Program Ancillary Genetic Study is provided. A computer program is freely available. CONCLUSIONS: Results from extensive simulation studies show that the proposed methods for testing the haplotype effects on quantitative traits have correct type I error rates and are more powerful than some existing methods.

On symmetric semiparametric two-sample problem.

We consider a two-sample problem where data come from symmetric distributions. Usual two-sample data with only magnitudes recorded, arising from case-control studies or logistic discriminant analyses, may constitute a symmetric two-sample problem. We propose a semiparametric model such that, in addition to symmetry, the log ratio of two unknown density functions is modeled in a known parametric form. The new semiparametric model, tailor-made for symmetric two-sample data, can also be viewed as a biased sampling model subject to symmetric constraint. A maximum empirical likelihood estimation approach is adopted to estimate the unknown model parameters, and the corresponding profile empirical likelihood ratio test is utilized to perform hypothesis testing regarding the two population distributions. Symmetry, however, comes with irregularity. It is shown that, under the null hypothesis of equal symmetric distributions, the maximum empirical likelihood estimator has degenerate Fisher information, and the test statistic has a mixture of χ2 -type asymptotic distribution. Extensive simulation studies have been conducted to demonstrate promising statistical powers under correct and misspecified models. We apply the proposed methods to two real examples.

Semiparametric frailty models for zero-inflated event count data in the presence of informative dropout.

Recurrent events data are commonly encountered in medical studies. In many applications, only the number of events during the follow-up period rather than the recurrent event times is available. Two important challenges arise in such studies: (a) a substantial portion of subjects may not experience the event, and (b) we may not observe the event count for the entire study period due to informative dropout. To address the first challenge, we assume that underlying population consists of two subpopulations: a subpopulation nonsusceptible to the event of interest and a subpopulation susceptible to the event of interest. In the susceptible subpopulation, the event count is assumed to follow a Poisson distribution given the follow-up time and the subject-specific characteristics. We then introduce a frailty to account for informative dropout. The proposed semiparametric frailty models consist of three submodels: (a) a logistic regression model for the probability such that a subject belongs to the nonsusceptible subpopulation; (b) a nonhomogeneous Poisson process model with an unspecified baseline rate function; and (c) a Cox model for the informative dropout time. We develop likelihood-based estimation and inference procedures. The maximum likelihood estimators are shown to be consistent. Additionally, the proposed estimators of the finite-dimensional parameters are asymptotically normal and the covariance matrix attains the semiparametric efficiency bound. Simulation studies demonstrate that the proposed methodologies perform well in practical situations. We apply the proposed methods to a clinical trial on patients with myelodysplastic syndromes.

Efficient methods for signal detection from correlated adverse events in clinical trials.

It is an important and yet challenging task to identify true signals from many adverse events that may be reported during the course of a clinical trial. One unique feature of drug safety data from clinical trials, unlike data from post-marketing spontaneous reporting, is that many types of adverse events are reported by only very few patients leading to rare events. Due to the limited study size, the p-values of testing whether the rate is higher in the treatment group across all types of adverse events are in general not uniformly distributed under the null hypothesis that there is no difference between the treatment group and the placebo group. A consequence is that typically fewer than 100α percent of the hypotheses are rejected under the null at the nominal significance level of α . The other challenge is multiplicity control. Adverse events from the same body system may be correlated. There may also be correlations between adverse events from different body systems. To tackle these challenging issues, we develop Monte-Carlo-based methods for the signal identification from patient-reported adverse events in clinical trials. The proposed methodologies account for the rare events and arbitrary correlation structures among adverse events within and/or between body systems. Extensive simulation studies demonstrate that the proposed method can accurately control the family-wise error rate and is more powerful than existing methods under many practical situations. Application to two real examples is provided.

Controlling false discovery proportion in identification of drug-related adverse events from multiple system organ classes.

Analyzing safety data from clinical trials to detect safety signals worth further examination involves testing multiple hypotheses, one for each observed adverse event (AE) type. There exists certain hierarchical structure for these hypotheses due to the classification of the AEs into system organ classes, and these AEs are also likely correlated. Many approaches have been proposed to identify safety signals under the multiple testing framework and tried to achieve control of false discovery rate (FDR). The FDR control concerns the expectation of the false discovery proportion (FDP). In practice, the control of the actual random variable FDP could be more relevant and has recently drawn much attention. In this paper, we proposed a two-stage procedure for safety signal detection with direct control of FDP, through a permutation-based approach for screening groups of AEs and a permutation-based approach of constructing simultaneous upper bounds for false discovery proportion. Our simulation studies showed that this new approach has controlled FDP. We demonstrate our approach using data sets derived from a drug clinical trial.

Quantifying time-varying cause-specific hazard and subdistribution hazard ratios with competing risks data.

Various non-proportional hazard models have been developed in the literature for competing risks data. The regression coefficients under these models, however, typically cannot be compared directly. We propose new methods to quantify the average of the time-varying cause-specific hazard ratios and subdistribution hazard ratios through two general classes of transformations and weight functions that are chosen to reflect the relative importance of the hazard ratios in different time periods. We further propose an L∞ -norm type of test statistic that incorporates the test statistics for all possible pairs of the transformation function and weight function under consideration. Extensive simulations are conducted under various settings of the hazards and demonstrate that the proposed test performs well under all settings. An application to a clinical trial in follicular lymphoma is examined in detail.

A class of semiparametric cure models with current status data.

Current status data occur in many biomedical studies where we only know whether the event of interest occurs before or after a particular time point. In practice, some subjects may never experience the event of interest, i.e., a certain fraction of the population is cured or is not susceptible to the event of interest. We consider a class of semiparametric transformation cure models for current status data with a survival fraction. This class includes both the proportional hazards and the proportional odds cure models as two special cases. We develop efficient likelihood-based estimation and inference procedures. We show that the maximum likelihood estimators for the regression coefficients are consistent, asymptotically normal, and asymptotically efficient. Simulation studies demonstrate that the proposed methods perform well in finite samples. For illustration, we provide an application of the models to a study on the calcification of the hydrogel intraocular lenses.

Robust Inference after Random Projections via Hellinger Distance for Location-Scale Family.

Big data and streaming data are encountered in a variety of contemporary applications in business and industry. In such cases, it is common to use random projections to reduce the dimension of the data yielding compressed data. These data however possess various anomalies such as heterogeneity, outliers, and round-off errors which are hard to detect due to volume and processing challenges. This paper describes a new robust and efficient methodology, using Hellinger distance, to analyze the compressed data. Using large sample methods and numerical experiments, it is demonstrated that a routine use of robust estimation procedure is feasible. The role of double limits in understanding the efficiency and robustness is brought out, which is of independent interest.

Biomarker threshold adaptive designs for survival endpoints.

Due to the importance of precision medicine, it is essential to identify the right patients for the right treatment. Biomarkers, which have been commonly used in clinical research as well as in clinical practice, can facilitate selection of patients with a good response to the treatment. In this paper, we describe a biomarker threshold adaptive design with survival endpoints. In the first stage, we determine subgroups for one or more biomarkers such that patients in these subgroups benefit the most from the new treatment. The analysis in this stage can be based on historical or pilot studies. In the second stage, we sample subjects from the subgroups determined in the first stage and randomly allocate them to the treatment or control group. Extensive simulation studies are conducted to examine the performance of the proposed design. Application to a real data example is provided for implementation of the first-stage algorithms.

Quantification of Muscle Tissue Properties by Modeling the Statistics of Ultrasound Image Intensities Using a Mixture of Gamma Distributions in Children With and Without Cerebral Palsy.

OBJECTIVES: To investigate whether quantitative ultrasound (US) imaging, based on the envelope statistics of the backscattered US signal, can describe muscle properties in typically developing children and those with cerebral palsy (CP). METHODS: Radiofrequency US data were acquired from the rectus femoris muscle of children with CP (n = 22) and an age-matched cohort without CP (n = 14) at rest and during maximal voluntary isometric contraction. A mixture of gamma distributions was used to model the histogram of the echo intensities within a region of interest in the muscle. RESULTS: Muscle in CP had a heterogeneous echo texture that was significantly different from that in healthy controls (P < .001), with larger deviations from Rayleigh scattering. A mixture of 2 gamma distributions showed an excellent fit to the US intensity, and the shape and rate parameters were significantly different between CP and control groups (P < .05). The rate parameters for both the single gamma distribution and mixture of gamma distributions were significantly higher for contracted muscles compared to resting muscles, but there was no significant interaction between these factors (CP and muscle contraction) for a mixed-model analysis of variance. CONCLUSIONS: Ultrasound tissue characterization indicates a more disorganized architecture and increased echogenicity in muscles in CP, consistent with previously documented increases in fibrous infiltration and connective tissue changes in this population. Our results indicate that quantitative US can be used to objectively differentiate muscle architecture and tissue properties.

A semiparametric method for comparing the discriminatory ability of biomarkers subject to limit of detection.

Receiver operating characteristic curves and the area under the curves (AUC) are often used to compare the discriminatory ability of potentially correlated biomarkers. Many biomarkers are subject to limit of detection due to the instrumental limitation in measurements and may not be normally distributed. Standard parametric methods assuming normality can lead to biased results when the normality assumption is violated. We propose new estimation and inference procedures for the AUCs of biomarkers subject to limit of detection by using the semiparametric transformation model allowing for heteroscedasticity. We obtain the nonparametric maximum likelihood estimators by maximizing the likelihood for the observed data with limit of detection. The proposed estimators are shown to be consistent, asymptotically normal, and asymptotically efficient. Additionally, we propose a Wald type test statistic to compare the AUCs of 2 potentially correlated biomarkers with limit of detection. Extensive simulation studies demonstrate that the proposed method is robust to nonnormality while performing as well as its parametric counterpart when the normality assumption is true. An application to an autism study is provided.