Genetic loci influencing kidney function and chronic kidney disease

Genetic Susceptibility to Chronic Kidney Disease - Some More Pieces for the Heritability Puzzle

Chronic kidney disease (CKD) is a major global health problem with an increasing prevalence partly driven by aging population structure. Both genomic and environmental factors contribute to this complex heterogeneous disease. CKD heritability is estimated to be high (30-75%). Genome-wide association studies (GWAS) and GWAS meta-analyses have identified several genetic loci associated with CKD, including variants in UMOD, SHROOM3, solute carriers, and E3 ubiquitin ligases. However, these genetic markers do not account for all the susceptibility to CKD, and the causal pathways remain incompletely understood; other factors must be contributing to the missing heritability. Less investigated biological factors such as telomere length; mitochondrial proteins, encoded by nuclear genes or specific mitochondrial DNA (mtDNA) encoded genes; structural variants, such as copy number variants (CNVs), insertions, deletions, inversions and translocations are poorly covered and may explain some of the missing heritability. The sex chromosomes, often excluded from GWAS studies, may also help explain gender imbalances in CKD. In this review, we outline recent findings on molecular biomarkers for CKD (telomeres, CNVs, mtDNA variants, sex chromosomes) that typically have received less attention than gene polymorphisms. Shorter telomere length has been associated with renal dysfunction and CKD progression, however, most publications report small numbers of subjects with conflicting findings. CNVs have been linked to congenital anomalies of the kidney and urinary tract, posterior urethral valves, nephronophthisis and immunoglobulin A nephropathy. Information on mtDNA biomarkers for CKD comes primarily from case reports, therefore the data are scarce and diverse. The most consistent finding is the A3243G mutation in the MT-TL1 gene, mainly associated with focal segmental glomerulosclerosis. Only one GWAS has found associations between X-chromosome and renal function (rs12845465 and rs5987107). No loci in the Y-chromosome have reached genome-wide significance. In conclusion, despite the efforts to find the genetic basis of CKD, it remains challenging to explain all of the heritability with currently available methods and datasets. Although additional biomarkers have been investigated in less common suspects such as telomeres, CNVs, mtDNA and sex chromosomes, hidden heritability in CKD remains elusive, and more comprehensive approaches, particularly through the integration of multiple -"omics" data, are needed.

The kidney transcriptome, from single cells to whole organs and back

PURPOSE OF REVIEW

Transcriptome analysis of human kidney samples provides an integrated output of genetic, physiological, or environmental inputs. This review summarizes recent findings including gene expression and genetic variation integration, bulk and single cell gene expression analysis, and describes how such studies have improved our understanding of kidney disease development.

RECENT FINDINGS

Bulk or whole tissue analysis of patient kidney samples identified a large number of genes, whose levels correlate with kidney function and/or structural damage. These genes were enriched for metabolic and immune functions. Using expression quantitative trait analysis, genetic variations-driven gene expression can be identified. Recent developments in single cell sequencing defined cell-type-specific gene expression changes and highlighted specific cell types for disease development.

SUMMARY

Recent advancement in whole tissue transcriptomics, specifically incorporating genotype information and single cell data have been powerful to identify kidney disease-associated genes, pathways, and cell types.

The mercapturomic profile of health and non-communicable diseases

The mercapturate pathway is a unique metabolic circuitry that detoxifies electrophiles upon adducts formation with glutathione. Since its discovery over a century ago, most of the knowledge on the mercapturate pathway has been provided from biomonitoring studies on environmental exposure to toxicants. However, the mercapturate pathway-related metabolites that is formed in humans—the mercapturomic profile—in health and disease is yet to be established. In this paper, we put forward the hypothesis that these metabolites are key pathophysiologic factors behind the onset and development of non-communicable chronic inflammatory diseases. This review goes from the evidence in the formation of endogenous metabolites undergoing the mercapturate pathway to the methodologies for their assessment and their association with cancer and respiratory, neurologic and cardiometabolic diseases.

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.

Uromodulin deficiency alters tubular injury and interstitial inflammation but not fibrosis in experimental obstructive nephropathy

Human GWAS and Mendelian genetic studies have linked polymorphic variants and mutations in the human uromodulin gene (UMOD) with chronic kidney disease. The primary function of this kidney‐specific and secreted protein remains elusive. This study investigated whether UMOD deficiency modified responses to unilateral ureteral obstruction (UUO)‐induced kidney injury. Kidneys harvested from groups of wild‐type (UMOD+/+) and knockout (UMOD−/−) male mice (n = 7–10 each) were studied on days 7, 14, and 21. Compared to sham kidneys, UMOD protein levels increased 9–13x after UUO and were associated with increased urinary UMOD levels. Kidney KIM‐1 protein levels were higher in the UMOD−/− groups at all time‐points (4–14x). The UMOD−/− groups also had higher KIM‐1 kidney‐to‐urine relative ratios (5–35x). In vitro studies using KIM‐1 expressing 769‐P cells showed lower KIM‐1 levels in the presence of UMOD protein. Levels of proapoptotic genes and the epithelial cell apoptotic protein marker M30 were significantly lower in the UMOD−/− groups. Both M30 and KIM‐1 colocalized with intraluminal UMOD protein deposits. Interstitial inflammation was less intense in the UMOD−/− groups. Renal fibrosis severity (kidney collagen mRNA and protein) was similar in both genotypic groups on days 7, 14, and 21. Our findings suggest a role for UMOD‐dependent inhibition of KIM‐1 expression and its apoptotic cell scavenging responses during chronic obstruction‐associated tubular injury.

Proteomic and metabolomic approaches in the search for biomarkers in chronic kidney disease

Chronic kidney disease (CKD) is an aging-related disorder that represents a major global public health burden. Current biochemical biomarkers, such as serum creatinine and urinary albumin, have important limitations when used to identify the earliest indication of CKD or in tracking the progression to more advanced CKD. These issues underline the importance of finding and testing new molecular biomarkers that are capable of successfully meeting this clinical need. The measurement of changes in nature and/or levels of proteins and metabolites in biological samples from patients provide insights into pathophysiological processes. Proteomic and metabolomic techniques provide opportunities to record dynamic chemical signatures in patients over time. This review article presents an overview of the recent developments in the fields of metabolomics and proteomics in relation to CKD. Among the many different proteomic biomarkers proposed, there is particular interest in the CKD273 classifier, a urinary proteome biomarker reported to predict CKD progression and with implementation potential. Other individual non-invasive peptidomic biomarkers that are potentially relevant for CKD detection include type 1 collagen, uromodulin and mucin-1. Despite the limited sample sizes and variability of the metabolomics studies, some metabolites such as trimethylamine N-oxide, kynurenine and citrulline stand out as potential biomarkers in CKD.

Multi-Imaging Investigation to Evaluate the Relationship between Serum Cystatin C and Features of Atherosclerosis in Non-ST-Segment Elevation Acute Coronary Syndrome

