Chen, Li

Abstract Background We performed a prospective study evaluating a baseline GA battery and clinical outcomes in fit older adults (age ≥60 years) with AML and normal cardiac and renal function in the recent E-A NCTN E2906 phase 3 study. A patient-reported outcomes (PRO) assessment of Depression, Activities of Daily Living (ADL), Cognitive Function, Social Support (SS) domains, and comorbidity was performed at study registration to evaluate the prevalence of vulnerabilities and their impact on overall survival (OS) and 30-day induction mortality rates (IM) in this cohort selected for fitness to receive IC. Methods E2906 study design and results have been presented previously [n=727, age ≥60 years, randomized 1:1 to ‘standard 7&3‘ & high dose cytarabine (Arm A) vs. single agent clofarabine (CLO, Arm B), as remission induction (Step 1) and consolidation (Step 2)]. There was no difference in composite complete remission (CCR, 50%) or IM (8.5%), and CLO was inferior for OS. A baseline GA was a key secondary protocol objective, offered to all patients, and n=532 (73%) participated. We assessed Geriatric Depression Scale (GDS, none vs. any depressive symptoms, n=524); cognitive screen using Mini-Mental Exam Section (MMES, score 8/8 vs. <8, n=516); ADL & Instrumental ADL (iADL), (need for assistances vs. not n=526); a Medical Outcomes Study SS Survey (MOS-SS, measured as per unit increase, n=527); a Comorbidity survey ( ≤4 vs. ≥5 patient-reported comorbidities, n=530), patient-reported prescription medications use (PM, 0-3 vs. ≥4, n=519); and a survey of Sexual Desire (n=491) and Sexual Activity (n=477) within 4 weeks of AML diagnosis. Association of GA domains with OS, IM, and achievement of CCR was evaluated using the Chi-squared, Fisher's exact, and Wilcoxon rank sum test. Survival analysis was performed using the Cox proportional hazards models, adjusted for treatment arm. All p-values are 2-sided. Results The median age among those participating in GA battery was 67 yrs (range 60-85), including 38% age ≥70. There were 42% females, and ECOG performance status was 0 (28%), 1 (54%), 2 (16%), or 3 (2%). Central Cytogenetic review was classified (ELN2017) as Favorable (3%), Intermediate (66%), and Adverse (31%) risk. There was a notable rate of GA scores outside ‘normal’ clinical cutoff, including Depression (GDS 31%), cognitive and functional impairment [MMSE (37%), ADL (25%) and iADL (32%)], ≥5 Comorbidities (23%), and ≥4 PM's (44%). 81% & 76%, respectively, reported no sexual activity and low/none sexual desire. MOS-SS tended to be high in this cohort [mean (SD) score 87 (17)]. We observed a significant association of GA domains with OS & IM. For OS, absence of depressive symptoms [GDS, Hazard Ratio (HR) 0.79, 95% confidence Intervals (CI) 0.64-0.97, p=0.025)], MMES 8/8 (HR 0.78, 95%CI 0.64-0.96, p=0.016), requiring no assistance with ADLs (HR 0.65, 95%CI 0.52-0.82, p<0.001) and iADLs (HR 0.69, 95%CI 0.57-0.85, p<0.001), and ≤4 Comorbidities (HR 0.68, 95%CI 0.54-0.84, p=0.001) were significant for superior OS; use of ≥4 PM's was associated with inferior OS (HR 1.32, 95%CI 1.09-1.61, p=0.005). There was no association of sexual activity/desire, BMI ≥30 kg/m2, or MOS-SS with OS. IM was significantly higher among those with any depression symptoms (GDS, p=0.005), any ADL (p=0.010) or iADL (p=0.013) impairment, ≥5 comorbidities (p=0.037) and ≥4 PM's (p=0.001); there was no association with MOS-SS or MMES. There was no significant association of GA scores with cytogenetic risk group or with CCR rates. Only ADL independence (p=0.048), and possibly Comorbidity ≤4 (p=0.071), were associated with proceeding to Step 2. Conclusions In this large prospective cohort of fit older adults with AML selected to receive IC, we observed high rates of PRO-assessed GA vulnerability, which was significantly associated with inferior OS and higher IM rates, and (for ADL impairment) lower rates of proceeding to consolidation therapy. This suggests that functional and cognitive impairments are prevalent and can significantly impact outcome even in fit older adults. Despite potential long-term OS advantages of IC previously observed in E2906, efforts to systematically recognize those at higher risk based on GA PROs are needed to help identify patients who may be more appropriate for study of lower intensity treatment options to mitigate impact on IM and OS. SImilarly we plan to evaluate the potential benefit of GA-driven supportive care to improve patients' outcomes.

