Genomic structural equation modeling reveals latent phenotypes in the human cortex with distinct genetic architecture.

Genetic contributions to human cortical structure manifest pervasive pleiotropy. This pleiotropy may be harnessed to identify unique genetically-informed parcellations of the cortex that are neurobiologically distinct from functional, cytoarchitectural, or other cortical parcellation schemes. We investigated genetic pleiotropy by applying genomic structural equation modeling (SEM) to map the genetic architecture of cortical surface area (SA) and cortical thickness (CT) for 34 brain regions recently reported in the ENIGMA cortical GWAS. Genomic SEM uses the empirical genetic covariance estimated from GWAS summary statistics with LD score regression (LDSC) to discover factors underlying genetic covariance, which we are denoting genetically informed brain networks (GIBNs). Genomic SEM can fit a multivariate GWAS from summary statistics for each of the GIBNs, which can subsequently be used for LD score regression (LDSC). We found the best-fitting model of cortical SA identified 6 GIBNs and CT identified 4 GIBNs, although sensitivity analyses indicated that other structures were plausible. The multivariate GWASs of the GIBNs identified 74 genome-wide significant (GWS) loci (p < 5 × 10-8), including many previously implicated in neuroimaging phenotypes, behavioral traits, and psychiatric conditions. LDSC of GIBN GWASs found that SA-derived GIBNs had a positive genetic correlation with bipolar disorder (BPD), and cannabis use disorder, indicating genetic predisposition to a larger SA in the specific GIBN is associated with greater genetic risk of these disorders. A negative genetic correlation was observed between attention deficit hyperactivity disorder (ADHD) and major depressive disorder (MDD). CT GIBNs displayed a negative genetic correlation with alcohol dependence. Even though we observed model instability in our application of genomic SEM to high-dimensional data, jointly modeling the genetic architecture of complex traits and investigating multivariate genetic links across neuroimaging phenotypes offers new insights into the genetics of cortical structure and relationships to psychopathology.

Genomic Structural Equation Modeling Reveals Latent Phenotypes in the Human Cortex with Distinct Genetic Architecture.

Genetic contributions to human cortical structure manifest pervasive pleiotropy. This pleiotropy may be harnessed to identify unique genetically-informed parcellations of the cortex that are neurobiologically distinct from functional, cytoarchitectural, or other cortical parcellation schemes. We investigated genetic pleiotropy by applying genomic structural equation modeling (SEM) to map the genetic architecture of cortical surface area (SA) and cortical thickness (CT) for the 34 brain regions recently reported in the ENIGMA cortical GWAS. Genomic SEM uses the empirical genetic covariance estimated from GWAS summary statistics with LD score regression (LDSC) to discover factors underlying genetic covariance, which we are denoting genetically informed brain networks (GIBNs). Genomic SEM can fit a multivariate GWAS from summary statistics for each of the GIBNs, which can subsequently be used for LD score regression (LDSC). We found the best-fitting model of cortical SA identified 6 GIBNs and CT identified 4 GIBNs. The multivariate GWASs of these GIBNs identified 74 genome-wide significant (GWS) loci (p<5×10-8), including many previously implicated in neuroimaging phenotypes, behavioral traits, and psychiatric conditions. LDSC of GIBN GWASs found that SA-derived GIBNs had a positive genetic correlation with bipolar disorder (BPD), and cannabis use disorder, indicating genetic predisposition to a larger SA in the specific GIBN is associated with greater genetic risk of these disorders. A negative genetic correlation was observed with attention deficit hyperactivity disorder (ADHD), major depressive disorder (MDD), and insomnia, indicating genetic predisposition to a larger SA in the specific GIBN is associated with lower genetic risk of these disorders. CT GIBNs displayed a negative genetic correlation with alcohol dependence. Jointly modeling the genetic architecture of complex traits and investigating multivariate genetic links across phenotypes offers a new vantage point for mapping the cortex into genetically informed networks.

A novel principal component based method for identifying differentially methylated regions in Illumina Infinium MethylationEPIC BeadChip data.