Objectives: High cystatin C(CysC) levels are associated with impaired cardiovascular outcome. Whether CysC levels are independently related to the atherosclerosis burden is still controversial. Methods: We enrolled 31 non-ST-segment elevation acute coronary syndrome patients undergoing percutaneous coronary intervention. Patients were divided into 2 groups on the basis of median value of serum CysC. Using the high CysC group as a dependent variable, univariable and multivariable analyses were used to evaluate the association between CysC and three different features of atherosclerosis: 1) coronary plaque vulnerability as assessed by optical coherence tomography (OCT), 2) coronary artery calcium (CAC) by means of computed tomography scan, and 3) aortic wall metabolic activity, as assessed using 18F-Fluorodeoxyglucose-positron emission tomography (18F-FDG-PET). Results: After univariable and multivariable analyses, 18F-FDG uptake in the descending aorta (DA) was independently associated with a low level of CysC [(Odds Ratio = 0.02; 95%CI 0.0004–0.89; p = 0.044; 18F-FDG uptake measured as averaged maximum target to blood ratio); (Odds Ratio = 0.89; 95%CI 0.82–0.98, p = 0.025; 18F-FDG uptake measured as number of active slices)]. No trend was found for the association between CysC and characteristics of OCT-assessed coronary plaque vulnerability or CAC score. Conclusions: In patients with non-ST-segment elevation acute coronary syndrome (NSTE-ACS), 18F-FDG uptake in the DA was associated with a low level of serum CysC. There was no relation between CysC levels and OCT-assessed coronary plaque vulnerability or CAC score. These findings suggest that high levels of CysC may not be considered as independent markers of atherosclerosis.

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.

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.

Functional Prediction of Chronic Kidney Disease Susceptibility Gene PRKAG2 by Comprehensively Bioinformatics Analysis

The genetic predisposition to chronic kidney disease (CKD) has been widely evaluated especially using the genome-wide association studies, which highlighted some novel genetic susceptibility variants in many genes, and estimated glomerular filtration rate to diagnose and stage CKD. Of these variants, rs7805747 in PRKAG2 was identified to be significantly associated with both serum creatinine and CKD with genome wide significance level. Until now, the potential mechanism by which rs7805747 affects CKD risk is still unclear. Here, we performed a functional analysis of rs7805747 variant using multiple bioinformatics software and databases. Using RegulomeDB and HaploReg (version 4.1), rs7805747 was predicated to locate in enhancer histone marks (Liver, Duodenum Mucosa, Fetal Intestine Large, Fetal Intestine Small, and Right Ventricle tissues). Using GWAS analysis in PhenoScanner, we showed that rs7805747 is not only associated with CKD, but also is significantly associated with other diseases or phenotypes. Using metabolite analysis in PhenoScanner, rs7805747 is identified to be significantly associated with not only the serum creatinine, but also with other 16 metabolites. Using eQTL analysis in PhenoScanner, rs7805747 is identified to be significantly associated with gene expression in multiple human tissues and multiple genes including PRKAG2. The gene expression analysis of PRKAG2 using 53 tissues from GTEx RNA-Seq of 8555 samples (570 donors) in GTEx showed that PRKAG2 had the highest median expression in Heart-Atrial Appendage. Using the gene expression profiles in human CKD, we further identified different expression of PRKAG2 gene in CKD cases compared with control samples. In summary, our findings provide new insight into the underlying susceptibility of PRKAG2 gene to CKD.

Modeling epigenetic modifications in renal development and disease with organoids and genome editing

Understanding epigenetic mechanisms is crucial to our comprehension of gene regulation in development and disease. In the past decades, different studies have shown the role of epigenetic modifications and modifiers in renal disease, especially during its progression towards chronic and end-stage renal disease. Thus, the identification of genetic variation associated with chronic kidney disease has resulted in better clinical management of patients. Despite the importance of these findings, the translation of genotype-phenotype data into gene-based medicine in chronic kidney disease populations still lacks faithful cellular or animal models that recapitulate the key aspects of the human kidney. The latest advances in the field of stem cells have shown that it is possible to emulate kidney development and function with organoids derived from human pluripotent stem cells. These have successfully recapitulated not only kidney differentiation, but also the specific phenotypical traits related to kidney function. The combination of this methodology with CRISPR/Cas9 genome editing has already helped researchers to model different genetic kidney disorders. Nowadays, CRISPR/Cas9-based approaches also allow epigenetic modifications, and thus represent an unprecedented tool for the screening of genetic variants, epigenetic modifications or even changes in chromatin structure that are altered in renal disease. In this Review, we discuss these technical advances in kidney modeling, and offer an overview of the role of epigenetic regulation in kidney development and disease.

Role of Alström syndrome 1 in the regulation of blood pressure and renal function

Elevated blood pressure (BP) and renal dysfunction are complex traits representing major global health problems. Single nucleotide polymorphisms identified by genome-wide association studies have identified the Alström syndrome 1 (ALMS1) gene locus to render susceptibility for renal dysfunction, hypertension, and chronic kidney disease (CKD). Mutations in the ALMS1 gene in humans causes Alström syndrome, characterized by progressive metabolic alterations including hypertension and CKD. Despite compelling genetic evidence, the underlying biological mechanism by which mutations in the ALMS1 gene lead to the above-mentioned pathophysiology is not understood. We modeled this effect in a KO rat model and showed that ALMS1 genetic deletion leads to hypertension. We demonstrate that the link between ALMS1 and hypertension involves the activation of the renal Na+/K+/2Cl- cotransporter NKCC2, mediated by regulation of its endocytosis. Our findings establish a link between the genetic susceptibility to hypertension, CKD, and the expression of ALMS1 through its role in a salt-reabsorbing tubular segment of the kidney. These data point to ALMS1 as a potentially novel gene involved in BP and renal function regulation.

Genomic approaches in the search for molecular biomarkers in chronic kidney disease

BACKGROUND

Chronic kidney disease (CKD) is recognised as a global public health problem, more prevalent in older persons and associated with multiple co-morbidities. Diabetes mellitus and hypertension are common aetiologies for CKD, but IgA glomerulonephritis, membranous glomerulonephritis, lupus nephritis and autosomal dominant polycystic kidney disease are also common causes of CKD.

MAIN BODY

Conventional biomarkers for CKD involving the use of estimated glomerular filtration rate (eGFR) derived from four variables (serum creatinine, age, gender and ethnicity) are recommended by clinical guidelines for the evaluation, classification, and stratification of CKD. However, these clinical biomarkers present some limitations, especially for early stages of CKD, elderly individuals, extreme body mass index values (serum creatinine), or are influenced by inflammation, steroid treatment and thyroid dysfunction (serum cystatin C). There is therefore a need to identify additional non-invasive biomarkers that are useful in clinical practice to help improve CKD diagnosis, inform prognosis and guide therapeutic management.

CONCLUSION

CKD is a multifactorial disease with associated genetic and environmental risk factors. Hence, many studies have employed genetic, epigenetic and transcriptomic approaches to identify biomarkers for kidney disease. In this review, we have summarised the most important studies in humans investigating genomic biomarkers for CKD in the last decade. Several genes, including UMOD, SHROOM3 and ELMO1 have been strongly associated with renal diseases, and some of their traits, such as eGFR and serum creatinine. The role of epigenetic and transcriptomic biomarkers in CKD and related diseases is still unclear. The combination of multiple biomarkers into classifiers, including genomic, and/or epigenomic, may give a more complete picture of kidney diseases.

A renal genetic risk score (GRS) is associated with kidney dysfunction in people with type 2 diabetes

This study aims to investigate whether renal and cardiovascular phenotypes in Italian patients with type 2 diabetes (T2D) could be influenced by a number of disease risk SNPs recently found in genome-wide association studies (GWAS). In 1591 Italian subjects with T2D: (1) 47 SNPs associated to kidney function and/or chronic kidney disease (CKD) and 49 SNPs associated to cardiovascular disease (CVD) risk were genotyped; (2) urinary albumin/creatinine (A/C) ratio, glomerular filtration rate (eGFR) and lipid profile were assessed; (3) a standard electrocardiogram was performed; (4) two genotype risk scores (GRS) were computed (a renal GRS calculated selecting 39 SNPs associated with intermediate traits of kidney damage and a cardiovascular GRS determined selecting 42 SNPs associated to CVD risk phenotypes). After correction for multiple comparisons, the renal GRS was not associated to A/C ratio (p = 0.33), but it was significantly related to decreased eGFR (p = 0.005). No association between the cardiovascular GRS and electrocardiogram was detected. Thus, in Italian patients with T2D a renal GRS might predict the decline in glomerular function, suggesting that the clock of diabetes associated CKD starts ticking long before hyperglycemia. Our data support the feasibility of gene-based prediction of complications in people with T2D.

