Longitudinal change in mitochondrial heteroplasmy exhibits positive selection for deleterious variants.

A common feature of human aging is the acquisition of somatic mutations, and mitochondria are particularly prone to mutation due to their inefficient DNA repair and close proximity to reactive oxygen species, leading to a state of mitochondrial DNA heteroplasmy1,2. Cross-sectional studies have demonstrated that detection of heteroplasmy increases with participant age3, a phenomenon that has been attributed to genetic drift4-7. In this first large-scale longitudinal study, we measured heteroplasmy in two prospective cohorts (combined n=1405) at two timepoints (mean time between visits, 8.6 years), demonstrating that deleterious heteroplasmies were more likely to increase in variant allele fraction (VAF). We further demonstrated that increase in VAF was associated with increased risk of overall mortality. These results challenge the claim that somatic mtDNA mutations arise mainly due to genetic drift, instead demonstrating positive selection for predicted deleterious mutations at the cellular level, despite an negative impact on overall mortality.

Clonal hematopoiesis of indeterminate potential is associated with acute kidney injury.

Age is a predominant risk factor for acute kidney injury (AKI), yet the biological mechanisms underlying this risk are largely unknown. Clonal hematopoiesis of indeterminate potential (CHIP) confers increased risk for several chronic diseases associated with aging. Here we sought to test whether CHIP increases the risk of AKI. In three population-based epidemiology cohorts, we found that CHIP was associated with a greater risk of incident AKI, which was more pronounced in patients with AKI requiring dialysis and in individuals with somatic mutations in genes other than DNMT3A, including mutations in TET2 and JAK2. Mendelian randomization analyses supported a causal role for CHIP in promoting AKI. Non-DNMT3A-CHIP was also associated with a nonresolving pattern of injury in patients with AKI. To gain mechanistic insight, we evaluated the role of Tet2-CHIP and Jak2V617F-CHIP in two mouse models of AKI. In both models, CHIP was associated with more severe AKI, greater renal proinflammatory macrophage infiltration and greater post-AKI kidney fibrosis. In summary, this work establishes CHIP as a genetic mechanism conferring impaired kidney function recovery after AKI via an aberrant inflammatory response mediated by renal macrophages.

Mitochondrial DNA copy number reduction via in vitro TFAM knockout remodels the nuclear epigenome and transcriptome.

Mitochondrial DNA copy number (mtDNA-CN) is associated with several age-related chronic diseases and is a predictor of all-cause mortality. Here, we examine site-specific differential nuclear DNA (nDNA) methylation and differential gene expression resulting from in vitro reduction of mtDNA-CN to uncover shared genes and biological pathways mediating the effect of mtDNA-CN on disease. Epigenome and transcriptome profiles were generated for three independent human embryonic kidney (HEK293T) cell lines harbouring a mitochondrial transcription factor A (TFAM) heterozygous knockout generated via CRISPR-Cas9, and matched control lines. We identified 4,242 differentially methylated sites, 228 differentially methylated regions, and 179 differentially expressed genes associated with mtDNA-CN. Integrated analysis uncovered 381 Gene-CpG pairs. GABAA receptor genes and related pathways, the neuroactive ligand receptor interaction pathway, ABCD1/2 gene activity, and cell signalling processes were overrepresented, providing insight into the underlying biological mechanisms facilitating these associations. We also report evidence implicating chromatin state regulatory mechanisms as modulators of mtDNA-CN effect on gene expression. We demonstrate that mitochondrial DNA variation signals to the nuclear DNA epigenome and transcriptome and may lead to nuclear remodelling relevant to development, aging, and complex disease.

Mitochondrial genetic variation and risk of chronic kidney disease and acute kidney injury in UK Biobank participants.