Differentially methylated regions (DMRs) are genomic regions with methylation patterns across multiple CpG sites that are associated with a phenotype. In this study, we proposed a Principal Component (PC) based DMR analysis method for use with data generated using the Illumina Infinium MethylationEPIC BeadChip (EPIC) array. We obtained methylation residuals by regressing the M-values of CpGs within a region on covariates, extracted PCs of the residuals, and then combined association information across PCs to obtain regional significance. Simulation-based genome-wide false positive (GFP) rates and true positive rates were estimated under a variety of conditions before determining the final version of our method, which we have named DMRPC. Then, DMRPC and another DMR method, coMethDMR, were used to perform epigenome-wide analyses of several phenotypes known to have multiple associated methylation loci (age, sex, and smoking) in a discovery and a replication cohort. Among regions that were analysed by both methods, DMRPC identified 50% more genome-wide significant age-associated DMRs than coMethDMR. The replication rate for the loci that were identified by only DMRPC was higher than the rate for those that were identified by only coMethDMR (90% for DMRPC vs. 76% for coMethDMR). Furthermore, DMRPC identified replicable associations in regions of moderate between-CpG correlation which are typically not analysed by coMethDMR. For the analyses of sex and smoking, the advantage of DMRPC was less clear. In conclusion, DMRPC is a new powerful DMR discovery tool that retains power in genomic regions with moderate correlation across CpGs.

An evaluation of the genome-wide false positive rates of common methods for identifying differentially methylated regions using illumina methylation arrays.

Differentially methylated regions (DMRs) are genomic regions with specific methylation patterns across multiple loci that are associated with a phenotype. We examined the genome-wide false positive (GFP) rates of five widely used DMR methods: comb-p, Bumphunter, DMRcate, mCSEA and coMethDMR using both Illumina HumanMethylation450 (450 K) and MethylationEPIC (EPIC) data and simulated continuous and dichotomous null phenotypes (i.e., generated independently of methylation data). coMethDMR provided well-controlled GFP rates (~5%) except when analysing skewed continuous phenotypes. DMRcate generally had well-controlled GFP rates when applied to 450 K data except for the skewed continuous phenotype and EPIC data only for the normally distributed continuous phenotype. GFP rates for mCSEA were at least 0.096 and comb-p yielded GFP rates above 0.34. Bumphunter had high GFP rates of at least 0.35 across conditions, reaching as high as 0.95. Analysis of the performance of these methods in specific regions of the genome found that regions with higher correlation across loci had higher regional false positive rates on average across methods. Based on the false positive rates, coMethDMR is the most recommended analysis method, and DMRcate had acceptable performance when analysing 450 K data. However, as both could display higher levels of FPs for skewed continuous distributions, a normalizing transformation of skewed continuous phenotypes is suggested. This study highlights the importance of genome-wide simulations when evaluating the performance of DMR-analysis methods.

Trauma and posttraumatic stress disorder modulate polygenic predictors of hippocampal and amygdala volume.

The volume of subcortical structures represents a reliable, quantitative, and objective phenotype that captures genetic effects, environmental effects such as trauma, and disease effects such as posttraumatic stress disorder (PTSD). Trauma and PTSD represent potent exposures that may interact with genetic markers to influence brain structure and function. Genetic variants, associated with subcortical volumes in two large normative discovery samples, were used to compute polygenic scores (PGS) for the volume of seven subcortical structures. These were applied to a target sample enriched for childhood trauma and PTSD. Subcortical volume PGS from the discovery sample were strongly associated in our trauma/PTSD enriched sample (n = 7580) with respective subcortical volumes of the hippocampus (p = 1.10 × 10-20), thalamus (p = 7.46 × 10-10), caudate (p = 1.97 × 10-18), putamen (p = 1.7 × 10-12), and nucleus accumbens (p = 1.99 × 10-7). We found a significant association between the hippocampal volume PGS and hippocampal volume in control subjects from our sample, but was absent in individuals with PTSD (GxE; (beta = -0.10, p = 0.027)). This significant GxE (PGS × PTSD) relationship persisted (p < 1 × 10-19) in four out of five threshold peaks (0.024, 0.133, 0.487, 0.730, and 0.889) used to calculate hippocampal volume PGSs. We detected similar GxE (G × ChildTrauma) relationships in the amygdala for exposure to childhood trauma (rs4702973; p = 2.16 × 10-7) or PTSD (rs10861272; p = 1.78 × 10-6) in the CHST11 gene. The hippocampus and amygdala are pivotal brain structures in mediating PTSD symptomatology. Trauma exposure and PTSD modulate the effect of polygenic markers on hippocampal volume (GxE) and the amygdala volume PGS is associated with PTSD risk, which supports the role of amygdala volume as a risk factor for PTSD.

Dynamic monitoring of circulating tumor DNA to predict prognosis and efficacy of adjuvant chemotherapy after resection of colorectal liver metastases.