Identification of 13 novel susceptibility loci for early-onset myocardial infarction, hypertension, or chronic kidney disease

Early-onset cardiovascular and renal diseases have a strong genetic component. In the present study, exome-wide association studies (EWASs) were performed to identify genetic variants that confer susceptibility to early-onset myocardial infarction (MI), hypertension, or chronic kidney disease (CKD) in Japanese individuals. A total of 8,093 individuals aged ≤65 years was enrolled in the study. The EWASs for MI, hypertension, and CKD were performed in 6,926 subjects (1,152 cases, 5,774 controls), 8,080 subjects (3,444 cases, 4,636 controls), and 2,556 subjects (1,051 cases, 1,505 controls), respectively. Genotyping of single nucleotide polymorphisms (SNPs) was performed with Illumina Human Exome-12 DNA Analysis BeadChip or Infinium Exome-24 BeadChip arrays. The associations of allele frequencies for 31,245, 31,276, or 31,514 SNPs that passed quality control to MI, hypertension, and CKD, respectively, was examined with Fisher's exact test. Bonferroni's correction for statistical significance of association was applied to compensate for multiple comparisons of genotypes with MI, hypertension, or CKD. The EWASs of allele frequencies revealed that 25, 11, and 11 SNPs were significantly associated with MI (P<1.60×10⁻⁶), hypertension (P<1.60×10⁻⁶), or CKD (P<1.59×10⁻⁶), respectively. Multivariable logistic regression analysis with adjustment for covariates showed that all 25, 11, and 11 SNPs were significantly associated with MI (P<0.0005), hypertension (P<0.0011), or CKD (P<0.0011), respectively. On examination of the results from previous genome-wide association studies and linkage disequilibrium of the identified SNPs, 11 loci (TMOD4, COL6A3, ADGRL3-CXCL8-MARCH1, OR52E4, TCHP-GIT2, CCDC63, 12q24.1, OAS3, PLCB2-VPS33B, GOSR2, ZNF77), six loci (MOB3C-TMOD4, COL6A3, COL6A5, CXCL8-MARCH1, NFKBIL1-6p21.3-NCR3, PLCB2-VPS33B), and seven loci (MOB3C-TMOD4, COL6A3, COL6A5, ADGRL3-CXCL8-MARCH1, MUC17, PLCB2-VPS33B, ZNF77) were identified as novel loci significantly associated with MI, hypertension, and CKD, respectively. Furthermore, six genes (TMOD4, COL6A3, CXCL8, MARCH1, PLCB2, VPS33B) were significantly associated with MI, hypertension and CKD; two genes (ADGRL3, ZNF77) with MI and CKD; and two genes (COL6A5, MOB3C) with hypertension and CKD. Therefore, 13 novel loci (MOB3C-TMOD4, COL6A3, ADGRL3-CXCL8-MARCH1, OR52E4, TCHP- GIT2, CCDC63, 12q24.1, OAS3, PLCB2-VPS33B, ZNF77, COL6A5, NFKBIL1-NCR3, MUC17) were identified that confer susceptibility to early-onset MI, hypertension, or CKD. The determination of genotypes for the SNPs at these loci may provide informative for assessment of the genetic risk for MI, hypertension, or CKD.

Adiponectin Receptor gene Polymorphisms are Associated with Kidney Function in Elderly Japanese Populations

Aim: Adiponectin exhibits its biological effects through adiponectin receptors (AdipoR1 and AdipoR2), which are distributed in the kidneys, and activation of those receptors could prevent or ameliorate diabetic nephropathy. This study aimed to evaluate the associations between AdipoR single nucleotide polymorphisms (SNPs) and kidney function in an elderly Japanese population.

Methods: A total of 271 elderly Japanese volunteers underwent anthropometric and laboratory tests (cystatin C-based eGFR and total and high molecular weight adiponectin levels at baseline and a follow-up visit). Genotype data were obtained for the selected 7 and 5 AdipoR1 and AdipoR2 SNPs, respectively.

Results: In a cross-sectional analysis at baseline, we found a significant association between the AdipoR2 SNP rs12230440 and kidney function; eGFRcys tended to increase as the number of carriers of T alleles increased after adjustment for covariates and Bonferroni correction, although the association of the SNP and annual eGFR decline could not be identified in the longitudinal data. Regarding the variants rs16850797, rs11061925, and rs10773983, each of the allele G, allele C, and allele G showed nominally significant associations with higher eGFRcys. However, this failed to reach significance after Bonferroni correction.

Conclusion: Here, an AdipoR2 SNP was associated with kidney function, suggesting that the effects of this polymorphism on adiponectin receptor may affect kidney function in the elderly Japanese population.

Cellular and molecular mechanisms of kidney fibrosis

Renal fibrosis is the final pathological process common to any ongoing, chronic kidney injury or maladaptive repair. It is considered as the underlying pathological process of chronic kidney disease (CKD), which affects more than 10% of world population and for which treatment options are limited. Renal fibrosis is defined by excessive deposition of extracellular matrix, which disrupts and replaces the functional parenchyma that leads to organ failure. Kidney's histological structure can be divided into three main compartments, all of which can be affected by fibrosis, specifically termed glomerulosclerosis in glomeruli, interstitial fibrosis in tubulointerstitium and arteriosclerosis and perivascular fibrosis in vasculature. In this review, we summarized the different appearance, cellular origin and major emerging processes and mediators of fibrosis in each compartment. We also depicted and discussed the challenges in translation of anti-fibrotic treatment to clinical practice and discuss possible solutions and future directions.

Assessing known chronic kidney disease associated genetic variants in Saudi Arabian populations

Background

Genome wide association studies of patients with European descent have identified common variants associated with risk of reduced estimated glomerular filtration rate (eGFR). A panel of eight variants were selected to evaluate their association and prevalence in a Saudi Arabian patient cohort with chronic kidney disease (CKD).

Methods

Eight genetic variants in four genes (SHROOM3, MYH9, SLC7A9, and CST3) were genotyped in 160 CKD patients and 189 ethnicity-matched healthy controls. Genetic variants were tested for association with the development of CKD (eGFR < 60 ml/min/1.73m²) and effects were compared with results obtained from 133,413 participants in the CKD genetics consortium. Multivariable regression was used to evaluate the role of these eight variants in improving prediction of CKD development.

Results

All eight variants were present in Saudi populations with minor allele frequency ranging from 16 to 46%. The risk variant in all four genes demonstrated the same direction of effect as observed in European populations. One variant, rs4821480, in MYH9 was significantly associated with increased risk of development of CKD (OR = 1.69, 95% CI 1.22–2.36, P = 0.002), but the additional variants were not statistically significant given our modest sample size.

Conclusions

CKD risk variants identified in European populations are present in Saudis. We did not find evidence to suggest heterogeneity of effect size compared to previously published estimates in European populations. Multivariable logistic regression analysis showed a statistically significant improvement in predicting the CKD using models with either FGF23 and vitamin D or FGF23, vitamin D level, and MYH9 genotypes (AUC = 0.93, 95% CI 0.90–0.95, P <  0.0001).

Electronic supplementary material

The online version of this article (10.1186/s12882-018-0890-9) contains supplementary material, which is available to authorized users.

Uromodulin rs4293393 T>C variation is associated with kidney disease in patients with type 2 diabetes

Background & objectives