Experimental models suggest an important role for mitochondrial dysfunction in the pathogenesis of chronic kidney disease (CKD) and acute kidney injury (AKI), but little is known regarding the impact of common mitochondrial genetic variation on kidney health. We sought to evaluate associations of inherited mitochondrial DNA (mtDNA) variation with risk of CKD and AKI in a large population-based cohort. We categorized UK Biobank participants who self-identified as white into eight distinct mtDNA haplotypes, which were previously identified based on their associations with phenotypes associated with mitochondrial DNA copy number, a measure of mitochondrial function. We used linear and logistic regression models to evaluate associations of these mtDNA haplotypes with estimated glomerular filtration rate by serum creatinine and cystatin C (eGFRCr-CysC, N = 362,802), prevalent (N = 416 cases) and incident (N = 405 cases) end-stage kidney disease (ESKD), AKI defined by diagnostic codes (N = 14,170 cases), and urine albumin/creatinine ratio (ACR, N = 114,662). The mean age was 57 ± 8 years and the mean eGFR was 90 ± 14 ml/min/1.73 m2. MtDNA haplotype was significantly associated with eGFR (p = 2.8E-12), but not with prevalent ESKD (p = 5.9E-2), incident ESKD (p = 0.93), AKI (p = 0.26), or urine ACR (p = 0.54). The association of mtDNA haplotype with eGFR remained significant after adjustment for diabetes mellitus and hypertension (p = 1.2E-10). When compared to the reference haplotype, mtDNA haplotypes I (ß = 0.402, standard error (SE) = 0.111; p = 2.7E-4), IV (ß = 0.430, SE = 0.073; p = 4.2E-9), and V (ß = 0.233, SE = 0.050; p = 2.7E-6) were each associated with higher eGFR. Among self-identified white UK Biobank participants, mtDNA haplotype was associated with eGFR, but not with ESKD, AKI or albuminuria.

Rare and common genetic determinants of mitochondrial function determine severity but not risk of amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease involving selective vulnerability of energy-intensive motor neurons (MNs). It has been unclear whether mitochondrial function is an upstream driver or a downstream modifier of neurotoxicity. We separated upstream genetic determinants of mitochondrial function, including genetic variation within the mitochondrial genome or autosomes; from downstream changeable factors including mitochondrial DNA copy number (mtCN). Across three cohorts including 6,437 ALS patients, we discovered that a set of mitochondrial haplotypes, chosen because they are linked to measurements of mitochondrial function, are a determinant of ALS survival following disease onset, but do not modify ALS risk. One particular haplotype appeared to be neuroprotective and was significantly over-represented in two cohorts of long-surviving ALS patients. Causal inference for mitochondrial function was achievable using mitochondrial haplotypes, but not autosomal SNPs in traditional Mendelian randomization (MR). Furthermore, rare loss-of-function genetic variants within, and reduced MN expression of, ACADM and DNA2 lead to ∼50 % shorter ALS survival; both proteins are implicated in mitochondrial function. Both mtCN and cellular vulnerability are linked to DNA2 function in ALS patient-derived neurons. Finally, MtCN responds dynamically to the onset of ALS independently of mitochondrial haplotype, and is correlated with disease severity. We conclude that, based on the genetic measures we have employed, mitochondrial function is a therapeutic target for amelioration of disease severity but not prevention of ALS.

Association analysis of mitochondrial DNA heteroplasmic variants: methods and application.

We rigorously assessed a comprehensive association testing framework for heteroplasmy, employing both simulated and real-world data. This framework employed a variant allele fraction (VAF) threshold and harnessed multiple gene-based tests for robust identification and association testing of heteroplasmy. Our simulation studies demonstrated that gene-based tests maintained an appropriate type I error rate at α=0.001. Notably, when 5% or more heteroplasmic variants within a target region were linked to an outcome, burden-extension tests (including the adaptive burden test, variable threshold burden test, and z-score weighting burden test) outperformed the sequence kernel association test (SKAT) and the original burden test. Applying this framework, we conducted association analyses on whole-blood derived heteroplasmy in 17,507 individuals of African and European ancestries (31% of African Ancestry, mean age of 62, with 58% women) with whole genome sequencing data. We performed both cohort- and ancestry-specific association analyses, followed by meta-analysis on both pooled samples and within each ancestry group. Our results suggest that mtDNA-encoded genes/regions are likely to exhibit varying rates in somatic aging, with the notably strong associations observed between heteroplasmy in the RNR1 and RNR2 genes (p<0.001) and advance aging by the Original Burden test. In contrast, SKAT identified significant associations (p<0.001) between diabetes and the aggregated effects of heteroplasmy in several protein-coding genes. Further research is warranted to validate these findings. In summary, our proposed statistical framework represents a valuable tool for facilitating association testing of heteroplasmy with disease traits in large human populations.