Rationale: Hepatectomy and adjuvant chemotherapy after resection of colorectal liver metastases (CRLM) may improve survival, however, patients which may benefit cannot currently be identified. Postoperative circulating tumor DNA (ctDNA) analysis can detect minimal residual disease (MRD) and predict the prognosis and efficacy of adjuvant chemotherapy. Our study aims to determine the impact of serial ctDNA analysis to predict the outcome among patients undergoing resection of CRLM. Methods: Between May 2018 and October 2019, 91 CRLM patients were prospectively enrolled. Whole exome sequencing was performed in 50 primary and 48 metastatic liver tissues. Targeted sequencing of 451 cancer relevant genes was performed in 50 baseline plasma to determine plasma-tissue concordance. We prospectively investigated changes in the amount and constitution of ctDNA in 271 serial plasma samples taken at different time points (baseline, pre-operation, post-operation, post-operative adjuvant chemotherapy (post-ACT) and recurrence) during the treatment of CRLM. Results: Detected molecular alterations were highly consistent among baseline ctDNA, primary and liver metastases tissue. Patients with a higher variant allele frequency (VAF) level at baseline ctDNA represent a higher tumor burden, and decreased ctDNA during pre-operative chemotherapy predicted better tumor response. Patients with detectable post-operative and post-ACT ctDNA were associated with significantly shorter recurrence-free survival (RFS). ROC analysis showed that post-ACT ctDNA status was superior to post-operative ctDNA status in predicting RFS with an AUROC of 0.79. A significant difference in overall recurrence rate was observed in patients with detectable vs undetectable levels of ctDNA after resection of CRLM (79.4% vs 41.7%) and after completion of adjuvant chemotherapy (77.3% vs 40.7%). During adjuvant chemotherapy, patients with decreased ctDNA VAF after adjuvant chemotherapy had a recurrence rate of 63.6%, compared to 92.3% in patients with increased ctDNA VAF. Conclusions: We envision that dynamic ctDNA analysis, especially in a post-ACT setting, might be used to not only reflect MRD but also to determine rational personalized adjuvant therapy after the resection of CRLM.

PTSD is associated with increased DNA methylation across regions of HLA-DPB1 and SPATC1L.

Posttraumatic stress disorder (PTSD) is characterized by intrusive thoughts, avoidance, negative alterations in cognitions and mood, and arousal symptoms that adversely affect mental and physical health. Recent evidence links changes in DNA methylation of CpG cites to PTSD. Since clusters of proximal CpGs share similar methylation signatures, identification of PTSD-associated differentially methylated regions (DMRs) may elucidate the pathways defining differential risk and resilience of PTSD. Here we aimed to identify epigenetic differences associated with PTSD. DNA methylation data profiled from blood samples using the MethylationEPIC BeadChip were used to perform a DMR analysis in 187 PTSD cases and 367 trauma-exposed controls from the Grady Trauma Project (GTP). DMRs were assessed with R package bumphunter. We identified two regions that associate with PTSD after multiple test correction. These regions were in the gene body of HLA-DPB1 and in the promoter of SPATC1L. The DMR in HLA-DPB1 was associated with PTSD in an independent cohort. Both DMRs included CpGs whose methylation associated with nearby sequence variation (meQTL) and that associated with expression of their respective genes (eQTM). This study supports an emerging literature linking PTSD risk to genetic and epigenetic variation in the HLA region.

Glycemic Control and Infections Among US Hemodialysis Patients With Diabetes Mellitus.

INTRODUCTION: Patients with diabetes mellitus (DM) on hemodialysis (HD) may be particularly vulnerable to infections. METHODS: We used merged data from the United States Renal Data System and electronic health records data from a large US dialysis provider to retrospectively examine the association between glycemic control and infections in these patients. Adult patients with DM aged ≥18 years who initiated in-center maintenance HD treatment from 2006 to 2011 and survived >90 days were included. Quarterly mean time-averaged hemoglobin A1c (HbA1c) values were categorized into <5.5%, 5.5 to <6.5%, 6.5 to <7.5%, 7.5 to <8.5%, and ≥8.5%. We used Medicare claims to ascertain infection-related outcomes and the ESRD Death Notification to identify death from infectious cause. We used Cox proportional hazards models to estimate multivariable-adjusted hazard ratios and 95% confidence intervals (CIs) for the associations between time-averaged HbA1c categories and infectious events. RESULTS: In a cohort of 33,753 eligible patients, those with higher HbA1c levels had higher rates of diabetic foot infections and skin and soft tissue infections, with patients with HbA1c ≥8.5% having 23% (95% CI, 5%, 45%) and 22% (95% CI, 5%, 42%) higher rates, respectively, compared with HbA1c 5.5 to <6.5%. Patients in the lower HbA1c categories had higher rates of infection-related and all-cause mortality (P-for-trend <0.001). CONCLUSION: This study highlights the need for greater attention to foot evaluation and skin and soft tissue infections among patients on HD with less than optimal diabetes control.

