Trans-ethnic Fine Mapping Highlights Kidney-Function Genes Linked to Salt Sensitivity

Hierarchical joint analysis of marginal summary statistics-Part I: Multipopulation fine mapping and credible set construction

Recent advancement in genome-wide association studies (GWAS) comes from not only increasingly larger sample sizes but also the shift in focus towards underrepresented populations. Multipopulation GWAS increase power to detect novel risk variants and improve fine-mapping resolution by leveraging evidence and differences in linkage disequilibrium (LD) from diverse populations. Here, we expand upon our previous approach for single-population fine-mapping through Joint Analysis of Marginal SNP Effects (JAM) to a multipopulation analysis (mJAM). Under the assumption that true causal variants are common across studies, we implement a hierarchical model framework that conditions on multiple SNPs while explicitly incorporating the different LD structures across populations. The mJAM framework can be used to first select index variants using the mJAM likelihood with different feature selection approaches. In addition, we present a novel approach leveraging the ideas of mediation to construct credible sets for these index variants. Construction of such credible sets can be performed given any existing index variants. We illustrate the implementation of the mJAM likelihood through two implementations: mJAM-SuSiE (a Bayesian approach) and mJAM-Forward selection. Through simulation studies based on realistic effect sizes and levels of LD, we demonstrated that mJAM performs well for constructing concise credible sets that include the underlying causal variants. In real data examples taken from the most recent multipopulation prostate cancer GWAS, we showed several practical advantages of mJAM over other existing multipopulation methods.

Genetic association and transferability for urinary albumin-creatinine ratio as a marker of kidney disease in four Sub-Saharan African populations and non-continental individuals of African ancestry

Background

Genome-wide association studies (GWAS) have predominantly focused on populations of European and Asian ancestry, limiting our understanding of genetic factors influencing kidney disease in Sub-Saharan African (SSA) populations. This study presents the largest GWAS for urinary albumin-to-creatinine ratio (UACR) in SSA individuals, including 8,970 participants living in different African regions and an additional 9,705 non-resident individuals of African ancestry from the UK Biobank and African American cohorts.

Methods

Urine biomarkers and genotype data were obtained from two SSA cohorts (AWI-Gen and ARK), and two non-resident African-ancestry studies (UK Biobank and CKD-Gen Consortium). Association testing and meta-analyses were conducted, with subsequent fine-mapping, conditional analyses, and replication studies. Polygenic scores (PGS) were assessed for transferability across populations.

Results

Two genome-wide significant (P < 5 × 10-8) UACR-associated loci were identified, one in the BMP6 region on chromosome 6, in the meta-analysis of resident African individuals, and another in the HBB region on chromosome 11 in the meta-analysis of non-resident SSA individuals, as well as the combined meta-analysis of all studies. Replication of previous significant results confirmed associations in known UACR-associated regions, including THB53, GATM, and ARL15. PGS estimated using previous studies from European ancestry, African ancestry, and multi-ancestry cohorts exhibited limited transferability of PGS across populations, with less than 1% of observed variance explained.

Conclusion

This study contributes novel insights into the genetic architecture of kidney disease in SSA populations, emphasizing the need for conducting genetic research in diverse cohorts. The identified loci provide a foundation for future investigations into the genetic susceptibility to chronic kidney disease in underrepresented African populations Additionally, there is a need to develop integrated scores using multi-omics data and risk factors specific to the African context to improve the accuracy of predicting disease outcomes.

Genetic Association and Transferability for Urinary Albumin-Creatinine Ratio as a Marker of Kidney Disease in four Sub-Saharan African Populations and non-continental Individuals of African Ancestry

BACKGROUND

Genome-wide association studies (GWAS) have predominantly focused on populations of European and Asian ancestry, limiting our understanding of genetic factors influencing kidney disease in Sub-Saharan African (SSA) populations. This study presents the largest GWAS for urinary albumin-to-creatinine ratio (UACR) in SSA individuals, including 8,970 participants living in different African regions and an additional 9,705 non-resident individuals of African ancestry from the UK Biobank and African American cohorts.

METHODS

Urine biomarkers and genotype data were obtained from two SSA cohorts (AWI-Gen and ARK), and two non-resident African-ancestry studies (UK Biobank and CKD-Gen Consortium). Association testing and meta-analyses were conducted, with subsequent fine-mapping, conditional analyses, and replication studies. Polygenic scores (PGS) were assessed for transferability across populations.

RESULTS

Two genome-wide significant (P<5x10-8) UACR-associated loci were identified, one in the BMP6 region on chromosome 6, in the meta-analysis of resident African individuals, and another in the HBB region on chromosome 11 in the meta-analysis of non-resident SSA individuals, as well as the combined meta-analysis of all studies. Replication of previous significant results confirmed associations in known UACR-associated regions, including THB53, GATM, and ARL15. PGS estimated using previous studies from European ancestry, African ancestry, and multi-ancestry cohorts exhibited limited transferability of PGS across populations, with less than 1% of observed variance explained.

CONCLUSION

This study contributes novel insights into the genetic architecture of kidney disease in SSA populations, emphasizing the need for conducting genetic research in diverse cohorts. The identified loci provide a foundation for future investigations into the genetic susceptibility to chronic kidney disease in underrepresented African populations Additionally, there is a need to develop integrated scores using multi-omics data and risk factors specific to the African context to improve the accuracy of predicting disease outcomes.

METHODS

Urine biomarkers and genotype data were obtained from two SSA cohorts (AWI-Gen and ARK), and two non-resident African-ancestry studies (UK Biobank and CKD-Gen Consortium). Association testing and meta-analyses were conducted, with subsequent fine-mapping, conditional analyses, and replication studies. Polygenic scores (PGS) were assessed for transferability across populations.

RESULTS

Two genome-wide significant (P<5x10-8) UACR-associated loci were identified, one in the BMP6 region on chromosome 6, in the meta-analysis of resident African individuals, and another in the HBB region on chromosome 11 in the meta-analysis of non-resident SSA individuals, as well as the combined meta-analysis of all studies. Replication of previous significant results confirmed associations in known UACR-associated regions, including THB53, GATM, and ARL15. PGS estimated using previous studies from European ancestry, African ancestry, and multi-ancestry cohorts exhibited limited transferability of PGS across populations, with less than 1% of observed variance explained.

CONCLUSION

This study contributes novel insights into the genetic architecture of kidney function in SSA populations, emphasizing the need for conducting genetic research in diverse cohorts. The identified loci provide a foundation for future investigations into the genetic susceptibility to chronic kidney disease in underrepresented African populations.

Multivariate canonical correlation analysis identifies additional genetic variants for chronic kidney disease

Chronic kidney diseases (CKD) have genetic associations with kidney function. Univariate genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with estimated glomerular filtration rate (eGFR) and blood urea nitrogen (BUN), two complementary kidney function markers. However, it is unknown whether additional SNPs for kidney function can be identified by multivariate statistical analysis. To address this, we applied canonical correlation analysis (CCA), a multivariate method, to two individual-level CKD genotype datasets, and metaCCA to two published GWAS summary statistics datasets. We identified SNPs previously associated with kidney function by published univariate GWASs with high replication rates, validating the metaCCA method. We then extended discovery and identified previously unreported lead SNPs for both kidney function markers, jointly. These showed expression quantitative trait loci (eQTL) colocalisation with genes having significant differential expression between CKD and healthy individuals. Several of these identified lead missense SNPs were predicted to have a functional impact, including in SLC14A2. We also identified previously unreported lead SNPs that showed significant correlation with both kidney function markers, jointly, in the European ancestry CKDGen, National Unified Renal Translational Research Enterprise (NURTuRE)-CKD and Salford Kidney Study (SKS) datasets. Of these, rs3094060 colocalised with FLOT1 gene expression and was significantly more common in CKD cases in both NURTURE-CKD and SKS, than in the general population. Overall, by using multivariate analysis by CCA, we identified additional SNPs and genes for both kidney function and CKD, that can be prioritised for further CKD analyses.

Genome-wide study investigating effector genes and polygenic prediction for kidney function in persons with ancestry from Africa and the Americas

Chronic kidney disease is a leading cause of death and disability globally and impacts individuals of African ancestry (AFR) or with ancestry in the Americas (AMS) who are under-represented in genome-wide association studies (GWASs) of kidney function. To address this bias, we conducted a large meta-analysis of GWASs of estimated glomerular filtration rate (eGFR) in 145,732 AFR and AMS individuals. We identified 41 loci at genome-wide significance (p < 5 × 10-8), of which two have not been previously reported in any ancestry group. We integrated fine-mapped loci with epigenomic and transcriptomic resources to highlight potential effector genes relevant to kidney physiology and disease, and reveal key regulatory elements and pathways involved in renal function and development. We demonstrate the varying but increased predictive power offered by a multi-ancestry polygenic score for eGFR and highlight the importance of population diversity in GWASs and multi-omics resources to enhance opportunities for clinical translation for all.

Evaluation of ceftriaxone pharmacokinetics in hospitalized Egyptian pediatric patients