Uromodulin, a UMOD gene encoded glycoprotein is synthesized exclusively in renal tubular cells and released into urine. Mutations lead to uromodulin misfolding and retention in the kidney, where it might stimulate cells of immune system to cause inflammation and progression of kidney disease. Genome-wide association studies (GWAS) have identified UMOD locus to be associated with hypertension and diabetic nephropathy (DN). In this study, we investigated the association between rs4293393 variation in UMOD gene and susceptibility to kidney disease in individuals with type 2 diabetes mellitus (T2DM).

Methods

A total of 646 individuals, 208 with T2DM without evidence of kidney disease (DM), 221 with DN and 217 healthy controls (HC) were genotyped for UMOD variant rs4293393T>C by restriction fragment length polymorphism. Serum uromodulin levels were quantified by enzyme-linked immunosorbent assay.

Results

A significant difference was found in genotype and allelic frequency among DM, DN and HC. TC+CC genotype and C allele were found more frequently in DN compared to HC (33.9 vs 23.0%, P=0.011 and 20.1 vs 12.9%, P=0.004, respectively). Compared to DM, C allele was found to be more frequent in individuals with DN (20.1 vs 14.7%, P=0.034). Those with DN had higher serum uromodulin levels compared to those with DM (P=0.001). Serum uromodulin levels showed a positive correlation with serum creatinine (r=0.431; P<0.001) and negative correlation with estimated glomerular filtration rate (r=-0.423; P<0.001).

Interpretation & conclusions

The frequency of UMOD rs4293393 variant with C allele was significantly higher in individuals with DN. UMOD rs4293393 T>C variation might have a bearing on susceptibility to nephropathy in north Indian individuals with type 2 diabetes.

Uromodulin associates with cardiorenal function in patients with hypertension and cardiovascular disease

OBJECTIVE

Common genetic variants in the gene encoding uromodulin (UMOD) have been associated with renal function, blood pressure (BP) and hypertension. We investigated the associations between an important single nucleotide polymorphism (SNP) in UMOD, that is rs12917707-G>T, and estimated glomerular filtration rate (eGFR), BP and cardiac organ damage as determined by echocardiography in patients with arterial hypertension.

METHODS

A cohort of 1218 treated high-risk patients (mean age 58.5 years, 83% men) with documented cardiovascular disease (81% with coronary heart disease) was analysed.

RESULTS

The mean values for 24-h SBP and DBP were 124.7 ± 14.7 and 73.9 ± 9.4 mmHg; mean eGFR was 77.5 ± 18.3 ml/min per 1.73 m, mean left ventricular ejection fraction was 59.3 ± 9.9% and mean left ventricular mass index in men and women was 53.9 ± 23.2 and 54.9 ± 23.7 g/m with 50.4% of patients having left ventricular hypertrophy. A significant association between rs12917707 and eGFR was observed with T-allele carriers showing significantly higher eGFR values (+2.6 ml/min per 1.73 m, P = 0.006) than noncarriers. This SNP associated also with left atrial diameter (P = 0.007); homozygous carriers of the T-allele had smaller left atrial diameter (-1.5 mm) than other genotype groups (P = 0.040). No significant associations between rs12917707 and other cardiac or BP phenotypes were observed.

CONCLUSIONS

These findings extend the previously documented role of UMOD for renal function also to treated high-risk patients with arterial hypertension and reveal a novel association with left atrial remodelling and thus a potential cardiorenal link modulated by UMOD.

Genome-wide association analysis identifies multiple loci associated with kidney disease-related traits in Korean populations

Chronic kidney disease (CKD) is an important social health problem characterized by a decrease in the kidney glomerular filtration rate (GFR). In this study, we analyzed genome-wide association studies for kidney disease-related traits using data from a Korean adult health screening cohort comprising 7,064 participants. Kidney disease-related traits analyzed include blood urea nitrogen (BUN), serum creatinine, estimated GFR, and uric acid levels. We detected two genetic loci (SLC14A2 and an intergenic region) and 8 single nucleotide polymorphisms (SNPs) associated with BUN, 3 genetic loci (BCAS3, C17orf82, ALDH2) and 6 SNPs associated with serum creatinine, 3 genetic loci (BCAS3, C17orf82/TBX2, LRP2) and 7 SNPs associated with GFR, and 14 genetic loci (3 in ABCG2/PKD2, 2 in SLC2A9, 3 in intergenic regions on chromosome 4; OTUB1, NRXN2/SLC22A12, CDC42BPG, RPS6KA4, SLC22A9, and MAP4K2 on chromosome 11) and 84 SNPs associated with uric acid levels. By comparing significant genetic loci associated with serum creatinine levels and GFR, rs9895661 in BCAS3 and rs757608 in C17orf82 were simultaneously associated with both traits. The SNPs rs11710227 in intergenic regions on chromosome 3 showing significant association with BUN is newly discovered. Genetic variations of multiple gene loci are associated with kidney disease-related traits, and differences in associations between kidney disease-related traits and genetic variation are dependent on the population. The meanings of the mutations identified in this study will need to be reaffirmed in other population groups in the future.

Characterization of Coding/Noncoding Variants for SHROOM3 in Patients with CKD

Background Interpreting genetic variants is one of the greatest challenges impeding analysis of rapidly increasing volumes of genomic data from patients. For example, SHROOM3 is an associated risk gene for CKD, yet causative mechanism(s) of SHROOM3 allele(s) are unknown.Methods We used our analytic pipeline that integrates genetic, computational, biochemical, CRISPR/Cas9 editing, molecular, and physiologic data to characterize coding and noncoding variants to study the human SHROOM3 risk locus for CKD.Results We identified a novel SHROOM3 transcriptional start site, which results in a shorter isoform lacking the PDZ domain and is regulated by a common noncoding sequence variant associated with CKD (rs17319721, allele frequency: 0.35). This variant disrupted allele binding to the transcription factor TCF7L2 in podocyte cell nuclear extracts and altered transcription levels of SHROOM3 in cultured cells, potentially through the loss of repressive looping between rs17319721 and the novel start site. Although common variant mechanisms are of high utility, sequencing is beginning to identify rare variants involved in disease; therefore, we used our biophysical tools to analyze an average of 112,849 individual human genome sequences for rare SHROOM3 missense variants, revealing 35 high-effect variants. The high-effect alleles include a coding variant (P1244L) previously associated with CKD (P=0.01, odds ratio=7.95; 95% CI, 1.53 to 41.46) that we find to be present in East Asian individuals at an allele frequency of 0.0027. We determined that P1244L attenuates the interaction of SHROOM3 with 14-3-3, suggesting alterations to the Hippo pathway, a known mediator of CKD.Conclusions These data demonstrate multiple new SHROOM3-dependent genetic/molecular mechanisms that likely affect CKD.

Identification of key microRNAs, transcription factors and genes associated with congenital obstructive nephropathy in a mouse model of megabladder

OBJECTIVE

The present study aimed to investigate the molecular mechanism underlying congenital obstructive nephropathy (CON).

METHODS

The microarray dataset GSE70879 was downloaded from the Gene Expression Omnibus, including 3 kidney samples of megabladder mice and 4 control kidneys. Using this dataset, differentially expressed miRNAs (DEMs) were identified between the kidney samples from megabladder mice and controls, followed by identification of the target genes for these DEMs and construction of a DEM and target gene interaction network. Additionally, the target genes were subjected to Gene Ontology and pathway enrichment analyses, and were used for construction of a protein-protein interaction (PPI) network. Finally, regulatory networks were constructed to analyze transcription factors for the key miRNAs.

RESULTS