An epigenome-wide association study of posttraumatic stress disorder in US veterans implicates several new DNA methylation loci.

BACKGROUND: Previous studies using candidate gene and genome-wide approaches have identified epigenetic changes in DNA methylation (DNAm) associated with posttraumatic stress disorder (PTSD). METHODS: In this study, we performed an EWAS of PTSD in a cohort of Veterans (n = 378 lifetime PTSD cases and 135 controls) from the Translational Research Center for TBI and Stress Disorders (TRACTS) cohort assessed using the Illumina EPIC Methylation BeadChip which assesses DNAm at more than 850,000 sites throughout the genome. Our model included covariates for ancestry, cell heterogeneity, sex, age, and a smoking score based on DNAm at 39 smoking-associated CpGs. We also examined in EPIC-based DNAm data generated from pre-frontal cortex (PFC) tissue from the National PTSD Brain Bank (n = 72). RESULTS: The analysis of blood samples yielded one genome-wide significant association with PTSD at cg19534438 in the gene G0S2 (p = 1.19 × 10-7, padj = 0.048). This association was replicated in an independent PGC-PTSD-EWAS consortium meta-analysis of military cohorts (p = 0.0024). We also observed association with the smoking-related locus cg05575921 in AHRR despite inclusion of a methylation-based smoking score covariate (p = 9.16 × 10-6), which replicates a previously observed PGC-PTSD-EWAS association (Smith et al. 2019), and yields evidence consistent with a smoking-independent effect. The top 100 EWAS loci were then examined in the PFC data. One of the blood-based PTSD loci, cg04130728 in CHST11, which was in the top 10 loci in blood, but which was not genome-wide significant, was significantly associated with PTSD in brain tissue (in blood p = 1.19 × 10-5, padj = 0.60, in brain, p = 0.00032 with the same direction of effect). Gene set enrichment analysis of the top 500 EWAS loci yielded several significant overlapping GO terms involved in pathogen response, including "Response to lipopolysaccharide" (p = 6.97 × 10-6, padj = 0.042). CONCLUSIONS: The cross replication observed in independent cohorts is evidence that DNA methylation in peripheral tissue can yield consistent and replicable PTSD associations, and our results also suggest that that some PTSD associations observed in peripheral tissue may mirror associations in the brain.

Screening Rates for Primary Aldosteronism in Resistant Hypertension: A Cohort Study.

Resistant hypertension is associated with higher rates of cardiovascular disease, kidney disease, and death than primary hypertension. Although clinical practice guidelines recommend screening for primary aldosteronism among persons with resistant hypertension, rates of screening are unknown. We identified 145 670 persons with hypertension and excluded persons with congestive heart failure or advanced chronic kidney disease. Among this cohort, we studied 4660 persons ages 18 to <90 from the years 2008 to 2014 with resistant hypertension and available laboratory tests within the following 24 months. The screening rate for primary aldosteronism in persons with resistant hypertension was 2.1%. Screened persons were younger (55.9±13.3 versus 65.5±11.6 years; P<0.0001) and had higher systolic (145.1±24.3 versus 139.6±20.5 mm Hg; P=0.04) and diastolic blood pressure (81.8±13.6 versus 74.4±13.8 mm Hg; P<0.0001), lower rates of coronary artery disease (5.2% versus 14.2%; P=0.01), and lower serum potassium concentrations (3.9±0.6 versus 4.1±0.5 mmol/L; P=0.04) than unscreened persons. Screened persons had significantly higher rates of prescription for calcium channel blockers, mixed α/ß-adrenergic receptor antagonists, sympatholytics, and vasodilators, and lower rates of prescription for loop, thiazide, and thiazide-type diuretics. The prescription of mineralocorticoid receptor antagonists or other potassium-sparing diuretics was not significantly different between groups (P=0.20). In conclusion, only 2.1% of eligible persons received a screening test within 2 years of meeting criteria for resistant hypertension. Low rates of screening were not due to the prescription of antihypertensive medications that may potentially interfere with interpretation of the screening test. Efforts to highlight guideline-recommended screening and targeted therapy are warranted.

Outcomes after left ventricular assist device implantation in patients with acute kidney injury.