Abstract Background: Population registries and the Alliance co-operative group database have highlighted racial disparities in AML survival attributed to differential prevalence and impact of oncogenic drivers on disease biology. However, the small number of Black patients with well-annotated genomic profiles limits the validity of these findings. The current study leverages ECOG-ACRIN cytogenetic and outcomes data from newly diagnosed adult AML patients treated on intensive chemotherapy-based clinical trials over three decades. Methods: Self-reported race was categorized as White, Black or Other. Descriptive statistics were used to summarize baseline characteristics using the Wilcoxon rank-sum and Fisher's exact tests. Multivariable Cox proportional hazards regression models evaluated the association between race and survival outcomes adjusting for age, hemoglobin, platelet count, white blood cell count, blast percentage, performance status and ELN cytogenetic risk score (RS). All Cox models were stratified by protocol. Subgroup analyses of overall survival (OS) and disease-free survival (DFS) were performed using Kaplan-Meier curves stratified by race using log-rank tests. All statistical tests were two-sided. Analyses included only patients with complete data for the variables of interest. Results: Data was collected across 10 phase II and III interventional trials for newly diagnosed adult AML that accrued patients between 1984 and 2019. The study population included 3469 White, 184 Black, and 156 Other patients with a Black representation of 1.1-11% across trials. Blacks were significantly younger at diagnosis with a median age of 47.9 years compared to 53.5 years in Whites (p<0.001). Racial groups were well matched for gender and performance status. Black race was an independent prognostic factor for inferior OS (HR 1.212, CI 1.01 – 1.453; p = 0.0383) and DFS (HR 1.313, CI 1.05-1.641; p = 0.017). Complete cytogenetic results were available for 117 Black and 2162 White patients and revealed no significant racial differences in the distribution of ELN 2017 cytogenetic RS, although core binding factor AML was non-significantly enriched in Black patients (12.8% vs 8.3%). No difference was observed between Blacks and Whites in the prevalence of complex karyotype (CK), CK with -5/del(5q), -7/del(7q), or -17/del(17p), and AML-MRC including +8 and del(20q). Notably, the ELN 2017 cytogenetic RS effectively predicted OS in Black (57.8 vs 14.2 vs 5.8 months in favorable, intermediate and adverse RS, p = 0.018), and White patients (60.1 vs 15.2 vs 7.1 months, p < 0.0004 for Whites) but only predicted DFS for White AML patients (p < 0.0001 for Whites, p=0.48 for Blacks). In patients with molecular data for FLT3-ITD (n=1218), CEBPA (n=576), TP53 (n=572), or NPM1(n=683), no racial difference in the prevalence of these mutations was observed. However, NPM1 mutations were associated with inferior OS in Blacks compared to Whites (19.1 vs 8.9 months; p = 0.0095). Black and White AML patients had similar responses with complete remission rates of 58.2% vs 57.5%, p=0.55, and treatment tolerability with 30-day mortality of 7.1% vs 8.9%, p=0.50. However, consolidation approaches differed significantly. Among patients that underwent stem cell transplant (SCT), a significantly higher percentage of White patients received allogenic SCT compared to Blacks (48.5% versus 37.1%, p <0.01). Conclusions: This study represents one of the largest race-based comparisons of cytogenetic abnormalities and clinical outcomes in prospective trials of intensive chemotherapy in AML . Our data confirms the independent prognostic relevance of race in AML survival and validates the predictive ability of the ELN 2017 cytogenetic RS for OS across racial groups. There is ancestry-based attenuation of the survival benefit of NPM1 mutations as previously reported, underscoring the need to tailor the current treatment approach to NPM1 mutated AML based on race. The availability of NGS data limited the genomic analysis to a subset of patients. Future studies should examine the depth of treatment response incorporating MRD biomarkers and explore the effect of race on co-occurring mutations and the leukemic stem cell transcriptome to better explain the observed disparities. Finally, persistent disparities in access to allo-SCT even within the context of clinical trials reinforces the need to address diversity in the stem cell donor pool and social barriers to transplant.

