Double-Negative Results Matter: A Reevaluation of Sensitivities for Detecting SARS-CoV-2 Infection Using Saliva Versus Nasopharyngeal Swabs.

In a recent systematic review, Bastos et al. (Ann Intern Med. 2021;174(4):501-510) compared the sensitivities of saliva sampling and nasopharyngeal swabs in the detection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection by assuming a composite reference standard defined as positive if either test is positive and negative if both tests are negative (double negative). Even under a perfect specificity assumption, this approach ignores the double-negative results and risks overestimating the sensitivities due to residual misclassification. In this article, we first illustrate the impact of double-negative results in the estimation of the sensitivities in a single study, and then propose a 2-step latent class meta-analysis method for reevaluating both sensitivities using the same published data set as that used in Bastos et al. by properly including the observed double-negative results. We also conduct extensive simulation studies to compare the performance of the proposed method with Bastos et al.'s method for varied levels of prevalence and between-study heterogeneity. The results demonstrate that the sensitivities are overestimated noticeably using Bastos et al.'s method, and the proposed method provides a more accurate evaluation with nearly no bias and close-to-nominal coverage probability. In conclusion, double-negative results can significantly impact the estimated sensitivities when a gold standard is absent, and thus they should be properly incorporated.

Influence of lack of blinding on the estimation of medication-related harms: a retrospective cohort study of randomized controlled trials.

BACKGROUND: Empirical evidence suggests that lack of blinding may be associated with biased estimates of treatment benefit in randomized controlled trials, but the influence on medication-related harms is not well-recognized. We aimed to investigate the association between blinding and clinical trial estimates of medication-related harms. METHODS: We searched PubMed from January 1, 2015, till January 1, 2020, for systematic reviews with meta-analyses of medication-related harms. Eligible meta-analyses must have contained trials both with and without blinding. Potential covariates that may confound effect estimates were addressed by restricting trials within the comparison or by hierarchical analysis of harmonized groups of meta-analyses (therefore harmonizing drug type, control, dosage, and registration status) across eligible meta-analyses. The weighted hierarchical linear regression was then used to estimate the differences in harm estimates (odds ratio, OR) between trials that lacked blinding and those that were blinded. The results were reported as the ratio of OR (ROR) with its 95% confidence interval (CI). RESULTS: We identified 629 meta-analyses of harms with 10,069 trials. We estimated a weighted average ROR of 0.68 (95% CI: 0.53 to 0.88, P < 0.01) among 82 trials in 20 meta-analyses where blinding of participants was lacking. With regard to lack of blinding of healthcare providers or outcomes assessors, the RORs were 0.68 (95% CI: 0.53 to 0.87, P < 0.01 from 81 trials in 22 meta-analyses) and 1.00 (95% CI: 0.94 to 1.07, P = 0.94 from 858 trials among 155 meta-analyses) respectively. Sensitivity analyses indicate that these findings are applicable to both objective and subjective outcomes. CONCLUSIONS: Lack of blinding of participants and health care providers in randomized controlled trials may underestimate medication-related harms. Adequate blinding in randomized trials, when feasible, may help safeguard against potential bias in estimating the effects of harms.

A Bayesian nonparametric meta-analysis model for estimating the reference interval.

A reference interval represents the normative range for measurements from a healthy population. It plays an important role in laboratory testing, as well as in differentiating healthy from diseased patients. The reference interval based on a single study might not be applicable to a broader population. Meta-analysis can provide a more generalizable reference interval based on the combined population by synthesizing results from multiple studies. However, the assumptions of normally distributed underlying study-specific means and equal within-study variances, which are commonly used in existing methods, are strong and may not hold in practice. We propose a Bayesian nonparametric model with more flexible assumptions to extend random effects meta-analysis for estimating reference intervals. We illustrate through simulation studies and two real data examples the performance of our proposed approach when the assumptions of normally distributed study means and equal within-study variances do not hold.

Association Between Maternal Depression and Lower Urinary Tract Symptoms in Their Primary School-Age Daughters: A Birth Cohort Study.