OBJECTIVE: The study objective was to compare outcomes for patients with and without acute kidney injury during hospitalizations when left ventricular assist devices are implanted. METHODS: By using the National Inpatient Sample from 2008 to 2013, we identified patients with an International Classification of Diseases, Ninth Revision procedure code for left ventricular assist device implantation (37.66). We ascertained the presence of acute kidney injury and acute kidney injury requiring dialysis using validated International Classification of Diseases, Ninth Revision codes. We used logistic regression to examine the association of nondialysis-requiring acute kidney injury and acute kidney injury requiring dialysis with mortality, procedural complications, and discharge destination. RESULTS: We identified 8362 patients who underwent left ventricular assist device implantation, of whom 3760 (45.0%) experienced nondialysis-requiring acute kidney injury and 426 (5.1%) experienced acute kidney injury requiring dialysis. In-hospital mortality was 3.9% for patients without acute kidney injury, 12.2% for patients with nondialysis-requiring acute kidney injury, and 47.4% for patients with acute kidney injury requiring dialysis. Patients with nondialysis-requiring acute kidney injury and acute kidney injury requiring dialysis had higher adjusted odds of mortality (3.24, 95% confidence interval [CI], 2.04-5.13 and 20.8, 95% CI, 9.7-44.2), major bleeding (1.38, 95% CI, 1.08-1.77 and 2.44, 95% CI, 1.47-4.04), sepsis (2.69, 95% CI, 1.93-3.75 and 5.75, 95% CI, 3.46-9.56), and discharge to a nursing facility (2.15, 95% CI, 1.51-3.07 and 5.89, 95% CI, 2.67-12.99). CONCLUSIONS: More than 1 in 10 patients with acute kidney injury and approximately 1 in 2 patients with acute kidney injury requiring dialysis died during their hospitalization, with only 30% of patients with acute kidney injury requiring dialysis discharged to home. This information is necessary to support shared decision-making for patients with advanced heart failure and acute kidney injury.

Early detection of chronic kidney disease in low-income and middle-income countries: development and validation of a point-of-care screening strategy for India.

INTRODUCTION: Although deaths due to chronic kidney disease (CKD) have doubled over the past two decades, few data exist to inform screening strategies for early detection of CKD in low-income and middle-income countries. METHODS: Using data from three population-based surveys in India, we developed a prediction model to identify a target population that could benefit from further CKD testing, after an initial screening implemented during home health visits. Using data from one urban survey (n=8698), we applied stepwise logistic regression to test three models: one comprised of demographics, self-reported medical history, anthropometry and point-of-care (urine dipstick or capillary glucose) tests; one with demographics and self-reported medical history and one with anthropometry and point-of-care tests. The 'gold-standard' definition of CKD was an estimated glomerular filtration rate <60 mL/min/1.73 m2 or urine albumin-to-creatinine ratio ≥30 mg/g. Models were internally validated via bootstrap. The most parsimonious model with comparable performance was externally validated on distinct urban (n=5365) and rural (n=6173) Indian cohorts. RESULTS: A model with age, sex, waist circumference, body mass index and urine dipstick had a c-statistic of 0.76 (95% CI 0.75 to 0.78) for predicting need for further CKD testing, with external validation c-statistics of 0.74 and 0.70 in the urban and rural cohorts, respectively. At a probability cut-point of 0.09, sensitivity was 71% (95% CI 68% to 74%) and specificity was 70% (95% CI 69% to 71%). The model captured 71% of persons with CKD and 90% of persons at highest risk of complications from untreated CKD (ie, CKD stage 3A2 and above). CONCLUSION: A point-of-care CKD screening strategy using three simple measures can accurately identify high-risk persons who require confirmatory kidney function testing.

Longitudinal changes in kidney function following heart transplantation: Stanford experience.

Many heart transplant recipients experience declining kidney function following transplantation. We aimed to quantify change in kidney function in heart transplant recipients stratified by pre-transplant kidney function. A total of 230 adult heart transplant recipients between May 1, 2008, and December 31, 2014, were evaluated for up to 5 years post-transplant (median 1 year). Using 19 398 total estimated glomerular filtration rate (eGFR) assessments, we evaluated trends in eGFR in recipients with normal/near-normal (eGFR ≥45 mL/min/1.73 m2 ) vs impaired (eGFR <45 mL/min/1.73 m2 ) kidney function and the likelihood of reaching an eGFR of 20 mL/min/1.73 m2 after heart transplant. Baseline characteristics were similar. Immediately following heart transplant, the impaired pre-transplant kidney function group showed a mean eGFR gain of 9.5 mL/min/1.73 m2 (n = 193) vs a mean decline of 4.9 mL/min/1.73 m2 (n = 37) in the normal/near-normal group. Subsequent rates of eGFR decline were 2.2 mL/min/1.73 m2 /y vs 2.9 mL/min/1.73 m2 /y, respectively. The probability of reaching an eGFR of 20 mL/min/1.73 m2 or less at 1, 5, and 10 years following heart transplant was 1%, 4%, and 30% in the impaired group, and <1%, <1%, and 10% in the normal/near-normal group. Estimates of expected recovery in kidney function and its decline over time will help inform decision making about kidney care after heart transplantation.