Abstract Introduction Nadir bone marrow biopsies (BM Bx) are routinely performed around day 14 to assess early response to AML intensive induction therapy. The presence of residual disease prompts initiation of reinduction therapy, often with another cycle of intensive chemotherapy. Whether there are characteristics of nadir BM Bx that predict for response to intensive reinduction is unknown. To evaluate this, we analyzed data from ECOG-ACRIN clinical trials in AML of intensive induction chemotherapy with nadir biopsy information available to determine predictors of outcomes after reinduction. Methods We included patients from ECOG-ACRIN E1900 (Fernandez NEJM 2009) and arm A of E2906 (Foran ASH 2015) who received induction therapy with daunorubicin and cytarabine (7+3), had residual disease (blasts ≥5% or cellularity ≥20%) at nadir BM Bx and received reinduction chemotherapy. Nadir BM Bx characteristics that were evaluated included blast percentage and cellularity as well as both absolute and relative (percentage change, e.g. 20% to 10% = -50%) differences calculated between diagnosis and nadir for each parameter. Complete remission (CR) after reinduction therapy, disease free survival (DFS) and overall survival (OS) were outcomes for all analyses. Logistic regression models were used to evaluate associations with CR, and Cox proportional hazards models were used for DFS and OS. For models evaluating changes in blast percentage or cellularity as predictors, baseline blast percentage or cellularity was included as a covariate, respectively. Results Of 1013 patients identified, 750 (74%) had a nadir biopsy, 233 (23%) had evidence of residual disease, received reinduction therapy, and were included in the primary analysis. The median age was 58 and 42% were female. Cytogenetic risk category (by ELN2017) was adverse in 31%, intermediate in 66% and favorable in 3%. The median nadir BM Bx blast was 30% (range 0 to +97%) and cellularity was 20% (0 to +100%). The median absolute difference in blast percentage between baseline and nadir BM Bx was -15% (-95% to +66%) with a median relative difference of -36% (-100% to +400%). After reinduction, 118 patients (52%) achieved a CR. Neither absolute difference in blast percentage (OR 0.99, 95%CI 0.98-1.00, p=0.20) nor relative difference (OR 1.00, 95%CI 0.99-1.00, p=0.33) predicted CR. There was similarly no association with DFS (HR = 1.00, 95%CI: 0.99-1.01, p=0.65) or OS (HR = 1.00, 95%CI 0.99-1.01, p=0.76). Blast reduction ≥50% similarly did not predict these outcomes. We also evaluated the impact of any reduction in blast percentage (responders) from baseline to nadir and again found no association with CR (OR 1.65, 95%CI 0.90-3.03, p=0.11). Similarly, there was no association with DFS (HR 0.74, 95%CI 0.51-1.09, p=0.13) or OS (HR 0.81, 95%CI 0.58-1.13, p=0.21) in responders vs nonresponders. No association with CR, DFS or OS were identified when performing the same analysis in intermediate and adverse cytogenetic subgroups. We found no association with CR when nadir cellularity was assessed as a continuous variable (OR 0.99, 95%CI 0.98-1.00, p=0.10) or by relative difference (OR 0.99, 95%CI 0.99-1.00, p=0.10). There was a trend towards association with absolute reduction in cellularity from baseline to nadir and CR (OR 0.99, 95%CI 0.98-1.00, p=0.06). No association was identified between relative difference in cellularity and DFS and OS, nor was there an association among cytogenetic risk subgroups with any outcome. Discussion In intensively treated patients with residual disease at nadir BM Bx (i.e. day 14) who received intensive reinduction, the likelihood of achieving a CR was not influenced by the magnitude of disease reduction nor the extent of disease at the time of nadir biopsy. Similarly, there were no nadir biopsy features that predicted DFS or OS after intensive reinduction therapy. Given alternative reinduction therapies such as hypomethylating agents and venetoclax are available, these results suggest that the decision to attempt a second intensive induction cycle following the unequivocal presence of residual disease should not be influenced by characteristics of the nadir BM Bx. Although patient-specific factors such as patient fitness may favor a switch to less intensive options, these results support consideration of intensive reinduction therapy in patients with residual disease at nadir biopsy, even without changes in blasts or cellularity.

