Genome-wide meta-analysis and omics integration identifies novel genes associated with diabetic kidney disease

Multi-omics characterization of type 2 diabetes associated genetic variation

Discerning the mechanisms driving type 2 diabetes (T2D) pathophysiology from genome-wide association studies (GWAS) remains a challenge. To this end, we integrated omics information from 16 multi-tissue and multi-ancestry expression, protein, and metabolite quantitative trait loci (QTL) studies and 46 multi-ancestry GWAS for T2D-related traits with the largest, most ancestrally diverse T2D GWAS to date. Of the 1,289 T2D GWAS index variants, 716 (56%) demonstrated strong evidence of colocalization with a molecular or T2D-related trait, implicating 657 cis-effector genes, 1,691 distal-effector genes, 731 metabolites, and 43 T2D-related traits. We identified 773 of these cis- and distal-effector genes using either expression QTL data from understudied ancestry groups or inclusion of T2D index variants enriched in underrepresented populations, emphasizing the value of increasing population diversity in functional mapping. Linking these variants, genes, metabolites, and traits into a network, we elucidated mechanisms through which T2D-associated variation may impact disease risk. Finally, we showed that drugs targeting effector proteins were enriched in those approved to treat T2D, highlighting the potential of these results to prioritize drug targets for T2D. These results represent a leap in the molecular characterization of T2D-associated genetic variation and will aid in translating genetic findings into novel therapeutic strategies.

Diabetes and Renal Complications: An Overview on Pathophysiology, Biomarkers and Therapeutic Interventions

Diabetic kidney disease (DKD) is a major microvascular complication of both type 1 and type 2 diabetes. DKD is characterised by injury to both glomerular and tubular compartments, leading to kidney dysfunction over time. It is one of the most common causes of chronic kidney disease (CKD) and end-stage renal disease (ESRD). Persistent high blood glucose levels can damage the small blood vessels in the kidneys, impairing their ability to filter waste and fluids from the blood effectively. Other factors like high blood pressure (hypertension), genetics, and lifestyle habits can also contribute to the development and progression of DKD. The key features of renal complications of diabetes include morphological and functional alterations to renal glomeruli and tubules leading to mesangial expansion, glomerulosclerosis, homogenous thickening of the glomerular basement membrane (GBM), albuminuria, tubulointerstitial fibrosis and progressive decline in renal function. In advanced stages, DKD may require treatments such as dialysis or kidney transplant to sustain life. Therefore, early detection and proactive management of diabetes and its complications are crucial in preventing DKD and preserving kidney function.

Novel genetic markers for chronic kidney disease in a geographically isolated population of Indigenous Australians: Individual and multiple phenotype genome-wide association study

BACKGROUND

Chronic kidney disease (CKD) is highly prevalent among Indigenous Australians, especially those in remote regions. The Tiwi population has been isolated from mainland Australia for millennia and exhibits unique genetic characteristics that distinguish them from other Indigenous and non-Indigenous populations. Notably, the rate of end-stage renal disease is up to 20 times greater in this population compared to non-Indigenous populations. Despite the identification of numerous genetic loci associated with kidney disease through GWAS, the Indigenous population such as Tiwi remains severely underrepresented and the increased prevalence of CKD in this population may be due to unique disease-causing alleles/genes.

METHODS

We used albumin-to-creatinine ratio (ACR) and estimated glomerular filtration rate (eGFR) to estimate the prevalence of kidney disease in the Tiwi population (N = 492) in comparison to the UK Biobank (UKBB) (N = 134,724) database. We then performed an exploratory factor analysis to identify correlations among 10 CKD-related phenotypes and identify new multi-phenotype factors. We subsequently conducted a genome-wide association study (GWAS) on all single and multiple phenotype factors using mixed linear regression models, adjusted for age, sex, population stratification, and genetic relatedness between individuals.

RESULTS

Based on ACR, 20.3% of the population was at severely increased risk of CKD progression and showed elevated levels of ACR compared to the UKBB population independent of HbA1c. A GWAS of ACR revealed novel association loci in the genes MEG3 (chr14:100812018:T:A), RAB36 (rs11704318), and TIAM2 (rs9689640). Additionally, multiple phenotypes GWAS of ACR, eGFR, urine albumin, and serum creatinine identified a novel variant that mapped to the gene MEIS2 (chr15:37218869:A:G). Most of the identified variants were found to be either absent or rare in the UKBB population.

CONCLUSIONS

Our study highlights the Tiwi population's predisposition towards elevated ACR, and the collection of novel genetic variants associated with kidney function. These associations may prove valuable in the early diagnosis and treatment of renal disease in this underrepresented population. Additionally, further research is needed to comprehensively validate the functions of the identified variants/genes.