The Relationship between Intradialytic Hypotension and Hospitalized Mesenteric Ischemia: A Case-Control Study.

BACKGROUND AND OBJECTIVES: Mesenteric ischemia is a rare but devastating condition caused by insufficient blood supply to meet the demands of intestinal metabolism. In patients with ESKD, it can be difficult to diagnose and has a >70% mortality rate. Patients on hemodialysis have a high prevalence of predisposing conditions for mesenteric ischemia, but the contribution of intradialytic hypotension, a potential modifiable risk factor, has not been well described. DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS: We used data from the US Renal Data System to identify 626 patients on hemodialysis with a hospitalized mesenteric ischemia event (cases). We selected 2428 controls in up to a 1:4 ratio matched by age, sex, black race, incident dialysis year, diabetes mellitus, coronary artery disease, and peripheral artery disease. We used six different definitions of intradialytic hypotension on the basis of prior studies, and categorized patients as having had intradialytic hypotension if ≥30% of hemodialysis sessions in the 30 days before the event met the specified definition. RESULTS: The proportion of patients with intradialytic hypotension varied depending on its definition: from 19% to 92% of cases and 11% to 94% of controls. Cases had a higher adjusted odds (1.82; 95% confidence interval, 1.47 to 2.26) of having had intradialytic hypotension in the preceding 30 days than controls when using nadir-based intradialytic hypotension definitions such as nadir systolic BP <90 mm Hg. To examine a potential dose-response association of intradialytic hypotension with hospitalized mesenteric ischemia, we categorized patients by the proportion of hemodialysis sessions having intradialytic hypotension, defined using the Nadir90 definition (0%, 1%-9%, 10%-29%, 30%-49%, and ≥50%), and found a direct association of proportion of intradialytic hypotension with hospitalized mesenteric ischemia (P-trend<0.001). CONCLUSIONS: Patients with hospitalized mesenteric ischemia had significantly higher odds of having had intradialytic hypotension in the preceding 30 days than controls, as defined by nadir-based definitions.

Trends in Rates of Lower Extremity Amputation Among Patients With End-stage Renal Disease Who Receive Dialysis.

Importance: Patients with end-stage renal disease (ESRD) who receive dialysis are at high risk of lower extremity amputation. Recent studies indicate decreasing rates of lower extremity amputation in non-ESRD populations, but contemporary data for patients with ESRD who receive dialysis are lacking. Objectives: To assess rates of lower extremity amputation among patients with ESRD who receive dialysis during a recent 15-year period; to analyze whether those rates differed by age, sex, diabetes, or geographic region; and to determine 1-year mortality rates in this population after lower extremity amputation. Design, Setting, and Participants: This retrospective study of 3 700 902 records obtained from a US national registry of patients with ESRD who receive dialysis assessed cross-sectional cohorts for each calendar year from 2000 through 2014. Adult patients with prevalent ESRD treated with hemodialysis or peritoneal dialysis covered by Medicare Part A and B on January 1 of each cohort year were included. Data analysis was conducted from August 2017 to April 2018. Exposures: Age, sex, diabetes, and hospital referral region. Main Outcomes and Measures: Annual rates per 100 person-years of nontraumatic major (above- or below-knee) and minor (below-ankle) amputations. Results: For each annual cohort, there were fewer women (47.5% in 2000, 46.2% in 2005, 44.9% in 2010, and 44.0% in 2014) than men, more than half the patients were white individuals (58.1% in 2000, 56.9% in 2005, 56.9% in 2010, and 56.7% in 2014), and a small proportion were employed (13.9% in 2000, 15.1% in 2005, 16.1% in 2010, and 16.5% in 2014). The rate of lower extremity amputations for patients with ESRD who receive dialysis decreased by 51.0% from 2000 to 2014, driven primarily by a decrease in the rate of major amputations (5.42 [95% CI, 5.28-5.56] in 2000 vs 2.66 [95% CI, 2.59-2.72] per 100 person-years in 2014). Patients with diabetes had amputation rates more than 5 times as high as patients without diabetes. Patients younger than 65 years had higher adjusted amputation rates than older patients, and men had consistently higher adjusted amputation rates than women. Adjusted 1-year mortality rates after lower extremity amputation for patients with ESRD who receive dialysis decreased from 52.2% (95% CI, 50.9%-53.4%) in 2000 to 43.6% (95% CI, 42.5%-44.8%) in 2013. In general, amputation rates decreased among all regions from 2000 to 2014, but regional variability persisted across time despite adjustment for differences in patient demographics and comorbid conditions. Conclusions and Relevance: Although rates of lower extremity amputations among US patients with ESRD who receive dialysis decreased by 51% during a recent 15-year period, mortality rates remained high, with nearly half of patients dying within a year after lower extremity amputation. Our results highlight the need for more research on ways to prevent lower extremity amputation in this extremely high-risk population.