Background: The addition of midostaurin (M), an oral multi-kinase inhibitor, to intensive chemotherapy (IC) prolongs survival in newly diagnosed (ND) FLT3 mutated (m) AML. It is not known if there is a benefit to more potent and selective FLT3 inhibitors. Among these is Gilteritinib (G), now approved as a single agent in relapsed/refractory FLT3m AML. Combined with IC, G also leads to high measurable residual disease (MRD)negative (-) composite complete remission (CRc) in FLT3 Internal Tandem Duplication (ITD)m AML. The aim of the PrECOG 0905 study was to compare FLT3m clearance and CRc rates of G vs M with IC in adults with ND FLT3m AML. Methods: We conducted a randomized, open-label, phase 2 trial comparing G to M in combination with IC during induction and consolidation in ND AML pts with centrally PCR detected FLT3 ITD or TKD mutations. Pts age 18-70 with non M3 AML and no known core binding factor rearrangement were eligible. Pts were stratified by mutation type (TKD only vs. ITD (+/- TKD)), NPM1 mutation status, and FLT3-ITD allelic ratio. Induction consisted of cytarabine 100 mg/m2 by continuous infusion daily, on day (d) 1-7 and daunorubicin 90 mg/m2 IV on d 1-3. G 120 mg daily or M 50 mg twice daily was given orally on d8-21. Consolidation (up to 4 cycles) began within 60 d of induction with Cytarabine (1.5-3 g/m2 IV) X 6 doses and G or M (as per randomization) on d8-21. Pts could proceed to transplant (HCT) at any time. Maintenance was not included. The primary outcome was FLT3m MRD (-) CRc after induction. MRD for TKD mutations was tested by PCR (sensitivity 10-2); ITD MRD was tested by amplicon-based NGS (sensitivity 10-4). Safety and efficacy were analyzed in all eligible pts who received treatment. Secondary endpoints included comparison of CRc, MRD by flow cytometry (sensitivity 10-3), and survival. Results: From 11/2019 to 11/2022, 722 pts were screened at 37 centers, 180 pts randomized, and 177 pts (58% female, 77% white, 8.5% black, 8.5% hispanic) treated on the G (n=90) or M (n=87) arm. Median age was 54, with 27.7% > age 60. FLT3-ITD was present in 80% in the G arm and 78.2% in the M arm. Of those with FLT3-ITD, NPM1 was detected in 56.9% pts in G arm and 58.8% in M arm. 24.4% G and 20.7% M pts had co-occuring FLT3-ITD, NPM1 and DNMT3A mutations. Adverse karyotype was seen in 3.4% G and 9.3% M while 10% G and 13.8% M were ELN 2022 adverse risk. For pts on G arm, 85.6% achieved CRc compared to 72.4% for M (p=0.042). 5 (5.6%) pts on G arm and 6 (6.9%) on M arm received 2 cycles in induction. Post induction, the FLT3m- CRc rate was 40% for G vs 47.1% for M (p=0.366). FLT3-TKD only mutated pts achieved CRc in 88.9 % with G and 63.2% with M and FLT3m- CRc in 72.2% and 63.2% respectively. FLT3-ITDm (+/- TKDm) pts achieved CRc/FLT3m- CRc in 84.7%/31.9% with G and 75%/42.6% with M, respectively. In pts with FLT3-ITD, NPM1 and DNMT3A mutations, 95.5% G and 72.2% M achieved CRc and 50% G vs 33.3% on M achieved FLT3m- CRc. Flow cytometry MRD negative CRc was documented in 64.4% on G arm and 59.8% on M arm after induction (p=0.539). On G arm 67 (74%) patients received at least one cycle of consolidation compared to 53 (61%) pts on M arm. On exploratory analysis, 54 (66%) G and 40 (46%) M were reported to proceed to HCT in first remission. FLT3m testing was done post consolidation cycle 1 in 27 pts in FLT3m+ CRc after induction and 15/18 (83%) on G and 4/9 (44.4%) on M converted to FLT3m-. On multivariate logistic regression analysis for FLT3m-CRc post-induction, FLT3-ITD allelic ratio, NPM1m, WBC and hemoglobin at baseline were significant at the 0.10 two-sided level, while treatment arm, and 2017 ELN risk were not. There were no deaths from treatment-related (TR) adverse events (TRAEs) while grade >=3 TRAE was reported in 73%G vs 70% M during induction and 79%G vs 73%M pts during consolidation. Conclusion: Induction therapy with daunorubicin 90 mg/m2 X 3 days, cytarabine 100 mg/m2 and gilteritinib results in an excellent CRc rate with no TR deaths in pts with ND FLT3m AML up to age 70. Compared to M, G increased the CRc but not the FLT3m- CRc rate after induction. More post induction FLT3m+CRc pts who received G became FLT3m- after consolidation cycle 1 and more G pts succeeded in proceeding to HCT. Future survival data will help evaluate the clinical impact of G and MRD assessments in ND FLT3mAML. Larger studies will be needed to establish most predictive timing of MRD and for definitive comparisons of these drugs in patients with specific mutation profiles ( ie TKD, ITD/NPM1/DNMT3).