PharmGWAS: a GWAS-based knowledgebase for drug repurposing

Leveraging genetics insights to promote drug repurposing has become a promising and active strategy in pharmacology. Indeed, among the 50 drugs approved by FDA in 2021, two-thirds have genetically supported evidence. In this regard, the increasing amount of widely available genome-wide association studies (GWAS) datasets have provided substantial opportunities for drug repurposing based on genetics discoveries. Here, we developed PharmGWAS, a comprehensive knowledgebase designed to identify candidate drugs through the integration of GWAS data. PharmGWAS focuses on novel connections between diseases and small-molecule compounds derived using a reverse relationship between the genetically-regulated expression signature and the drug-induced signature. Specifically, we collected and processed 1929 GWAS datasets across a diverse spectrum of diseases and 724 485 perturbation signatures pertaining to a substantial 33609 molecular compounds. To obtain reliable and robust predictions for the reverse connections, we implemented six distinct connectivity methods. In the current version, PharmGWAS deposits a total of 740 227 genetically-informed disease-drug pairs derived from drug-perturbation signatures, presenting a valuable and comprehensive catalog. Further equipped with its user-friendly web design, PharmGWAS is expected to greatly aid the discovery of novel drugs, the exploration of drug combination therapies and the identification of drug resistance or side effects. PharmGWAS is available at https://ngdc.cncb.ac.cn/pharmgwas.

Methylation analysis by targeted bisulfite sequencing in large for gestational age (LGA) newborns: the LARGAN cohort

BACKGROUND

In 1990, David Barker proposed that prenatal nutrition is directly linked to adult cardiovascular disease. Since then, the relationship between adult cardiovascular risk, metabolic syndrome and birth weight has been widely documented. Here, we used the TruSeq Methyl Capture EPIC platform to compare the methylation patterns in cord blood from large for gestational age (LGA) vs adequate for gestational age (AGA) newborns from the LARGAN cohort.

RESULTS

We found 1672 differentially methylated CpGs (DMCs) with a nominal p < 0.05 and 48 differentially methylated regions (DMRs) with a corrected p < 0.05 between the LGA and AGA groups. A systems biology approach identified several biological processes significantly enriched with genes in association with DMCs with FDR < 0.05, including regulation of transcription, regulation of epinephrine secretion, norepinephrine biosynthesis, receptor transactivation, forebrain regionalization and several terms related to kidney and cardiovascular development. Gene ontology analysis of the genes in association with the 48 DMRs identified several significantly enriched biological processes related to kidney development, including mesonephric duct development and nephron tubule development. Furthermore, our dataset identified several DNA methylation markers enriched in gene networks involved in biological pathways and rare diseases of the cardiovascular system, kidneys, and metabolism.

CONCLUSIONS

Our study identified several DMCs/DMRs in association with fetal overgrowth. The use of cord blood as a material for the identification of DNA methylation biomarkers gives us the possibility to perform follow-up studies on the same patients as they grow. These studies will not only help us understand how the methylome responds to continuum postnatal growth but also link early alterations of the DNA methylome with later clinical markers of growth and metabolic fitness.

Kidney fibrosis: Emerging diagnostic and therapeutic strategies

An increasing number of patients worldwide suffers from chronic kidney disease (CKD). CKD is accompanied by kidney fibrosis, which affects all compartments of the kidney, i.e., the glomeruli, tubulointerstitium, and vasculature. Fibrosis is the best predictor of progression of kidney diseases. Currently, there is no specific anti-fibrotic therapy for kidney patients and invasive renal biopsy remains the only option for specific detection and quantification of kidney fibrosis. Here we review emerging diagnostic approaches and potential therapeutic options for fibrosis. We discuss how translational research could help to establish fibrosis-specific endpoints for clinical trials, leading to improved patient stratification and potentially companion diagnostics, and facilitating and optimizing development of novel anti-fibrotic therapies for kidney patients.

Circulating epigenomic biomarkers correspond with kidney disease susceptibility in high-risk populations with type 2 diabetes mellitus

AIMS

To investigate epigenomic indices of diabetic kidney disease (DKD) susceptibility among high-risk populations with type 2 diabetes mellitus.

METHODS

KDIGO (Kidney Disease: Improving Global Outcomes) clinical guidelines were used to classify people living with or without DKD. Differential gene methylation of DKD was then assessed in a discovery Aboriginal Diabetes Study cohort (PROPHECY, 89 people) and an external independent study from Thailand (THEPTARIN, 128 people). Corresponding mRNA levels were also measured and linked to levels of albuminuria and eGFR.

RESULTS

Increased DKD risk was associated with reduced methylation and elevated gene expression in the PROPHECY discovery cohort of Aboriginal Australians and these findings were externally validated in the THEPTARIN diabetes registry of Thai people living with type 2 diabetes mellitus.