Acute Kidney Injury Due to Diarrheal Illness Requiring Hospitalization: Data from the National Inpatient Sample.

BACKGROUND: Diarrheal illness is a major reason for hospitalization, but data on consequent acute kidney injury (AKI) are sparse. OBJECTIVE: To determine the incidence of AKI in infectious and non-infectious diarrheal illness requiring hospitalization and to identify correlates and outcomes of diarrhea-associated AKI. DESIGN: Using data from the 2012 National Inpatient Sample (NIS), we created a cohort of patients with a primary diagnosis of diarrheal illness. Diarrheal illness, disease correlates, and AKI were defined by ICD-9 diagnosis codes. We used logistic regression with backward variable selection to determine factors independently associated with AKI in infectious and non-infectious diarrheal illness, as well as to determine the association of AKI with in-hospital mortality. We used generalized linear models to assess differences in length of stay and costs of hospitalization. MAIN MEASURES: The primary outcome was AKI in hospitalized diarrheal illness. Secondary outcomes were in-hospital mortality, length of stay, and cost of hospitalization associated with AKI. KEY RESULTS: One in ten adults hospitalized with diarrheal illness experienced AKI, with higher incidence rates in older adults. Chronic kidney disease (CKD) and hypertension were associated with increased odds of AKI (all diarrhea OR 4.81, 95% CI 4.52 to 5.12 and OR 1.33, 95% CI 1.27 to 1.40, respectively). AKI in diarrheal illness was associated with substantial increase in mortality (OR 5.05, 95% CI 4.47 to 5.72), length of stay (mean increase 1.7 days [95% CI 1.6 to 1.8]), and cost of hospitalization (mean increase $4411 [95% CI 4023 to 4800]). CONCLUSION: Acute kidney injury is common and consequential among patients hospitalized for diarrheal illness. Persons with CKD and hypertension are the most susceptible, possibly due to diminished renal reserve and exacerbating effects of treatment with diuretics and renin-angiotensin-aldosterone system blockers. Proactive management of these unique pharmacologic and physiologic factors is necessary to prevent AKI in this vulnerable population.

Market Competition and Health Outcomes in Hemodialysis.

OBJECTIVE: To examine whether market competition is associated with improved health outcomes in hemodialysis. DATA SOURCES: Secondary analysis of data from a national dialysis registry between 2001 and 2011. STUDY DESIGN: We conducted one- and two-part linear regression models, using each hospital service area (HSA) as its own control, to examine the independent associations among market concentration and health outcomes. DATA COLLECTION: We selected cohorts of patients receiving in-center hemodialysis in the United States at the start of each calendar year. We used information about dialysis facility ownership and the location where patients received dialysis to measure an index of market concentration-the Hirschman-Herfindahl Index (HHI)-for HSA and year, which ranges from near zero (perfect competition) to one (monopoly). PRINCIPAL FINDINGS: An average reduction in HHI by 0.2 (one standard deviation in 2011) was associated with 2.9 fewer hospitalizations per 100 patient-years (95 percent CI, 0.4 to 5.4). If these findings were generalized to the entire in-center hemodialysis population, this would translate to 8,100 (95 percent CI 1,200 to 15,000) fewer hospitalizations in 2011. There was no association between change in market competition and mortality. CONCLUSIONS: Market competition in dialysis may lead to improved health outcomes.

Infrequent Provision of Palliative Care to Patients with Dialysis-Requiring AKI.

