Plasma proteins associated with plant-based diets: Results from the Atherosclerosis Risk in Communities (ARIC) study and Framingham Heart Study (FHS)

BACKGROUND & AIMS

Plant-based diets are associated with a lower risk of chronic diseases. Large-scale proteomics can identify objective biomarkers of plant-based diets, and improve our understanding of the pathways that link plant-based diets to health outcomes. This study investigated the plasma proteome of four different plant-based diets [overall plant-based diet (PDI), provegetarian diet, healthful plant-based diet (hPDI), and unhealthful plant-based diet (uPDI)] in the Atherosclerosis Risk in Communities (ARIC) Study and replicated the findings in the Framingham Heart Study (FHS) Offspring cohort.

METHODS

ARIC Study participants at visit 3 (1993-1995) with completed food frequency questionnaire (FFQ) data and proteomics data were divided into internal discovery (n = 7690) and replication (n = 2543) data sets. Multivariable linear regression was used to examine associations between plant-based diet indices (PDIs) and 4955 individual proteins in the discovery sample. Then, proteins that were internally replicated in the ARIC Study were tested for external replication in FHS (n = 1358). Pathway overrepresentation analysis was conducted for diet-related proteins. C-statistics were used to predict if the proteins improved prediction of plant-based diet indices beyond participant characteristics.

RESULTS

In ARIC discovery, a total of 837 diet-protein associations (PDI = 233; provegetarian = 182; hPDI = 406; uPDI = 16) were observed at false discovery rate (FDR) < 0.05. Of these, 453 diet-protein associations (PDI = 132; provegetarian = 104; hPDI = 208; uPDI = 9) were internally replicated. In FHS, 167/453 diet-protein associations were available for external replication, of which 8 proteins (PDI = 1; provegetarian = 0; hPDI = 8; uPDI = 0) replicated. Complement and coagulation cascades, cell adhesion molecules, and retinol metabolism were over-represented. C-C motif chemokine 25 for PDI and 8 proteins for hPDI modestly but significantly improved the prediction of these indices individually and collectively (P value for difference in C-statistics<0.05 for all tests).

CONCLUSIONS

Using large-scale proteomics, we identified potential candidate biomarkers of plant-based diets, and pathways that may partially explain the associations between plant-based diets and chronic conditions.

Plasma proteomic comparisons change as coverage expands for SomaLogic and Olink

The number of assays on highly-multiplexed proteomic platforms has grown ten-fold over the past 15 years from less than 1,000 to >11,000. The leading aptamer-based and antibody-based platforms have different strengths. For example, Eldjarn et al 1 demonstrated that the aptamer-based SomaScan 5k (4,907 assays, assessed in the Icelandic 36K) and the antibody-based Olink Explore 3072 (2,931 assays, assessed in the UK BioBank) had a similar number of cis -pQTLs among all targets (2,120 vs. 2,101) but Olink had a greater number of cis -pQTLs among the overlapping targets (1,164 vs. 1,467). Analysis of split plasma measures showed the SomaScan assays to be more precise: median coefficient of variation (CV) of 9.9% vs. 16.5% for Olink. 1 Precision of the newest versions of the platforms-SomaScan 11k (>11,000 assays, released in December 2023) and Olink Explore HT (>5,400 assays, released in July 2023)-has not yet been established. We assessed the reproducibility of the SomaScan 11k and Olink Explore HT using split plasma samples from 102 Atherosclerosis Risk in Communities (ARIC) Study participants. We found that the SomaScan 11k assays had a median CV of 6.8% (vs 6.6% for the subset of assays available on the SomaScan 5k) and the Olink Explore HT assays had a median CV of 35.7% (vs 19.8% for the subset of assays available on the Olink Explore 3072). Across Olink assays, the CVs were strongly negatively correlated with protein detectability, i.e., percent of samples above the limit of detection (LOD). For the 4,443 overlapping assays, the distribution of between-platform correlations was bimodal with a peak at r ∼0 and with another smaller peak at r ∼0.8. These findings on precision are consistent with the updated results by Eldjarn et al 1 but indicate that precision of these two leading platforms in human plasma has diverged as the number of included proteins has increased.

Association of Integrated Proteomic and Metabolomic Modules with Risk of Kidney Disease Progression

Key Points

Integrated analysis of proteome and metabolome identifies modules associated with CKD progression and kidney failure. Ephrin transmembrane proteins and podocyte-expressed CRIM1 and NPNT emerged as central components and warrant experimental and clinical investigation.

Background

Proteins and metabolites play crucial roles in various biological functions and are frequently interconnected through enzymatic or transport processes.

Methods

We present an integrated analysis of 4091 proteins and 630 metabolites in the Chronic Renal Insufficiency Cohort study (N =1708; average follow-up for kidney failure, 9.5 years, with 537 events). Proteins and metabolites were integrated using an unsupervised clustering method, and we assessed associations between clusters and CKD progression and kidney failure using Cox proportional hazards models. Analyses were adjusted for demographics and risk factors, including the eGFR and urine protein–creatinine ratio. Associations were identified in a discovery sample (random two thirds, n =1139) and then evaluated in a replication sample (one third, n =569).

Results

We identified 139 modules of correlated proteins and metabolites, which were represented by their principal components. Modules and principal component loadings were projected onto the replication sample, which demonstrated a consistent network structure. Two modules, representing a total of 236 proteins and 82 metabolites, were robustly associated with both CKD progression and kidney failure in both discovery and validation samples. Using gene set enrichment, several transmembrane-related terms were identified as overrepresented in these modules. Transmembrane–ephrin receptor activity displayed the largest odds (odds ratio=13.2, P value = 5.5×10−5). A module containing CRIM1 and NPNT expressed in podocytes demonstrated particularly strong associations with kidney failure (P value = 2.6×10−5).