PURPOSE: Although maternal depression is associated with adverse outcomes in women and children, its relationship with lower urinary tract symptoms (LUTS) in offspring is less well-characterized. We examined the association between prenatal and postpartum maternal depression and LUTS in primary school-age daughters. DESIGN: Observational cohort study. SUBJECTS AND SETTING: The sample comprised 7148 mother-daughter dyads from the Avon Longitudinal Study of Parents and Children. METHOD: Mothers completed questionnaires about depressive symptoms at 18 and 32 weeks' gestation and 21 months postpartum and their children's LUTS (urinary urgency, nocturia, and daytime and nighttime wetting) at 6, 7, and 9 years of age. Multivariable logistic regression models were used to estimate the association between maternal depression and LUTS in daughters. RESULTS: Compared to daughters of mothers without depression, those born to mothers with prenatal and postpartum depression had higher odds of LUTS, including urinary urgency (adjusted odds ratio [aOR] range = 1.99-2.50) and nocturia (aOR range = 1.67-1.97) at 6, 7, and 9 years of age. Additionally, daughters born to mothers with prenatal and postpartum depression had higher odds of daytime wetting (aOR range = 1.81-1.99) and nighttime wetting (aOR range = 1.63-1.95) at 6 and 7 years of age. Less consistent associations were observed for depression limited to the prenatal or postpartum periods only. CONCLUSIONS: Exposure to maternal depression in the prenatal and postpartum periods was associated with an increased likelihood of LUTS in daughters. This association may be an important opportunity for childhood LUTS prevention. Prevention strategies should reflect an understanding of potential biological and environmental mechanisms through which maternal depression may influence childhood LUTS.

Assessing the robustness of results from clinical trials and meta-analyses with the fragility index.

The fragility index has been increasingly used to assess the robustness of the results of clinical trials since 2014. It aims at finding the smallest number of event changes that could alter originally statistically significant results. Despite its popularity, some researchers have expressed several concerns about the validity and usefulness of the fragility index. It offers a comprehensive review of the fragility index's rationale, calculation, software, and interpretation, with emphasis on application to studies in obstetrics and gynecology. This article presents the fragility index in the settings of individual clinical trials, standard pairwise meta-analyses, and network meta-analyses. Moreover, this article provides worked examples to demonstrate how the fragility index can be appropriately calculated and interpreted. In addition, the limitations of the traditional fragility index and some solutions proposed in the literature to address these limitations were reviewed. In summary, the fragility index is recommended to be used as a supplemental measure in the reporting of clinical trials and a tool to communicate the robustness of trial results to clinicians. Other considerations that can aid in the fragility index's interpretation include the loss to follow-up and the likelihood of data modifications that achieve the loss of statistical significance.

Accounting for post-randomization variables in meta-analysis: A joint meta-regression approach.

Meta-regression is widely used in systematic reviews to investigate sources of heterogeneity and the association of study-level covariates with treatment effectiveness. Existing meta-regression approaches are successful in adjusting for baseline covariates, which include real study-level covariates (e.g., publication year) that are invariant within a study and aggregated baseline covariates (e.g., mean age) that differ for each participant but are measured before randomization within a study. However, these methods have several limitations in adjusting for post-randomization variables. Although post-randomization variables share a handful of similarities with baseline covariates, they differ in several aspects. First, baseline covariates can be aggregated at the study level presumably because they are assumed to be balanced by the randomization, while post-randomization variables are not balanced across arms within a study and are commonly aggregated at the arm level. Second, post-randomization variables may interact dynamically with the primary outcome. Third, unlike baseline covariates, post-randomization variables are themselves often important outcomes under investigation. In light of these differences, we propose a Bayesian joint meta-regression approach adjusting for post-randomization variables. The proposed method simultaneously estimates the treatment effect on the primary outcome and on the post-randomization variables. It takes into consideration both between- and within-study variability in post-randomization variables. Studies with missing data in either the primary outcome or the post-randomization variables are included in the joint model to improve estimation. Our method is evaluated by simulations and a real meta-analysis of major depression disorder treatments.

Bayesian hierarchical models incorporating study-level covariates for multivariate meta-analysis of diagnostic tests without a gold standard with application to COVID-19.

When evaluating a diagnostic test, it is common that a gold standard may not be available. One example is the diagnosis of SARS-CoV-2 infection using saliva sampling or nasopharyngeal swabs. Without a gold standard, a pragmatic approach is to postulate a "reference standard," defined as positive if either test is positive, or negative if both are negative. However, this pragmatic approach may overestimate sensitivities because subjects infected with SARS-CoV-2 may still have double-negative test results even when both tests exhibit perfect specificity. To address this limitation, we propose a Bayesian hierarchical model for simultaneously estimating sensitivity, specificity, and disease prevalence in the absence of a gold standard. The proposed model allows adjusting for study-level covariates. We evaluate the model performance using an example based on a recently published meta-analysis on the diagnosis of SARS-CoV-2 infection and extensive simulations. Compared with the pragmatic reference standard approach, we demonstrate that the proposed Bayesian method provides a more accurate evaluation of prevalence, specificity, and sensitivity in a meta-analytic framework.