This study aimed to evaluate ceftriaxone pharmacokinetics that affects the achievement of targets in the treatment of critically ill children (meningitis, pneumonia, urinary tract infection, peritonitis, and infective endocarditis( who were admitted to Zagazig University Pediatric hospital in Egypt to monitor for the drug adverse effects.Blood samples were obtained from 24 hospitalized pediatric patients (ages ranging from 2.5 months to 12 years) after administering the calculated dose of ceftriaxone via intravenous bolus route. Then, ceftriaxone plasma concentrations were measured using a validated HPLC method with ultraviolet detection. The pharmacokinetic analysis was conducted using Phoenix Winnonlin Program® software.Data for total and free ceftriaxone best fitted on a one-compartment model with the first-order elimination process. Clearance of ceftriaxone is reduced for patients with reduced kidney function and increased with those with augmented renal clearance. The volume of distribution and the free fraction are increased in these patients, especially those with hypoalbuminemia with a shorter half-life time were detected. A slight increase in total bilirubin and liver enzymes has been observed after treatment with ceftriaxone in these patients.   Conclusion: In most critically ill pediatric patients, the current ceftriaxone treatment regimen (50 to 100 mg/kg) offers adequate pathogenic coverage. The clearance of free ceftriaxone in all patients correlates well with their renal function (eGFR), with r2 = 0.7252. During therapy with ceftriaxone at all doses ranging from 50 to 100 mg/kg, a rise in total bilirubin was observed in these patients. Moreover, liver enzymes (ALT and AST) increased moderately (p 0.0001). So, it is recommended to monitor total bilirubin and liver enzymes during the treatment with ceftriaxone, especially for a long duration (more than 5 days) or use another agent in patients with high baseline values. What is Known: • The dosing regimen of ceftriaxone (50 to 100 mg/kg) provided optimum therapeutic outcomes. • Some studies show data for total and free Ceftriaxone best fitted on a one-compartment model while other studies show data for total and free Ceftriaxone best fitted on a two-compartment model. What is New: • Up to my knowledge this is the first study ,considering individual pharmacokinetic analysis, conducted on hospitalized Egyptian pediatric population most of them with reduced kidney function with ages ranging from 2.5 months to 12 years. Data for total and free Ceftriaxone best fitted on a one-compartment model with linear clearance of the free ceftriaxone. • In all patients, total bilirubin and liver function tests were mildly increased, making them at risk for cholestasis or ceftriaxone-induced cholestatic hepatitis.

Genome-Wide Association Study of CKD Progression

SIGNIFICANCE STATEMENT

Rapid progression of CKD is associated with poor clinical outcomes. Most previous studies looking for genetic factors associated with low eGFR have used cross-sectional data. The authors conducted a meta-analysis of genome-wide association studies of eGFR decline among 116,870 participants with CKD, focusing on longitudinal data. They identified three loci (two of them novel) associated with longitudinal eGFR decline. In addition to the known UMOD/PDILT locus, variants within BICC1 were associated with significant differences in longitudinal eGFR slope. Variants within HEATR4 also were associated with differences in eGFR decline, but only among Black/African American individuals without diabetes. These findings help characterize molecular mechanisms of eGFR decline in CKD and may inform new therapeutic approaches for progressive kidney disease.

BACKGROUND

Rapid progression of CKD is associated with poor clinical outcomes. Despite extensive study of the genetics of cross-sectional eGFR, only a few loci associated with eGFR decline over time have been identified.

METHODS

We performed a meta-analysis of genome-wide association studies of eGFR decline among 116,870 participants with CKD-defined by two outpatient eGFR measurements of <60 ml/min per 1.73 m 2 , obtained 90-365 days apart-from the Million Veteran Program and Vanderbilt University Medical Center's DNA biobank. The primary outcome was the annualized relative slope in outpatient eGFR. Analyses were stratified by ethnicity and diabetes status and meta-analyzed thereafter.

RESULTS

In cross-ancestry meta-analysis, the strongest association was rs77924615, near UMOD / PDILT ; each copy of the G allele was associated with a 0.30%/yr faster eGFR decline ( P = 4.9×10 -27 ). We also observed an association within BICC1 (rs11592748), where every additional minor allele was associated with a 0.13%/yr slower eGFR decline ( P = 5.6×10 -9 ). Among participants without diabetes, the strongest association was the UMOD/PDILT variant rs36060036, associated with a 0.27%/yr faster eGFR decline per copy of the C allele ( P = 1.9×10 -17 ). Among Black participants, a significantly faster eGFR decline was associated with variant rs16996674 near APOL1 (R 2 =0.29 with the G1 high-risk genotype); among Black participants with diabetes, lead variant rs11624911 near HEATR4 also was associated with a significantly faster eGFR decline. We also nominally replicated loci with known associations with eGFR decline, near PRKAG2, FGF5, and C15ORF54.

CONCLUSIONS

Three loci were significantly associated with longitudinal eGFR change at genome-wide significance. These findings help characterize molecular mechanisms of eGFR decline and may contribute to the development of new therapeutic approaches for progressive CKD.

Albuminuria-Related Genetic Biomarkers: Replication and Predictive Evaluation in Individuals with and without Diabetes from the UK Biobank

Increased albuminuria indicates underlying glomerular pathology and is associated with worse renal disease outcomes, especially in diabetic kidney disease. Many single nucleotide polymorphisms (SNPs), associated with albuminuria, could be potentially useful to construct polygenic risk scores (PRSs) for kidney disease. We investigated the diagnostic accuracy of SNPs, previously associated with albuminuria-related traits, on albuminuria and renal injury in the UK Biobank population, with a particular interest in diabetes. Multivariable logistic regression was used to evaluate the influence of 91 SNPs on urine albumin-to-creatinine ratio (UACR)-related traits and kidney damage (any pathology indicating renal injury), stratifying by diabetes. Weighted PRSs for microalbuminuria and UACR from previous studies were used to calculate the area under the receiver operating characteristic curve (AUROC). CUBN-rs1801239 and DDR1-rs116772905 were associated with all the UACR-derived phenotypes, in both the overall and non-diabetic cohorts, but not with kidney damage. Several SNPs demonstrated different effects in individuals with diabetes compared to those without. SNPs did not improve the AUROC over currently used clinical variables. Many SNPs are associated with UACR or renal injury, suggesting a role in kidney dysfunction, dependent on the presence of diabetes in some cases. However, individual SNPs or PRSs did not improve the diagnostic accuracy for albuminuria or renal injury compared to standard clinical variables.

Genome-wide association study of the risk of chronic kidney disease and kidney-related traits in the Japanese population: J-Kidney-Biobank

Chronic kidney disease (CKD) is a syndrome characterized by a gradual loss of kidney function with decreased estimated glomerular filtration rate (eGFR), which may be accompanied by an increase in the urine albumin-to-creatinine ratio (UACR). Although trans-ethnic genome-wide association studies (GWASs) have been conducted for kidney-related traits, there have been few analyses in the Japanese population, especially for the UACR trait. In this study, we conducted a GWAS to identify loci related to multiple kidney-related traits in Japanese individuals. First, to detect loci associated with CKD, eGFR, and UACR, we performed separate GWASs with the following two datasets: 475 cases of CKD diagnosed at seven university hospitals and 3471 healthy subjects (dataset 1) and 3664 cases of CKD-suspected individuals with eGFR <60 ml/min/1.73 m² or urinary protein ≥ 1+ and 5952 healthy subjects (dataset 2). Second, we performed a meta-analysis between these two datasets and detected the following associated loci: 10 loci for CKD, 9 loci for eGFR, and 22 loci for UACR. Among the loci detected, 22 have never been reported previously. Half of the significant loci for CKD were shared with those for eGFR, whereas most of the loci associated with UACR were different from those associated with CKD or eGFR. The GWAS of the Japanese population identified novel genetic components that were not previously detected. The results also suggest that the group primarily characterized by increased UACR possessed genetically different features from the group characterized by decreased eGFR.

Proteome-wide Mendelian randomization in global biobank meta-analysis reveals multi-ancestry drug targets for common diseases

Proteome-wide Mendelian randomization (MR) shows value in prioritizing drug targets in Europeans but with limited evidence in other ancestries. Here, we present a multi-ancestry proteome-wide MR analysis based on cross-population data from the Global Biobank Meta-analysis Initiative (GBMI). We estimated the putative causal effects of 1,545 proteins on eight diseases in African (32,658) and European (1,219,993) ancestries and identified 45 and 7 protein-disease pairs with MR and genetic colocalization evidence in the two ancestries, respectively. A multi-ancestry MR comparison identified two protein-disease pairs with MR evidence in both ancestries and seven pairs with specific effects in the two ancestries separately. Integrating these MR signals with clinical trial evidence, we prioritized 16 pairs for investigation in future drug trials. Our results highlight the value of proteome-wide MR in informing the generalizability of drug targets for disease prevention across ancestries and illustrate the value of meta-analysis of biobanks in drug development.

Towards a global view of multiple sclerosis genetics

Multiple sclerosis (MS) is a neuroimmunological disorder of the CNS with a strong heritable component. The genetic architecture of MS susceptibility is well understood in populations of European ancestry. However, the extent to which this architecture explains MS susceptibility in populations of non-European ancestry remains unclear. In this Perspective article, we outline the scientific arguments for studying MS genetics in ancestrally diverse populations. We argue that this approach is likely to yield insights that could benefit individuals with MS from all ancestral groups. We explore the logistical and theoretical challenges that have held back this field to date and conclude that, despite these challenges, inclusion of participants of non-European ancestry in MS genetics studies will ultimately be of value to all patients with MS worldwide.

Phenome-Wide Association Study of UMOD Gene Variants and Differential Associations With Clinical Outcomes Across Populations in the Million Veteran Program a Multiethnic Biobank

Introduction

Common variants in the UMOD gene are considered an evolutionary adaptation against urinary tract infections (UTIs) and have been implicated in kidney stone formation, chronic kidney disease (CKD), and hypertension. However, differences in UMOD variant-phenotype associations across population groups are unclear.

Methods

We tested associations between UMOD/PDILT variants and up to 1528 clinical diagnosis codes mapped to phenotype groups in the Million Veteran Program (MVP), using published phenome-wide association study (PheWAS) methodology. Associations were tested using logistic regression adjusted for age, sex, and 10 principal components of ancestry. Bonferroni correction for multiple comparisons was applied.

Results

Among 648,593 veterans, mean (SD) age was 62 (14) years; 9% were female, 19% Black, and 8% Hispanic. In White patients, the rs4293393 UMOD risk variant associated with increased uromodulin was associated with increased odds of CKD (odds ratio [OR]: 1.22, 95% CI: 1.20-1.24, P = 5.90 × 10-111), end-stage kidney disease (OR: 1.17, 95% CI: 1.11-1.24, P = 2.40 × 10-09), and hypertension (OR: 1.03, 95% CI: 1.05-1.05, P = 2.11 × 10-06) and significantly lower odds of UTIs (OR: 0.94, 95% CI: 0.92-0.96, P = 1.21 × 10-10) and kidney calculus (OR: 0.85, 95% CI: 0.83-0.86, P = 4.27 × 10-69). Similar findings were observed across UMOD/PDILT variants. The rs77924615 PDILT variant had stronger associations with acute cystitis in White female (OR: 0.73, 95% CI: 0.59-0.91, P = 4.98 × 10-03) versus male (OR: 0.99, 95% CI: 0.89-1.11, P = 8.80 × 10-01) (P interaction = 0.01) patients. In Black patients, the rs77924615 PDILT variant was significantly associated with pyelonephritis (OR: 0.65, 95% CI: 0.54-0.79, P = 1.05 × 10-05), whereas associations with UMOD promoter variants were attenuated.

Conclusion

Robust associations were observed between UMOD/PDILT variants linked with increased uromodulin expression and lower odds of UTIs and calculus and increased odds of CKD and hypertension. However, these associations varied significantly across ancestry groups and sex.

RGS14 regulates PTH- and FGF23-sensitive NPT2A-mediated renal phosphate uptake via binding to the NHERF1 scaffolding protein

Phosphate homeostasis, mediated by dietary intake, renal absorption, and bone deposition, is incompletely understood because of the uncharacterized roles of numerous implicated protein factors. Here, we identified a novel role for one such element, regulator of G protein signaling 14 (RGS14), suggested by genome-wide association studies to associate with dysregulated Pi levels. We show that human RGS14 possesses a carboxy-terminal PDZ ligand required for sodium phosphate cotransporter 2a (NPT2A) and sodium hydrogen exchanger regulatory factor-1 (NHERF1)-mediated renal Pi transport. In addition, we found using isotope uptake measurements combined with bioluminescence resonance energy transfer assays, siRNA knockdown, pull-down and overlay assays, and molecular modeling that secreted proteins parathyroid hormone (PTH) and fibroblast growth factor 23 inhibited Pi uptake by inducing dissociation of the NPT2A-NHERF1 complex. PTH failed to affect Pi transport in cells expressing RGS14, suggesting that it suppresses hormone-sensitive but not basal Pi uptake. Interestingly, RGS14 did not affect PTH-directed G protein activation or cAMP formation, implying a postreceptor site of action. Further pull-down experiments and direct binding assays indicated that NPT2A and RGS14 bind distinct PDZ domains on NHERF1. We showed that RGS14 expression in human renal proximal tubule epithelial cells blocked the effects of PTH and fibroblast growth factor 23 and stabilized the NPT2A-NHERF1 complex. In contrast, RGS14 genetic variants bearing mutations in the PDZ ligand disrupted RGS14 binding to NHERF1 and subsequent PTH-sensitive Pi transport. In conclusion, these findings identify RGS14 as a novel regulator of hormone-sensitive Pi transport. The results suggest that changes in RGS14 function or abundance may contribute to the hormone resistance and hyperphosphatemia observed in kidney diseases.

Contribution of the innate and adaptive immune systems to aortic dilation in murine mucopolysaccharidosis type I

BACKGROUND

Adult immunocompetent male C57Bl/6 mucopolysaccharidosis, type I (MPSI) mice develop aortic insufficiency (AI), dilated ascending aortas and decreased cardiac function, findings not observed in immune incompetent adult male NSG MPSI mice. We sought to determine why.

METHODS

Cardiac ultrasound measurements of ascending aorta and left ventricular dimensions and Doppler interrogation for AI were performed in 6-month-old male B6 MPSI (N = 12), WT (N = 6), NSG MPSI (N = 8), NSG (N = 6) mice. Urinary glycosaminoglycans, RNA sequencing with quantitative PCR were performed and aortic pathology assessed by routine and immunohistochemical staining on subsets of murine aortas.

RESULTS

Ascending aortic diameters were significantly greater, left ventricular function significantly decreased, and AI significantly more frequent in B6 MPSI mice compared to NSG MPSI mice (p < 0.0001, p = 0.008 and p = 0.02, respectively); NSG and B6 WT mice showed no changes. Urinary glycosaminoglycans were significantly greater in B6 and NSG MPSI mice and both were significantly elevated compared to WT controls (p = 0.003 and p < 0.0001, respectively). By RNA sequencing, all 11 components of the inflammasome pathway were upregulated in B6 MUT, but only Aim2 and Ctsb in NSG MUT mice and none in WT controls. Both B6 and NSG MUT mice demonstrated variably-severe intramural inflammation, vacuolated cells, elastin fragmentation and disarray, and intense glycosaminoglycans on histological staining. B6 MPSI mice demonstrated numerous medial MAC2+ macrophages and adventitial CD3+ T-cells while MAC2+ macrophages were sparse and CD3+ T-cells absent in NSG MPSI mice.

CONCLUSIONS

Aortic dilation, AI and decreased cardiac function occur in immunocompetent B6 MPSI male mice but not in immune incompetent NSG MPSI mice, unrelated to GAG excretion, upregulation of Ctsb, or routine histologic appearance. Upregulation of all components of the inflammasome pathway in B6 MUT, but not NSG MUT mice, and abundant medial MAC2 and adventitial CD3 infiltrates in B6, but not NSG, MPSI aortas differentiated the two strains. These results suggest that the innate and adaptive immune systems play a role in these cardiac findings which may be relevant to human MPSI.

Genome-Wide Admixture Mapping of Estimated Glomerular Filtration Rate and Chronic Kidney Disease Identifies European and African Ancestry-of-Origin Loci in Hispanic and Latino Individuals in the United States

BACKGROUND

Admixture mapping is a powerful approach for gene mapping of complex traits that leverages the diverse genetic ancestry in populations with recent admixture, such as Hispanic or Latino individuals in the United States. These individuals have an increased risk of CKD.

METHODS

We performed genome-wide admixture mapping for both CKD and eGFR in a sample of 12,601 participants from the Hispanic Community Health Study/Study of Latinos, with validation in a sample of 8191 Black participants from the Women's Health Initiative (WHI). We also compared the findings with those from a conventional genome-wide association study.

RESULTS

Three novel ancestry-of-origin loci were identified on chromosomes 2, 14, and 15 for CKD and eGFR. The chromosome 2 locus comprises two European ancestry regions encompassing the FSHR and NRXN1 genes, with European ancestry at this locus associated with increased CKD risk. The chromosome 14 locus, found within the DLK1-DIO3 imprinted domain, was associated with lower eGFR and driven by European ancestry. The eGFR-associated locus on chromosome 15 included intronic variants of RYR3 and was within an African-specific genomic region associated with higher eGFR. The genome-wide association study failed to identify significant associations in these regions. We validated the chromosome 14 and 15 loci associated with eGFR in the WHI Black participants.

CONCLUSIONS

This study provides evidence of shared ancestry-specific genomic regions influencing eGFR in Hispanic or Latino individuals and Black individuals and illustrates the potential for leveraging genetic ancestry in recently admixed populations for the discovery of novel candidate loci for kidney phenotypes.

Discovery of novel eGFR-associated multiple independent signals using a quasi-adaptive method

A decreased estimated glomerular filtration rate (eGFR) leading to chronic kidney disease is a significant public health problem. Kidney function is a heritable trait, and recent application of genome-wide association studies (GWAS) successfully identified multiple eGFR-associated genetic loci. To increase statistical power for detecting independent associations in GWAS loci, we improved our recently developed quasi-adaptive method estimating SNP-specific alpha levels for the conditional analysis, and applied it to the GWAS meta-analysis results of eGFR among 783,978 European-ancestry individuals. Among known eGFR loci, we revealed 19 new independent association signals that were subsequently replicated in the United Kingdom Biobank (n = 408,608). These associations have remained undetected by conditional analysis using the established conservative genome-wide significance level of 5 × 10^-8. Functional characterization of known index SNPs and novel independent signals using colocalization of conditional eGFR association results and gene expression in cis across 51 human tissues identified two potentially causal genes across kidney tissues: TSPAN33 and TFDP2, and three candidate genes across other tissues: SLC22A2, LRP2, and CDKN1C. These colocalizations were not identified in the original GWAS. By applying our improved quasi-adaptive method, we successfully identified additional genetic variants associated with eGFR. Considering these signals in colocalization analyses can increase the precision of revealing potentially functional genes of GWAS loci.

RGS14 Regulation of Post-Synaptic Signaling and Spine Plasticity in Brain

The regulator of G-protein signaling 14 (RGS14) is a multifunctional signaling protein that regulates post synaptic plasticity in neurons. RGS14 is expressed in the brain regions essential for learning, memory, emotion, and stimulus-induced behaviors, including the basal ganglia, limbic system, and cortex. Behaviorally, RGS14 regulates spatial and object memory, female-specific responses to cued fear conditioning, and environmental- and psychostimulant-induced locomotion. At the cellular level, RGS14 acts as a scaffolding protein that integrates G protein, Ras/ERK, and calcium/calmodulin signaling pathways essential for spine plasticity and cell signaling, allowing RGS14 to naturally suppress long-term potentiation (LTP) and structural plasticity in hippocampal area CA2 pyramidal cells. Recent proteomics findings indicate that RGS14 also engages the actomyosin system in the brain, perhaps to impact spine morphogenesis. Of note, RGS14 is also a nucleocytoplasmic shuttling protein, where its role in the nucleus remains uncertain. Balanced nuclear import/export and dendritic spine localization are likely essential for RGS14 neuronal functions as a regulator of synaptic plasticity. Supporting this idea, human genetic variants disrupting RGS14 localization also disrupt RGS14’s effects on plasticity. This review will focus on the known and unexplored roles of RGS14 in cell signaling, physiology, disease and behavior.

Mendelian pathway analysis of laboratory traits reveals distinct roles for ciliary subcompartments in common disease pathogenesis

Rare monogenic disorders of the primary cilium, termed ciliopathies, are characterized by extreme presentations of otherwise common diseases, such as diabetes, hepatic fibrosis, and kidney failure. However, despite a recent revolution in our understanding of the cilium's role in rare disease pathogenesis, the organelle's contribution to common disease remains largely unknown. Hypothesizing that common genetic variants within Mendelian ciliopathy genes might contribute to common complex diseases pathogenesis, we performed association studies of 16,874 common genetic variants across 122 ciliary genes with 12 quantitative laboratory traits characteristic of ciliopathy syndromes in 452,593 individuals in the UK Biobank. We incorporated tissue-specific gene expression analysis, expression quantitative trait loci, and Mendelian disease phenotype information into our analysis and replicated our findings in meta-analysis. 101 statistically significant associations were identified across 42 of the 122 examined ciliary genes (including eight novel replicating associations). These ciliary genes were widely expressed in tissues relevant to the phenotypes being studied, and eQTL analysis revealed strong evidence for correlation between ciliary gene expression levels and laboratory traits. Perhaps most interestingly, our analysis identified different ciliary subcompartments as being specifically associated with distinct sets of phenotypes. Taken together, our data demonstrate the utility of a Mendelian pathway-based approach to genomic association studies, challenge the widely held belief that the cilium is an organelle important mainly in development and in rare syndromic disease pathogenesis, and provide a framework for the continued integration of common and rare disease genetics to provide insight into the pathophysiology of human diseases of immense public health burden.

Calcineurin and Systemic Lupus Erythematosus: The Rationale for Using Calcineurin Inhibitors in the Treatment of Lupus Nephritis

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease with a broad spectrum of clinical presentations that can affect almost all organ systems. Lupus nephritis (LN) is a severe complication that affects approximately half of the systemic erythematosus lupus (SLE) patients, which significantly increases the morbidity and the mortality risk. LN is characterized by the accumulation of immune complexes, ultimately leading to renal failure. Aberrant activation of T cells plays a critical role in the pathogenesis of both SLE and LN and is involved in the production of inflammatory cytokines, the recruitment of inflammatory cells to the affected tissues and the co-stimulation of B cells. Calcineurin is a serine-threonine phosphatase that, as a consequence of the T cell hyperactivation, induces the production of inflammatory mediators. Moreover, calcineurin is also involved in the alterations of the podocyte phenotype, which contribute to proteinuria and kidney damage observed in LN patients. Therefore, calcineurin inhibitors have been postulated as a potential treatment strategy in LN, since they reduce T cell activation and promote podocyte cytoskeleton stabilization, both being key aspects in the development of LN. Here, we review the role of calcineurin in SLE and the latest findings about calcineurin inhibitors and their mechanisms of action in the treatment of LN.

Large genome-wide association study identifies three novel risk variants for restless legs syndrome

Restless legs syndrome (RLS) is a common neurological sensorimotor disorder often described as an unpleasant sensation associated with an urge to move the legs. Here we report findings from a meta-analysis of genome-wide association studies of RLS including 480,982 Caucasians (cases = 10,257) and a follow up sample of 24,977 (cases = 6,651). We confirm 19 of the 20 previously reported RLS sequence variants at 19 loci and report three novel RLS associations; rs112716420-G (OR = 1.25, P = 1.5 × 10-18), rs10068599-T (OR = 1.09, P = 6.9 × 10-10) and rs10769894-A (OR = 0.90, P = 9.4 × 10-14). At four of the 22 RLS loci, cis-eQTL analysis indicates a causal impact on gene expression. Through polygenic risk score for RLS we extended prior epidemiological findings implicating obesity, smoking and high alcohol intake as risk factors for RLS. To improve our understanding, with the purpose of seeking better treatments, more genetics studies yielding deeper insights into the disease biology are needed.

Genome-Wide Association Studies of CKD and Related Traits

The past few years have seen major advances in genome-wide association studies (GWAS) of CKD and kidney function-related traits in several areas: increases in sample size from >100,000 to >1 million, enabling the discovery of >250 associated genetic loci that are highly reproducible; the inclusion of participants not only of European but also of non-European ancestries; and the use of advanced computational methods to integrate additional genomic and other unbiased, high-dimensional data to characterize the underlying genetic architecture and prioritize potentially causal genes and variants. Together with other large-scale biobank and genetic association studies of complex traits, these GWAS of kidney function-related traits have also provided novel insight into the relationship of kidney function to other diseases with respect to their genetic associations, genetic correlation, and directional relationships. A number of studies also included functional experiments using model organisms or cell lines to validate prioritized potentially causal genes and/or variants. In this review article, we will summarize these recent GWAS of CKD and kidney function-related traits, explain approaches for downstream characterization of associated genetic loci and the value of such computational follow-up analyses, and discuss related challenges along with potential solutions to ultimately enable improved treatment and prevention of kidney diseases through genetics.

Regulator of G protein signaling 10: Structure, expression and functions in cellular physiology and diseases

Regulator of G protein signaling 10 (RGS10) belongs to the superfamily of RGS proteins, defined by the presence of a conserved RGS domain that canonically binds and deactivates heterotrimeric G-proteins. RGS proteins act as GTPase activating proteins (GAPs), which accelerate GTP hydrolysis on the G-protein α subunits and result in termination of signaling pathways downstream of G protein-coupled receptors. RGS10 is the smallest protein of the D/R12 subfamily and selectively interacts with Gαi proteins. It is widely expressed in many cells and tissues, with the highest expression found in the brain and immune cells. RGS10 expression is transcriptionally regulated via epigenetic mechanisms. Although RGS10 lacks multiple of the defined regulatory domains found in other RGS proteins, RGS10 contains post-translational modification sites regulating its expression, localization, and function. Additionally, RGS10 is a critical protein in the regulation of physiological processes in multiple cells, where dysregulation of its expression has been implicated in various diseases including Parkinson's disease, multiple sclerosis, osteopetrosis, chemoresistant ovarian cancer and cardiac hypertrophy. This review summarizes RGS10 features and its regulatory mechanisms, and discusses the known functions of RGS10 in cellular physiology and pathogenesis of several diseases.

Identification of potential key genes and key pathways related to clear cell renal cell carcinoma through bioinformatics analysis

Clear cell renal cell carcinoma (ccRCC) is a common malignancy of the genitourinary system and is associated with high mortality rates. However, the molecular mechanism of ccRCC pathogenesis is still unclear, which translates to few effective diagnostic and prognostic biomarkers. In this study, we conducted a bioinformatics analysis on three Gene Expression Omnibus datasets and identified 437 differentially expressed genes (DEGs) related to ccRCC development and prognosis, of which 311 and 126 genes are respectively down-regulated and up-regulated. The protein-protein interaction network of these DEGs consists of 395 nodes and 1872 interactions and 2 prominent modules. The Staphylococcus aureus infection and complement and coagulation cascades are significantly enriched in module 1 and are likely involved in ccRCC progression. Forty-two hub genes were screened, of which von Willebrand factor, TIMP metallopeptidase inhibitor 1, plasminogen, formimidoyltransferase cyclodeaminase, solute carrier family 34 member 1, hydroxyacid oxidase 2, alanine-glyoxylate aminotransferase 2, phosphoenolpyruvate carboxykinase 1, and 3-hydroxy-3-methylglutaryl-CoA synthase 2 are possibly related to the prognosis of ccRCC. The differential expression of all nine genes was confirmed by quantitative real-time polymerase chain reaction analysis of the ccRCC and normal renal tissues. These key genes are potential biomarkers for the diagnosis and prognosis of ccRCC and warrant further investigation.

Phosphate Transport in Epithelial and Nonepithelial Tissue

Phosphate is an essential nutrient for life and is a critical component of bone formation, a major signaling molecule, and structural component of cell walls. Phosphate is also a component of high-energy compounds (i.e., AMP, ADP, and ATP) and essential for nucleic acid helical structure (i.e., RNA and DNA). Phosphate plays a central role in the process of mineralization, normal serum levels being associated with appropriate bone mineralization, while high and low serum levels are associated with soft tissue calcification. The serum concentration of phosphate and the total body content of phosphate are highly regulated, a process that is accomplished by the coordinated effort of two families of sodium-dependent transporter proteins. The three isoforms of the SLC34 family (SLC34A1-A3) show very restricted tissue expression and regulate intestinal absorption and renal excretion of phosphate. SLC34A2 also regulates the phosphate concentration in multiple lumen fluids including milk, saliva, pancreatic fluid, and surfactant. Both isoforms of the SLC20 family exhibit ubiquitous expression (with some variation as to which one or both are expressed), are regulated by ambient phosphate, and likely serve the phosphate needs of the individual cell. These proteins exhibit similarities to phosphate transporters in nonmammalian organisms. The proteins are nonredundant as mutations in each yield unique clinical presentations. Further research is essential to understand the function, regulation, and coordination of the various phosphate transporters, both the ones described in this review and the phosphate transporters involved in intracellular transport.

Sodium-hydrogen exchanger regulatory factor-1 (NHERF1) confers salt sensitivity in both male and female models of hypertension in aging

Hypertension is a risk factor for premature death and roughly 50% of hypertensive patients are salt-sensitive. The incidence of salt-sensitive hypertension increases with age. However, the mechanisms of salt-sensitive hypertension are not well understood. We had demonstrated decreased renal sodium‑hydrogen exchanger regulatory factor 1 (NHERF1) expression in old salt-resistant F344 rats. Based on those studies we hypothesized that NHERF1 expression is required for the development of some forms of salt-sensitive hypertension. To address this hypothesis, we measured blood pressure in NHERF1 expressing salt-sensitive 4-mo and 24-mo-old male and female Fischer Brown Norway (FBN) rats male and female 18-mo-old NHERF1 knock-out (NHERF1-/-) mice and wild-type (WT) littermates on C57BL/6J background after feeding high salt (8% NaCl) diet for 7 days. Our data demonstrate that 8% salt diet increased blood pressure in both male and female 24-mo-old FBN rats but not in 4-mo-old FBN rats and in 18-mo-old male and female WT mice but not in NHERF1-/- mice. Renal dopamine 1 receptor (D1R) expression was decreased in 24-mo-old rats, compared with 4-mo-old FBN rats. However, sodium chloride cotransporter (NCC) expression increased in 24-mo-old FBN rats. In FBN rats, age had no effect on NaK ATPase α1 and NKCC2 expression. By contrast, high salt diet increased the renal expressions of NKCC2, and NCC in 24-mo-old FBN rats. High salt diet also increased NKCC2 and NCC expression in WT mice but not NHERF1-/- mice. Our data suggest that renal NHERF1 expression confers salt sensitivity with aging, associated with increased expression of sodium transporters.

Brain scans from 21,297 individuals reveal the genetic architecture of hippocampal subfield volumes

The hippocampus is a heterogeneous structure, comprising histologically distinguishable subfields. These subfields are differentially involved in memory consolidation, spatial navigation and pattern separation, complex functions often impaired in individuals with brain disorders characterized by reduced hippocampal volume, including Alzheimer's disease (AD) and schizophrenia. Given the structural and functional heterogeneity of the hippocampal formation, we sought to characterize the subfields' genetic architecture. T1-weighted brain scans (n = 21,297, 16 cohorts) were processed with the hippocampal subfields algorithm in FreeSurfer v6.0. We ran a genome-wide association analysis on each subfield, co-varying for whole hippocampal volume. We further calculated the single-nucleotide polymorphism (SNP)-based heritability of 12 subfields, as well as their genetic correlation with each other, with other structural brain features and with AD and schizophrenia. All outcome measures were corrected for age, sex and intracranial volume. We found 15 unique genome-wide significant loci across six subfields, of which eight had not been previously linked to the hippocampus. Top SNPs were mapped to genes associated with neuronal differentiation, locomotor behaviour, schizophrenia and AD. The volumes of all the subfields were estimated to be heritable (h2 from 0.14 to 0.27, all p < 1 × 10^-16) and clustered together based on their genetic correlations compared with other structural brain features. There was also evidence of genetic overlap of subicular subfield volumes with schizophrenia. We conclude that hippocampal subfields have partly distinct genetic determinants associated with specific biological processes and traits. Taking into account this specificity may increase our understanding of hippocampal neurobiology and associated pathologies.

A catalog of genetic loci associated with kidney function from analyses of a million individuals

Chronic kidney disease (CKD) is responsible for a public health burden with multi-systemic complications. Through trans-ancestry meta-analysis of genome-wide association studies of estimated glomerular filtration rate (eGFR) and independent replication (n = 1,046,070), we identified 264 associated loci (166 new). Of these, 147 were likely to be relevant for kidney function on the basis of associations with the alternative kidney function marker blood urea nitrogen (n = 416,178). Pathway and enrichment analyses, including mouse models with renal phenotypes, support the kidney as the main target organ. A genetic risk score for lower eGFR was associated with clinically diagnosed CKD in 452,264 independent individuals. Colocalization analyses of associations with eGFR among 783,978 European-ancestry individuals and gene expression across 46 human tissues, including tubulo-interstitial and glomerular kidney compartments, identified 17 genes differentially expressed in kidney. Fine-mapping highlighted missense driver variants in 11 genes and kidney-specific regulatory variants. These results provide a comprehensive priority list of molecular targets for translational research.

Integrated analysis of genomics, longitudinal metabolomics, and Alzheimer's risk factors among 1,111 cohort participants

Although Alzheimer's disease (AD) is highly heritable, genetic variants are known to be associated with AD only explain a small proportion of its heritability. Genetic factors may only convey disease risk in individuals with certain environmental exposures, suggesting that a multiomics approach could reveal underlying mechanisms contributing to complex traits, such as AD. We developed an integrated network to investigate relationships between metabolomics, genomics, and AD risk factors using Wisconsin Registry for Alzheimer's Prevention participants. Analyses included 1,111 non-Hispanic Caucasian participants with whole blood expression for 11,376 genes (imputed from dense genome-wide genotyping), 1,097 fasting plasma metabolites, and 17 AD risk factors. A subset of 155 individuals also had 364 fastings cerebral spinal fluid (CSF) metabolites. After adjusting each of these 12,854 variables for potential confounders, we developed an undirected graphical network, representing all significant pairwise correlations upon adjusting for multiple testing. There were many instances of genes being indirectly linked to AD risk factors through metabolites, suggesting that genes may influence AD risk through particular metabolites. Follow-up analyses suggested that glycine mediates the relationship between carbamoyl-phosphate synthase 1 and measures of cardiovascular and diabetes risk, including body mass index, waist-hip ratio, inflammation, and insulin resistance. Further, 38 CSF metabolites explained more than 60% of the variance of CSF levels of tau, a detrimental protein that accumulates in the brain of AD patients and is necessary for its diagnosis. These results further our understanding of underlying mechanisms contributing to AD risk while demonstrating the utility of generating and integrating multiple omics data types.

Genome-wide association study identifies novel loci for type 2 diabetes-attributed end-stage kidney disease in African Americans

BACKGROUND

End-stage kidney disease (ESKD) is a significant public health concern disproportionately affecting African Americans (AAs). Type 2 diabetes (T2D) is the leading cause of ESKD in the USA, and efforts to uncover genetic susceptibility to diabetic kidney disease (DKD) have had limited success. A prior genome-wide association study (GWAS) in AAs with T2D-ESKD was expanded with additional AA cases and controls and genotypes imputed to the higher density 1000 Genomes reference panel. The discovery analysis included 3432 T2D-ESKD cases and 6977 non-diabetic non-nephropathy controls (N = 10,409), followed by a discrimination analysis in 2756 T2D non-nephropathy controls to exclude T2D-associated variants.

RESULTS

Six independent variants located in or near RND3/RBM43, SLITRK3, ENPP7, GNG7, and APOL1 achieved genome-wide significant association (P < 5 × 10-8) with T2D-ESKD. Following extension analyses in 1910 non-diabetic ESKD cases and 908 non-diabetic non-nephropathy controls, a meta-analysis of 5342 AA all-cause ESKD cases and 6977 AA non-diabetic non-nephropathy controls revealed an additional novel all-cause ESKD locus at EFNB2 (rs77113398; P = 9.84 × 10-9; OR = 1.94). Exclusion of APOL1 renal-risk genotype carriers identified two additional genome-wide significant T2D-ESKD-associated loci at GRAMD3 and MGAT4C. A second variant at GNG7 (rs373971520; P = 2.17 × 10-8, OR = 1.46) remained associated with all-cause ESKD in the APOL1-negative analysis.

CONCLUSIONS

Findings provide further evidence for genetic factors associated with advanced kidney disease in AAs with T2D.

Sex-specific and pleiotropic effects underlying kidney function identified from GWAS meta-analysis

Chronic kidney disease (CKD) is a growing health burden currently affecting 10–15% of adults worldwide. Estimated glomerular filtration rate (eGFR) as a marker of kidney function is commonly used to diagnose CKD. We analyze eGFR data from the Nord-Trøndelag Health Study and Michigan Genomics Initiative and perform a GWAS meta-analysis with public summary statistics, more than doubling the sample size of previous meta-analyses. We identify 147 loci (53 novel) associated with eGFR, including genes involved in transcriptional regulation, kidney development, cellular signaling, metabolism, and solute transport. Additionally, sex-stratified analysis identifies one locus with more significant effects in women than men. Using genetic risk scores constructed from these eGFR meta-analysis results, we show that associated variants are generally predictive of CKD with only modest improvements in detection compared with other known clinical risk factors. Collectively, these results yield additional insight into the genetic factors underlying kidney function and progression to CKD.

Estimated glomerular filtration rate (eGFR) is a measure of kidney function and used to characterize chronic kidney disease. Here, Graham et al. identify 53 novel loci for eGFR in a GWAS meta-analysis, a subset of which are associated with other common diseases, such as diabetes and hypertension, based on PheWAS.

Multi-ancestry genome-wide gene-smoking interaction study of 387,272 individuals identifies new loci associated with serum lipids

The concentrations of high- and low-density-lipoprotein cholesterol and triglycerides are influenced by smoking, but it is unknown whether genetic associations with lipids may be modified by smoking. We conducted a multi-ancestry genome-wide gene-smoking interaction study in 133,805 individuals with follow-up in an additional 253,467 individuals. Combined meta-analyses identified 13 new loci associated with lipids, some of which were detected only because association differed by smoking status. Additionally, we demonstrate the importance of including diverse populations, particularly in studies of interactions with lifestyle factors, where genomic and lifestyle differences by ancestry may contribute to novel findings.

Renal phosphate handling and inherited disorders of phosphate reabsorption: an update

Renal phosphate handling critically determines plasma phosphate and whole body phosphate levels. Filtered phosphate is mostly reabsorbed by Na⁺-dependent phosphate transporters located in the brush border membrane of the proximal tubule: NaPi-IIa (SLC34A1), NaPi-IIc (SLC34A3), and Pit-2 (SLC20A2). Here we review new evidence for the role and relevance of these transporters in inherited disorders of renal phosphate handling. The importance of NaPi-IIa and NaPi-IIc for renal phosphate reabsorption and mineral homeostasis has been highlighted by the identification of mutations in these transporters in a subset of patients with infantile idiopathic hypercalcemia and patients with hereditary hypophosphatemic rickets with hypercalciuria. Both diseases are characterized by disturbed calcium homeostasis secondary to elevated 1,25-(OH)₂ vitamin D₃ as a consequence of hypophosphatemia. In vitro analysis of mutated NaPi-IIa or NaPi-IIc transporters suggests defective trafficking underlying disease in most cases. Monoallelic pathogenic mutations in both SLC34A1 and SLC34A3 appear to be very frequent in the general population and have been associated with kidney stones. Consistent with these findings, results from genome-wide association studies indicate that variants in SLC34A1 are associated with a higher risk to develop kidney stones and chronic kidney disease, but underlying mechanisms have not been addressed to date.

Evidence for Recent Polygenic Selection on Educational Attainment and Intelligence Inferred from Gwas Hits: A Replication of Previous Findings Using Recent Data

Genetic variants identified by three large genome-wide association studies (GWAS) of educational attainment (EA) were used to test a polygenic selection model. Weighted and unweighted polygenic scores (PGS) were calculated and compared across populations using data from the 1000 Genomes (n = 26), HGDP-CEPH (n = 52) and gnomAD (n = 8) datasets. The PGS from the largest EA GWAS was highly correlated to two previously published PGSs (r = 0.96–0.97, N = 26). These factors are both highly predictive of average population IQ (r = 0.9, N = 23) and Learning index (r = 0.8, N = 22) and are robust to tests of spatial autocorrelation. Monte Carlo simulations yielded highly significant p values. In the gnomAD samples, the correlation between PGS and IQ was almost perfect (r = 0.98, N = 8), and ANOVA showed significant population differences in allele frequencies with positive effect. Socioeconomic variables slightly improved the prediction accuracy of the model (from 78–80% to 85–89%), but the PGS explained twice as much of the variance in IQ compared to socioeconomic variables. In both 1000 Genomes and gnomAD, there was a weak trend for lower GWAS significance SNPs to be less predictive of population IQ. Additionally, a subset of SNPs were found in the HGDP-CEPH sample (N = 127). The analysis of this sample yielded a positive correlation with latitude and a low negative correlation with distance from East Africa. This study provides robust results after accounting for spatial autocorrelation with Fst distances and random noise via an empirical Monte Carlo simulation using null SNPs.

Association of Candidate Gene Polymorphisms With Chronic Kidney Disease: Results of a Case-Control Analysis in the Nefrona Cohort

Chronic kidney disease (CKD) is a major risk factor for end-stage renal disease, cardiovascular disease and premature death. Despite classical clinical risk factors for CKD and some genetic risk factors have been identified, the residual risk observed in prediction models is still high. Therefore, new risk factors need to be identified in order to better predict the risk of CKD in the population. Here, we analyzed the genetic association of 79 SNPs of proteins associated with mineral metabolism disturbances with CKD in a cohort that includes 2,445 CKD cases and 559 controls. Genotyping was performed with matrix assisted laser desorption ionization–time of flight mass spectrometry. We used logistic regression models considering different genetic inheritance models to assess the association of the SNPs with the prevalence of CKD, adjusting for known risk factors. Eight SNPs (rs1126616, rs35068180, rs2238135, rs1800247, rs385564, rs4236, rs2248359, and rs1564858) were associated with CKD even after adjusting by sex, age and race. A model containing five of these SNPs (rs1126616, rs35068180, rs1800247, rs4236, and rs2248359), diabetes and hypertension showed better performance than models considering only clinical risk factors, significantly increasing the area under the curve of the model without polymorphisms. Furthermore, one of the SNPs (the rs2248359) showed an interaction with hypertension, being the risk genotype affecting only hypertensive patients. We conclude that 5 SNPs related to proteins implicated in mineral metabolism disturbances (Osteopontin, osteocalcin, matrix gla protein, matrix metalloprotease 3 and 24 hydroxylase) are associated to an increased risk of suffering CKD.

Single-cell transcriptomics reveals gene expression dynamics of human fetal kidney development

The current understanding of mammalian kidney development is largely based on mouse models. Recent landmark studies revealed pervasive differences in renal embryogenesis between mouse and human. The scarcity of detailed gene expression data in humans therefore hampers a thorough understanding of human kidney development and the possible developmental origin of kidney diseases. In this paper, we present a single-cell transcriptomics study of the human fetal kidney. We identified 22 cell types and a host of marker genes. Comparison of samples from different developmental ages revealed continuous gene expression changes in podocytes. To demonstrate the usefulness of our data set, we explored the heterogeneity of the nephrogenic niche, localized podocyte precursors, and confirmed disease-associated marker genes. With close to 18,000 renal cells from five different developmental ages, this study provides a rich resource for the elucidation of human kidney development, easily accessible through an interactive web application.

Genome-Wide Association Studies of Metabolite Concentrations (mGWAS): Relevance for Nephrology

Metabolites are small molecules that are intermediates or products of metabolism, many of which are freely filtered by the kidneys. In addition, the kidneys have a central role in metabolite anabolism and catabolism, as well as in active metabolite reabsorption and/or secretion during tubular passage. This review article illustrates how the coupling of genomics and metabolomics in genome-wide association analyses of metabolites can be used to illuminate mechanisms underlying human metabolism, with a special focus on insights relevant to nephrology. First, genetic susceptibility loci for reduced kidney function and chronic kidney disease (CKD) were reviewed systematically for their associations with metabolite concentrations in metabolomics studies of blood and urine. Second, kidney function and CKD-associated metabolites reported from observational studies were interrogated for metabolite-associated genetic variants to generate and discuss complementary insights. Finally, insights originating from the simultaneous study of both blood and urine or by modeling intermetabolite relationships are summarized. We also discuss methodologic questions related to the study of metabolite concentrations in urine as well as among CKD patients. In summary, genome-wide association analyses of metabolites using metabolite concentrations quantified from blood and/or urine are a promising avenue of research to illuminate physiological and pathophysiological functions of the kidney.

Trans-ethnic kidney function association study reveals putative causal genes and effects on kidney-specific disease aetiologies

Chronic kidney disease (CKD) affects ~10% of the global population, with considerable ethnic differences in prevalence and aetiology. We assemble genome-wide association studies of estimated glomerular filtration rate (eGFR), a measure of kidney function that defines CKD, in 312,468 individuals of diverse ancestry. We identify 127 distinct association signals with homogeneous effects on eGFR across ancestries and enrichment in genomic annotations including kidney-specific histone modifications. Fine-mapping reveals 40 high-confidence variants driving eGFR associations and highlights putative causal genes with cell-type specific expression in glomerulus, and in proximal and distal nephron. Mendelian randomisation supports causal effects of eGFR on overall and cause-specific CKD, kidney stone formation, diastolic blood pressure and hypertension. These results define novel molecular mechanisms and putative causal genes for eGFR, offering insight into clinical outcomes and routes to CKD treatment development.

Estimated glomerular filtration rate (eGFR) is a measure of kidney function used to define chronic kidney disease. Here, Morris et al. perform trans-ethnic genome-wide meta-analyses for eGFR in 312,468 individuals and identify novel loci and downstream putative causal genes.

Clinical and genetic associations of renal function and diabetic kidney disease in the United Arab Emirates: a cross-sectional study

OBJECTIVES

Within the Emirati population, risk factors and genetic predisposition to diabetic kidney disease (DKD) have not yet been investigated. The aim of this research was to determine potential clinical, laboratory and reported genetic loci as risk factors for DKD.

RESEARCH DESIGN AND METHODS

Four hundred and ninety unrelated Emirati nationals with type 2 diabetes mellitus (T2DM) were recruited with and without DKD, and clinical and laboratory data were obtained. Following adjustments for possible confounders, a logistic regression model was developed to test the associations of 63 single nucleotide polymorphisms (SNPs) in 43 genetic loci with DKD (145 patients with DKD and 265 without DKD). Linear regression models, adjusted for age and gender, were then used to study the genetic associations of five renal function traits, including 83 SNPs with albumin-to-creatinine ratio, 92 SNPs with vitamin D (25-OH cholecalciferol), 288 SNPs with estimated glomerular filtration rate (eGFR), 363 SNPs with serum creatinine and 73 SNPs with blood urea.

RESULTS

Patients with DKD, as compared with those without the disease, were mostly men (52%vs38% for controls), older (67vs59 years) and had significant rates of hypertension and dyslipidaemia. Furthermore, patients with DKD had T2DM for a longer duration of time (16vs10 years), which in an additive manner was the single factor that significantly contributed to the development of DKD (p=0.02, OR=3.12, 95% CI 1.21 to 8.02). Among the replicated associations of the genetic loci with different renal function traits, the most notable included SHROOM3 with levels of serum creatinine, eGFR and DKD (P_adjusted=0.04, OR=1.46); CASR, GC and CYP2R1 with vitamin D levels; as well as WDR72 with serum creatinine and eGFR levels.

CONCLUSIONS

Associations were found between several genetic loci and risk markers for DKD, which may influence kidney function traits and DKD in a population of Arab ancestry.

Genetic Analysis of Rare Human Variants of Regulators of G Protein Signaling Proteins and Their Role in Human Physiology and Disease

Regulators of G protein signaling (RGS) proteins modulate the physiologic actions of many neurotransmitters, hormones, and other signaling molecules. Human RGS proteins comprise a family of 20 canonical proteins that bind directly to G protein-coupled receptors/G protein complexes to limit the lifetime of their signaling events, which regulate all aspects of cell and organ physiology. Genetic variations account for diverse human traits and individual predispositions to disease. RGS proteins contribute to many complex polygenic human traits and pathologies such as hypertension, atherosclerosis, schizophrenia, depression, addiction, cancers, and many others. Recent analysis indicates that most human diseases are due to extremely rare genetic variants. In this study, we summarize physiologic roles for RGS proteins and links to human diseases/traits and report rare variants found within each human RGS protein exome sequence derived from global population studies. Each RGS sequence is analyzed using recently described bioinformatics and proteomic tools for measures of missense tolerance ratio paired with combined annotation-dependent depletion scores, and protein post-translational modification (PTM) alignment cluster analysis. We highlight selected variants within the well-studied RGS domain that likely disrupt RGS protein functions and provide comprehensive variant and PTM data for each RGS protein for future study. We propose that rare variants in functionally sensitive regions of RGS proteins confer profound change-of-function phenotypes that may contribute, in newly appreciated ways, to complex human diseases and/or traits. This information provides investigators with a valuable database to explore variation in RGS protein function, and for targeting RGS proteins as future therapeutic targets.

Genetic Variants Associated with Circulating Fibroblast Growth Factor 23

BACKGROUND

Fibroblast growth factor 23 (FGF23), a bone-derived hormone that regulates phosphorus and vitamin D metabolism, contributes to the pathogenesis of mineral and bone disorders in CKD and is an emerging cardiovascular risk factor. Central elements of FGF23 regulation remain incompletely understood; genetic variation may help explain interindividual differences.

METHODS

We performed a meta-analysis of genome-wide association studies of circulating FGF23 concentrations among 16,624 participants of European ancestry from seven cohort studies, excluding participants with eGFR<30 ml/min per 1.73 m2 to focus on FGF23 under normal conditions. We evaluated the association of single-nucleotide polymorphisms (SNPs) with natural log-transformed FGF23 concentration, adjusted for age, sex, study site, and principal components of ancestry. A second model additionally adjusted for BMI and eGFR.

RESULTS

We discovered 154 SNPs from five independent regions associated with FGF23 concentration. The SNP with the strongest association, rs17216707 (P=3.0×10-24), lies upstream of CYP24A1, which encodes the primary catabolic enzyme for 1,25-dihydroxyvitamin D and 25-hydroxyvitamin D. Each additional copy of the T allele at this locus is associated with 5% higher FGF23 concentration. Another locus strongly associated with variations in FGF23 concentration is rs11741640, within RGS14 and upstream of SLC34A1 (a gene involved in renal phosphate transport). Additional adjustment for BMI and eGFR did not materially alter the magnitude of these associations. Another top locus (within ABO, the ABO blood group transferase gene) was no longer statistically significant at the genome-wide level.

CONCLUSIONS

Common genetic variants located near genes involved in vitamin D metabolism and renal phosphate transport are associated with differences in circulating FGF23 concentrations.

Genome-wide association studies and CRISPR/Cas9-mediated gene editing identify regulatory variants influencing eyebrow thickness in humans

Hair plays an important role in primates and is clearly subject to adaptive selection. While humans have lost most facial hair, eyebrows are a notable exception. Eyebrow thickness is heritable and widely believed to be subject to sexual selection. Nevertheless, few genomic studies have explored its genetic basis. Here, we performed a genome-wide scan for eyebrow thickness in 2961 Han Chinese. We identified two new loci of genome-wide significance, at 3q26.33 near SOX2 (rs1345417: P = 6.51×10(-10)) and at 5q13.2 near FOXD1 (rs12651896: P = 1.73×10(-8)). We further replicated our findings in the Uyghurs, a population from China characterized by East Asian-European admixture (N = 721), the CANDELA cohort from five Latin American countries (N = 2301), and the Rotterdam Study cohort of Dutch Europeans (N = 4411). A meta-analysis combining the full GWAS results from the three cohorts of full or partial Asian descent (Han Chinese, Uyghur and Latin Americans, N = 5983) highlighted a third signal of genome-wide significance at 2q12.3 (rs1866188: P = 5.81×10(-11)) near EDAR. We performed fine-mapping and prioritized four variants for further experimental verification. CRISPR/Cas9-mediated gene editing provided evidence that rs1345417 and rs12651896 affect the transcriptional activity of the nearby SOX2 and FOXD1 genes, which are both involved in hair development. Finally, suitable statistical analyses revealed that none of the associated variants showed clear signals of selection in any of the populations tested. Contrary to popular speculation, we found no evidence that eyebrow thickness is subject to strong selective pressure.

Genetic contributions to lupus nephritis in a multi-ethnic cohort of systemic lupus erythematous patients

OBJECTIVE

African Americans, East Asians, and Hispanics with systemic lupus erythematous (SLE) are more likely to develop lupus nephritis (LN) than are SLE patients of European descent. The etiology of this difference is not clear, and this study was undertaken to investigate how genetic variants might explain this effect.

METHODS

In this cross-sectional study, 1244 SLE patients from multiethnic case collections were genotyped for 817,810 single-nucleotide polymorphisms (SNPs) across the genome. Continental genetic ancestry was estimated utilizing the program ADMIXTURE. Gene-based testing and pathway analysis was performed within each ethnic group and meta-analyzed across ethnicities. We also performed candidate SNP association tests with SNPs previously established as risk alleles for SLE, LN, and chronic kidney disease (CKD). Association testing and logistic regression models were performed with LN as the outcome, adjusted for continental ancestries, sex, disease duration, and age.

RESULTS

We studied 255 North European, 263 South European, 238 Hispanic, 224 African American and 264 East Asian SLE patients, of whom 606 had LN (48.7%). In genome-wide gene-based and candidate SNP analyses, we found distinct genes, pathways and established risk SNPs associated with LN for each ethnic group. Gene-based analyses showed significant associations between variation in ZNF546 (p = 1.0E-06), TRIM15 (p = 1.0E-06), and TRIMI0 (p = 1.0E-06) and LN among South Europeans, and TTC34 (p = 8.0E-06) was significantly associated with LN among Hispanics. The SNP rs8091180 in NFATC1 was associated with LN (OR 1.43, p = 3.3E-04) in the candidate SNP meta-analysis with the highest OR among African-Americans (OR 2.17, p = 0.0035).

CONCLUSION

Distinct genetic factors are associated with the risk of LN in SLE patients of different ethnicities. CKD risk alleles may play a role in the development of LN in addition to SLE-associated risk variants. These findings may further explain the clinical heterogeneity of LN risk and response to therapy observed between different ethnic groups.

Genome-Wide Association Study of Renal Function Traits: Results from the Japan Multi-Institutional Collaborative Cohort Study

BACKGROUND

Chronic kidney disease (CKD) is a rapidly growing, worldwide public health problem. Recent advances in genome-wide-association studies (GWAS) revealed several genetic loci associated with renal function traits worldwide.

METHODS

We investigated the association of genetic factors with the levels of serum creatinine (SCr) and the estimated glomerular filtration rate (eGFR) in Japanese population-based cohorts analyzing the GWAS imputed data with 11,221 subjects and 12,617,569 variants, and replicated the findings with the 148,829 hospital-based Japanese subjects.

RESULTS

In the discovery phase, 28 variants within 4 loci (chromosome [chr] 2 with 8 variants including rs3770636 in the LDL receptor related protein 2 gene locus, on chr 5 with 2 variants including rs270184, chr 17 with 15 variants including rs3785837 in the BCAS3 gene locus, and chr 18 with 3 variants including rs74183647 in the nuclear factor of -activated T-cells 1 gene locus) reached the suggestive level of p < 1 × 10-6 in association with eGFR and SCr, and 2 variants on chr 4 (including rs78351985 in the microsomal triglyceride transfer protein gene locus) fulfilled the suggestive level in association with the risk of CKD. In the replication phase, 25 variants within 3 loci (chr 2 with 7 variants, chr 17 with 15 variants and chr 18 with 3 variants) in association with eGFR and SCr, and 2 variants on chr 4 associated with the risk of CKD became nominally statistically significant after Bonferroni correction, among which 15 variants on chr 17 and 3 variants on chr 18 reached genome-wide significance of p < 5 × 10-8 in the combined study meta-analysis. The associations of the loci on chr 2 and 18 with eGFR and SCr as well as that on chr 4 with CKD risk have not been previously reported in the Japanese and East Asian populations.

CONCLUSION

Although the present GWAS of renal function traits included the largest sample of Japanese participants to date, we did not identify novel loci for renal traits. However, we identified the novel associations of the genetic loci on chr 2, 4, and 18 with renal function traits in the Japanese population, suggesting these are transethnic loci. Further investigations of these associations are expected to further validate our findings for the potential establishment of personalized prevention of renal disease in the Japanese and East Asian populations.

Improved detection of genetic loci in estimated glomerular filtration rate and type 2 diabetes using a pleiotropic cFDR method

Genome-wide association studies (GWAS) have been shown to have the potential of explaining more of the "missing heritability" of complex human phenotypes by improving statistical approaches. Here, we applied a genetic-pleiotropy-informed conditional false discovery rate (cFDR) to capture additional polygenic effects associated with estimated glomerular filtration rate (creatinine) (eGFRcrea) and type 2 diabetes (T2D). The cFDR analysis improves the identification of pleiotropic variants by incorporating potentially shared genetic mechanisms between two related traits. The Q-Q and fold-enrichment plots were used to assess the enrichment of SNPs associated with eGFRcrea or T2D, and Manhattan plots were used for showing chromosomal locations of the significant loci detected. By applying the cFDR method, we newly identified 74 loci for eGFRcrea and 3 loci for T2D with the cFDR criterion of 0.05 compared with previous related GWAS studies. Four shared SNPs were detected to be associated with both eGFRcrea and T2D at the significant conjunction cFDR level of 0.05, and these shared SNPs had not been reported in previous studies. In addition, we used DAVID analysis to perform functional analysis of the shared SNPs' annotated genes and found their potential hidden associations with eGFRcrea and T2D. In this study, the cFDR method shows the feasibility to detect more genetic variants underlying the heritability of eGFRcrea and T2D, and the overlapping SNPs identified could be regarded as candidate loci that provide a thread of genetic mechanisms between eGFRcrea and T2D in future research.

Trans-ethnic meta-regression of genome-wide association studies accounting for ancestry increases power for discovery and improves fine-mapping resolution

Trans-ethnic meta-analysis of genome-wide association studies (GWAS) across diverse populations can increase power to detect complex trait loci when the underlying causal variants are shared between ancestry groups. However, heterogeneity in allelic effects between GWAS at these loci can occur that is correlated with ancestry. Here, a novel approach is presented to detect SNP association and quantify the extent of heterogeneity in allelic effects that is correlated with ancestry. We employ trans-ethnic meta-regression to model allelic effects as a function of axes of genetic variation, derived from a matrix of mean pairwise allele frequency differences between GWAS, and implemented in the MR-MEGA software. Through detailed simulations, we demonstrate increased power to detect association for MR-MEGA over fixed- and random-effects meta-analysis across a range of scenarios of heterogeneity in allelic effects between ethnic groups. We also demonstrate improved fine-mapping resolution, in loci containing a single causal variant, compared to these meta-analysis approaches and PAINTOR, and equivalent performance to MANTRA at reduced computational cost. Application of MR-MEGA to trans-ethnic GWAS of kidney function in 71,461 individuals indicates stronger signals of association than fixed-effects meta-analysis when heterogeneity in allelic effects is correlated with ancestry. Application of MR-MEGA to fine-mapping four type 2 diabetes susceptibility loci in 22,086 cases and 42,539 controls highlights: (i) strong evidence for heterogeneity in allelic effects that is correlated with ancestry only at the index SNP for the association signal at the CDKAL1 locus; and (ii) 99% credible sets with six or fewer variants for five distinct association signals.

Lessons from CKD-Related Genetic Association Studies-Moving Forward

Over the past decade, genetic association studies have uncovered numerous determinants of kidney function in the general, diabetic, hypertensive, CKD, ESRD, and GN-based study populations (e.g., IgA nephropathy, membranous nephropathy, FSGS). These studies have led to numerous novel and unanticipated findings, which are helping improve our understanding of factors and pathways affecting both normal and pathologic kidney function. In this review, we report on major discoveries and advances resulting from this rapidly progressing research domain. We also predict some of the next steps the nephrology community should embrace to accelerate the identification of genetic and molecular processes leading to kidney dysfunction, pathophysiologically based disease subgroups, and specific therapeutic targets, as we attempt to transition toward a more precision-based medicine approach.

The UMOD Locus: Insights into the Pathogenesis and Prognosis of Kidney Disease

The identification of genetic factors associated with kidney disease has the potential to provide critical insights into disease mechanisms. Genome-wide association studies have uncovered genomic regions associated with renal function metrics and risk of CKD. UMOD is among the most outstanding loci associated with CKD in the general population, because it has a large effect on eGFR and CKD risk that is consistent across different ethnic groups. The relevance of UMOD for CKD is clear, because the encoded protein, uromodulin (Tamm-Horsfall protein), is exclusively produced by the kidney tubule and has specific biochemical properties that mediate important functions in the kidney and urine. Rare mutations in UMOD are the major cause of autosomal dominant tubulointerstitial kidney disease, a condition that leads to CKD and ESRD. In this brief review, we use the UMOD paradigm to describe how population genetic studies can yield insight into the pathogenesis and prognosis of kidney diseases.