CONCLUSIONS

Novel epigenomic scores can improve diagnostic performance over clinical modelling using albuminuria and GFR alone and can distinguish DKD susceptibility.

Genetics of diabetes-associated microvascular complications

Diabetes is associated with excess morbidity and mortality due to both micro- and macrovascular complications, as well as a range of non-classical comorbidities. Diabetes-associated microvascular complications are those considered most closely related to hyperglycaemia in a causal manner. However, some individuals with hyperglycaemia (even those with severe hyperglycaemia) do not develop microvascular diseases, which, together with evidence of co-occurrence of microvascular diseases in families, suggests a role for genetics. While genome-wide association studies (GWASs) produced firm evidence of multiple genetic variants underlying differential susceptibility to type 1 and type 2 diabetes, genetic determinants of microvascular complications are mostly suggestive. Identified susceptibility variants of diabetic kidney disease (DKD) in type 2 diabetes mirror variants underlying chronic kidney disease (CKD) in individuals without diabetes. As for retinopathy and neuropathy, reported risk variants currently lack large-scale replication. The reported associations between type 2 diabetes risk variants and microvascular complications may be explained by hyperglycaemia. More extensive phenotyping, along with adjustments for unmeasured confounding, including both early (fetal) and late-life (hyperglycaemia, hypertension, etc.) environmental factors, are urgently needed to understand the genetics of microvascular complications. Finally, genetic variants associated with reduced glycolysis, mitochondrial dysfunction and DNA damage and sustained cell regeneration may protect against microvascular complications, illustrating the utility of studies in individuals who have escaped these complications.

The Microbiome, Epigenome, and Diet in Adults with Obesity during Behavioral Weight Loss

Obesity has been linked to the gut microbiome, epigenome, and diet, yet these factors have not been studied together during obesity treatment. Our objective was to evaluate associations among gut microbiota (MB), DNA methylation (DNAme), and diet prior to and during a behavioral weight loss intervention. Adults (n = 47, age 40.9 ± 9.7 years, body mass index (BMI) 33.5 ± 4.5 kg/m2, 77% female) with data collected at baseline (BL) and 3 months (3 m) were included. Fecal MB was assessed via 16S sequencing and whole blood DNAme via the Infinium EPIC array. Food group and nutrient intakes and Healthy Eating Index (HEI) scores were calculated from 7-day diet records. Linear models were used to test for the effect of taxa relative abundance on DNAme and diet cross-sectionally at each time point, adjusting for confounders and a false discovery rate of 5%. Mean weight loss was 6.2 ± 3.9% at 3 m. At BL, one MB taxon, Ruminiclostridium, was associated with DNAme of the genes COL20A1 (r = 0.651, p = 0.029), COL18A1 (r = 0.578, p = 0.044), and NT5E (r = 0.365, p = 0.043). At 3 m, there were 14 unique MB:DNAme associations, such as Akkermansia with DNAme of GUSB (r = -0.585, p = 0.003), CRYL1 (r = -0.419, p = 0.007), C9 (r = -0.439, p = 0.019), and GMDS (r = -0.559, p = 0.046). Among taxa associated with DNAme, no significant relationships were seen with dietary intakes of relevant nutrients, food groups, or HEI scores. Our findings indicate that microbes linked to mucin degradation, short-chain fatty acid production, and body weight are associated with DNAme of phenotypically relevant genes. These relationships offer an initial understanding of the possible routes by which alterations in gut MB may influence metabolism during weight loss.

Genetic and epigenetic background of diabetic kidney disease

Diabetic kidney disease (DKD) is a severe diabetic complication that affects up to half of the individuals with diabetes. Elevated blood glucose levels are a key underlying cause of DKD, but DKD is a complex multifactorial disease, which takes years to develop. Family studies have shown that inherited factors also contribute to the risk of the disease. During the last decade, genome-wide association studies (GWASs) have emerged as a powerful tool to identify genetic risk factors for DKD. In recent years, the GWASs have acquired larger number of participants, leading to increased statistical power to detect more genetic risk factors. In addition, whole-exome and whole-genome sequencing studies are emerging, aiming to identify rare genetic risk factors for DKD, as well as epigenome-wide association studies, investigating DNA methylation in relation to DKD. This article aims to review the identified genetic and epigenetic risk factors for DKD.

Multi-omic single cell sequencing: Overview and opportunities for kidney disease therapeutic development