From 17 DEMs identified between kidney samples of megabladder mice and controls, 3 key miRNAs were screened, including mmu-miR-150-5p, mmu-miR-374b-5p and mmu-miR-126a-5p. The regulatory networks identified vascular endothelial growth factor A (Vegfa) as the common target gene of mmu-miR-150-5p and five transcription factors, including nuclear receptor subfamily 4, group A, member 2 (Nr4a2), Jun dimerisation protein 2 (Jdp2), Kruppel-like factor 6 (Klf6), Neurexophilin-3 (Nxph3) and RNA binding motif protein 17 (Rbm17). The gene encoding phosphatase and tensin homolog (Pten) was found to be co-regulated by mmu-miR-374b-5p and high mobility group protein A1 (Hmga1), whereas the kirsten rat sarcoma viral oncogene (Kras) was identified as a common target gene of mmu-miR-126a-5p and paired box 6 (Pax6).

CONCLUSIONS

In summary, the above-listed key miRNAs, transcription factors and key genes may be involved in the development of CON.

Advances in Computational Analysis of Metabolomic NMR Data

In this chapter we discuss some of the more recent developments in preprocessing and statistical analysis of NMR spectra in metabolomics. Bayesian methods for analyzing NMR spectra are summarized and we describe one particular approach, BATMAN, in more detail. We consider techniques based on statistical associations, such as correlation spectroscopy (e.g. STOCSY and recent variants), as well as approaches that model the associations as a network and how these change under different biological conditions. The link between metabolism and genotype is explored by looking at metabolic GWAS and related techniques. Finally, we describe the relevance and current status of data standards for NMR metabolomics.

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.

RNA-seq of serial kidney biopsies obtained during progression of chronic kidney disease from dogs with X-linked hereditary nephropathy

Dogs with X-linked hereditary nephropathy (XLHN) have a glomerular basement membrane defect that leads to progressive juvenile-onset renal failure. Their disease is analogous to Alport syndrome in humans, and they also serve as a good model of progressive chronic kidney disease (CKD). However, the gene expression profile that affects progression in this disease has only been partially characterized. To help fill this gap, we used RNA sequencing to identify differentially expressed genes (DEGs), over-represented pathways, and upstream regulators that contribute to kidney disease progression. Total RNA from kidney biopsies was isolated at 3 clinical time points from 3 males with rapidly-progressing CKD, 3 males with slowly-progressing CKD, and 2 age-matched controls. We identified 70 DEGs by comparing rapid and slow groups at specific time points. Based on time course analysis, 1,947 DEGs were identified over the 3 time points revealing upregulation of inflammatory pathways: integrin signaling, T cell activation, and chemokine and cytokine signaling pathways. T cell infiltration was verified by immunohistochemistry. TGF-β1 was identified as the primary upstream regulator. These results provide new insights into the underlying molecular mechanisms of disease progression in XLHN, and the identified DEGs can be potential biomarkers and therapeutic targets translatable to all CKDs.

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.

Genome-Wide Analysis Studies and Chronic Kidney Disease

In recent years, the very high worldwide prevalence of chronic kidney disease (CKD) has led some authors to talk of an "epidemic." The progression of CKD varies considerably among individuals despite similar aetiologies, optimal blood pressure, and glycaemic control. Over the last decade, through genome-wide association studies (GWAS), more than 50 genetic loci have been identified in association with CKD. Understanding the genetic basis of CKD could provide a better knowledge of the biology of the involved pathways, thus potentially leading to novel tools for the diagnosis, prevention, and therapy of CKD. In this review, we will analyse the role of GWAS in the study of CKD.

Identification of CDC42BPG as a novel susceptibility locus for hyperuricemia in a Japanese population

Chronic kidney disease and hyperuricemia are serious global health problems. Recent genome-wide association studies have identified various genetic variants related to these disorders. However, most studies have been conducted in a cross-sectional manner. To identify novel susceptibility loci for chronic kidney disease or hyperuricemia, we performed longitudinal exome-wide association studies (EWASs), using ~ 244,000 genetic variants and clinical data of Japanese individuals who had undergone annual health checkups for several years. After establishing quality controls, the association of renal function-related traits in 5648 subjects (excluding patients with dialysis and population outliers) with 24,579 single nucleotide variants (SNVs) for three genetic models (P < 3.39 × 10^− 7) was tested using generalized estimating equation models. The longitudinal EWASs revealed novel relations of five SNVs to renal function-related traits. Cross-sectional data for renal function-related traits in 7699 Japanese subjects were examined in a replication study. Among the five SNVs, rs55975541 in CDC42BPG was significantly (P < 4.90 × 10^− 4) related to the serum concentration of uric acid in the replication cohort. We also examined the SNVs detected in our longitudinal EWASs with the information on P values in GKDGEN meta-analysis data. Four SNVs in SLC15A2 were significantly associated with the estimated glomerular filtration rate in European ancestry populations, although these SNVs were related to the serum concentration of uric acid with borderline significance in our longitudinal EWASs. Our findings indicate that CDC42BPG may be a novel susceptibility locus for hyperuricemia.

Chromatin Conformation Links Distal Target Genes to CKD Loci

Genome-wide association studies (GWASs) have identified many genetic risk factors for CKD. However, linking common variants to genes that are causal for CKD etiology remains challenging. By adapting self-transcribing active regulatory region sequencing, we evaluated the effect of genetic variation on DNA regulatory elements (DREs). Variants in linkage with the CKD-associated single-nucleotide polymorphism rs11959928 were shown to affect DRE function, illustrating that genes regulated by DREs colocalizing with CKD-associated variation can be dysregulated and therefore, considered as CKD candidate genes. To identify target genes of these DREs, we used circular chromosome conformation capture (4C) sequencing on glomerular endothelial cells and renal tubular epithelial cells. Our 4C analyses revealed interactions of CKD-associated susceptibility regions with the transcriptional start sites of 304 target genes. Overlap with multiple databases confirmed that many of these target genes are involved in kidney homeostasis. Expression quantitative trait loci analysis revealed that mRNA levels of many target genes are genotype dependent. Pathway analyses showed that target genes were enriched in processes crucial for renal function, identifying dysregulated geranylgeranyl diphosphate biosynthesis as a potential disease mechanism. Overall, our data annotated multiple genes to previously reported CKD-associated single-nucleotide polymorphisms and provided evidence for interaction between these loci and target genes. This pipeline provides a novel technique for hypothesis generation and complements classic GWAS interpretation. Future studies are required to specify the implications of our dataset and further reveal the complex roles that common variants have in complex diseases, such as CKD.

Genetic epidemiology in kidney disease

Familial aggregation of chronic kidney disease and its component phenotypes-reduced glomerular filtration rate, proteinuria and renal histologic changes-has long been recognized. Rates of severe kidney disease are also known to differ markedly between populations based on ancestry. These epidemiologic observations support the existence of nephropathy susceptibility genes. Several molecular genetic technologies are now available to identify causative loci. The present article summarizes available strategies useful for identifying nephropathy susceptibility genes, including candidate gene association, family-based linkage, genome-wide association and admixture mapping (mapping by admixture linkage disequilibrium) approaches. Examples of loci detected using these techniques are provided. Epigenetic studies and future directions are also discussed. The identification of nephropathy susceptibility genes, coupled with modifiable environmental triggers impacting their function, is likely to improve risk prediction and transform care. Development of novel therapies to prevent progression of kidney disease will follow.

Uromodulin: from physiology to rare and complex kidney disorders