Conclusions

This study demonstrates that integration of the proteome and metabolome can identify functions of pathophysiologic importance in kidney disease.

Plasma Proteins Associated with Chronic Histopathologic Lesions on Kidney Biopsy

Key Points

Proteomic profiling identified 35 blood proteins associated with chronic histopathologic lesions in the kidney. Testican-2 was expressed in the glomerulus, released by the kidney into circulation, and inversely associated with glomerulosclerosis severity. NELL1 was expressed in tubular epithelial cells, released by the kidney into circulation, and inversely associated with interstitial fibrosis and tubular atrophy severity.

Background

The severity of chronic histopathologic lesions on kidney biopsy is independently associated with higher risk of progressive CKD. Because kidney biopsies are invasive, identification of blood markers that report on underlying kidney histopathology has the potential to enhance CKD care.

Methods

We examined the association between 6592 plasma protein levels measured by aptamers and the severity of interstitial fibrosis and tubular atrophy (IFTA), glomerulosclerosis, arteriolar sclerosis, and arterial sclerosis among 434 participants of the Boston Kidney Biopsy Cohort. For proteins significantly associated with at least one histologic lesion, we assessed renal arteriovenous protein gradients among 21 individuals who had undergone invasive catheterization and assessed the expression of the cognate gene among 47 individuals with single-cell RNA sequencing data in the Kidney Precision Medicine Project.

Results

In models adjusted for eGFR, proteinuria, and demographic factors, we identified 35 proteins associated with one or more chronic histologic lesions, including 20 specific for IFTA, eight specific for glomerulosclerosis, and one specific for arteriolar sclerosis. In general, higher levels of these proteins were associated with more severe histologic score and lower eGFR. Exceptions included testican-2 and NELL1, which were associated with less glomerulosclerosis and IFTA, respectively, and higher eGFR; notably, both of these proteins demonstrated significantly higher levels from artery to renal vein, demonstrating net kidney release. In the Kidney Precision Medicine Project, 13 of the 35 protein hits had cognate gene expression enriched in one or more cell types in the kidney, including podocyte expression of select glomerulosclerosis markers (including testican-2) and tubular expression of several IFTA markers (including NELL1).

Conclusions

Proteomic analysis identified circulating proteins associated with chronic histopathologic lesions, some of which had concordant site-specific expression within the kidney.

Extreme Humid-Heat Exposure and Mortality Among Patients Receiving Dialysis

RATIONALE & OBJECTIVE

Exposure to extreme heat events has been linked to increased morbidity and mortality in the general population. Patients receiving maintenance dialysis may be vulnerable to greater risks from these events, but this is not well understood. We characterized the association of extreme heat events and the risk of death among patients receiving dialysis in the United States.

STUDY DESIGN

Retrospective cohort study.

SETTING & PARTICIPANTS

Data from the US Renal Data System were used to identify adults living in US urban settlements prone to extreme heat who initiated maintenance dialysis between 1997 and 2016.

EXPOSURE

An extreme heat event, defined as a time-updated heat index (a humid-heat metric) exceeding 40.6°C for≥2 days or 46.1°C for≥1day.

OUTCOME

Death.

ANALYTICAL APPROACH

Cox proportional hazards regression to estimate the elevation in risk of death during a humid-heat event adjusted for age, sex, year of dialysis initiation, dialysis modality, poverty level, and climate region. Interactions between humid-heat and these same factors were explored.

RESULTS

Among 945,251 adults in 245 urban settlements, the mean age was 63 years, and 44% were female. During a median follow-up period of 3.6 years, 498,049 adults were exposed to at least 1 of 7,154 extreme humid-heat events, and 500,025 deaths occurred. In adjusted models, there was an increased risk of death (hazard ratio 1.18 [95% CI, 1.15-1.20]) during extreme humid-heat exposure. The relative mortality risk was higher among patients living in the Southeast (P<0.001) compared with the Southwest.

LIMITATIONS

Possibility of exposure misclassification, did not account for land use and air pollution co-exposures.

CONCLUSIONS

This study suggests that patients receiving dialysis face an increased risk of death during extreme humid-heat exposure.

PLAIN-LANGUAGE SUMMARY

Patients who receive dialysis are vulnerable to extreme weather events, and rising global temperatures may bring more frequent extreme heat events. We sought to determine whether extreme heat exposure was associated with an increased risk of death in urban-dwelling patients receiving dialysis across the United States. We found that people receiving dialysis were more likely to die during extreme humid-heat events, defined by a heat index exceeding 40.6°C (105°F) for≥2 days or 46.1°C (115°F) for≥1day. These findings inform the nephrology community about the potential importance of protecting patients receiving maintenance dialysis from the risks associated with extreme heat.

Mid-life plasma proteins associated with late-life prefrailty and frailty: a proteomic analysis