The association of sensitivity and specificity with disease prevalence: analysis of 6909 studies of diagnostic test accuracy.

BACKGROUND: Sensitivity and specificity are characteristics of a diagnostic test and are not expected to change as the prevalence of the target condition changes. We sought to evaluate the association between prevalence and changes in sensitivity and specificity. METHODS: We retrieved data from meta-analyses of diagnostic test accuracy published in the Cochrane Database of Systematic Reviews (2003-2020). We used mixed-effects random-intercept linear regression models to evaluate the association between prevalence and logit-transformed sensitivity and specificity. The model evaluated all meta-analyses as nested within each systematic review. RESULTS: We analyzed 6909 diagnostic test accuracy studies from 552 meta-analyses that were included in 92 systematic reviews. For sensitivity, compared with the lowest quartile of prevalence, the second, third and fourth quartiles were associated with significantly higher odds of identifying a true positive case (odds ratio [OR] 1.17, 95% confidence interval [CI] 1.09-1.26; OR 1.32, 95% CI 1.23-1.41; OR 1.47, 95% CI 1.37-1.58; respectively). For specificity, compared with the lowest quartile of prevalence, the second, third and fourth quartiles were associated with significantly lower odds of identifying a true negative case (OR 0.74, 95% CI 0.69-0.80; OR 0.65, 95% CI 0.60-0.70; OR 0.47, 95% CI 0.44-0.51; respectively). Pooled regression coefficients from bivariate models conducted within each meta-analysis showed that prevalence was positively associated with sensitivity and negatively associated with specificity. Findings were consistent across subgroups. INTERPRETATION: In this large sample of diagnostic studies, higher prevalence was associated with higher estimated sensitivity and lower estimated specificity. Clinicians should consider the implications of disease prevalence and spectrum when interpreting the results from studies of diagnostic test accuracy.

A Guide to Estimating the Reference Range From a Meta-Analysis Using Aggregate or Individual Participant Data.

Clinicians frequently must decide whether a patient's measurement reflects that of a healthy "normal" individual. Thus, the reference range is defined as the interval in which some proportion (frequently 95%) of measurements from a healthy population is expected to fall. One can estimate it from a single study or preferably from a meta-analysis of multiple studies to increase generalizability. This range differs from the confidence interval for the pooled mean and the prediction interval for a new study mean in a meta-analysis, which do not capture natural variation across healthy individuals. Methods for estimating the reference range from a meta-analysis of aggregate data that incorporates both within- and between-study variations were recently proposed. In this guide, we present 3 approaches for estimating the reference range: one frequentist, one Bayesian, and one empirical. Each method can be applied to either aggregate or individual-participant data meta-analysis, with the latter being the gold standard when available. We illustrate the application of these approaches to data from a previously published individual-participant data meta-analysis of studies measuring liver stiffness by transient elastography in healthy individuals between 2006 and 2016.

Estimating the Complier Average Causal Effect in a Meta-Analysis of Randomized Clinical Trials With Binary Outcomes Accounting for Noncompliance: A Generalized Linear Latent and Mixed Model Approach.

Noncompliance, a common problem in randomized clinical trials (RCTs), can bias estimation of the effect of treatment receipt using a standard intention-to-treat analysis. The complier average causal effect (CACE) measures the effect of an intervention in the latent subpopulation that would comply with their assigned treatment. Although several methods have been developed to estimate the CACE in analyzing a single RCT, methods for estimating the CACE in a meta-analysis of RCTs with noncompliance await further development. This article reviews the assumptions needed to estimate the CACE in a single RCT and proposes a frequentist alternative for estimating the CACE in a meta-analysis, using a generalized linear latent and mixed model with SAS software (SAS Institute, Inc.). The method accounts for between-study heterogeneity using random effects. We implement the methods and describe an illustrative example of a meta-analysis of 10 RCTs evaluating the effect of receiving epidural analgesia in labor on cesarean delivery, where noncompliance varies dramatically between studies. Simulation studies are used to evaluate the performance of the proposed method.

Empirical Comparisons of 12 Meta-analysis Methods for Synthesizing Proportions of Binary Outcomes.