The progression of autosomal dominant polycystic kidney disease (ADPKD), an inherited kidney disease, is associated with renal interstitial inflammation and fibrosis. CD74 has been known not only as a receptor of macrophage migration inhibitory factor (MIF) it can also have MIF independent functions. In this study, we report unknown roles and function of CD74 in ADPKD. We show that knockout of CD74 delays cyst growth in Pkd1 mutant kidneys. Knockout and knockdown of CD74 (1) normalize PKD associated signaling pathways, including ERK, mTOR and Rb to decrease Pkd1 mutant renal epithelial cell proliferation, (2) decrease the activation of NF-κB and the expression of MCP-1 and TNF-alpha (TNF-α) which decreases the recruitment of macrophages in Pkd1 mutant kidneys, and (3) decrease renal fibrosis in Pkd1 mutant kidneys. We show for the first time that CD74 functions as a transcriptional factor to regulate the expression of fibrotic markers, including collagen I (Col I), fibronectin, and α-smooth muscle actin (α-SMA), through binding on their promoters. Interestingly, CD74 also regulates the transcription of MIF to form a positive feedback loop in that MIF binds with its receptor CD74 to regulate the activity of intracellular signaling pathways and CD74 increases the expression of MIF in ADPKD kidneys during cyst progression. We further show that knockout of MIF and targeting MIF with its inhibitor ISO-1 not only delay cyst growth but also ameliorate renal fibrosis through blocking the activation of renal fibroblasts and CD74 mediated the activation of TGF-β-Smad3 signaling, supporting the idea that CD74 is a key and novel upstream regulator of cyst growth and interstitial fibrosis. Thus, targeting MIF-CD74 axis is a novel therapeutic strategy for ADPKD treatment.

