Federico Mattiello

High-throughput sequencing technologies allow easy characterization of the human microbiome, but the statistical methods to analyze microbiome data are still in their infancy. Differential abundance methods aim at detecting associations between the abundances of bacterial species and subject grouping factors. The results of such methods are important to identify the microbiome as a prognostic or diagnostic biomarker or to demonstrate efficacy of prodrug or antibiotic drugs. Because of a lack of benchmarking studies in the microbiome field, no consensus exists on the performance of the statistical methods. We have compared a large number of popular methods through extensive parametric and nonparametric simulation as well as real data shuffling algorithms. The results are consistent over the different approaches and all point to an alarming excess of false discoveries. This raises great doubts about the reliability of discoveries in past studies and imperils reproducibility of microbiome experiments. To further improve method benchmarking, we introduce a new simulation tool that allows to generate correlated count data following any univariate count distribution; the correlation structure may be inferred from real data. Most simulation studies discard the correlation between species, but our results indicate that this correlation can negatively affect the performance of statistical methods.

Purpose To perform an updated analysis of the randomized phase III GADOLIN trial in patients with rituximab-refractory indolent non-Hodgkin lymphoma treated with obinutuzumab (GA101; G) and bendamustine (B). Patients and Methods Patients with histologically documented, rituximab-refractory CD20 + indolent non-Hodgkin lymphoma received G 1,000 mg (days 1, 8, and 15, cycle 1; day 1, cycles 2 to 6) plus B 90 mg/m 2 /d (days 1 and 2, all cycles) or B 120 mg/m 2 /d monotherapy. Patients who did not experience disease progression with G-B received G maintenance (1,000 mg every 2 months) for up to 2 years. The primary end point was progression-free survival (PFS). Results Of 413 randomly assigned patients (intention-to-treat [ITT]: G-B, n = 204; B monotherapy, n = 209), 335 had follicular lymphoma (FL; G-B, n = 164; B monotherapy, n = 171). After a median follow-up of 31.8 months, median PFS in ITT patients was 25.8 months (G-B) and 14.1 months (B monotherapy; hazard ratio [HR], 0.57; 95% CI, 0.44 to 0.73; P < .001). Overall survival (OS) also was prolonged (HR, 0.67; 95% CI, 0.47 to 0.96; P = .027). PFS and OS benefits were similar in patients with FL. Grade 3 to 5 adverse events (AEs) were reported by 148 (72.5%) and 133 (65.5%) patients in the G-B and B monotherapy arms, respectively, most commonly neutropenia (G-B, 34.8%; B monotherapy, 27.1%), thrombocytopenia (10.8% and 15.8%), anemia (7.4% and 10.8%), and infusion-related reactions (9.3% and 3.4%). Serious AEs occurred in 89 G-B patients (43.6%) and 75 B monotherapy patients (36.9%); fatal AEs occurred in 16 (7.8%) and 13 (6.4%), respectively. Conclusion This updated analysis confirms the PFS benefit for G-B shown in the primary analysis. A substantial OS benefit also was demonstrated in the ITT population and in patients with FL. Toxicity was similar for both treatments.