BACKGROUND: Meta-analysis is increasingly used to synthesize proportions (e.g., disease prevalence). It can be implemented with widely used two-step methods or one-step methods, such as generalized linear mixed models (GLMMs). Existing simulation studies have shown that GLMMs outperform the two-step methods in some settings. It is, however, unclear whether these simulation settings are common in the real world. We aim to compare the real-world performance of various meta-analysis methods for synthesizing proportions. METHODS: We extracted datasets of proportions from the Cochrane Library and applied 12 two-step and one-step methods to each dataset. We used Spearman's ρ and the Bland-Altman plot to assess their results' correlation and agreement. The GLMM with the logit link was chosen as the reference method. We calculated the absolute difference and fold change (ratio of estimates) of the overall proportion estimates produced by each method vs. the reference method. RESULTS: We obtained a total of 43,644 datasets. The various methods generally had high correlations (ρ > 0.9) and agreements. GLMMs had computational issues more frequently than the two-step methods. However, the two-step methods generally produced large absolute differences from the GLMM with the logit link for small total sample sizes (< 50) and crude event rates within 10-20% and 90-95%, and large fold changes for small total event counts (< 10) and low crude event rates (< 20%). CONCLUSIONS: Although different methods produced similar overall proportion estimates in most datasets, one-step methods should be considered in the presence of small total event counts or sample sizes and very low or high event rates.

A penalization approach to random-effects meta-analysis.

Systematic reviews and meta-analyses are principal tools to synthesize evidence from multiple independent sources in many research fields. The assessment of heterogeneity among collected studies is a critical step when performing a meta-analysis, given its influence on model selection and conclusions about treatment effects. A common-effect (CE) model is conventionally used when the studies are deemed homogeneous, while a random-effects (RE) model is used for heterogeneous studies. However, both models have limitations. For example, the CE model produces excessively conservative confidence intervals with low coverage probabilities when the collected studies have heterogeneous treatment effects. The RE model, on the other hand, assigns higher weights to small studies compared to the CE model. In the presence of small-study effects or publication bias, the over-weighted small studies from a RE model can lead to substantially biased overall treatment effect estimates. In addition, outlying studies may exaggerate between-study heterogeneity. This article introduces penalization methods as a compromise between the CE and RE models. The proposed methods are motivated by the penalized likelihood approach, which is widely used in the current literature to control model complexity and reduce variances of parameter estimates. We compare the existing and proposed methods with simulated data and several case studies to illustrate the benefits of the penalization methods.

Accounting for publication bias using a bivariate trim and fill meta-analysis procedure.

In research synthesis, publication bias (PB) refers to the phenomenon that the publication of a study is associated with the direction and statistical significance of its results. Consequently, it may lead to biased (commonly optimistic) estimates of treatment effects. Visualization tools such as funnel plots have been widely used to investigate PB in univariate meta-analyses. The trim and fill procedure is a nonparametric method to identify and adjust for PB. It is popular among applied scientists due to its simplicity. However, most visualization tools and PB correction methods focus on univariate outcomes. For a meta-analysis with multiple outcomes, the conventional univariate trim and fill method can only account for different outcomes separately and thus may lead to inconsistent conclusions. In this article, we propose a bivariate trim and fill procedure to simultaneously account for PB in the presence of two outcomes that are possibly associated. Based on a recently developed galaxy plot for bivariate meta-analysis, the proposed procedure uses a data-driven imputation algorithm to detect and adjust PB. The method relies on the symmetry of the galaxy plot and assumes that some studies are suppressed based on a linear combination of outcomes. The method projects bivariate outcomes along a particular direction, uses the univariate trim and fill method to estimate the number of trimmed and filled studies, and yields consistent conclusions about PB. The proposed approach is validated using simulated data and is applied to a meta-analysis of the efficacy and safety of antidepressant drugs.

A Bayesian hierarchical model for individual participant data meta-analysis of demand curves.

Individual participant data meta-analysis is a frequently used method to combine and contrast data from multiple independent studies. Bayesian hierarchical models are increasingly used to appropriately take into account potential heterogeneity between studies. In this paper, we propose a Bayesian hierarchical model for individual participant data generated from the Cigarette Purchase Task (CPT). Data from the CPT details how demand for cigarettes varies as a function of price, which is usually described as an exponential demand curve. As opposed to the conventional random-effects meta-analysis methods, Bayesian hierarchical models are able to estimate both the study-specific and population-level parameters simultaneously without relying on the normality assumptions. We applied the proposed model to a meta-analysis with baseline CPT data from six studies and compared the results from the proposed model and a two-step conventional random-effects meta-analysis approach. We conducted extensive simulation studies to investigate the performance of the proposed approach and discussed the benefits of using the Bayesian hierarchical model for individual participant data meta-analysis of demand curves.