Background There is limited prospective or comparative data evaluating the HRQoL impact of intensive chemotherapy (IC) in older adults (age ≥60 years) with AML. In the recent EA E2906 phase 3 study of '7&3' (daunorubicin & cytarabine) vs. clofarabine (CLO), a putative lower intensity therapy, we incorporated a prospective patient-reported outcomes (PRO) assessment of HRQoL and fatigue to understand treatment impact from the patients' (pts) perspective. Methods E2906 study design and results have been presented previously [n=727, med. age 68 years (range 60-86), randomized 1:1 to ‘standard’ ‘7&3‘ & high dose cytarabine (Arm A) vs. single agent CLO (Arm B), as remission induction (Step 1) and consolidation (Step 2), respectively]. There was no difference in composite complete remission (CCR, 50%) or 30-day mortality rates (8.5%), and in primary analysis CLO was inferior for Overall Survival. HRQoL was a key secondary protocol endpoint. PROs assessed HRQoL and leukemia-specific symptoms and concerns using Functional Assessment of Cancer Therapy-Leukemia scale (FACT-Leu), & Fatigue using FACT-Fatigue subscale. PROs were administered at prespecified time points: (1) Baseline; (2) Cycle 1 at nadir (approx. Day 15); (3) End of Step 1 restaging (approx. D 40); (4) Beginning of Step 2 (1 st consolidation, approx. D 45); and (5) 2 nd cycle consolidation/End Step 2 (approx. D 60). The FACT-Leu Trial Outcome Index ( TOI) was calculated by summing physical well-being, functional well-being, and leukemia-specific concerns subscales to create a single composite score. Higher TOI scores correspond to better HRQoL. The primary endpoint was to evaluate TOI difference from randomization to day 30 after induction therapy between arms, evaluated using Wilcoxon rank sum test. We also compared Clinically Meaningful differences (CMiD) in TOI (defined as patient change ≥½ standard deviation from baseline TOI) between arms, using the X2 test. Pts not completing at least 2 assessment forms (including baseline) were excluded. Testing of constructs of Physical Well-being (PWB), Functional, Leukemia & Fatigue scales demonstrated high Cronbach's a values (a >0.7) except for PWB (a=0.672) at a single time point (End of Step 2). Results Baseline TOI was available for n=489, and 364 pts (n=182 each Arm) were evaluable in this planned analysis. Evaluating CMiD during Cycle 1 induction, pts in Arm B were less likely to experience ‘Worse’ TOI score (32% vs. 40.7%, Arm A), more likely to maintain stable TOI scores (49.1% vs. 39%, Arm A), and had a similar proportion of ‘Better’ TOI scores (20.3%, A vs. 18.9%, B) although differences were not significant (p=0.1389). Conversely in cycle 1 consolidation, more patients in Arm A experienced Better TOI (45.1% vs. 29.2%, Arm B), fewer had stable TOI (45.1% vs. 58.3%, Arm B), and rates of ‘Worse’ TOI were similar (9.8%, A vs 12.5%, Arm B), again differences were not significant (p=0.2616). There was no difference in Fatigue scores between Arms (p=0.49). We evaluated whether TOI and Fatigue scores could predict the ability to complete protocol therapy. In univariate analysis baseline TOI (p=0.0059) and Fatigue (p=0.0057) were associated with Step 1 completion, however on multivariate analysis (adjusting for age, sex, and Arm) were no longer significant [TOI, Odds Ratio (OR) 1.03 (0.98-1.08, p=0.27); Fatigue, OR 1.03 (0.96-1.11, p=0.42)]. There was a borderline association of TOI with Step 2 completion [TOI OR 1.10 (1.00-1.23, p=0.06)]. Using a multivariable linear mixed effects model, we observed similar and significant changes in TOI scores in both Arms, demonstrating the effect of IC & CLO on HRQoL (Figure). Comparing to baseline, TOI scores decreased significantly at C1 induction nadir (-4.0873, p<0.0001), and increased by 5.0762 (p=0.0004) at end Step 1, by 6.9047 at cycle 1 consolidation (p<0.0001), and by 11.4601 (p<0.0001) at end Step 2. Differences were mainly apparent in those who achieved CCR, and there was significant attrition in PRO responses in non-CCR pts. Conclusions Older adults experienced significant improvement in HRQoL after nadir regardless of therapy, and those achieving CCR report significantly higher scores. We did not observe differences in leukemia specific HRQoL or fatigue between treatment Arms (7&3 vs. CLO). Patient-reported HRQoL & Fatigue was not predictive of completion of protocol therapy. These results serve as a benchmark for expectations with IC in older adults.

Normalization is the first critical step in microbiome sequencing data analysis used to account for variable library sizes. Current RNA-Seq based normalization methods that have been adapted for microbiome data fail to consider the unique characteristics of microbiome data, which contain a vast number of zeros due to the physical absence or under-sampling of the microbes. Normalization methods that specifically address the zero-inflation remain largely undeveloped. Here we propose geometric mean of pairwise ratios—a simple but effective normalization method—for zero-inflated sequencing data such as microbiome data. Simulation studies and real datasets analyses demonstrate that the proposed method is more robust than competing methods, leading to more powerful detection of differentially abundant taxa and higher reproducibility of the relative abundances of taxa.

Normalization is the first critical step in microbiome sequencing data analysis used to account for variable library sizes. Current RNA-Seq based normalization methods that have been adapted for microbiome data fail to consider the unique characteristics of microbiome data, which contain a vast number of zeros due to the physical absence or under-sampling of the microbes. Normalization methods that specifically address the zero inflation remain largely undeveloped. Here we propose GMPR - a simple but effective normalization method - for zero-inflated sequencing data such as microbiome data. Simulation studies and real datasets analyses demonstrate that the proposed method is more robust than competing methods, leading to more powerful detection of differentially abundant taxa and higher reproducibility of the relative abundances of taxa.

Normalization is the first critical step in microbiome sequencing data analysis used to account for variable library sizes. Current RNA-Seq based normalization methods that have been adapted for microbiome data fail to consider the unique characteristics of microbiome data, which contain a vast number of zeros due to the physical absence or under-sampling of the microbes. Normalization methods that specifically address the zero inflation remain largely undeveloped. Here we propose GMPR - a simple but effective normalization method - for zero-inflated sequencing data such as microbiome data. Simulation studies and real datasets analyses demonstrate that the proposed method is more robust than competing methods, leading to more powerful detection of differentially abundant taxa and higher reproducibility of the relative abundances of taxa.