Risk Prediction in Male Adolescents With Congenital Long QT Syndrome: Implications for Sex-Specific Risk Stratification in Potassium Channel-Mediated Long QT Syndrome.

BACKGROUND: Sex-specific risk management may improve outcomes in congenital long QT syndrome (LQTS). We recently developed a prediction score for cardiac events (CEs) and life-threatening events (LTEs) in postadolescent women with LQTS. In the present study, we aimed to develop personalized risk estimates for the burden of CEs and LTEs in male adolescents with potassium channel-mediated LQTS. METHODS AND RESULTS: The prognostic model was derived from the LQTS Registry headquartered in Rochester, NY, comprising 611 LQT1 or LQT2 male adolescents from age 10 through 20 years, using the following variables: genotype/mutation location, QTc-specific thresholds, history of syncope, and ß-blocker therapy. Anderson-Gill modeling was performed for the end point of CE burden (total number of syncope, aborted cardiac arrest, and appropriate defibrillator shocks). The applicability of the CE prediction model was tested for the end point of the first LTE (excluding syncope and adding sudden cardiac death) using Cox modeling. A total of 270 CEs occurred during follow-up. The genotype-phenotype risk prediction model identified low-, intermediate-, and high-risk groups, comprising 74%, 14%, and 12% of the study population, respectively. Compared with the low-risk group, high-risk male subjects experienced a pronounced 5.2-fold increased risk of recurrent CEs (P<0.001), whereas intermediate-risk patients had a 2.1-fold (P=0.004) increased risk . At age 20 years, the low-, intermediate-, and high-risk adolescent male patients had on average 0.3, 0.6, and 1.4 CEs per person, respectively. Corresponding 10-year adjusted probabilities for a first LTE were 2%, 6%, and 8%. CONCLUSIONS: Personalized genotype-phenotype risk estimates can be used to guide sex-specific management in male adolescents with potassium channel-mediated LQTS.

Polygenic Risk Score Associates With Atherosclerotic Plaque Characteristics at Autopsy.

BACKGROUND: Polygenic risk scores (PRSs) for coronary artery disease (CAD) potentially improve cardiovascular risk prediction. However, their relationship with histopathologic features of CAD has never been examined systematically. METHODS: From 4327 subjects referred to CVPath by the State of Maryland Office Chief Medical Examiner for sudden death between 1994 and 2015, 2455 cases were randomly selected for genotyping. We generated PRS from 291 known CAD risk loci. Detailed histopathologic examination of the coronary arteries was performed in all subjects. The primary study outcome measurements were histopathologic plaque features determining severity of atherosclerosis, including %stenosis, calcification, thin-cap fibroatheromas, and thrombotic CAD. RESULTS: After exclusion of cases with insufficient DNA sample quality or with missing data, 954 cases (mean age, 48.8±14.7 years; 75.7% men) remained in the final study cohort. Subjects in the highest PRS quintile exhibited more severe atherosclerosis compared with subjects in the lowest quintile, with greater %stenosis (80.3%±27.0% versus 50.4%±38.7%; adjusted P<0.001) and a higher frequency of calcification (69.6% versus 35.8%; adjusted P=0.004) and thin-cap fibroatheroma (26.7% versus 9.5%; adjusted P=0.007). Even after adjustment for traditional CAD risk factors, subjects within the highest PRS quintile had higher odds of severe atherosclerosis (ie, ≥75% stenosis; adjusted odds ratio, 3.77 [95% CI, 2.10-6.78]; P<0.001) and plaque rupture (adjusted odds ratio, 4.05 [95% CI, 2.26-7.24]; P<0.001). Moreover, subjects within the highest quintile had higher odds of CAD-associated cause of death, especially among those aged ≤50 years (adjusted odds ratio, 4.08 [95% CI, 2.01-8.30]; P<0.001). No statistically significant associations were observed with plaque erosion after adjusting for covariates. CONCLUSIONS: This is the first autopsy study investigating associations between PRS and atherosclerosis severity at the histopathologic level in subjects with sudden death. Our pathological analysis suggests PRS correlates with plaque burden and features of advanced atherosclerosis and may be useful as a method for CAD risk stratification, especially in younger subjects.

Mitochondrial DNA Copy Number Variation in Asthma Risk, Severity, and Exacerbations.