Postextraction ridge preservation by using dense PTFE membranes: A systematic review and meta-analysis.

STATEMENT OF PROBLEM: The use of dense polytetrafluoroethylene (dPTFE) membranes in alveolar ridge preservation may help reduce the risk of bacterial contamination and infection, maintaining the soft-tissue anatomy. However, systematic reviews on their efficacy in postextraction sites are lacking. PURPOSE: The purpose of this systematic review and meta-analysis was to assess the efficacy of alveolar ridge preservation with dPTFE membranes when used alone or in combination with bone grafting materials in postextraction sites. MATERIAL AND METHODS: An electronic search up to February 2021 was conducted by using PubMed, Embase, and the Cochrane library to detect studies using dPTFE membranes in postextraction sites. An additional manual search was performed in relevant journals. Clinical and radiographic dimensional changes of the alveolar ridge, histomorphometric, microcomputed tomography, implant-related findings, and rate of complications were recorded. One-dimensional meta-analysis was performed to calculate the overall means and 95% confidence intervals (α=.05). RESULTS: A total of 23 studies, 14 randomized controlled trials, 4 retrospective cohort studies, 3 case series, and 2 prospective nonrandomized clinical trials, met the inclusion criteria. Five studies were included in the quantitative analysis. The meta-analysis revealed that the use of dPTFE membranes resulted in a statistically significant (P=.042) increase in clinical keratinized tissue of 3.49 mm (95% confidence interval [CI]: 0.16, 6.83) when compared with extraction alone. Metaregression showed that the difference of 1.10 mm (95% CI: -0.14, 2.35) in the radiographic horizontal measurements was not significant (P=.082), but the difference of 1.06 mm (95% CI: 0.51, 1.62) in the radiographic vertical dimensional change between dPTFE membranes+allograft and extraction alone was statistically significant (P<.001). CONCLUSIONS: The use of dPTFE membranes was better than extraction alone in terms of keratinized tissue width and radiographic vertical bone loss.

Effectiveness of Messenger RNA Coronavirus Disease 2019 (COVID-19) Vaccines Against Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Infection in a Cohort of Healthcare Personnel.

In a large cohort of United States healthcare personnel without prior coronavirus disease 2019 (COVID-19) infection, 94 382 doses of messenger RNA (mRNA) COVID-19 vaccine were administered to 49 220 individuals. The adjusted vaccine effectiveness following 2 doses of each of the 2 available brands of mRNA vaccine exceeded 96%.

Meta-Analysis and Sparse-Data Bias.

Meta-analyses are undertaken to combine information from a set of studies, often in settings where some of the individual study-specific estimates are based on relatively small study samples. Finite sample bias may occur when maximum likelihood estimates of associations are obtained by fitting logistic regression models to sparse data sets. Here we show that combining information from small studies by undertaking a meta-analytical summary of logistic regression estimates can propagate such sparse-data bias. In simulations, we illustrate 2 challenges encountered in meta-analyses of logistic regression results in settings of sparse data: 1) bias in the summary meta-analytical result and 2) confidence interval coverage that can worsen rather than improve, in terms of being less than nominal, as the number of studies in the meta-analysis increases.

Identification of novel candidate genes in rosacea by bioinformatic methods.

BACKGROUND: Rosacea is a chronic inflammatory skin disease whose psychological consequences severely affect patient's quality of life. OBJECTIVE: To identify candidate genes of rosacea for potential development of new target therapies. METHODS: Gene Expression Omnibus datasets were retrieved to obtain differentially expressed genes (DEGs) between rosacea patients and healthy controls. Gene ontology (GO) analyses were used to identify functions of candidate genes. Related signaling pathways of DEGs were analyzed using Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene set enrichment analysis. Protein-protein interaction (PPI) networks were applied using search tools for the retrieval of interacting genes/proteins and modulations involving PPI networks were evaluated with use of the MCODE app. RESULTS: Samples from 19 rosacea patients and 10 healthy controls of dataset GSE65914 were enrolled. A total of 215 DEGs, 115 GO terms and 6 KEGG pathways were identified. A total of 182 nodes and 456 edges were enriched in PPI networks. Maximal clusters showed 15 central nodes and 96 edges. The toll-like receptor (TLR) signaling pathway was the most significant pathway detected and 5 DEGs were identified as candidate genes which included TLR2, C-C motif chemokine (CCL) 5, C-X-C motif chemokine ligand (CXCL) 9, CXCL10 and CXCL11. The results were verified in rosacea patients with use of real-time polymerase chain reaction and immunohistochemistry. Cell-type enrichment analysis revealed 8 lymphocytes that were enriched in rosacea patients. CONCLUSIONS: The results suggest that both innate and adaptive immune responses were involved in the etiology of rosacea. Five DEGs in the TLR signaling pathway may serve as potential therapeutic target genes.