Physical frailty is a syndrome that typically manifests in later life, although the pathogenic process causing physical frailty likely begins decades earlier. To date, few studies have examined the biological signatures in mid-life associated with physical frailty later in life. Among 4,189 middle-aged participants (57.8 ± 5.0 years, 55.8% women) from the Atherosclerosis Risk in Community (ARIC) study, we evaluated the associations of 4,955 plasma proteins (log 2-transformed and standardized) measured using the SomaScan platform with their frailty status approximately 20 years later. Using multinomial logistic regression models adjusting for demographics, health behaviors, kidney function, total cholesterol, and comorbidities, 12 and 221 proteins were associated with prefrailty and frailty in later life, respectively (FDR p < 0.05). Top frailty-associated proteins included neurocan core protein (NCAN, OR = 0.66), fatty acid-binding protein heart (FABP3, OR = 1.62) and adipocyte (FABP4, OR = 1.65), as well proteins involved in the contactin-1 (CNTN1), toll-like receptor 5 (TLR5), and neurogenic locus notch homolog protein 1 (NOTCH1) signaling pathway relevant to skeletal muscle regeneration, myelination, and inflammation. Pathway analyses suggest midlife dysregulation of inflammation, metabolism, extracellular matrix, angiogenesis, and lysosomal autophagy among those at risk for late-life frailty. After further adjusting for midlife body mass index (BMI) - an established frailty risk factor - only CNTN1 (OR = 0.75) remained significantly associated with frailty. Post-hoc analyses demonstrated that the top 41 midlife frailty-associated proteins mediate 32% of the association between mid-life BMI and late-life frailty. Our findings provide new insights into frailty etiology earlier in the life course, enhancing the potential for prevention.

Serum and Urine Metabolites and Kidney Function

Key Points

Background

Metabolites represent a read-out of cellular processes underlying states of health and disease.

Methods

We evaluated cross-sectional and longitudinal associations between 1255 serum and 1398 urine known and unknown (denoted with “X” in name) metabolites (Metabolon HD4, 721 detected in both biofluids) and kidney function in 1612 participants of the Atherosclerosis Risk in Communities study. All analyses were adjusted for clinical and demographic covariates, including for baseline eGFR and urine albumin-creatinine ratio (UACR) in longitudinal analyses.

Results

At visit 5 of the Atherosclerosis Risk in Communities study, the mean age of participants was 76 years (SD 6); 56% were women, mean eGFR was 62 ml/min per 1.73 m2 (SD 20), and median UACR level was 13 mg/g (interquartile range, 25). In cross-sectional analysis, 675 serum and 542 urine metabolites were associated with eGFR (Bonferroni-corrected P < 4.0E-5 for serum analyses and P < 3.6E-5 for urine analyses), including 248 metabolites shared across biofluids. Fewer metabolites (75 serum and 91 urine metabolites, including seven metabolites shared across biofluids) were cross-sectionally associated with albuminuria. Guanidinosuccinate; N2,N2-dimethylguanosine; hydroxy-N6,N6,N6-trimethyllysine; X-13844; and X-25422 were significantly associated with both eGFR and albuminuria. Over a mean follow-up of 6.6 years, serum mannose (hazard ratio [HR], 2.3 [1.6–3.2], P = 2.7E-5) and urine X-12117 (HR, 1.7 [1.3–2.2], P = 1.9E-5) were risk factors of UACR doubling, whereas urine sebacate (HR, 0.86 [0.80–0.92], P = 1.9E-5) was inversely associated. Compared with clinical characteristics alone, including the top five endogenous metabolites in serum and urine associated with longitudinal outcomes improved the outcome prediction (area under the receiver operating characteristic curves for eGFR decline: clinical model=0.79, clinical+metabolites model=0.87, P = 8.1E-6; for UACR doubling: clinical model=0.66, clinical+metabolites model=0.73, P = 2.9E-5).

Conclusions

Metabolomic profiling in different biofluids provided distinct and potentially complementary insights into the biology and prognosis of kidney diseases.

Cigarette Smoking and Risk of Hospitalization With Acute Kidney Injury: The Atherosclerosis Risk in Communities (ARIC) Study

RATIONALE & OBJECTIVE

Smoking is a modifiable risk factor for various adverse events. However, little is known about the association of smoking with the incidence of acute kidney injury (AKI) in the general population. This study investigated the association of cigarette smoking with the risk of AKI.

STUDY DESIGN

Prospective observational study.

SETTING & PARTICIPANTS

14,571 participants (mean age 55±6 years, 55% women, and 25% Black participants) from the ARIC study visit 1 (1987-1989) followed through December 31, 2019.

EXPOSURE

Smoking parameters (status, duration, pack-years, intensity, and years since cessation).

OUTCOME

Incident hospitalization with AKI, defined by a hospital discharge with a diagnostic code relevant to AKI.

ANALYTICAL APPROACH

Multivariable Cox regression models.

RESULTS

Over a median follow-up period of 26.3 years, 2,984 participants had an incident hospitalization with AKI. Current and former smokers had a significantly higher risk of AKI compared to never smokers after adjusting for potential confounders (HR, 2.22 [95% CI, 2.02-2.45] and 1.12 [1.02-1.23], respectively). A dose-response association was consistently seen for each of smoking duration, pack-years, and intensity with AKI (eg, HR, 1.19 [95% CI, 1.16-1.22] per 10 years of smoking). When years since cessation were considered as a time-varying exposure, the risk of AKI associated with smoking compared with current smokers began to decrease after 10 years, and became similar to never smokers at 30 years (HR for≥30 years, 1.07 [95% CI, 0.97-1.20] vs never smokers).

LIMITATIONS

Self-reported smoking measurements and missing outpatient AKI cases.

CONCLUSIONS

In a community-based cohort, all smoking parameters were robustly associated with the risk of AKI. Smoking cessation was associated with decreased risk of AKI, although the excess risk lasted up to 30 years. Our study supports the importance of preventing smoking initiation and promoting smoking cessation for the risk of AKI.

PLAIN-LANGUAGE SUMMARY

Smoking is a behavior that is associated with many negative health effects. It is not well understood how smoking relates to the occurrence of acute kidney injury (AKI) in the community. In this study, we looked at data from a group of 14,571 adults who were followed for 26 years to see how different aspects of smoking (such as whether someone smoked, how long they smoked for, how many cigarettes they smoked per day, and whether they quit smoking) were related to AKI. We found that smoking was strongly linked to an increased risk of AKI. This risk decreased after 5-10 years of quitting smoking, but the excess risk lasted up to 30 years. This study shows the importance of preventing people from starting smoking and to encourage smokers to quit to reduce their risk of AKI.

