Whole-exome and whole-genome sequencing of 1064 individuals with type 1 diabetes reveals novel genes for diabetic kidney disease.

AIMS/HYPOTHESIS: Diabetic kidney disease (DKD) is a severe diabetic complication that affects one third of individuals with type 1 diabetes. Although several genes and common variants have been shown to be associated with DKD, much of the predicted inheritance remains unexplained. Here, we performed next-generation sequencing to assess whether low-frequency variants, extending to a minor allele frequency (MAF) ≤10% (single or aggregated) contribute to the missing heritability in DKD. METHODS: We performed whole-exome sequencing (WES) of 498 individuals and whole-genome sequencing (WGS) of 599 individuals with type 1 diabetes. After quality control, next-generation sequencing data were available for a total of 1064 individuals, of whom 541 had developed either severe albuminuria or end-stage kidney disease, and 523 had retained normal albumin excretion despite a long duration of type 1 diabetes. Single-variant and gene-aggregate tests for protein-altering variants (PAV) and protein-truncating variants (PTV) were performed separately for WES and WGS data and combined in a meta-analysis. We also performed genome-wide aggregate analyses on genomic windows (sliding window), promoters and enhancers using the WGS dataset. RESULTS: In the single-variant meta-analysis, no variant reached genome-wide significance, but a suggestively associated common THAP7 rs369250 variant (p=1.50 × 10-5, MAF=49%) was replicated in the FinnGen general population genome-wide association study (GWAS) data for chronic kidney disease and DKD phenotypes. The gene-aggregate meta-analysis provided suggestive evidence (p<4.0 × 10-4) at four genes for DKD, of which NAT16 (MAFPAV≤10%) and LTA (also known as TNFß, MAFPAV≤5%) are replicated in the FinnGen general population GWAS data. The LTA rs2229092 C allele was associated with significantly lower TNFR1, TNFR2 and TNFR3 serum levels in a subset of FinnDiane participants. Of the intergenic regions suggestively associated with DKD, the enhancer on chromosome 18q12.3 (p=3.94 × 10-5, MAFvariants≤5%) showed interaction with the METTL4 gene; the lead variant was replicated, and predicted to alter binding of the MafB transcription factor. CONCLUSIONS/INTERPRETATION: Our sequencing-based meta-analysis revealed multiple genes, variants and regulatory regions that were suggestively associated with DKD. However, as no variant or gene reached genome-wide significance, further studies are needed to validate the findings.

A methylation risk score for chronic kidney disease: a HyperGEN study.

Chronic kidney disease (CKD) impacts about 1 in 7 adults in the United States, but African Americans (AAs) carry a disproportionately higher burden of disease. Epigenetic modifications, such as DNA methylation at cytosine-phosphate-guanine (CpG) sites, have been linked to kidney function and may have clinical utility in predicting the risk of CKD. Given the dynamic relationship between the epigenome, environment, and disease, AAs may be especially sensitive to environment-driven methylation alterations. Moreover, risk models incorporating CpG methylation have been shown to predict disease across multiple racial groups. In this study, we developed a methylation risk score (MRS) for CKD in cohorts of AAs. We selected nine CpG sites that were previously reported to be associated with estimated glomerular filtration rate (eGFR) in epigenome-wide association studies to construct a MRS in the Hypertension Genetic Epidemiology Network (HyperGEN). In logistic mixed models, the MRS was significantly associated with prevalent CKD and was robust to multiple sensitivity analyses, including CKD risk factors. There was modest replication in validation cohorts. In summary, we demonstrated that an eGFR-based CpG score is an independent predictor of prevalent CKD, suggesting that MRS should be further investigated for clinical utility in evaluating CKD risk and progression.

Advancing Genetic Testing in Kidney Diseases: Report From a National Kidney Foundation Working Group.

About 37 million people in the United States have chronic kidney disease, a disease that encompasses diseases of multiple causes. About 10% or more of kidney diseases in adults and about 70% of selected chronic kidney diseases in children are expected to be explained by genetic causes. Despite the advances in genetic testing and an increasing understanding of the genetic bases of certain kidney diseases, genetic testing in nephrology lags behind other medical fields. More understanding of the benefits and logistics of genetic testing is needed to advance the implementation of genetic testing in chronic kidney diseases. Accordingly, the National Kidney Foundation convened a Working Group of experts with diverse expertise in genetics, nephrology, and allied fields to develop recommendations for genetic testing for monogenic disorders and to identify genetic risk factors for oligogenic and polygenic causes of kidney diseases. Algorithms for clinical decision making on genetic testing and a road map for advancing genetic testing in kidney diseases were generated. An important aspect of this initiative was the use of a modified Delphi process to reach group consensus on the recommendations. The recommendations and resources described herein provide support to nephrologists and allied health professionals to advance the use of genetic testing for diagnosis and screening of kidney diseases.