Prior Choices of Between-Study Heterogeneity in Contemporary Bayesian Network Meta-analyses: an Empirical Study.

BACKGROUND: Network meta-analysis (NMA) is a popular tool to compare multiple treatments in medical research. It is frequently implemented via Bayesian methods. The prior choice of between-study heterogeneity is critical in Bayesian NMAs. This study evaluates the impact of different priors for heterogeneity on NMA results. METHODS: We identified all NMAs with binary outcomes published in The BMJ, JAMA, and The Lancet during 2010-2018, and extracted information about their prior choices for heterogeneity. Our primary analyses focused on those with publicly available full data. We re-analyzed the NMAs using 3 commonly-used non-informative priors and empirical informative log-normal priors. We obtained the posterior median odds ratios and 95% credible intervals of all comparisons, assessed the correlation among different priors, and used Bland-Altman plots to evaluate their agreement. The kappa statistic was also used to evaluate the agreement among these priors regarding statistical significance. RESULTS: Among the selected Bayesian NMAs, 52.3% did not specify the prior choice for heterogeneity, and 84.1% did not provide rationales. We re-analyzed 19 NMAs with full data available, involving 894 studies, 173 treatments, and 395,429 patients. The correlation among posterior median (log) odds ratios using different priors were generally very strong for NMAs with over 20 studies. The informative priors produced substantially narrower credible intervals than non-informative priors, especially for NMAs with few studies. Bland-Altman plots and kappa statistics indicated strong overall agreement, but this was not always the case for a specific NMA. CONCLUSIONS: Priors should be routinely reported in Bayesian NMAs. Sensitivity analyses are recommended to examine the impact of priors, especially for NMAs with relatively small sample sizes. Informative priors may produce substantially narrower credible intervals for such NMAs.

A Plant-Centered Diet and Markers of Early Chronic Kidney Disease during Young to Middle Adulthood: Findings from the Coronary Artery Risk Development in Young Adults (CARDIA) Cohort.

BACKGROUND: Few studies have evaluated whether plant-centered diets prevent progression of early stage chronic kidney disease (CKD). OBJECTIVES: We examined the association between plant-centered diet quality and early CKD markers. METHODS: We prospectively examined 2869 black and white men and women in the Coronary Artery Risk Development in Young Adults Study free of diagnosed kidney failure in 2005-2006 [examination year 20 (Y20); mean age: 45.3 ± 3.6  y]. CKD marker changes from Y20 to 2015-2016 (Y30) were considered, including estimated glomerular filtration rate (eGFR; serum creatinine), urinary albumin-to-creatinine ratio (ACR), and both. Diet was assessed through interviewer-administered diet histories at Y0, Y7, and Y20, and plant-centered diet quality was quantified with the A Priori Diet Quality Score (APDQS). Linear regression models were used to examine the association of APDQS and subsequent 10-y changes in CKD markers. RESULTS: After adjustment for sociodemographic, behavioral, and diet factors, we found that higher APDQS was related to less adverse changes in CKD markers in the subsequent 10-y period. Compared with the lowest APDQS quintile, the highest quintile was associated with an attenuated increase in lnACR (-0.25 mg/g; 95% CI: -0.37, -0.13 mg/g; P-trend < 0.001), whereas the highest quintile was associated with an attenuated decrease in eGFR (4.45 mL·min-1·1.73 m-2; 95% CI: 2.46, 6.43 mL·min-1·1.73 m-2). There was a 0.50 lower increase in combined CKD markers [ln(ACR) z score - eGFR z score] when comparing the extreme quintiles. Associations remained similar after further adjustment for hypertension, diabetes, and obesity as potential mediating factors. The attenuated worsening CKD marker changes associated with higher APDQS strengthened across increasing initial CKD category; those with the best diet and microalbuminuria in Y10-Y20 returned to high normal albuminuria (all P-interaction < 0.001). CONCLUSIONS: Individuals who consumed plant-centered, high-quality diets were less likely to experience deterioration of kidney function through midlife, especially among participants with initial stage characterized as mild CKD.