Pneumococcal vaccination effectiveness (PCV13 and PPSV23) in individuals with and without reduced kidney function: a test-negative design study

Background

Streptococcus pneumoniae vaccination effectiveness (VE) in individuals with reduced kidney function is unknown. We estimated pneumococcal conjugate vaccine (PCV13), pneumococcal polysaccharide vaccine (PPSV23), and combined PCV13 and PPSV23 effectiveness against pneumococcal disease in individuals with and without reduced estimated glomerular filtration rate (eGFR).

Methods

All eligible individuals (case and controls) were adults (aged ≥18 years) hospitalized within the Geisinger Health System and required to have S. pneumoniae urinary antigen testing (i.e. test-negative design). Vaccination records were obtained from the electronic health record and statewide vaccination registry. After controlling for the probability of receiving a pneumococcal vaccine, we used multivariable logistic regression models to estimate the odds ratios (ORs) of vaccination between those who did and did not meet the S. pneumoniae case definition. VE was calculated as (1 - OR) × 100%.

Results

There were 180 cases and 3889 controls (mean age 69 years, female 48%, white 97%, mean eGFR 71 mL/min/1.73 m2). The adjusted population PCV13 VE was 39% (95% CI 13%-58%), and combination PCV13 and PPSV23 was 39% (95% CI 12%-58%). PPSV23 VE was -3.7% (95% CI -57% to 32%). Stratified by eGFR, adjusted PCV13 VE was consistent in eGFR ≥60 [VE 38% (95% CI 2.9%-61%)] and 30-59 [VE 61% (95% CI 24%-80%)] without significant interaction. VE was not calculable for eGFR <30 due to small sample size.

Conclusion

PCV13 vaccination was associated with reduced risk of S. pneumoniae hospitalization in individuals with a reduced eGFR (30-59 mL/min/1.73 m2).

Consistency of metabolite associations with measured glomerular filtration rate in children and adults

Background

There is interest in identifying novel filtration markers that lead to more accurate GFR estimates than current markers (creatinine and cystatin C) and are more consistent across demographic groups. We hypothesize that large-scale metabolomics can identify serum metabolites that are strongly influenced by glomerular filtration rate (GFR) and are more consistent across demographic variables than creatinine, which would be promising filtration markers for future investigation.

Methods

We evaluated the consistency of associations between measured GFR (mGFR) and 887 common, known metabolites quantified by an untargeted chromatography- and spectroscopy-based metabolomics platform (Metabolon) performed on frozen blood samples from 580 participants in Chronic Kidney Disease in Children (CKiD), 674 participants in Modification of Diet in Renal Disease (MDRD) Study and 962 participants in African American Study of Kidney Disease and Hypertension (AASK). We evaluated metabolite-mGFR correlation association with metabolite class, molecular weight, assay platform and measurement coefficient of variation (CV). Among metabolites with strong negative correlations with mGFR (r < -0.5), we assessed additional variation by age (height in children), sex, race and body mass index (BMI).

Results

A total of 561 metabolites (63%) were negatively correlated with mGFR. Correlations with mGFR were highly consistent across study, sex, race and BMI categories (correlation of metabolite-mGFR correlations between 0.88 and 0.95). Amino acids, carbohydrates and nucleotides were more often negatively correlated with mGFR compared with lipids, but there was no association with metabolite molecular weight, liquid chromatography/mass spectrometry platform and measurement CV. Among 114 metabolites with strong negative associations with mGFR (r < -0.5), 27 were consistently not associated with age (height in children), sex or race.

Conclusions

The majority of metabolite-mGFR correlations were negative and consistent across sex, race, BMI and study. Metabolites with consistent strong negative correlations with mGFR and non-association with demographic variables may represent candidate markers to improve estimation of GFR.

Associations of Urine and Plasma Metabolites with Kidney Failure and Death in a CKD Cohort

RATIONALE & OBJECTIVE

Biomarkers that enable better identification of persons with chronic kidney disease (CKD) who are at higher risk for disease progression and adverse events are needed. This study sought to identify urine and plasma metabolites associated with progression of kidney disease.

STUDY DESIGN

Prospective metabolome-wide association study.

SETTING & PARTICIPANTS

Persons with CKD enrolled in the German CKD Study (GCKD) with metabolite measurements; with external validation within the Atherosclerosis Risk in Communities Study.

EXPOSURES

1,513 urine and 1,416 plasma metabolites (Metabolon, Inc.) measured at study entry using untargeted mass spectrometry.

OUTCOMES

Main endpoints were kidney failure (KF), and a composite endpoint of KF, eGFR <15 mL/min/1.73m2, or 40% decline in eGFR (CKE). Death from any cause was a secondary endpoint. After a median of 6.5 years follow-up, 500 persons experienced KF, 1,083 experienced CKE and 680 died.

ANALYTICAL APPROACH

Time-to-event analyses using multivariable proportional hazard regression models in a discovery-replication design, with external validation.

RESULTS

5,088 GCKD participants were included in analyses of urine metabolites and 5,144 in analyses of plasma metabolites. Among 182 unique metabolites, 30 were significantly associated with KF, 49 with CKE, and 163 with death. The strongest association with KF was observed for plasma hydroxyasparagine (hazard ratio: 1.95, 95% confidence interval: 1.68-2.25). An unnamed metabolite measured in plasma and urine was significantly associated with KF, CKE, and death. External validation of the identified associations of metabolites with KF or CKE revealed direction-consistency for 88% of observed associations. Selected associations of 18 metabolites with study outcomes have not been previously reported.