Uromodulin (also known as Tamm-Horsfall protein) is exclusively produced in the kidney and is the most abundant protein in normal urine. The function of uromodulin remains elusive, but the available data suggest that this protein might regulate salt transport, protect against urinary tract infection and kidney stones, and have roles in kidney injury and innate immunity. Interest in uromodulin was boosted by genetic studies that reported involvement of the UMOD gene, which encodes uromodulin, in a spectrum of rare and common kidney diseases. Rare mutations in UMOD cause autosomal dominant tubulointerstitial kidney disease (ADTKD), which leads to chronic kidney disease (CKD). Moreover, genome-wide association studies have identified common variants in UMOD that are strongly associated with risk of CKD and also with hypertension and kidney stones in the general population. These findings have opened up a new field of kidney research. In this Review we summarize biochemical, physiological, genetic and pathological insights into the roles of uromodulin; the mechanisms by which UMOD mutations cause ADTKD, and the association of common UMOD variants with complex disorders.

Identification of C21orf59 and ATG2A as novel determinants of renal function-related traits in Japanese by exome-wide association studies

We have performed exome-wide association studies to identify genetic variants that influence renal function-related traits or confer susceptibility to chronic kidney disease or hyperuricemia in Japanese. Exome-wide association studies for estimated glomerular filtration rate and the serum concentration of creatinine were performed with 12,565 individuals, that for the serum concentration of uric acid with 9934 individuals, and those for chronic kidney disease or hyperuricemia with 5161 individuals (3270 cases, 1891 controls) or 11,686 individuals (2045 cases, 9641 controls), respectively. The relation of genotypes of single nucleotide polymorphisms to estimated glomerular filtration rate or the serum concentrations of creatinine or uric acid was examined by linear regression analysis, and that of allele frequencies of single nucleotide polymorphisms to chronic kidney disease or hyperuricemia was examined with Fisher's exact test. The exome-wide association studies revealed that 25, seven, and six single nucleotide polymorphisms were significantly (P <1.21 × 10-6) associated with estimated glomerular filtration rate or the serum concentrations of creatinine or uric acid, respectively, and that 49 and 35 polymorphisms were significantly associated with chronic kidney disease or hyperuricemia, respectively. Subsequent multivariable logistic regression analysis with adjustment for covariates revealed that four and three single nucleotide polymorphisms were related (P < 0.05) to chronic kidney disease or hyperuricemia, respectively. Among polymorphisms identified in the present study, rs76974938 [C/T (D67N)] of C21orf59 and rs188780113 [G/A (R478C)] of ATG2A may be novel determinants of estimated glomerular filtration rate and chronic kidney disease or of the serum concentration of uric acid, respectively.

Cystatin C and Cardiovascular Disease: A Mendelian Randomization Study

BACKGROUND

Epidemiological studies show that high circulating cystatin C is associated with risk of cardiovascular disease (CVD), independent of creatinine-based renal function measurements. It is unclear whether this relationship is causal, arises from residual confounding, and/or is a consequence of reverse causation.

OBJECTIVES

The aim of this study was to use Mendelian randomization to investigate whether cystatin C is causally related to CVD in the general population.

METHODS

We incorporated participant data from 16 prospective cohorts (n = 76,481) with 37,126 measures of cystatin C and added genetic data from 43 studies (n = 252,216) with 63,292 CVD events. We used the common variant rs911119 in CST3 as an instrumental variable to investigate the causal role of cystatin C in CVD, including coronary heart disease, ischemic stroke, and heart failure.

RESULTS

Cystatin C concentrations were associated with CVD risk after adjusting for age, sex, and traditional risk factors (relative risk: 1.82 per doubling of cystatin C; 95% confidence interval [CI]: 1.56 to 2.13; p = 2.12 × 10⁻¹⁴). The minor allele of rs911119 was associated with decreased serum cystatin C (6.13% per allele; 95% CI: 5.75 to 6.50; p = 5.95 × 10⁻²¹¹), explaining 2.8% of the observed variation in cystatin C. Mendelian randomization analysis did not provide evidence for a causal role of cystatin C, with a causal relative risk for CVD of 1.00 per doubling cystatin C (95% CI: 0.82 to 1.22; p = 0.994), which was statistically different from the observational estimate (p = 1.6 × 10⁻⁵). A causal effect of cystatin C was not detected for any individual component of CVD.

CONCLUSIONS

Mendelian randomization analyses did not support a causal role of cystatin C in the etiology of CVD. As such, therapeutics targeted at lowering circulating cystatin C are unlikely to be effective in preventing CVD.

Predicting acute kidney injury: current status and future challenges

Acute kidney injury (AKI) is characterized by an acute decline in renal function and is associated to increased mortality rate, hospitalization time, and total health-related costs. The severity of this ‘fearsome’ clinical complication might depend on, or even be worsened by, the late detection of AKI, when the diagnosis is based on the elevation of serum creatinine (SCr). For these reasons, in recent years a great number of new tools, biomarkers and predictive models have been proposed to clinicians in order to improve diagnosis and prevent the development of AKI. The purpose of this narrative paper is to review the current state of the art in prediction and early detection of AKI and outline future challenges.

Development of next generation sequencing panel for UMOD and association with kidney disease

Chronic kidney disease (CKD) has a prevalence of approximately 10% in adult populations. CKD can progress to end-stage renal disease (ESRD) and this is usually fatal unless some form of renal replacement therapy (chronic dialysis or renal transplantation) is provided. There is an inherited predisposition to CKD with several genetic risk markers now identified. The UMOD gene has been associated with CKD of varying aetiologies. An AmpliSeq next generation sequencing panel was developed to facilitate comprehensive sequencing of the UMOD gene, covering exonic and regulatory regions. SNPs and CpG sites in the genomic region encompassing UMOD were evaluated for association with CKD in two studies; the UK Wellcome Trust Case-Control 3 Renal Transplant Dysfunction Study (n = 1088) and UK-ROI GENIE GWAS (n = 1726). A technological comparison of two Ion Torrent machines revealed 100% allele call concordance between S5 XL™ and PGM™ machines. One SNP (rs183962941), located in a non-coding region of UMOD, was nominally associated with ESRD (p = 0.008). No association was identified between UMOD variants and estimated glomerular filtration rate. Analysis of methylation data for over 480,000 CpG sites revealed differential methylation patterns within UMOD, the most significant of these was cg03140788 p = 3.7 x 10-10.

1000 Genomes-based meta-analysis identifies 10 novel loci for kidney function

HapMap imputed genome-wide association studies (GWAS) have revealed >50 loci at which common variants with minor allele frequency >5% are associated with kidney function. GWAS using more complete reference sets for imputation, such as those from The 1000 Genomes project, promise to identify novel loci that have been missed by previous efforts. To investigate the value of such a more complete variant catalog, we conducted a GWAS meta-analysis of kidney function based on the estimated glomerular filtration rate (eGFR) in 110,517 European ancestry participants using 1000 Genomes imputed data. We identified 10 novel loci with p-value < 5 × 10^-8 previously missed by HapMap-based GWAS. Six of these loci (HOXD8, ARL15, PIK3R1, EYA4, ASTN2, and EPB41L3) are tagged by common SNPs unique to the 1000 Genomes reference panel. Using pathway analysis, we identified 39 significant (FDR < 0.05) genes and 127 significantly (FDR < 0.05) enriched gene sets, which were missed by our previous analyses. Among those, the 10 identified novel genes are part of pathways of kidney development, carbohydrate metabolism, cardiac septum development and glucose metabolism. These results highlight the utility of re-imputing from denser reference panels, until whole-genome sequencing becomes feasible in large samples.

Effect of UMOD genotype on long-term graft survival after kidney transplantation in patients treated with cyclosporine-based therapy