Single cell sequencing technologies have rapidly advanced in the last decade and are increasingly applied to gain unprecedented insights by deconstructing complex biology to its fundamental unit, the individual cell. First developed for measurement of gene expression, single cell sequencing approaches have evolved to allow simultaneous profiling of multiple additional features, including chromatin accessibility within the nucleus and protein expression at the cell surface. These multi-omic approaches can now further be applied to cells in situ, capturing the spatial context within which their biology occurs. To extract insights from these complex datasets, new computational tools have facilitated the integration of information across different data types and the use of machine learning approaches. Here, we summarize current experimental and computational methods for generation and integration of single cell multi-omic datasets. We focus on opportunities for multi-omic single cell sequencing to augment therapeutic development for kidney disease, including applications for biomarkers, disease stratification and target identification.

Progress in genetics of type 2 diabetes and diabetic complications

Type 2 diabetes results from a complex interaction between genetic and environmental factors. Precision medicine for type 2 diabetes using genetic data is expected to predict the risk of developing diabetes and complications and to predict the effects of medications and life-style intervention more accurately for individuals. Genome-wide association studies (GWAS) have been conducted in European and Asian populations and new genetic loci have been identified that modulate the risk of developing type 2 diabetes. Novel loci were discovered by GWAS in diabetic complications with increasing sample sizes. Large-scale genome-wide association analysis and polygenic risk scores using biobank information is making it possible to predict the development of type 2 diabetes. In the ADVANCE clinical trial of type 2 diabetes, a multi-polygenic risk score was useful to predict diabetic complications and their response to treatment. Proteomics and metabolomics studies have been conducted and have revealed the associations between type 2 diabetes and inflammatory signals and amino acid synthesis. Using multi-omics analysis, comprehensive molecular mechanisms have been elucidated to guide the development of targeted therapy for type 2 diabetes and diabetic complications.

A genome-wide association study identifies a possible role for cannabinoid signalling in the pathogenesis of diabetic kidney disease

Diabetic kidney disease (DKD), also known as diabetic nephropathy, is the leading cause of renal impairment and end-stage renal disease. Patients with diabetes are at risk for DKD because of poor control of their blood glucose, as well as nonmodifiable risk factors including age, ethnicity, and genetics. This genome-wide association study (GWAS) was conducted for the first time in the Emirati population to investigate possible genetic factors associated with the development and progression of DKD. We included data on 7,921,925 single nucleotide polymorphism (SNPs) in 258 cases of type 2 diabetes mellitus (T2DM) who developed DKD and 938 control subjects with T2DM who did not develop DKD. GWAS suggestive results (P < 1 × 10-5) were further replicated using summary statistics from three cohorts with T2DM-induced DKD (Bio Bank Japan data, UK Biobank, and FinnGen Project data) and T1DM-induced DKD (UK-ROI cohort data from Belfast, UK). When conducting a multiple linear regression model for gene-set analyses, the CNR2 gene demonstrated genome-wide significance at 1.46 × 10-6. SNPs in CNR2 gene, encodes cannabinoid receptor 2 or CB2, were replicated in Japanese samples with the leading SNP rs2501391 showing a Pcombined = 9.3 × 10-7, and odds ratio = 0.67 in association with DKD associated with T2DM, but not with T1DM, without any significant association with T2DM itself. The allele frequencies of our cohort and those of the replication cohorts were in most cases markedly different. In addition, we replicated the association between rs1564939 in the GLRA3 gene and DKD in T2DM (P = 0.016, odds ratio = 0.54 per allele C). Our findings suggest evidence that cannabinoid signalling may be involved in the development of DKD through CB2, which is expressed in different kidney regions and known to be involved in insulin resistance, inflammation, and the development of kidney fibrosis.

Precision Medicine in Type 1 Diabetes

Type 1 diabetes is a complex, chronic disease in which the insulin-producing beta cells in the pancreas are sufficiently altered or impaired to result in requirement of exogenous insulin for survival. The development of type 1 diabetes is thought to be an autoimmune process, in which an environmental (unknown) trigger initiates a T cell-mediated immune response in genetically susceptible individuals. The presence of islet autoantibodies in the blood are signs of type 1 diabetes development, and risk of progressing to clinical type 1 diabetes is correlated with the presence of multiple islet autoantibodies. Currently, a "staging" model of type 1 diabetes proposes discrete components consisting of normal blood glucose but at least two islet autoantibodies (Stage 1), abnormal blood glucose with at least two islet autoantibodies (Stage 2), and clinical diagnosis (Stage 3). While these stages may, in fact, not be discrete and vary by individual, the format suggests important applications of precision medicine to diagnosis, prevention, prognosis, treatment and monitoring. In this paper, applications of precision medicine in type 1 diabetes are discussed, with both opportunities and barriers to global implementation highlighted. Several groups have implemented components of precision medicine, yet the integration of the necessary steps to achieve both short- and long-term solutions will need to involve researchers, patients, families, and healthcare providers to fully impact and reduce the burden of type 1 diabetes.