LIMITATIONS

Use of observational data and semi-quantitative metabolite measurements at a single time point.

CONCLUSIONS

The observed associations between metabolites and KF, CKE or death in persons with CKD confirmed previously reported findings and also revealed several associations not previously described. These findings warrant confirmatory research in other study cohorts.

Glucagon-like peptide-1 receptor agonists and the risk of atrial fibrillation in adults with diabetes: a real-world study

BACKGROUND

Glucagon-like peptide-1 receptor agonists (GLP-1RA) have cardiovascular benefits in type 2 diabetes, but none of the cardiovascular trials studied atrial fibrillation/atrial flutter (AF) as a primary endpoint. Data from post-marketing surveillance studies remains sparse.

OBJECTIVE

To examine the real-world risk of AF comparing GLP-1RA with other non-insulin glucose-lowering agents.

DESIGN

Cohort study using de-identified electronic health record data from the Optum Labs Data Warehouse.

PARTICIPANTS

Adult patients with diabetes who were newly prescribed add-on non-insulin glucose-lowering agents and were on metformin between 2005-2020.

EXPOSURES

New users of GLP-1RA were separately compared with new users of dipeptidyl peptidase-4 inhibitors (DPP4i) and sodium-glucose cotransporter 2 inhibitors (SGLT2i), using 1:1 propensity score matching to adjust for differences in patient characteristics.

MAIN MEASURES

The primary outcome was incident AF, defined and captured by diagnosis code for AF. Incidence rate difference (IRD) and hazard ratio (HR) were estimated in the matched cohorts.

KEY RESULTS

In the matched cohort of 14,566 pairs of GLP-1RA and DPP4i followed for a median of 3.8 years, GLP-1RA use was associated with a lower risk of AF (IRD, -1.0; 95% CI, -1.8 to -0.2 per 1000 person-years; HR, 0.82; 95% CI, 0.70 to 0.96). In the matched cohort of 9,424 pairs of patients on GLP-1RA and SGLT2i with a median follow-up of 2.9 years, there was no difference in the risk for AF (IRD, 0.4; 95% CI -0.7 to 1.5 per 1000 person-years; HR, 1.12; 95% CI, 0.89 to 1.42).

CONCLUSIONS

In this real-word study, GLP-1RA was associated with a lower risk of AF compared with DPP4i, but no difference compared with SGLT2i, suggesting that cardiovascular benefits of GLP-1RA use may extend to prevention for AF in patients with diabetes. Our findings call for future randomized controlled trials to focus on the effects of GLP-1RA on AF prevention.

Dihydropyridine Calcium Channel Blockers and Kidney Outcomes

BACKGROUND

Early trials of dihydropyridine calcium channel blockers (DCCBs) suggest a detrimental effect on intraglomerular pressure and an association with albuminuria.

OBJECTIVE

We sought to evaluate the associations of DCCB initiation with albuminuria and kidney failure with replacement therapy (KFRT) and to determine whether renin-angiotensin system (RAS) blockade modified these associations.

DESIGN

We conducted a target trial emulation study using a new user, active comparator design and electronic health record data from Geisinger Health.

PARTICIPANTS

We included patients without severe albuminuria or KFRT who were initiated on a DCCB or thiazide (active comparator) between January 1, 2004, and December 31, 2019.

MAIN MEASURES

Using inverse probability of treatment weighting, we performed doubly robust Cox proportional hazards regression to estimate the association of DCCB initiation with incident severe albuminuria (urine albumin to creatinine ratio > 300 mg/g) and KFRT, overall and stratified by RAS blocker use.

KEY RESULTS

There were 11,747 and 26,758 eligible patients initiating a DCCB and thiazide, respectively, with a weighted baseline mean age of 60 years, systolic blood pressure of 143 mm Hg, and eGFR of 86 mL/min/1.73 m2, and with a mean follow-up of 8 years. Compared with thiazides, DCCBs were significantly associated with the development of severe albuminuria (hazard ratio [HR], 1.29; 95% confidence interval [CI], 1.16-1.43), with attenuation of risk in the presence of RAS blockade (P for interaction < 0.001). The risk of KFRT was increased among patients without RAS blockade (HR, 1.66; 95% CI, 1.19-2.31), but not with RAS blockade (P for interaction = 0.005).

CONCLUSIONS

DCCBs were associated with increased risk of albuminuria and, in the absence of RAS blockade, KFRT. These findings suggest coupling DCCB therapy with RAS blockade may mitigate adverse kidney outcomes.

Serum Metabolomic Markers of Protein-Rich Foods and Incident CKD: Results From the Atherosclerosis Risk in Communities Study

Rationale & Objective

While urine excretion of nitrogen estimates the total protein intake, biomarkers of specific dietary protein sources have been sparsely studied. Using untargeted metabolomics, this study aimed to identify serum metabolomic markers of 6 protein-rich foods and to examine whether dietary protein-related metabolites are associated with incident chronic kidney disease (CKD).

Study Design

Prospective cohort study.

Setting & Participants

A total of 3,726 participants from the Atherosclerosis Risk in Communities study without CKD at baseline.

Exposures

Dietary intake of 6 protein-rich foods (fish, nuts, legumes, red and processed meat, eggs, and poultry), serum metabolites.

Outcomes

Incident CKD (estimated glomerular filtration rate < 60 mL/min/1.73 m2 with ≥25% estimated glomerular filtration rate decline relative to visit 1, hospitalization or death related to CKD, or end-stage kidney disease).