The genetic rs12917707-G>T variant in uromodulin (UMOD) has been associated with renal function, chronic kidney disease and hypertension with the minor T-allele showing a protective effect. Hypertension and nephrotoxicity are adverse effects of chronic cyclosporine treatment. We tested whether UMOD rs12917707-T in donor kidneys associates with long-term graft survival in 393 Caucasian patients with stable graft function for more than 10 weeks after kidney transplantation treated with a cyclosporine-based maintenance therapy (mean graft survival 9 years). Presence of the donor T-allele had no effect on blood pressure, serum creatinine 1 year after transplantation, and on number of acute graft rejections during the first year. No significant effect on overall graft survival was observed in Kaplan-Meier analysis (P=0.65). In death-censored adjusted multivariate analysis, presence of donor T-allele associated with a significant lower hazard ratio of 0.67 (95% confidence interval: 0.46-0.97, P=0.05) for graft loss. This protective effect of the donor T-allele on graft loss observed in multivariate adjusted analysis justifies further investigations including patients treated with similar or other immunosuppressive regimens.

Missing heritability: is the gap closing? An analysis of 32 complex traits in the Lifelines Cohort Study

Despite the recent explosive rise in number of genetic markers for complex disease traits identified in genome-wide association studies, there is still a large gap between the known heritability of these traits and the part explained by these markers. To gauge whether this 'heritability gap' is closing, we first identified genome-wide significant SNPs from the literature and performed replication analyses for 32 highly relevant traits from five broad disease areas in 13 436 subjects of the Lifelines Cohort. Next, we calculated the variance explained by multi-SNP genetic risk scores (GRSs) for each trait, and compared it to their broad- and narrow-sense heritabilities captured by all common SNPs. The majority of all previously-associated SNPs (median=75%) were significantly associated with their respective traits. All GRSs were significant, with unweighted GRSs generally explaining less phenotypic variance than weighted GRSs, for which the explained variance was highest for height (15.5%) and varied between 0.02 and 6.7% for the other traits. Broad-sense common-SNP heritability estimates were significant for all traits, with the additive effect of common SNPs explaining 48.9% of the variance for height and between 5.6 and 39.2% for the other traits. Dominance effects were uniformly small (0-1.5%) and not significant. On average, the variance explained by the weighted GRSs accounted for only 10.7% of the common-SNP heritability of the 32 traits. These results indicate that GRSs may not yet be ready for accurate personalized prediction of complex disease traits limiting widespread adoption in clinical practice.

Genetic variation underlying renal uric acid excretion in Hispanic children: the Viva La Familia Study

Background

Reduced renal excretion of uric acid plays a significant role in the development of hyperuricemia and gout in adults. Hyperuricemia has been associated with chronic kidney disease and cardiovascular disease in children and adults. There are limited genome-wide association studies associating genetic polymorphisms with renal urate excretion measures. Therefore, we investigated the genetic factors that influence the excretion of uric acid and related indices in 768 Hispanic children of the Viva La Familia Study.

Methods

We performed a genome-wide association analysis for 24-h urinary excretion measures such as urinary uric acid/urinary creatinine ratio, uric acid clearance, fractional excretion of uric acid, and glomerular load of uric acid in SOLAR, while accounting for non-independence among family members.

Results

All renal urate excretion measures were significantly heritable (p <2 × 10⁻⁶) and ranged from 0.41 to 0.74. Empirical threshold for genome-wide significance was set at p <1 × 10⁻⁷. We observed a strong association (p < 8 × 10⁻⁸) of uric acid clearance with a single nucleotide polymorphism (SNP) in zinc finger protein 446 (ZNF446) (rs2033711 (A/G), MAF: 0.30). The minor allele (G) was associated with increased uric acid clearance. Also, we found suggestive associations of uric acid clearance with SNPs in ZNF324, ZNF584, and ZNF132 (in a 72 kb region of 19q13; p <1 × 10⁻⁶, MAFs: 0.28–0.31).

Conclusion

For the first time, we showed the importance of 19q13 region in the regulation of renal urate excretion in Hispanic children. Our findings indicate differences in inherent genetic architecture and shared environmental risk factors between our cohort and other pediatric and adult populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12881-016-0366-3) contains supplementary material, which is available to authorized users.

A single nucleotide polymorphism in the UMOD promoter is associated with end stage renal disease

BACKGROUND

Several genome-wide association studies revealed that several variants of UMOD gene were related to the estimated glomerular filtration rate (eGFR), CKD or hypertension. In this study, we investigated the association between a common variant rs13333226 in the promoter region of UMOD gene and end stage renal disease (ESRD).

METHODS

Variant rs13333226 of UMOD gene was genotyped by using the ABI Real time TaqMan allelic discrimination assay in a case-control study including 638 unrelated patients with ESRD and 366 controls.

RESULTS

The frequency of UMOD SNP rs13333226 GG/GA genotype was significantly higher (36.83% vs. 20.22%, P = 4.02 × 10^-8) and the frequency of G allele was much higher (19.04% vs. 11.20%, P = 4.00 × 10^-6) in the patients with ESRD than in the controls. The G allele was associated with an increased risk of ESRD (odds ratio 2.30, 95% confidence interval 1.70-3.11, P = 6.10 × 10^-8). And G allele (odds ratio 2.33, 95% confidence interval 1.32-4.13, P = 3.65 × 10^-3) was associated independently with ESRD.

CONCLUSIONS

A common variation rs13333226 in the promoter region of UMOD gene was independently associated with ESRD in Han Chinese.

Pharmacogenetics of lithium effects on glomerular function in bipolar disorder patients under chronic lithium treatment: a pilot study

Bipolar disorder (BD) is a psychiatric disease characterized by alternating episodes of mania and depression. Lithium (Li) represents the mainstay treatment for BD, although a significant proportion of patients shows insufficient or no response. Li is also associated with potentially severe side effects, including renal effects. Several studies reported that Li may induce reduction of glomerular filtration rate (GFR) in patients under long-term treatment. The biological systems and the genetic factors involved in susceptibility to Li-induced renal-side effects have been scarcely explored. The aim of our study was to test the contribution of putatively risk genetic variants in Li-induced reduction in estimated GFR (eGFR) in BD patients under long-term Li treatment. Tag SNPs, located in genes previously shown to be associated with kidney dysfunction or Li mechanism of action, were selected and genotyped in a sample of 70 BD patients of Sardinian origin. SNP rs378448, located in Acid Sensing Ion Channel Neurona-1 (ACCN1) gene, showed a significant interaction with duration of Li treatment on eGFR (F₂=3.623, p=0.033). Our preliminary findings suggest that rs378448 could predispose BD subjects to a detrimental effect of chronic Li treatment on kidney functioning.

Single nucleotide polymorphisms in the intergenic region between metformin transporter OCT2 and OCT3 coding genes are associated with short-term response to metformin monotherapy in type 2 diabetes mellitus patients

OBJECTIVES

High variability in clinical response to metformin is often observed in type 2 diabetes (T2D) patients, and it highlights the need for identification of genetic components affecting the efficiency of metformin therapy. Aim of this observational study is to evaluate the role of tagSNPs (tagging single nucleotide polymorphisms) from genomic regions coding for six metformin transporter genes with respect to the short-term efficiency.

DESIGN

102 tagSNPs in 6 genes coding for metformin transporters were genotyped in the group of 102 T2D patients treated with metformin for 3 months.

METHODS

Most significant hits were analyzed in the group of 131 T2D patients from Slovakia. Pharmacokinetic study in 25 healthy nondiabetic volunteers was conducted to investigate the effects of identified polymorphisms.

RESULTS

In the discovery group of 102 patients, minor alleles of rs3119309, rs7757336 and rs2481030 were significantly nominally associated with metformin inefficiency (P = 1.9 × 10^-6 to 8.1 × 10^-6). Effects of rs2481030 and rs7757336 did not replicate in the group of 131 T2DM patients from Slovakia alone, whereas rs7757336 was significantly associated with a reduced metformin response in combined group. In pharmacokinetic study, group of individuals harboring risk alleles of rs7757336 and rs2481030 displayed significantly reduced AUC_∞ of metformin in plasma.