Rationale: Although airway oxidative stress and inflammation are central to asthma pathogenesis, there is limited knowledge of the relationship of asthma risk, severity, or exacerbations to mitochondrial dysfunction, which is pivotal to oxidant generation and inflammation. Objectives: We investigated whether mitochondrial DNA copy number (mtDNA-CN) as a measure of mitochondrial function is associated with asthma diagnosis, severity, oxidative stress, and exacerbations. Methods: We measured mtDNA-CN in blood in two cohorts. In the UK Biobank (UKB), we compared mtDNA-CN in mild and moderate-severe asthmatics to non-asthmatics. In the Severe Asthma Research Program (SARP), we evaluated mtDNA-CN in relation to asthma severity, biomarkers of oxidative stress and inflammation, and exacerbations. Measures and Main Results: In UK Biobank, asthmatics (n = 29,768) have lower mtDNA-CN compared to non-asthmatics (n = 239,158) (beta, -0.026 [95% CI, -0.038 to -0.014], P = 2.46×10-5). While lower mtDNA-CN is associated with asthma, mtDNA-CN did not differ by asthma severity in either UKB or SARP. Biomarkers of inflammation show that asthmatics have higher white blood cells (WBC), neutrophils, eosinophils, fraction exhaled nitric oxide (FENO), and lower superoxide dismutase (SOD) than non-asthmatics, confirming greater oxidative stress in asthma. In one year follow-up in SARP, higher mtDNA-CN is associated with reduced risk of three or more exacerbations in the subsequent year (OR 0.352 [95% CI, 0.164 to 0.753], P = 0.007). Conclusions: Asthma is characterized by mitochondrial dysfunction. Higher mtDNA-CN identifies an exacerbation-resistant asthma phenotype, suggesting mitochondrial function is important in exacerbation risk.

Association Between Whole Blood-Derived Mitochondrial DNA Copy Number, Low-Density Lipoprotein Cholesterol, and Cardiovascular Disease Risk.

Background The relationship between mitochondrial DNA copy number (mtDNA CN) and cardiovascular disease remains elusive. Methods and Results We performed cross-sectional and prospective association analyses of blood-derived mtDNA CN and cardiovascular disease outcomes in 27 316 participants in 8 cohorts of multiple racial and ethnic groups with whole-genome sequencing. We also performed Mendelian randomization to explore causal relationships of mtDNA CN with coronary heart disease (CHD) and cardiometabolic risk factors (obesity, diabetes, hypertension, and hyperlipidemia). P<0.01 was used for significance. We validated most of the previously reported associations between mtDNA CN and cardiovascular disease outcomes. For example, 1-SD unit lower level of mtDNA CN was associated with 1.08 (95% CI, 1.04-1.12; P<0.001) times the hazard for developing incident CHD, adjusting for covariates. Mendelian randomization analyses showed no causal effect from a lower level of mtDNA CN to a higher CHD risk (ß=0.091; P=0.11) or in the reverse direction (ß=-0.012; P=0.076). Additional bidirectional Mendelian randomization analyses revealed that low-density lipoprotein cholesterol had a causal effect on mtDNA CN (ß=-0.084; P<0.001), but the reverse direction was not significant (P=0.059). No causal associations were observed between mtDNA CN and obesity, diabetes, and hypertension, in either direction. Multivariable Mendelian randomization analyses showed no causal effect of CHD on mtDNA CN, controlling for low-density lipoprotein cholesterol level (P=0.52), whereas there was a strong direct causal effect of higher low-density lipoprotein cholesterol on lower mtDNA CN, adjusting for CHD status (ß=-0.092; P<0.001). Conclusions Our findings indicate that high low-density lipoprotein cholesterol may underlie the complex relationships between mtDNA CN and vascular atherosclerosis.

The functional impact of rare variation across the regulatory cascade.

Each human genome has tens of thousands of rare genetic variants; however, identifying impactful rare variants remains a major challenge. We demonstrate how use of personal multi-omics can enable identification of impactful rare variants by using the Multi-Ethnic Study of Atherosclerosis, which included several hundred individuals, with whole-genome sequencing, transcriptomes, methylomes, and proteomes collected across two time points, 10 years apart. We evaluated each multi-omics phenotype's ability to separately and jointly inform functional rare variation. By combining expression and protein data, we observed rare stop variants 62 times and rare frameshift variants 216 times as frequently as controls, compared to 13-27 times as frequently for expression or protein effects alone. We extended a Bayesian hierarchical model, "Watershed," to prioritize specific rare variants underlying multi-omics signals across the regulatory cascade. With this approach, we identified rare variants that exhibited large effect sizes on multiple complex traits including height, schizophrenia, and Alzheimer's disease.