Analytical Approach

Multivariable linear regression models estimated cross-sectional associations between protein-rich foods and serum metabolites. C statistics assessed the ability of the metabolites to improve the discrimination of highest versus lower 3 quartiles of intake of protein-rich foods beyond covariates (demographics, clinical factors, health behaviors, and the intake of nonprotein food groups). Cox regression models identified prospective associations between protein-related metabolites and incident CKD.

Results

Thirty significant associations were identified between protein-rich foods and serum metabolites (fish, n = 8; nuts, n = 5; legumes, n = 0; red and processed meat, n = 5; eggs, n = 3; and poultry, n = 9). Metabolites collectively and significantly improved the discrimination of high intake of protein-rich foods compared with covariates alone (difference in C statistics = 0.033, 0.051, 0.003, 0.024, and 0.025 for fish, nuts, red and processed meat, eggs, and poultry-related metabolites, respectively; P < 1.00 × 10-16 for all). Dietary intake of fish was positively associated with 1-docosahexaenoylglycerophosphocholine (22:6n3), which was inversely associated with incident CKD (HR, 0.82; 95% CI, 0.75-0.89; P = 7.81 × 10-6).

Limitations

Residual confounding and sample-storage duration.

Conclusions

We identified candidate biomarkers of fish, nuts, red and processed meat, eggs, and poultry. A fish-related metabolite, 1-docosahexaenoylglycerophosphocholine (22:6n3), was associated with a lower risk of CKD.

Association of Sevelamer Initiation with Gastrointestinal Bleeding Hospitalization in Individuals Requiring Hemodialysis

Introduction Case reports have suggested a causative role between sevelamer use and subsequent gastrointestinal bleeding (GIB), but no large observational studies have evaluated this association. Methods Using the United States Renal Data System database from 2015 to 2019, we examined the association between initiation of sevelamer (versus non-sevelamer containing phosphate binders) and GIB hospitalization as well as all-cause mortality among individuals on hemodialysis. We emulated a target trial using Cox regression models and inverse probability of treatment weights to estimate the adjusted hazard ratios (HR) across outcomes and subgroups. Results Among 21,354 new users of phosphate binders (11,276 sevelamer and 10,078 non-sevelamer) with baseline lab data (calcium, phosphorus, hemoglobin, and albumin), there were 2,811 GIB hospitalizations and 5,920 deaths after a median follow-up of 1.3 years. Compared with the initiation of non-sevelamer binders, sevelamer was not associated with an increased risk of GIB hospitalization (89 vs. 90 events per 1000 person-years; IPTW-HR 0.98, 95% CI 0.91 - 1.06) or all-cause mortality (220 vs. 224 events per 1000 person-years; IPTW-HR 0.98 95% CI 0.93 - 1.03). Subgroup analyses (such as diabetes and anti-coagulation use) were generally consistent, and there was no association between sevelamer dose and GIB hospitalization. Conclusion Among patients requiring hemodialysis, sevelamer (vs non-sevelamer) containing phosphate binders was not associated with increased risk of GIB hospitalization.

Accuracy of GFR estimating equations based on creatinine, cystatin C or both in routine care

BACKGROUND

New equations to estimate glomerular filtration rate based on creatinine (eGFRcr), cystatin C (eGFRcys) or both (eGFRcr-cys) have been developed by the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) and the European Kidney Function Consortium (EKFC). There is a need to evaluate the performance of these equations in diverse European settings to inform implementation decisions, especially among people with key comorbid conditions.

METHODS

We performed a cross-sectional study including 6174 adults referred for single-point plasma clearance of iohexol in Stockholm, Sweden, with 9579 concurrent measurements of creatinine and cystatin C. We assessed the performance of the CKD-EPI 2009/2012/2021, EKFC 2021/2023, revised Lund-Malmö (RLM) 2011 and Caucasian, Asian, Pediatric and Adult (CAPA) 2014 equations against measured GFR (mGFR).

RESULTS

Mean age was 56 years, median mGFR was 62 mL/min/1.73 m2 and 40% were female. Comorbid conditions were common: cardiovascular disease (30%), liver disease (28%), diabetes (26%) and cancer (26%). All eGFRcr-cys equations had small bias and P30 (the percentage of estimated values within 30% of mGFR) close to 90%, and performed better than eGFRcr or eGFRcys equations. Among eGFRcr equations, CKD-EPI 2009 and CKD-EPI 2021 showed larger bias and lower P30 than EKFC 2021 and RLM. There were no meaningful differences in performance across eGFRcys equations. Findings were consistent across comorbid conditions, and eGFRcr-cys equations showed good performance in patients with liver disease, cancer and heart failure.

CONCLUSIONS

In conclusion, eGFRcr-cys equations performed best, with minimal variation among equations in this Swedish cohort. The lower performance of CKD-EPI eGFRcr equations compared with EKFC and RLM may reflect differences in population characteristics and mGFR methods. Implementing eGFRcr equations will require a trade-off between accuracy and uniformity across regions.

Association of Low Glomerular Filtration Rate With Adverse Outcomes at Older Age in a Large Population With Routinely Measured Cystatin C

BACKGROUND

The commonly accepted threshold of glomerular filtration rate (GFR) to define chronic kidney disease (CKD) is less than 60 mL/min/1.73 m2. This threshold is based partly on associations between estimated GFR (eGFR) and the frequency of adverse outcomes. The association is weaker in older adults, which has created disagreement about the appropriateness of the threshold for these persons. In addition, the studies measuring these associations included relatively few outcomes and estimated GFR on the basis of creatinine level (eGFRcr), which may be less accurate in older adults.

OBJECTIVE