Admixture Mapping of Chronic Kidney Disease and Risk Factors in Hispanic/Latino Individuals From Central America Country of Origin.

BACKGROUND: Chronic kidney disease (CKD) is highly prevalent in Central America, and genetic factors may contribute to CKD risk. To understand the influences of genetic admixture on CKD susceptibility, we conducted an admixture mapping screening of CKD traits and risk factors in US Hispanic and Latino individuals from Central America country of origin. METHODS: We analyzed 1023 participants of HCHS/SOL (Hispanic Community Health Study/Study of Latinos) who reported 4 grandparents originating from the same Central America country. Ancestry admixture findings were validated on 8191 African Americans from WHI (Women's Health Initiative), 3141 American Indians from SHS (Strong Heart Study), and over 1.1 million European individuals from a multistudy meta-analysis. RESULTS: We identified 3 novel genomic regions for albuminuria (chromosome 14q24.2), CKD (chromosome 6q25.3), and type 2 diabetes (chromosome 3q22.2). The 14q24.2 locus driven by a Native American ancestry had a protective effect on albuminuria and consisted of 2 nearby regions spanning the RGS6 gene. Variants at this locus were validated in American Indians. The 6q25.3 African ancestry-derived locus, encompassing the ARID1B gene, was associated with increased risk for CKD and replicated in African Americans through admixture mapping. The European ancestry type 2 diabetes locus at 3q22.2, encompassing the EPHB1 and KY genes, was validated in European individuals through variant association. CONCLUSIONS: US Hispanic/Latino populations are culturally and genetically diverse. This study focusing on Central America grandparent country of origin provides new loci discovery and insights into the ancestry-of-origin influences on CKD and risk factors in US Hispanic and Latino individuals.

Understanding the genetic complexity of puberty timing across the allele frequency spectrum.

Pubertal timing varies considerably and is associated with later health outcomes. We performed multi-ancestry genetic analyses on ~800,000 women, identifying 1,080 signals for age at menarche. Collectively, these explained 11% of trait variance in an independent sample. Women at the top and bottom 1% of polygenic risk exhibited ~11 and ~14-fold higher risks of delayed and precocious puberty, respectively. We identified several genes harboring rare loss-of-function variants in ~200,000 women, including variants in ZNF483, which abolished the impact of polygenic risk. Variant-to-gene mapping approaches and mouse gonadotropin-releasing hormone neuron RNA sequencing implicated 665 genes, including an uncharacterized G-protein-coupled receptor, GPR83, which amplified the signaling of MC3R, a key nutritional sensor. Shared signals with menopause timing at genes involved in DNA damage response suggest that the ovarian reserve might signal centrally to trigger puberty. We also highlight body size-dependent and independent mechanisms that potentially link reproductive timing to later life disease.

Nutritional Strategies against Diabetic Nephropathy: Insights from Animal Studies and Human Trials.

Diabetic nephropathy (DN), defined as continuously elevated urinary albumin and a diminished estimated glomerular filtration rate, is a serious complication of both type 1 diabetes and type 2 diabetes and is the main cause of end-stage kidney disease. Patients with end-stage renal disease require chronic kidney dialysis and/or a kidney transplantation. Research highlights the role of diet in modulating specific signaling pathways that are instrumental in the progression of DN. Nutrient-sensitive pathways, affected by nutritional compounds and dietary components, offer a novel perspective on the management of DN by influencing inflammation, oxidative stress, and nutrient metabolism. Animal models have identified signaling pathways related to glucose metabolism, inflammation responses, autophagy, and lipid metabolism, while human population studies have contributed to the clinical significance of designing medical and nutritional therapies to attenuate DN progression. Here, we will update recent progress in research into the renoprotective or therapeutic effects of nutritional compounds, and potential nutrition-modulated pathways.

Determinants of mosaic chromosomal alteration fitness.

Clonal hematopoiesis (CH) is characterized by the acquisition of a somatic mutation in a hematopoietic stem cell that results in a clonal expansion. These driver mutations can be single nucleotide variants in cancer driver genes or larger structural rearrangements called mosaic chromosomal alterations (mCAs). The factors that influence the variations in mCA fitness and ultimately result in different clonal expansion rates are not well understood. We used the Passenger-Approximated Clonal Expansion Rate (PACER) method to estimate clonal expansion rate as PACER scores for 6,381 individuals in the NHLBI TOPMed cohort with gain, loss, and copy-neutral loss of heterozygosity mCAs. Our mCA fitness estimates, derived by aggregating per-individual PACER scores, were correlated (R2 = 0.49) with an alternative approach that estimated fitness of mCAs in the UK Biobank using population-level distributions of clonal fraction. Among individuals with JAK2 V617F clonal hematopoiesis of indeterminate potential or mCAs affecting the JAK2 gene on chromosome 9, PACER score was strongly correlated with erythrocyte count. In a cross-sectional analysis, genome-wide association study of estimates of mCA expansion rate identified a TCL1A locus variant associated with mCA clonal expansion rate, with suggestive variants in NRIP1 and TERT.

African ancestry-derived APOL1 risk genotypes show proximal epigenetic associations.

Apolipoprotein L1 (APOL1) coding variants, termed G1 and G2, are established genetic risk factors for a growing spectrum of diseases, including kidney disease, in individuals of African ancestry. Evidence suggests that the risk variants, which show a recessive mode of inheritance, lead to toxic gain-of-function changes of the APOL1 protein. Disease occurrence and presentation vary, likely due to modifiers or second hits. To understand the role of the epigenetic landscape in relation to APOL1 risk variants, we performed methylation quantitative trait locus (meQTL) analysis to identify differentially methylated CpGs influenced by APOL1 risk variants in 611 African American individuals. We identified five CpGs that were significantly associated with APOL1 risk alleles in discovery and replication studies, and one CpG-APOL1 association was independent of other genomic variants. Our study highlights proximal DNA methylation alterations that may help explain the variable disease risk and clinical manifestation of APOL1 variants.

Genetic variant rs1205 is associated with COVID-19 outcomes: The Strong Heart Study and Strong Heart Family Study.

BACKGROUND: Although COVID-19 infection has been associated with a number of clinical and environmental risk factors, host genetic variation has also been associated with the incidence and morbidity of infection. The CRP gene codes for a critical component of the innate immune system and CRP variants have been reported associated with infectious disease and vaccination outcomes. We investigated possible associations between COVID-19 outcome and a limited number of candidate gene variants including rs1205. METHODOLOGY/PRINCIPAL FINDINGS: The Strong Heart and Strong Heart Family studies have accumulated detailed genetic, cardiovascular risk and event data in geographically dispersed American Indian communities since 1988. Genotypic data and 91 COVID-19 adjudicated deaths or hospitalizations from 2/1/20 through 3/1/23 were identified among 3,780 participants in two subsets. Among 21 candidate variants including genes in the interferon response pathway, APOE, TMPRSS2, TLR3, the HLA complex and the ABO blood group, only rs1205, a 3' untranslated region variant in the CRP gene, showed nominally significant association in T-dominant model analyses (odds ratio 1.859, 95%CI 1.001-3.453, p = 0.049) after adjustment for age, sex, center, body mass index, and a history of cardiovascular disease. Within the younger subset, association with the rs1205 T-Dom genotype was stronger, both in the same adjusted logistic model and in the SOLAR analysis also adjusting for other genetic relatedness. CONCLUSION: A T-dominant genotype of rs1205 in the CRP gene is associated with COVID-19 death or hospitalization, even after adjustment for relevant clinical factors and potential participant relatedness. Additional study of other populations and genetic variants of this gene are warranted.

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.

Multiple cardiovascular risk factor care in 55 low- and middle-income countries: A cross-sectional analysis of nationally-representative, individual-level data from 280,783 adults.

The prevalence of multiple age-related cardiovascular disease (CVD) risk factors is high among individuals living in low- and middle-income countries. We described receipt of healthcare services for and management of hypertension and diabetes among individuals living with these conditions using individual-level data from 55 nationally representative population-based surveys (2009-2019) with measured blood pressure (BP) and diabetes biomarker. We restricted our analysis to non-pregnant individuals aged 40-69 years and defined three mutually exclusive groups (i.e., hypertension only, diabetes only, and both hypertension-diabetes) to compare individuals living with concurrent hypertension and diabetes to individuals with each condition separately. We included 90,086 individuals who lived with hypertension only, 11,975 with diabetes only, and 16,228 with hypertension-diabetes. We estimated the percentage of individuals who were aware of their diagnosis, used pharmacological therapy, or achieved appropriate hypertension and diabetes management. A greater percentage of individuals with hypertension-diabetes were fully diagnosed (64.1% [95% CI: 61.8-66.4]) than those with hypertension only (47.4% [45.3-49.6]) or diabetes only (46.7% [44.1-49.2]). Among the hypertension-diabetes group, pharmacological treatment was higher for individual conditions (38.3% [95% CI: 34.8-41.8] using antihypertensive and 42.3% [95% CI: 39.4-45.2] using glucose-lowering medications) than for both conditions jointly (24.6% [95% CI: 22.1-27.2]).The percentage of individuals achieving appropriate management was highest in the hypertension group (17.6% [16.4-18.8]), followed by diabetes (13.3% [10.7-15.8]) and hypertension-diabetes (6.6% [5.4-7.8]) groups. Although health systems in LMICs are reaching a larger share of individuals living with both hypertension and diabetes than those living with just one of these conditions, only seven percent achieved both BP and blood glucose treatment targets. Implementation of cost-effective population-level interventions that shift clinical care paradigm from disease-specific to comprehensive CVD care are urgently needed for all three groups, especially for those with multiple CVD risk factors.

Genome-Wide Interaction Analysis With DASH Diet Score Identified Novel Loci for Systolic Blood Pressure.

BACKGROUND: The Dietary Approaches to Stop Hypertension (DASH) diet score lowers blood pressure (BP). We examined interactions between genotype and the DASH diet score in relation to systolic BP. METHODS: We analyzed up to 9 420 585 single nucleotide polymorphisms in up to 127 282 individuals of 6 population groups (91% of European population) from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (n=35 660) and UK Biobank (n=91 622) and performed European population-specific and cross-population meta-analyses. RESULTS: We identified 3 loci in European-specific analyses and an additional 4 loci in cross-population analyses at Pinteraction<5e-8. We observed a consistent interaction between rs117878928 at 15q25.1 (minor allele frequency, 0.03) and the DASH diet score (Pinteraction=4e-8; P for heterogeneity, 0.35) in European population, where the interaction effect size was 0.42±0.09 mm Hg (Pinteraction=9.4e-7) and 0.20±0.06 mm Hg (Pinteraction=0.001) in Cohorts for Heart and Aging Research in Genomic Epidemiology and the UK Biobank, respectively. The 1 Mb region surrounding rs117878928 was enriched with cis-expression quantitative trait loci (eQTL) variants (P=4e-273) and cis-DNA methylation quantitative trait loci variants (P=1e-300). Although the closest gene for rs117878928 is MTHFS, the highest narrow sense heritability accounted by single nucleotide polymorphisms potentially interacting with the DASH diet score in this locus was for gene ST20 at 15q25.1. CONCLUSIONS: We demonstrated gene-DASH diet score interaction effects on systolic BP in several loci. Studies with larger diverse populations are needed to validate our findings.

FORGEdb: a tool for identifying candidate functional variants and uncovering target genes and mechanisms for complex diseases.

The majority of disease-associated variants identified through genome-wide association studies are located outside of protein-coding regions. Prioritizing candidate regulatory variants and gene targets to identify potential biological mechanisms for further functional experiments can be challenging. To address this challenge, we developed FORGEdb ( https://forgedb.cancer.gov/ ; https://forge2.altiusinstitute.org/files/forgedb.html ; and https://doi.org/10.5281/zenodo.10067458 ), a standalone and web-based tool that integrates multiple datasets, delivering information on associated regulatory elements, transcription factor binding sites, and target genes for over 37 million variants. FORGEdb scores provide researchers with a quantitative assessment of the relative importance of each variant for targeted functional experiments.

Mendelian randomization analyses suggest a causal role for circulating GIP and IL-1RA levels in homeostatic model assessment-derived measures of ß-cell function and insulin sensitivity in Africans without type 2 diabetes.

BACKGROUND: In vitro and in vivo studies have shown that certain cytokines and hormones may play a role in the development and progression of type 2 diabetes (T2D). However, studies on their role in T2D in humans are scarce. We evaluated associations between 11 circulating cytokines and hormones with T2D among a population of sub-Saharan Africans and tested for causal relationships using Mendelian randomization (MR) analyses. METHODS: We used logistic regression analysis adjusted for age, sex, body mass index, and recruitment country to regress levels of 11 cytokines and hormones (adipsin, leptin, visfatin, PAI-1, GIP, GLP-1, ghrelin, resistin, IL-6, IL-10, IL-1RA) on T2D among Ghanaians, Nigerians, and Kenyans from the Africa America Diabetes Mellitus study including 2276 individuals with T2D and 2790 non-T2D individuals. Similar linear regression models were fitted with homeostatic modelling assessments of insulin sensitivity (HOMA-S) and ß-cell function (HOMA-B) as dependent variables among non-T2D individuals (n = 2790). We used 35 genetic variants previously associated with at least one of these 11 cytokines and hormones among non-T2D individuals as instrumental variables in univariable and multivariable MR analyses. Statistical significance was set at 0.0045 (0.05/11 cytokines and hormones). RESULTS: Circulating GIP and IL-1RA levels were associated with T2D. Nine of the 11 cytokines and hormones (exceptions GLP-1 and IL-6) were associated with HOMA-S, HOMA-B, or both among non-T2D individuals. Two-stage least squares MR analysis provided evidence for a causal effect of GIP and IL-RA on HOMA-S and HOMA-B in multivariable analyses (GIP ~ HOMA-S ß = - 0.67, P-value = 1.88 × 10-6 and HOMA-B ß = 0.59, P-value = 1.88 × 10-5; IL-1RA ~ HOMA-S ß = - 0.51, P-value = 8.49 × 10-5 and HOMA-B ß = 0.48, P-value = 5.71 × 10-4). IL-RA was partly mediated via BMI (30-34%), but GIP was not. Inverse variance weighted MR analysis provided evidence for a causal effect of adipsin on T2D (multivariable OR = 1.83, P-value = 9.79 × 10-6), though these associations were not consistent in all sensitivity analyses. CONCLUSIONS: The findings of this comprehensive MR analysis indicate that circulating GIP and IL-1RA levels are causal for reduced insulin sensitivity and increased ß-cell function. GIP's effect being independent of BMI suggests that circulating levels of GIP could be a promising early biomarker for T2D risk. Our MR analyses do not provide conclusive evidence for a causal role of other circulating cytokines in T2D among sub-Saharan Africans.

A statistical framework for powerful multi-trait rare variant analysis in large-scale whole-genome sequencing studies.

Large-scale whole-genome sequencing (WGS) studies have improved our understanding of the contributions of coding and noncoding rare variants to complex human traits. Leveraging association effect sizes across multiple traits in WGS rare variant association analysis can improve statistical power over single-trait analysis, and also detect pleiotropic genes and regions. Existing multi-trait methods have limited ability to perform rare variant analysis of large-scale WGS data. We propose MultiSTAAR, a statistical framework and computationally-scalable analytical pipeline for functionally-informed multi-trait rare variant analysis in large-scale WGS studies. MultiSTAAR accounts for relatedness, population structure and correlation among phenotypes by jointly analyzing multiple traits, and further empowers rare variant association analysis by incorporating multiple functional annotations. We applied MultiSTAAR to jointly analyze three lipid traits (low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides) in 61,861 multi-ethnic samples from the Trans-Omics for Precision Medicine (TOPMed) Program. We discovered new associations with lipid traits missed by single-trait analysis, including rare variants within an enhancer of NIPSNAP3A and an intergenic region on chromosome 1.

Genome-Wide Interaction Analysis with DASH Diet Score Identified Novel Loci for Systolic Blood Pressure.

Objective: We examined interactions between genotype and a Dietary Approaches to Stop Hypertension (DASH) diet score in relation to systolic blood pressure (SBP). Methods: We analyzed up to 9,420,585 biallelic imputed single nucleotide polymorphisms (SNPs) in up to 127,282 individuals of six population groups (91% of European population) from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (CHARGE; n=35,660) and UK Biobank (n=91,622) and performed European population-specific and cross-population meta-analyses. Results: We identified three loci in European-specific analyses and an additional four loci in cross-population analyses at P for interaction < 5e-8. We observed a consistent interaction between rs117878928 at 15q25.1 (minor allele frequency = 0.03) and the DASH diet score (P for interaction = 4e-8; P for heterogeneity = 0.35) in European population, where the interaction effect size was 0.42±0.09 mm Hg (P for interaction = 9.4e-7) and 0.20±0.06 mm Hg (P for interaction = 0.001) in CHARGE and the UK Biobank, respectively. The 1 Mb region surrounding rs117878928 was enriched with cis-expression quantitative trait loci (eQTL) variants (P = 4e-273) and cis-DNA methylation quantitative trait loci (mQTL) variants (P = 1e-300). While the closest gene for rs117878928 is MTHFS, the highest narrow sense heritability accounted by SNPs potentially interacting with the DASH diet score in this locus was for gene ST20 at 15q25.1. Conclusion: We demonstrated gene-DASH diet score interaction effects on SBP in several loci. Studies with larger diverse populations are needed to validate our findings.

Determinants of mosaic chromosomal alteration fitness.

Clonal hematopoiesis (CH) is characterized by the acquisition of a somatic mutation in a hematopoietic stem cell that results in a clonal expansion. These driver mutations can be single nucleotide variants in cancer driver genes or larger structural rearrangements called mosaic chromosomal alterations (mCAs). The factors that influence the variations in mCA fitness and ultimately result in different clonal expansion rates are not well-understood. We used the Passenger-Approximated Clonal Expansion Rate (PACER) method to estimate clonal expansion rate for 6,381 individuals in the NHLBI TOPMed cohort with gain, loss, and copy-neutral loss of heterozygosity mCAs. Our estimates of mCA fitness were correlated (R 2 = 0.49) with an alternative approach that estimated fitness of mCAs in the UK Biobank using a theoretical probability distribution. Individuals with lymphoid-associated mCAs had a significantly higher white blood cell count and faster clonal expansion rate. In a cross-sectional analysis, genome-wide association study of estimates of mCA expansion rate identified TCL1A , NRIP1 , and TERT locus variants as modulators of mCA clonal expansion rate.

Interplay Between Chronic Kidney Disease, Hypertension, and Stroke: Insights From a Multivariable Mendelian Randomization Analysis.

BACKGROUND AND OBJECTIVES: Chronic kidney disease (CKD) increases the risk of stroke, but the extent through which this association is mediated by hypertension is unknown. We leveraged large-scale genetic data to explore causal relationships between CKD, hypertension, and cerebrovascular disease phenotypes. METHODS: We used data from genome-wide association studies of European ancestry to identify genetic proxies for kidney function (CKD diagnosis, estimated glomerular filtration rate [eGFR], and urinary albumin-to-creatinine ratio [UACR]), systolic blood pressure (SBP), and cerebrovascular disease (ischemic stroke and its subtypes and intracerebral hemorrhage). We then conducted univariable, multivariable, and mediation Mendelian randomization (MR) analyses to investigate the effect of kidney function on stroke risk and the proportion of this effect mediated through hypertension. RESULTS: Univariable MR revealed associations between genetically determined lower eGFR and risk of all stroke (odds ratio [OR] per 1-log decrement in eGFR, 1.77; 95% CI 1.31-2.40; p < 0.001), ischemic stroke (OR 1.81; 95% CI 1.31-2.51; p < 0.001), and most strongly with large artery stroke (LAS) (OR 3.00; 95% CI 1.33-6.75; p = 0.008). These associations remained significant in the multivariable MR analysis, controlling for SBP (OR 1.98; 95% CI 1.39-2.82; p < 0.001 for all stroke; OR 2.16; 95% CI 1.48-3.17; p < 0.001 for ischemic stroke; OR 4.35; 95% CI 1.84-10.27; p = 0.001 for LAS), with only a small proportion of the total effects mediated by SBP (6.5% [0.7%-16.8%], 6.6% [0.8%-18.3%], and 7.2% [0.5%-24.8%], respectively). Total, direct and indirect effect estimates were similar across a number of sensitivity analyses (weighted median, MR-Egger regression). DISCUSSION: Our results demonstrate an independent causal effect of impaired kidney function, as assessed by decreased eGFR, on stroke risk, particularly LAS, even when controlled for SBP. Targeted prevention of kidney disease could lower atherosclerotic stroke risk independent of hypertension.