Abstract Central nervous system (CNS) relapse carries a poor prognosis in diffuse large B-cell lymphoma (DLBCL). Integrating biomarkers into the CNS–International Prognostic Index (CNS-IPI) risk model may improve identification of patients at high risk for developing secondary CNS disease. CNS relapse was analyzed in 1418 DLBCL patients treated with obinutuzumab or rituximab plus cyclophosphamide, doxorubicin, vincristine, prednisone chemotherapy in the phase 3 GOYA study. Cell of origin (COO) was assessed using gene-expression profiling. BCL2 and MYC protein expression was analyzed by immunohistochemistry. The impact of CNS-IPI, COO, and BCL2/MYC dual-expression status on CNS relapse was assessed using a multivariate Cox regression model (data available in n = 1418, n = 933, and n = 688, respectively). High CNS-IPI score (hazard ratio [HR], 4.0; 95% confidence interval [CI], 1.3-12.3; P = .02) and activated B-cell‒like (ABC) (HR, 5.2; 95% CI, 2.1-12.9; P = .0004) or unclassified COO subtypes (HR, 4.2; 95% CI, 1.5-11.7; P = .006) were independently associated with CNS relapse. BCL2/MYC dual-expression status did not impact CNS relapse risk. Three risk subgroups were identified based on the presence of high CNS-IPI score and/or ABC/unclassified COO (CNS-IPI-C model): low risk (no risk factors, n = 450 [48.2%]), intermediate risk (1 factor, n = 408 [43.7%]), and high risk (both factors, n = 75 [8.0%]). Two-year CNS relapse rates were 0.5%, 4.4%, and 15.2% in the respective risk subgroups. Combining high CNS-IPI and ABC/unclassified COO improved CNS relapse prediction and identified a patient subgroup at high risk for developing CNS relapse. The study was registered at www.clinicaltrials.gov as #NCT01287741.

UNLABELLED: : When designing a case-control study to investigate differences in microbial composition, it is fundamental to assess the sample sizes needed to detect an hypothesized difference with sufficient statistical power. Our application includes power calculation for (i) a recoded version of the two-sample generalized Wald test of the 'HMP' R-package for comparing community composition, and (ii) the Wilcoxon-Mann-Whitney test for comparing operational taxonomic unit-specific abundances between two samples (optional). The simulation-based power calculations make use of the Dirichlet-Multinomial model to describe and generate abundances. The web interface allows for easy specification of sample and effect sizes. As an illustration of our application, we compared the statistical power of the two tests, with and without stratification of samples. We observed that statistical power increases considerably when stratification is employed, meaning that less samples are needed to detect the same effect size with the same power. AVAILABILITY AND IMPLEMENTATION: The web interface is written in R code using Shiny (RStudio Inc., 2016) and it is available at https://fedematt.shinyapps.io/shinyMB The R code for the recoded generalized Wald test can be found at https://github.com/mafed/msWaldHMP CONTACT: Federico.Mattiello@UGent.be.

GOYA was a randomized phase 3 study comparing obinutuzumab plus cyclophosphamide, doxorubicin, vincristine, and prednisone (CHOP) vs standard-of-care rituximab plus CHOP in patients with previously untreated diffuse large B-cell lymphoma (DLBCL). This retrospective analysis of GOYA aimed to assess the association between progression-free survival (PFS) and overall survival (OS) with positron emission tomography (PET)-based complete response (CR) status. Overall, 1418 patients were randomly assigned to receive 8 21-day cycles of obinutuzumab (n = 706) or rituximab (n = 712) plus 6 or 8 cycles of CHOP. Patients received a mandatory fluoro-2-deoxy-d-glucose-PET/computed tomography scan at baseline and end of treatment. After a median follow-up of 29 months, the numbers of independent review committee-assessed PFS and OS events in the entire cohort were 416 (29.3%) and 252 (17.8%), respectively. End-of-treatment PET CR was highly prognostic for PFS and OS according to Lugano 2014 criteria (PFS: hazard ratio [HR], 0.26; 95% confidence interval [CI], 0.19-0.38; P < .0001; OS: HR, 0.12; 95% CI, 0.08-0.17; P < .0001), irrespective of international prognostic index score and cell of origin. In conclusion, the results from this prospectively acquired large cohort corroborated previously published data from smaller sample sizes showing that end-of-treatment PET CR is an independent predictor of PFS and OS and a promising prognostic marker in DLBCL. Long-term survival analysis confirmed the robustness of these data over time. Additional meta-analyses including other prospective studies are necessary to support the substitution of PET CR for PFS as an effective and practical surrogate end point. This trial was registered at www.clinicaltrials.gov as #NCT01287741.