To evaluate associations in older adults between eGFRcr versus eGFR based on creatinine and cystatin C levels (eGFRcr-cys) and 8 outcomes.

DESIGN

Population-based cohort study.

SETTING

Stockholm, Sweden, 2010 to 2019.

PARTICIPANTS

82 154 participants aged 65 years or older with outpatient creatinine and cystatin C testing.

MEASUREMENTS

Hazard ratios for all-cause mortality, cardiovascular mortality, and kidney failure with replacement therapy (KFRT); incidence rate ratios for recurrent hospitalizations, infection, myocardial infarction or stroke, heart failure, and acute kidney injury.

RESULTS

The associations between eGFRcr-cys and outcomes were monotonic, but most associations for eGFRcr were U-shaped. In addition, eGFRcr-cys was more strongly associated with outcomes than eGFRcr. For example, the adjusted hazard ratios for 60 versus 80 mL/min/1.73 m2 for all-cause mortality were 1.2 (95% CI, 1.1 to 1.3) for eGFRcr-cys and 1.0 (CI, 0.9 to 1.0) for eGFRcr, and for KFRT they were 2.6 (CI, 1.2 to 5.8) and 1.4 (CI, 0.7 to 2.8), respectively. Similar findings were observed in subgroups, including those with a urinary albumin-creatinine ratio below 30 mg/g.

LIMITATION

No GFR measurements.

CONCLUSION

Compared with low eGFRcr in older patients, low eGFRcr-cys was more strongly associated with adverse outcomes and the associations were more uniform.

PRIMARY FUNDING SOURCE

Swedish Research Council, National Institutes of Health, and Dutch Kidney Foundation.

Association of eGFR and Albuminuria with Venous Thromboembolism

BACKGROUND

CKD has been implicated as a risk factor of venous thromboembolism, but the evidence is limited to relatively healthy populations. The objective of this study was to discern whether parameters of kidney function and damage are associated with the occurrence of venous thromboembolism after hospitalization.

METHODS

We conducted a retrospective study including 23,899 and 11,552 adult individuals hospitalized within Geisinger Health System and New York University (NYU) Langone Health from 2004 to 2019 and 2012 to 2022, respectively. A Poisson model was used to evaluate adjusted incidence rates of venous thromboembolism according to eGFR and albuminuria categories in each cohort. Cox proportional hazards models were used to analyze associations of eGFR and urinary albumin-to-creatinine ratio (UACR) with venous thromboembolism, and hazard ratios (HRs) were meta-analyzed across cohorts.

RESULTS

Both lower eGFR and higher UACR were associated with higher risks of venous thromboembolism. In the Geisinger cohort, the incidence of venous thromboembolism after hospital discharge ranged from 10.7 (95% confidence interval [CI], 9.2 to 12.6) events per 1000 person-years in individuals in G1A1 (eGFR >90 ml/min per 1.73 m 2 and UACR <30 mg/g) to 27.7 (95% CI, 20.6 to 37.2) events per 1000 person-years in individuals with G4-5A3 (eGFR <30 ml/min per 1.73 m 2 and UACR >300 mg/g). A similar pattern was observed in the NYU cohort. Meta-analyses of the two cohorts showed that every 10 ml/min per 1.73 m 2 reduction in eGFR below 60 ml/min per 1.73 m 2 was associated with a 6% higher risk of venous thromboembolism (HR 1.06 [95% CI, 1.02 to 1.11], P = 0.01), and each two-fold higher UACR was associated with a 5% higher risk of venous thromboembolism (HR 1.05 [95% CI, 1.03 to 1.07], P < 0.001).

CONCLUSIONS

Both eGFR and UACR were independently associated with higher risk of venous thromboembolism after hospitalization. The incidence rate was higher with greater severity of CKD.

PODCAST

This article contains a podcast at https://dts.podtrac.com/redirect.mp3/www.asn-online.org/media/podcast/CJASN/2023_12_14_CJN0000000000000352.mp3.

Development and Validation of the American Heart Association's PREVENT Equations

BACKGROUND

Multivariable equations are recommended by primary prevention guidelines to assess absolute risk of cardiovascular disease (CVD). However, current equations have several limitations. Therefore, we developed and validated the American Heart Association Predicting Risk of CVD EVENTs (PREVENT) equations among US adults 30 to 79 years of age without known CVD.

METHODS

The derivation sample included individual-level participant data from 25 data sets (N=3 281 919) between 1992 and 2017. The primary outcome was CVD (atherosclerotic CVD and heart failure). Predictors included traditional risk factors (smoking status, systolic blood pressure, cholesterol, antihypertensive or statin use, and diabetes) and estimated glomerular filtration rate. Models were sex-specific, race-free, developed on the age scale, and adjusted for competing risk of non-CVD death. Analyses were conducted in each data set and meta-analyzed. Discrimination was assessed using the Harrell C-statistic. Calibration was calculated as the slope of the observed versus predicted risk by decile. Additional equations to predict each CVD subtype (atherosclerotic CVD and heart failure) and include optional predictors (urine albumin-to-creatinine ratio and hemoglobin A1c), and social deprivation index were also developed. External validation was performed in 3 330 085 participants from 21 additional data sets.

RESULTS