BACKGROUND AND OBJECTIVES: The use of palliative care in AKI is not well described. We sought to better understand palliative care practice patterns for hospitalized patients with AKI requiring dialysis in the United States. DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS: Using the 2012 National Inpatient Sample, we identified patients with AKI and palliative care encounters using validated International Classification of Diseases, Ninth Revision, Clinical Modification codes. We compared palliative care encounters in patients with AKI requiring dialysis, patients with AKI not requiring dialysis, and patients without AKI. We described the provision of palliative care in patients with AKI requiring dialysis and compared the frequency of palliative care encounters for patients with AKI requiring dialysis with that for patients with other illnesses with similarly poor prognoses. We used logistic regression to determine factors associated with the provision of palliative care, adjusting for demographics, hospital-level variables, and patient comorbidities. RESULTS: We identified 3,031,036 patients with AKI, of whom 91,850 (3%) received dialysis. We observed significant patient- and hospital-level differences in the provision of palliative care for patients with AKI requiring dialysis; adjusted odds were 26% (95% confidence interval, 12% to 38%) lower in blacks and 23% (95% confidence interval, 3% to 39%) lower in Hispanics relative to whites. Lower provision of palliative care was observed for rural and urban nonteaching hospitals relative to urban teaching hospitals, small and medium hospitals relative to large hospitals, and hospitals in the Northeast compared with the South. After adjusting for age and sex, there was low utilization of palliative care services for patients with AKI requiring dialysis (8%)-comparable with rates of utilization by patients with other illnesses with poor prognosis, including cardiogenic shock (9%), intracranial hemorrhage (10%), and acute respiratory distress syndrome (10%). CONCLUSIONS: The provision of palliative care varied widely by patient and facility characteristics. Palliative care was infrequently used in hospitalized patients with AKI requiring dialysis, despite its poor prognosis and the regular application of life-sustaining therapy.

Do attributes of persons with chronic kidney disease differ in low-income and middle-income countries compared with high-income countries? Evidence from population-based data in six countries.

Kidney biopsies to elucidate the cause of chronic kidney disease (CKD) are performed in a minority of persons with CKD living in high-income countries, since associated conditions-that is, diabetes mellitus, vascular disease or obesity with pre-diabetes, prehypertension or dyslipidaemia-can inform management targeted at slowing CKD progression in a majority. However, attributes of CKD may differ substantially among persons living in low-income and middle-income countries (LMICs). We used data from population or community-based studies from five LMICs (China, urban India, Moldova, Nepal and Nigeria) to determine what proportion of persons with CKD living in diverse regions fit one of the three major clinical profiles, with data from the US National Health Nutrition and Examination Survey as reference. In the USA, urban India and Moldova, 79.0%-83.9%; in China and Nepal, 62.4%-66.7% and in Nigeria, 51.6% persons with CKD fit one of three established risk profiles. Diabetes was most common in urban India and vascular disease in Moldova (50.7% and 33.2% of persons with CKD in urban India and Moldova, respectively). In Nigeria, 17.8% of persons with CKD without established risk factors had albuminuria ≥300 mg/g, the highest proportion in any country. While the majority of persons with CKD in LMICs fit into one of three established risk profiles, the proportion of persons who have CKD without established risk factors is higher than in the USA. These findings can inform tailored CKD detection and management systems and highlight the importance of studying potential causes and outcomes of CKD without established risk factors in LMICs.

Hospitalizations and Nursing Facility Stays During the Transition from CKD to ESRD on Dialysis: An Observational Study.

BACKGROUND: There is little information on hospital and nursing facility stays during the transition from pre-dialysis kidney disease to end-stage renal disease treated with dialysis. OBJECTIVES: To examine hospital and nursing facility stays in the years pre- and post-dialysis initiation, and to develop a novel method for visualizing these data. DESIGN: Observational study of patients in the US Renal Data System initiating dialysis from October 2011 to October 2012. PARTICIPANTS: Patients aged ≥67 years with Medicare Part A/B coverage for 1 year pre-dialysis initiation. MAIN MEASURES: Proportion of patients with ≥1 facility day, and among these, the mean number of days and the mean proportion of time spent in a facility in the first year post-dialysis initiation. We created "heat maps" to represent data visually. KEY RESULTS: Among 28,049 patients, > 60% initiated dialysis in the hospital. Patients with at least 1 facility day spent 37-42 days in a facility in the year pre-dialysis initiation and 59-67 facility days in the year post-dialysis initiation. The duration of facility stay varied by age: patients aged 67-70 years spent 60 (95% CI 57-62) days or 25.8% of the first year post-dialysis initiation in a facility, while patients aged >80 years spent 67 (CI 65-69) days or 36.8% of the first year post-dialysis initiation in a facility. Patterns varied depending on the presence or absence of certain comorbid conditions, with dementia having a particularly large effect: patients with dementia spent approximately 50% of the first year post-dialysis initiation in a facility, regardless of age. CONCLUSIONS: Older patients, particularly octogenarians and patients with dementia or other comorbidities, spend a large proportion of time in a facility during the first year after dialysis initiation. Our heat maps provide a novel and concise visual representation of a large amount of quantitative data regarding expected outcomes after initiation of dialysis.