Mitochondrial DNA copy number is associated with incident chronic kidney disease and proteinuria in the AIDS linked to the intravenous experience cohort.

We evaluated the prospective association of mitochondrial DNA copy number (mtDNA CN) with markers of kidney function among a cohort of persons who inject drugs (PWID). This is a Prospective cohort study nested in the AIDS linked to the intravenous experience cohort (community-based cohort of PWID in Baltimore, MD). mtDNA CN was measured at two time-points 5 years apart using a real-time polymerase chain reaction. Kidney function (estimated glomerular filtration rate [eGFR], serum creatinine, urine protein) was measured annually. We used linear mixed effects models to evaluate kidney function trajectories (N = 946) and Cox regression models to assess hazard of incident CKD (eGFR < 60 at two consecutive visits, N = 739) and proteinuria (urine protein:creatinine ratio > 200, N = 573) by level of mtDNA CN (Low [lowest quartile], vs high [other three quartiles]. Models were adjusted for demographic and behavioral characteristics, HIV and/or HCV infection, and comorbidity burden. Low mtDNA CN was independently associated with higher hazard of incident CKD (aHR: 2.33, 95% CI 1.42, 3.80) and proteinuria (aHR: 1.42, 95% CI 1.04, 1.96). Participants with low mtDNA CN had greater declines in eGFR and greater increases in serum creatinine over time. Low mtDNA CN is associated with more rapid kidney function decline and risk of incident CKD and proteinuria.

Deleterious heteroplasmic mitochondrial mutations are associated with an increased risk of overall and cancer-specific mortality.

Mitochondria carry their own circular genome and disruption of the mitochondrial genome is associated with various aging-related diseases. Unlike the nuclear genome, mitochondrial DNA (mtDNA) can be present at 1000 s to 10,000 s copies in somatic cells and variants may exist in a state of heteroplasmy, where only a fraction of the DNA molecules harbors a particular variant. We quantify mtDNA heteroplasmy in 194,871 participants in the UK Biobank and find that heteroplasmy is associated with a 1.5-fold increased risk of all-cause mortality. Additionally, we functionally characterize mtDNA single nucleotide variants (SNVs) using a constraint-based score, mitochondrial local constraint score sum (MSS) and find it associated with all-cause mortality, and with the prevalence and incidence of cancer and cancer-related mortality, particularly leukemia. These results indicate that mitochondria may have a functional role in certain cancers, and mitochondrial heteroplasmic SNVs may serve as a prognostic marker for cancer, especially for leukemia.

Genetic insights into resting heart rate and its role in cardiovascular disease.

Resting heart rate is associated with cardiovascular diseases and mortality in observational and Mendelian randomization studies. The aims of this study are to extend the number of resting heart rate associated genetic variants and to obtain further insights in resting heart rate biology and its clinical consequences. A genome-wide meta-analysis of 100 studies in up to 835,465 individuals reveals 493 independent genetic variants in 352 loci, including 68 genetic variants outside previously identified resting heart rate associated loci. We prioritize 670 genes and in silico annotations point to their enrichment in cardiomyocytes and provide insights in their ECG signature. Two-sample Mendelian randomization analyses indicate that higher genetically predicted resting heart rate increases risk of dilated cardiomyopathy, but decreases risk of developing atrial fibrillation, ischemic stroke, and cardio-embolic stroke. We do not find evidence for a linear or non-linear genetic association between resting heart rate and all-cause mortality in contrast to our previous Mendelian randomization study. Systematic alteration of key differences between the current and previous Mendelian randomization study indicates that the most likely cause of the discrepancy between these studies arises from false positive findings in previous one-sample MR analyses caused by weak-instrument bias at lower P-value thresholds. The results extend our understanding of resting heart rate biology and give additional insights in its role in cardiovascular disease development.

Genetic variant annotation scores in congenital long QT syndrome.

BACKGROUND: Congenital Long QT Syndrome (LQTS) is a hereditary arrhythmic disorder. We aimed to assess the performance of current genetic variant annotation scores among LQTS patients and their predictive impact. METHODS: We evaluated 2025 patients with unique mutations for LQT1-LQT3. A patient-specific score was calculated for each of four established genetic variant annotation algorithms: CADD, SIFT, REVEL, and PolyPhen-2. The scores were tested for the identification of LQTS and their predictive performance for cardiac events (CE) and life-threatening events (LTE) and then compared with the predictive performance of LQTS categorization based on mutation location/function. Score performance was tested using Harrell's C-index. RESULTS: A total of 917 subjects were classified as LQT1, 838 as LQT2, and 270 as LQT3. The identification of a pathogenic variant occurred in 99% with CADD, 92% with SIFT, 100% with REVEL, and 86% with PolyPhen-2. However, none of the genetic scores correlated with the risk of CE (Harrell's C-index: CADD = 0.50, SIFT = 0.51, REVEL = 0.50, and PolyPhen-2 = 0.52) or LTE (Harrell's C-index: CADD = 0.50, SIFT = 0.53, REVEL = 0.54, and PolyPhen-2 = 0.52). In contrast, high-risk mutation categorization based on location/function was a powerful independent predictor of CE (HR = 1.88; p < .001) and LTE (HR = 1.89, p < .001). CONCLUSION: In congenital LQTS patients, well-established algorithms (CADD, SIFT, REVEL, and PolyPhen-2) were able to identify the majority of the causal variants as pathogenic. However, the scores did not predict clinical outcomes. These results indicate that mutation location/functional assays are essential for accurate interpretation of the risk associated with LQTS mutations.

Polygenic Risk Score Associates with Atherosclerotic Plaque Characteristics at Autopsy.

Background: Polygenic risk scores (PRS) for coronary artery disease (CAD) potentially improve cardiovascular risk prediction. However, their relationship with histopathologic features of CAD has never been examined systematically. Methods: From 4,327 subjects referred to CVPath by the State of Maryland Office Chief Medical Examiner (OCME) for sudden death between 1994 and 2015, 2,455 cases were randomly selected for genotyping. We generated PRS from 291 known CAD risk loci. Detailed histopathologic examination of the coronary arteries was performed in all subjects. The primary study outcome measurements were histopathologic plaque features determining severity of atherosclerosis, including %stenosis, calcification, thin-cap fibroatheromas (TCFA), and thrombotic CAD. Results: After exclusion of cases with insufficient DNA sample quality or with missing data, 954 cases (mean age 48.8±14.7; 75.7% men) remained in the final study cohort. Subjects in the highest PRS quintile exhibited more severe atherosclerosis compared to subjects in the lowest quintile, with greater %stenosis (80.3%±27.0% vs. 50.4%±38.7%; adjusted p<0.001) and a higher frequency of calcification (69.6% vs. 35.8%; adjusted p=0.004) and TCFAs (26.7% vs. 9.5%; adjusted p=0.007). Even after adjustment for traditional CAD risk factors subjects within the highest PRS quintile had higher odds of severe atherosclerosis (i.e., ≥75% stenosis; adjusted OR 3.77; 95%CI 2.10-6.78; p<0.001) and plaque rupture (adjusted OR 4.05; 95%CI 2.26-7.24; p<0.001). Moreover, subjects within the highest quintile had higher odds of CAD-associated cause of death, especially among those aged 50 years and younger (adjusted OR 4.08; 95%CI 2.01-8.30; p<0.001). No associations were observed with plaque erosion. Conclusions: This is the first autopsy study investigating associations between PRS and atherosclerosis severity at the histopathologic level in subjects with sudden death. Our pathological analysis suggests PRS correlates with plaque burden and features of advanced atherosclerosis and may be useful as a method for CAD risk stratification, especially in younger subjects. Highlights: In this autopsy study including 954 subjects within the CVPath Sudden Death Registry, high PRS correlated with plaque burden and atherosclerosis severity.The PRS showed differential associations with plaque rupture and plaque erosion, suggesting different etiologies to these two causes of thrombotic CAD.PRS may be useful for risk stratification, particularly in the young. Further examination of individual risk loci and their association with plaque morphology may help understand molecular mechanisms of atherosclerosis, potentially revealing new therapy targets of CAD. Graphic Abstract: A polygenic risk score, generated from 291 known CAD risk loci, was assessed in 954 subjects within the CVPath Sudden Death Registry. Histopathologic examination of the coronary arteries was performed in all subjects. Subjects in the highest PRS quintile exhibited more severe atherosclerosis as compared to subjects in the lowest quintile, with a greater plaque burden, more calcification, and a higher frequency of plaque rupture.

Clonal Hematopoiesis of Indeterminate Potential is Associated with Acute Kidney Injury.

Age is a predominant risk factor for acute kidney injury (AKI), yet the biological mechanisms underlying this risk are largely unknown and to date no genetic mechanisms for AKI have been established. Clonal hematopoiesis of indeterminate potential (CHIP) is a recently recognized biological mechanism conferring risk of several chronic aging diseases including cardiovascular disease, pulmonary disease and liver disease. In CHIP, blood stem cells acquire mutations in myeloid cancer driver genes such as DNMT3A, TET2, ASXL1 and JAK2 and the myeloid progeny of these mutated cells contribute to end-organ damage through inflammatory dysregulation. We sought to establish whether CHIP causes acute kidney injury (AKI). To address this question, we first evaluated associations with incident AKI events in three population-based epidemiology cohorts (N = 442,153). We found that CHIP was associated with a greater risk of AKI (adjusted HR 1.26, 95% CI: 1.19-1.34, p<0.0001), which was more pronounced in patients with AKI requiring dialysis (adjusted HR 1.65, 95% CI: 1.24-2.20, p=0.001). The risk was particularly high in the subset of individuals where CHIP was driven by mutations in genes other than DNMT3A (HR: 1.49, 95% CI: 1.37-1.61, p<0.0001). We then examined the association between CHIP and recovery from AKI in the ASSESS-AKI cohort and identified that non-DNMT3A CHIP was more common among those with a non-resolving pattern of injury (HR 2.3, 95% CI: 1.14-4.64, p = 0.03). To gain mechanistic insight, we evaluated the role of Tet2-CHIP to AKI in ischemia-reperfusion injury (IRI) and unilateral ureteral obstruction (UUO) mouse models. In both models, we observed more severe AKI and greater post-AKI kidney fibrosis in Tet2-CHIP mice. Kidney macrophage infiltration was markedly increased in Tet2-CHIP mice and Tet2-CHIP mutant renal macrophages displayed greater proinflammatory responses. In summary, this work establishes CHIP as a genetic mechanism conferring risk of AKI and impaired kidney function recovery following AKI via an aberrant inflammatory response in CHIP derived renal macrophages.

Transcriptome- and proteome-wide association studies nominate determinants of kidney function and damage.

BACKGROUND: The pathophysiological causes of kidney disease are not fully understood. Here we show that the integration of genome-wide genetic, transcriptomic, and proteomic association studies can nominate causal determinants of kidney function and damage. RESULTS: Through transcriptome-wide association studies (TWAS) in kidney cortex, kidney tubule, liver, and whole blood and proteome-wide association studies (PWAS) in plasma, we assess for effects of 12,893 genes and 1342 proteins on kidney filtration (glomerular filtration rate (GFR) estimated by creatinine; GFR estimated by cystatin C; and blood urea nitrogen) and kidney damage (albuminuria). We find 1561 associations distributed among 260 genomic regions that are supported as putatively causal. We then prioritize 153 of these genomic regions using additional colocalization analyses. Our genome-wide findings are supported by existing knowledge (animal models for MANBA, DACH1, SH3YL1, INHBB), exceed the underlying GWAS signals (28 region-trait combinations without significant GWAS hit), identify independent gene/protein-trait associations within the same genomic region (INHBC, SPRYD4), nominate tissues underlying the associations (tubule expression of NRBP1), and distinguish markers of kidney filtration from those with a role in creatinine and cystatin C metabolism. Furthermore, we follow up on members of the TGF-beta superfamily of proteins and find a prognostic value of INHBC for kidney disease progression even after adjustment for measured glomerular filtration rate (GFR). CONCLUSION: In summary, this study combines multimodal, genome-wide association studies to generate a catalog of putatively causal target genes and proteins relevant to kidney function and damage which can guide follow-up studies in physiology, basic science, and clinical medicine.