CONCLUSIONS

For the first time, we have identified an association between the lack of metformin response and SNPs rs3119309 and rs7757336 located in the 5' flanking region of the genes coding for Organic cation transporter 2 and rs2481030 located in the 5' flanking region of Organic cation transporter 3 that was supported by the results of a pharmacokinetic study on 25 healthy volunteers.

A principal component meta-analysis on multiple anthropometric traits identifies novel loci for body shape

Large consortia have revealed hundreds of genetic loci associated with anthropometric traits, one trait at a time. We examined whether genetic variants affect body shape as a composite phenotype that is represented by a combination of anthropometric traits. We developed an approach that calculates averaged PCs (AvPCs) representing body shape derived from six anthropometric traits (body mass index, height, weight, waist and hip circumference, waist-to-hip ratio). The first four AvPCs explain >99% of the variability, are heritable, and associate with cardiometabolic outcomes. We performed genome-wide association analyses for each body shape composite phenotype across 65 studies and meta-analysed summary statistics. We identify six novel loci: LEMD2 and CD47 for AvPC1, RPS6KA5/C14orf159 and GANAB for AvPC3, and ARL15 and ANP32 for AvPC4. Our findings highlight the value of using multiple traits to define complex phenotypes for discovery, which are not captured by single-trait analyses, and may shed light onto new pathways.

Familial aggregation of albuminuria and arterial hypertension in an Aboriginal Australian community and the contribution of variants in ACE and TP53

BACKGROUND

Aboriginal Australians are at high risk of cardiovascular, metabolic and renal diseases, resulting in a marked reduction in life expectancy when compared to the rest of the Australian population. This is partly due to recognized environmental and lifestyle risk factors, but a contribution of genetic susceptibility is also likely.

METHODS

Using results from a comprehensive survey of one community (N = 1350 examined individuals), we have tested for familial aggregation of plasma glucose, arterial blood pressure, albuminuria (measured as urinary albumin to creatinine ratio, UACR) and estimated glomerular filtration rate (eGFR), and quantified the contribution of variation at four candidate genes (ACE; TP53; ENOS3; MTHFR).

RESULTS

In the subsample of 357 individuals with complete genotype and phenotype data we showed that both UACR (h2 = 64%) and blood pressure (sBP h2 = 29%, dBP, h2 = 11%) were significantly heritable. The ACE insertion-deletion (P = 0.0009) and TP53 codon72 polymorphisms (P = 0.003) together contributed approximately 15% of the total heritability of UACR, with an effect of ACE genotype on BP also clearly evident.

CONCLUSIONS

While the effects of the ACE insertion-deletion on risk of renal disease (especially in the setting of diabetes) are well recognized, this is only the second study to implicate p53 genotype as a risk factor for albuminuria - the other being an earlier study we performed in a different Aboriginal community (McDonald et al., J Am Soc Nephrol 13: 677-83, 2002). We conclude that there are significant genetic contributions to the high prevalence of chronic diseases observed in this population.

Meta-analysis identifies common and rare variants influencing blood pressure and overlapping with metabolic trait loci

Meta-analyses of association results for blood pressure using exome-centric single-variant and gene-based tests identified 31 new loci in a discovery stage among 146,562 individuals, with follow-up and meta-analysis in 180,726 additional individuals (total n = 327,288). These blood pressure-associated loci are enriched for known variants for cardiometabolic traits. Associations were also observed for the aggregation of rare and low-frequency missense variants in three genes, NPR1, DBH, and PTPMT1. In addition, blood pressure associations at 39 previously reported loci were confirmed. The identified variants implicate biological pathways related to cardiometabolic traits, vascular function, and development. Several new variants are inferred to have roles in transcription or as hubs in protein-protein interaction networks. Genetic risk scores constructed from the identified variants were strongly associated with coronary disease and myocardial infarction. This large collection of blood pressure-associated loci suggests new therapeutic strategies for hypertension, emphasizing a link with cardiometabolic risk.

Meta-analysis of rare and common exome chip variants identifies S1PR4 and other loci influencing blood cell traits

Hematologic measures such as hematocrit and white blood cell (WBC) count are heritable and clinically relevant. We analyzed erythrocyte and WBC phenotypes in 52,531 individuals (37,775 of European ancestry, 11,589 African Americans, and 3,167 Hispanic Americans) from 16 population-based cohorts with Illumina HumanExome BeadChip genotypes. We then performed replication analyses of new discoveries in 18,018 European-American women and 5,261 Han Chinese. We identified and replicated four new erythrocyte trait-locus associations (CEP89, SHROOM3, FADS2, and APOE) and six new WBC loci for neutrophil count (S1PR4), monocyte count (BTBD8, NLRP12, and IL17RA), eosinophil count (IRF1), and total WBC count (MYB). The association of a rare missense variant in S1PR4 supports the role of sphingosine-1-phosphate signaling in leukocyte trafficking and circulating neutrophil counts. Loss-of-function experiments for S1pr4 in mouse and s1pr4 in zebrafish demonstrated phenotypes consistent with the association observed in humans and altered kinetics of neutrophil recruitment and resolution in response to tissue injury.

Insights into kidney diseases from genome-wide association studies

Over the past decade, genome-wide association studies (GWAS) have considerably improved our understanding of the genetic basis of kidney function and disease. Population-based studies, used to investigate traits that define chronic kidney disease (CKD), have identified >50 genomic regions in which common genetic variants associate with estimated glomerular filtration rate or urinary albumin-to-creatinine ratio. Case-control studies, used to study specific CKD aetiologies, have yielded risk loci for specific kidney diseases such as IgA nephropathy and membranous nephropathy. In this Review, we summarize important findings from GWAS and clinical and experimental follow-up studies. We also compare risk allele frequency, effect sizes, and specificity in GWAS of CKD-defining traits and GWAS of specific CKD aetiologies and the implications for study design. Genomic regions identified in GWAS of CKD-defining traits can contain causal genes for monogenic kidney diseases. Population-based research on kidney function traits can therefore generate insights into more severe forms of kidney diseases. Experimental follow-up studies have begun to identify causal genes and variants, which are potential therapeutic targets, and suggest mechanisms underlying the high allele frequency of causal variants. GWAS are thus a useful approach to advance knowledge in nephrology.

Influence of common genetic variants on childhood kidney outcomes

BACKGROUND

Kidney measures in early life are associated with kidney disease in later life. We hypothesized that these associations are partly explained by common genetic variants that lead to both smaller kidneys with lower kidney function in early childhood and kidney disease in adulthood.

METHODS

We examined in a population-based prospective cohort study among 4,119 children the associations of a weighted genetic risk score combining 20 previously identified common genetic variants related to adult eGFRcreat with kidney outcomes in children aged 6.0 years (95% range 5.7-7.8). Childhood kidney outcomes included combined kidney volume, glomerular filtration rate (eGFR) based on creatinine levels, and microalbuminuria based on albumin and creatinine urine levels.

RESULTS

We observed that the genetic risk score based on variants related to impaired kidney function in adults was associated with a smaller combined kidney volume (P value 3.0 × 10(-3)) and with a lower eGFR (P value 4.0 × 10(-4)) in children. The genetic risk score was not associated with microalbuminuria.

CONCLUSION

Common genetic variants related to impaired kidney function in adults already lead to subclinical changes in childhood kidney outcomes. The well-known associations of kidney measures in early life with kidney disease in later life may at least be partly explained by common genetic variants.