Among 6 612 004 adults included, mean±SD age was 53±12 years, and 56% were women. Over a mean±SD follow-up of 4.8±3.1 years, there were 211 515 incident total CVD events. The median C-statistics in external validation for CVD were 0.794 (interquartile interval, 0.763-0.809) in female and 0.757 (0.727-0.778) in male participants. The calibration slopes were 1.03 (interquartile interval, 0.81-1.16) and 0.94 (0.81-1.13) among female and male participants, respectively. Similar estimates for discrimination and calibration were observed for atherosclerotic CVD- and heart failure-specific models. The improvement in discrimination was small but statistically significant when urine albumin-to-creatinine ratio, hemoglobin A1c, and social deprivation index were added together to the base model to total CVD (ΔC-statistic [interquartile interval] 0.004 [0.004-0.005] and 0.005 [0.004-0.007] among female and male participants, respectively). Calibration improved significantly when the urine albumin-to-creatinine ratio was added to the base model among those with marked albuminuria (>300 mg/g; 1.05 [0.84-1.20] versus 1.39 [1.14-1.65]; P=0.01).

CONCLUSIONS

PREVENT equations accurately and precisely predicted risk for incident CVD and CVD subtypes in a large, diverse, and contemporary sample of US adults by using routinely available clinical variables.

Associations of Baseline and Longitudinal Serum Uromodulin With Kidney Failure and Mortality: Results From the African American Study of Kidney Disease and Hypertension (AASK) Trial

RATIONALE & OBJECTIVE

Uromodulin (UMOD) is the most abundant protein found in urine and has emerged as a promising biomarker of tubule health. Circulating UMOD is also detectable, but at lower levels. We evaluated whether serum UMOD levels were associated with the risks of incident kidney failure with replacement therapy (KFRT) and mortality.

STUDY DESIGN

Prospective cohort.

SETTING & PARTICIPANTS

Participants in AASK (the African American Study of Kidney Disease and Hypertension) with available stored serum samples from the 0-, 12-, and 24-month visits for biomarker measurement.

PREDICTORS

Baseline log-transformed UMOD and change in UMOD over 2 years.

OUTCOMES

KFRT and mortality.

ANALYTICAL APPROACH

Cox proportional hazards and mixed-effects models.

RESULTS

Among 500 participants with baseline serum UMOD levels (mean age, 54y; 37% female), 161 KFRT events occurred during a median of 8.5 years. After adjusting for baseline demographic factors, clinical factors, glomerular filtration rate, log-transformed urine protein-creatinine ratio, and randomized treatment groups, a 50% lower baseline UMOD level was independently associated with a 35% higher risk of KFRT (adjusted HR, 1.35; 95% CI, 1.07-1.70). For annual UMOD change, each 1-standard deviation lower change was associated with a 67% higher risk of KFRT (adjusted HR, 1.67; 95% CI, 1.41-1.99). Baseline UMOD and UMOD change were not associated with mortality. UMOD levels declined more steeply for metoprolol versus ramipril (P<0.001) as well as for intensive versus standard blood pressure goals (P = 0.002).

LIMITATIONS

Small sample size and limited generalizability.

CONCLUSIONS

Lower UMOD levels at baseline and steeper declines in UMOD over time were associated with a higher risk of subsequent KFRT in a cohort of African American adults with chronic kidney disease and hypertension.

PLAIN-LANGUAGE SUMMARY

Prior studies of uromodulin (UMOD), the most abundant protein in urine, and kidney disease have focused primarily on urinary UMOD levels. The present study evaluated associations of serum UMOD levels with the risks of kidney failure with replacement therapy (KFRT) and mortality in a cohort of African American adults with hypertension and chronic kidney disease. It found that participants with lower levels of UMOD at baseline were more likely to experience KFRT even after accounting for baseline kidney measures. Similarly, participants who experienced steeper annual declines in UMOD also had a heightened risk of kidney failure. Neither baseline nor annual change in UMOD was associated with mortality. Serum UMOD is a promising biomarker of kidney health.

Eicosanoids and Related Metabolites Associated with ESKD in a Community-Based Cohort

Key Points

High-throughput eicosanoid profiling can identify metabolites that may play a protective role in the development of kidney disease. In contrast to many other nonlipid metabolites, eicosanoid levels are minimally related with kidney filtration cross-sectionally.

Background

Eicosanoids are derivatives of polyunsaturated fatty acids and participate in the inflammatory response and the maintenance of endothelial function. Specific eicosanoids have been linked to various diseases, including hypertension and asthma, and may also reduce renal blood flow. A systematic investigation of eicosanoid-related metabolites and adverse kidney outcomes could identify key mediators of kidney disease and inform ongoing work in drug development.

Methods

Profiling of eicosanoid-related metabolites was performed in 9650 participants in the Atherosclerosis Risk in Communities Study (visit 2; mean age, 57 years). The associations between metabolite levels and the development of ESKD was investigated using Cox proportional hazards regression (n =256 events; median follow-up, 25.5 years). Metabolites with statistically significant associations with ESKD were evaluated for a potential causal role using bidirectional Mendelian randomization techniques, linking genetic instruments for eicosanoid levels to genomewide association study summary statistics of eGFR.

Results

The 223 eicosanoid-related metabolites that were profiled and passed quality control (QC) were generally uncorrelated with eGFR in cross-sectional analyses (median Spearman correlation, −0.03; IQR, −0.05 to 0.002). In models adjusted for multiple covariates, including baseline eGFR, three metabolites had statistically significant associations with ESKD (P value < 0.05/223). These included a hydroxyoctadecenoic acid, a dihydroxydocosapentaenoic acid, and arachidonic acid, with higher levels of the former two protective against ESKD and higher levels of arachidonic acid having a positive association with risk of ESKD. Mendelian randomization analyses suggested a causal role for the hydroxyoctadecenoic and arachidonic acid in determining eGFR. Spectral analysis identified the former metabolite as either 11-hydroxy-9-octadecenoic acid or 10-hydroxy-11-octadecenoic acid.

Conclusions

High-throughput eicosanoid profiling can identify metabolites that may play a protective role in the development of kidney disease.