metaGWASmanager: a toolbox for an automated workflow from phenotypes to meta-analysis in GWAS consortia.

SUMMARY: This article introduces the metaGWASmanager, which streamlines genome-wide association studies within large-scale meta-analysis consortia. It is a toolbox for both the central consortium analysis group and participating studies to generate homogeneous phenotypes, minimize unwanted variability from inconsistent methodologies, ensure high-quality association results, and implement time-efficient quality control workflows. The toolbox features a plug-in-based approach for customization of association testing. RESULTS: The metaGWASmanager toolbox has been successfully deployed in both the CKDGen and MetalGWAS Initiative consortia across hundreds of participating studies, demonstrating its effectiveness in GWAS analysis optimization by automating routine tasks and ensuring the value and reliability of association results, thus, ultimately promoting scientific discovery. We provide a simulated data set with examples for script customization so that readers can reproduce the pipeline at their convenience. AVAILABILITY AND IMPLEMENTATION: GitHub: https://github.com/genepi-freiburg/metaGWASmanager.

KidneyGPS: a user-friendly web application to help prioritize kidney function genes and variants based on evidence from genome-wide association studies.

BACKGROUND: Genome-wide association studies (GWAS) have identified hundreds of genetic loci associated with kidney function. By combining these findings with post-GWAS information (e.g., statistical fine-mapping to identify independent association signals and to narrow down signals to causal variants; or different sources of annotation data), new hypotheses regarding physiology and disease aetiology can be obtained. These hypotheses need to be tested in laboratory experiments, for example, to identify new therapeutic targets. For this purpose, the evidence obtained from GWAS and post-GWAS analyses must be processed and presented in a way that they are easily accessible to kidney researchers without specific GWAS expertise. MAIN: Here we present KidneyGPS, a user-friendly web-application that combines genetic variant association for estimated glomerular filtration rate (eGFR) from the Chronic Kidney Disease Genetics consortium with annotation of (i) genetic variants with functional or regulatory effects ("SNP-to-gene" mapping), (ii) genes with kidney phenotypes in mice or human ("gene-to-phenotype"), and (iii) drugability of genes (to support re-purposing). KidneyGPS adopts a comprehensive approach summarizing evidence for all 5906 genes in the 424 GWAS loci for eGFR identified previously and the 35,885 variants in the 99% credible sets of 594 independent signals. KidneyGPS enables user-friendly access to the abundance of information by search functions for genes, variants, and regions. KidneyGPS also provides a function ("GPS tab") to generate lists of genes with specific characteristics thus enabling customizable Gene Prioritisation (GPS). These specific characteristics can be as broad as any gene in the 424 loci with a known kidney phenotype in mice or human; or they can be highly focussed on genes mapping to genetic variants or signals with particularly with high statistical support. KidneyGPS is implemented with RShiny in a modularized fashion to facilitate update of input data ( https://kidneygps.ur.de/gps/ ). CONCLUSION: With the focus on kidney function related evidence, KidneyGPS fills a gap between large general platforms for accessing GWAS and post-GWAS results and the specific needs of the kidney research community. This makes KidneyGPS an important platform for kidney researchers to help translate in silico research results into in vitro or in vivo research.

Genetic loci and prioritization of genes for kidney function decline derived from a meta-analysis of 62 longitudinal genome-wide association studies.

Estimated glomerular filtration rate (eGFR) reflects kidney function. Progressive eGFR-decline can lead to kidney failure, necessitating dialysis or transplantation. Hundreds of loci from genome-wide association studies (GWAS) for eGFR help explain population cross section variability. Since the contribution of these or other loci to eGFR-decline remains largely unknown, we derived GWAS for annual eGFR-decline and meta-analyzed 62 longitudinal studies with eGFR assessed twice over time in all 343,339 individuals and in high-risk groups. We also explored different covariate adjustment. Twelve genome-wide significant independent variants for eGFR-decline unadjusted or adjusted for eGFR-baseline (11 novel, one known for this phenotype), including nine variants robustly associated across models were identified. All loci for eGFR-decline were known for cross-sectional eGFR and thus distinguished a subgroup of eGFR loci. Seven of the nine variants showed variant-by-age interaction on eGFR cross section (further about 350,000 individuals), which linked genetic associations for eGFR-decline with age-dependency of genetic cross-section associations. Clinically important were two to four-fold greater genetic effects on eGFR-decline in high-risk subgroups. Five variants associated also with chronic kidney disease progression mapped to genes with functional in-silico evidence (UMOD, SPATA7, GALNTL5, TPPP). An unfavorable versus favorable nine-variant genetic profile showed increased risk odds ratios of 1.35 for kidney failure (95% confidence intervals 1.03-1.77) and 1.27 for acute kidney injury (95% confidence intervals 1.08-1.50) in over 2000 cases each, with matched controls). Thus, we provide a large data resource, genetic loci, and prioritized genes for kidney function decline, which help inform drug development pipelines revealing important insights into the age-dependency of kidney function genetics.

Meta-analysis uncovers genome-wide significant variants for rapid kidney function decline.

Rapid decline of glomerular filtration rate estimated from creatinine (eGFRcrea) is associated with severe clinical endpoints. In contrast to cross-sectionally assessed eGFRcrea, the genetic basis for rapid eGFRcrea decline is largely unknown. To help define this, we meta-analyzed 42 genome-wide association studies from the Chronic Kidney Diseases Genetics Consortium and United Kingdom Biobank to identify genetic loci for rapid eGFRcrea decline. Two definitions of eGFRcrea decline were used: 3 mL/min/1.73m2/year or more ("Rapid3"; encompassing 34,874 cases, 107,090 controls) and eGFRcrea decline 25% or more and eGFRcrea under 60 mL/min/1.73m2 at follow-up among those with eGFRcrea 60 mL/min/1.73m2 or more at baseline ("CKDi25"; encompassing 19,901 cases, 175,244 controls). Seven independent variants were identified across six loci for Rapid3 and/or CKDi25: consisting of five variants at four loci with genome-wide significance (near UMOD-PDILT (2), PRKAG2, WDR72, OR2S2) and two variants among 265 known eGFRcrea variants (near GATM, LARP4B). All these loci were novel for Rapid3 and/or CKDi25 and our bioinformatic follow-up prioritized variants and genes underneath these loci. The OR2S2 locus is novel for any eGFRcrea trait including interesting candidates. For the five genome-wide significant lead variants, we found supporting effects for annual change in blood urea nitrogen or cystatin-based eGFR, but not for GATM or LARP4B. Individuals at high compared to those at low genetic risk (8-14 vs. 0-5 adverse alleles) had a 1.20-fold increased risk of acute kidney injury (95% confidence interval 1.08-1.33). Thus, our identified loci for rapid kidney function decline may help prioritize therapeutic targets and identify mechanisms and individuals at risk for sustained deterioration of kidney function.

On the differences between mega- and meta-imputation and analysis exemplified on the genetics of age-related macular degeneration.

While current genome-wide association analyses often rely on meta-analysis of study-specific summary statistics, individual participant data (IPD) from multiple studies increase options for modeling. When multistudy IPD is available, however, it is unclear whether this data is to be imputed and modeled across all participants (mega-imputation and mega-analysis) or study-specifically (meta-imputation and meta-analysis). Here, we investigated different approaches toward imputation and analysis using 52,189 subjects from 25 studies of the International Age-related Macular Degeneration (AMD) Genomics Consortium including, 16,144 AMD cases and 17,832 controls for association analysis. From 27,448,454 genetic variants after 1,000-Genomes-based imputation, mega-imputation yielded ~400,000 more variants with high imputation quality (mostly rare variants) compared to meta-imputation. For AMD signal detection (P < 5 × 10-8 ) in mega-imputed data, most loci were detected with mega-analysis without adjusting for study membership (40 loci, including 34 known); we considered these loci genuine, since genetic effects and P-values were comparable across analyses. In meta-imputed data, we found 31 additional signals, mostly near chromosome tails or reference panel gaps, which disappeared after accounting for interaction of whole-genome amplification (WGA) with study membership or after excluding studies with WGA-participants. For signal detection with multistudy IPD, we recommend mega-imputation and mega-analysis, with meta-imputation followed by meta-analysis being a computationally appealing alternative.

NFAT5 and SLC4A10 Loci Associate with Plasma Osmolality.

Disorders of water balance, an excess or deficit of total body water relative to body electrolyte content, are common and ascertained by plasma hypo- or hypernatremia, respectively. We performed a two-stage genome-wide association study meta-analysis on plasma sodium concentration in 45,889 individuals of European descent (stage 1 discovery) and 17,637 additional individuals of European descent (stage 2 replication), and a transethnic meta-analysis of replicated single-nucleotide polymorphisms in 79,506 individuals (63,526 individuals of European descent, 8765 individuals of Asian Indian descent, and 7215 individuals of African descent). In stage 1, we identified eight loci associated with plasma sodium concentration at P<5.0 × 10-6 Of these, rs9980 at NFAT5 replicated in stage 2 meta-analysis (P=3.1 × 10-5), with combined stages 1 and 2 genome-wide significance of P=5.6 × 10-10 Transethnic meta-analysis further supported the association at rs9980 (P=5.9 × 10-12). Additionally, rs16846053 at SLC4A10 showed nominally, but not genome-wide, significant association in combined stages 1 and 2 meta-analysis (P=6.7 × 10-8). NFAT5 encodes a ubiquitously expressed transcription factor that coordinates the intracellular response to hypertonic stress but was not previously implicated in the regulation of systemic water balance. SLC4A10 encodes a sodium bicarbonate transporter with a brain-restricted expression pattern, and variant rs16846053 affects a putative intronic NFAT5 DNA binding motif. The lead variants for NFAT5 and SLC4A10 are cis expression quantitative trait loci in tissues of the central nervous system and relevant to transcriptional regulation. Thus, genetic variation in NFAT5 and SLC4A10 expression and function in the central nervous system may affect the regulation of systemic water balance.

SOS2 and ACP1 Loci Identified through Large-Scale Exome Chip Analysis Regulate Kidney Development and Function.

Genome-wide association studies have identified >50 common variants associated with kidney function, but these variants do not fully explain the variation in eGFR. We performed a two-stage meta-analysis of associations between genotypes from the Illumina exome array and eGFR on the basis of serum creatinine (eGFRcrea) among participants of European ancestry from the CKDGen Consortium (nStage1: 111,666; nStage2: 48,343). In single-variant analyses, we identified single nucleotide polymorphisms at seven new loci associated with eGFRcrea (PPM1J, EDEM3, ACP1, SPEG, EYA4, CYP1A1, and ATXN2L; PStage1<3.7×10-7), of which most were common and annotated as nonsynonymous variants. Gene-based analysis identified associations of functional rare variants in three genes with eGFRcrea, including a novel association with the SOS Ras/Rho guanine nucleotide exchange factor 2 gene, SOS2 (P=5.4×10-8 by sequence kernel association test). Experimental follow-up in zebrafish embryos revealed changes in glomerular gene expression and renal tubule morphology in the embryonic kidney of acp1- and sos2-knockdowns. These developmental abnormalities associated with altered blood clearance rate and heightened prevalence of edema. This study expands the number of loci associated with kidney function and identifies novel genes with potential roles in kidney formation.

Gene-centric meta-analyses for central adiposity traits in up to 57 412 individuals of European descent confirm known loci and reveal several novel associations.

Waist circumference (WC) and waist-to-hip ratio (WHR) are surrogate measures of central adiposity that are associated with adverse cardiovascular events, type 2 diabetes and cancer independent of body mass index (BMI). WC and WHR are highly heritable with multiple susceptibility loci identified to date. We assessed the association between SNPs and BMI-adjusted WC and WHR and unadjusted WC in up to 57 412 individuals of European descent from 22 cohorts collaborating with the NHLBI's Candidate Gene Association Resource (CARe) project. The study population consisted of women and men aged 20-80 years. Study participants were genotyped using the ITMAT/Broad/CARE array, which includes ∼50 000 cosmopolitan tagged SNPs across ∼2100 cardiovascular-related genes. Each trait was modeled as a function of age, study site and principal components to control for population stratification, and we conducted a fixed-effects meta-analysis. No new loci for WC were observed. For WHR analyses, three novel loci were significantly associated (P < 2.4 × 10(-6)). Previously unreported rs2811337-G near TMCC1 was associated with increased WHR (ß ± SE, 0.048 ± 0.008, P = 7.7 × 10(-9)) as was rs7302703-G in HOXC10 (ß = 0.044 ± 0.008, P = 2.9 × 10(-7)) and rs936108-C in PEMT (ß = 0.035 ± 0.007, P = 1.9 × 10(-6)). Sex-stratified analyses revealed two additional novel signals among females only, rs12076073-A in SHC1 (ß = 0.10 ± 0.02, P = 1.9 × 10(-6)) and rs1037575-A in ATBDB4 (ß = 0.046 ± 0.01, P = 2.2 × 10(-6)), supporting an already established sexual dimorphism of central adiposity-related genetic variants. Functional analysis using ENCODE and eQTL databases revealed that several of these loci are in regulatory regions or regions with differential expression in adipose tissue.

EasyStrata: evaluation and visualization of stratified genome-wide association meta-analysis data.

UNLABELLED: The R package EasyStrata facilitates the evaluation and visualization of stratified genome-wide association meta-analyses (GWAMAs) results. It provides (i) statistical methods to test and account for between-strata difference as a means to tackle gene-strata interaction effects and (ii) extended graphical features tailored for stratified GWAMA results. The software provides further features also suitable for general GWAMAs including functions to annotate, exclude or highlight specific loci in plots or to extract independent subsets of loci from genome-wide datasets. It is freely available and includes a user-friendly scripting interface that simplifies data handling and allows for combining statistical and graphical functions in a flexible fashion. AVAILABILITY: EasyStrata is available for free (under the GNU General Public License v3) from our Web site www.genepi-regensburg.de/easystrata and from the CRAN R package repository cran.r-project.org/web/packages/EasyStrata/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Harmonization of study and reference data by PhaseLift: saving time when imputing study data.

Genome-wide association studies are usually accompanied by imputation techniques to complement genome-wide SNP chip genotypes. Current imputation approaches separate the phasing of study data from imputing, which makes the phasing independent from the reference data. The two-step approach allows for updating the imputation for a new reference panel without repeating the tedious phasing step. This advantage, however, does no longer hold, when the build of the study data differs from the build of the reference data. In this case, the current approach is to harmonize the study data annotation with the reference data (prephasing lift-over), requiring rephasing and re-imputing. As a novel approach, we propose to harmonize study haplotypes with reference haplotypes (postphasing lift-over). This allows for updating imputed study data for new reference panels without requiring rephasing. With continuously updated reference panels, our approach can save considerable computing time of up to 1 month per re-imputation. We evaluated the rephasing and postphasing lift-over approaches by using data from 1,644 unrelated individuals imputed by both approaches and comparing it with directly typed genotypes. On average, both approaches perform equally well with mean concordances of 93% between imputed and typed genotypes for both approaches. Also, imputation qualities are similar (mean difference in RSQ < 0.1%). We demonstrate that our novel postphasing lift-over approach is a practical and time-saving alternative to the prephasing lift-over. This might encourage study partners to accommodate updated reference builds and ultimately improve the information content of study data. Our novel approach is implemented in the software PhaseLift.

Genome-wide association study of kidney function decline in individuals of European descent.

Genome-wide association studies (GWASs) have identified multiple loci associated with cross-sectional eGFR, but a systematic genetic analysis of kidney function decline over time is missing. Here we conducted a GWAS meta-analysis among 63,558 participants of European descent, initially from 16 cohorts with serial kidney function measurements within the CKDGen Consortium, followed by independent replication among additional participants from 13 cohorts. In stage 1 GWAS meta-analysis, single-nucleotide polymorphisms (SNPs) at MEOX2, GALNT11, IL1RAP, NPPA, HPCAL1, and CDH23 showed the strongest associations for at least one trait, in addition to the known UMOD locus, which showed genome-wide significance with an annual change in eGFR. In stage 2 meta-analysis, the significant association at UMOD was replicated. Associations at GALNT11 with Rapid Decline (annual eGFR decline of 3 ml/min per 1.73 m(2) or more), and CDH23 with eGFR change among those with CKD showed significant suggestive evidence of replication. Combined stage 1 and 2 meta-analyses showed significance for UMOD, GALNT11, and CDH23. Morpholino knockdowns of galnt11 and cdh23 in zebrafish embryos each had signs of severe edema 72 h after gentamicin treatment compared with controls, but no gross morphological renal abnormalities before gentamicin administration. Thus, our results suggest a role in the deterioration of kidney function for the loci GALNT11 and CDH23, and show that the UMOD locus is significantly associated with kidney function decline.

Genome-wide association and functional follow-up reveals new loci for kidney function.

Chronic kidney disease (CKD) is an important public health problem with a genetic component. We performed genome-wide association studies in up to 130,600 European ancestry participants overall, and stratified for key CKD risk factors. We uncovered 6 new loci in association with estimated glomerular filtration rate (eGFR), the primary clinical measure of CKD, in or near MPPED2, DDX1, SLC47A1, CDK12, CASP9, and INO80. Morpholino knockdown of mpped2 and casp9 in zebrafish embryos revealed podocyte and tubular abnormalities with altered dextran clearance, suggesting a role for these genes in renal function. By providing new insights into genes that regulate renal function, these results could further our understanding of the pathogenesis of CKD.

Integration of genome-wide association studies with biological knowledge identifies six novel genes related to kidney function.

In conducting genome-wide association studies (GWAS), analytical approaches leveraging biological information may further understanding of the pathophysiology of clinical traits. To discover novel associations with estimated glomerular filtration rate (eGFR), a measure of kidney function, we developed a strategy for integrating prior biological knowledge into the existing GWAS data for eGFR from the CKDGen Consortium. Our strategy focuses on single nucleotide polymorphism (SNPs) in genes that are connected by functional evidence, determined by literature mining and gene ontology (GO) hierarchies, to genes near previously validated eGFR associations. It then requires association thresholds consistent with multiple testing, and finally evaluates novel candidates by independent replication. Among the samples of European ancestry, we identified a genome-wide significant SNP in FBXL20 (P = 5.6 × 10(-9)) in meta-analysis of all available data, and additional SNPs at the INHBC, LRP2, PLEKHA1, SLC3A2 and SLC7A6 genes meeting multiple-testing corrected significance for replication and overall P-values of 4.5 × 10(-4)-2.2 × 10(-7). Neither the novel PLEKHA1 nor FBXL20 associations, both further supported by association with eGFR among African Americans and with transcript abundance, would have been implicated by eGFR candidate gene approaches. LRP2, encoding the megalin receptor, was identified through connection with the previously known eGFR gene DAB2 and extends understanding of the megalin system in kidney function. These findings highlight integration of existing genome-wide association data with independent biological knowledge to uncover novel candidate eGFR associations, including candidates lacking known connections to kidney-specific pathways. The strategy may also be applicable to other clinical phenotypes, although more testing will be needed to assess its potential for discovery in general.

Association of eGFR-Related Loci Identified by GWAS with Incident CKD and ESRD.

Family studies suggest a genetic component to the etiology of chronic kidney disease (CKD) and end stage renal disease (ESRD). Previously, we identified 16 loci for eGFR in genome-wide association studies, but the associations of these single nucleotide polymorphisms (SNPs) for incident CKD or ESRD are unknown. We thus investigated the association of these loci with incident CKD in 26,308 individuals of European ancestry free of CKD at baseline drawn from eight population-based cohorts followed for a median of 7.2 years (including 2,122 incident CKD cases defined as eGFR <60ml/min/1.73m(2) at follow-up) and with ESRD in four case-control studies in subjects of European ancestry (3,775 cases, 4,577 controls). SNPs at 11 of the 16 loci (UMOD, PRKAG2, ANXA9, DAB2, SHROOM3, DACH1, STC1, SLC34A1, ALMS1/NAT8, UBE2Q2, and GCKR) were associated with incident CKD; p-values ranged from p = 4.1e-9 in UMOD to p = 0.03 in GCKR. After adjusting for baseline eGFR, six of these loci remained significantly associated with incident CKD (UMOD, PRKAG2, ANXA9, DAB2, DACH1, and STC1). SNPs in UMOD (OR = 0.92, p = 0.04) and GCKR (OR = 0.93, p = 0.03) were nominally associated with ESRD. In summary, the majority of eGFR-related loci are either associated or show a strong trend towards association with incident CKD, but have modest associations with ESRD in individuals of European descent. Additional work is required to characterize the association of genetic determinants of CKD and ESRD at different stages of disease progression.

LMX1B is essential for the maintenance of differentiated podocytes in adult kidneys.

Mutations of the LMX1B gene cause nail-patella syndrome, a rare autosomal-dominant disorder affecting the development of the limbs, eyes, brain, and kidneys. The characterization of conventional Lmx1b knockout mice has shown that LMX1B regulates the development of podocyte foot processes and slit diaphragms, but studies using podocyte-specific Lmx1b knockout mice have yielded conflicting results regarding the importance of LMX1B for maintaining podocyte structures. In order to address this question, we generated inducible podocyte-specific Lmx1b knockout mice. One week of Lmx1b inactivation in adult mice resulted in proteinuria with only minimal foot process effacement. Notably, expression levels of slit diaphragm and basement membrane proteins remained stable at this time point, and basement membrane charge properties also did not change, suggesting that alternative mechanisms mediate the development of proteinuria in these mice. Cell biological and biophysical experiments with primary podocytes isolated after 1 week of Lmx1b inactivation indicated dysregulation of actin cytoskeleton organization, and time-resolved DNA microarray analysis identified the genes encoding actin cytoskeleton-associated proteins, including Abra and Arl4c, as putative LMX1B targets. Chromatin immunoprecipitation experiments in conditionally immortalized human podocytes and gel shift assays showed that LMX1B recognizes AT-rich binding sites (FLAT elements) in the promoter regions of ABRA and ARL4C, and knockdown experiments in zebrafish support a model in which LMX1B and ABRA act in a common pathway during pronephros development. Our report establishes the importance of LMX1B in fully differentiated podocytes and argues that LMX1B is essential for the maintenance of an appropriately structured actin cytoskeleton in podocytes.

Common variants in Mendelian kidney disease genes and their association with renal function.

Many common genetic variants identified by genome-wide association studies for complex traits map to genes previously linked to rare inherited Mendelian disorders. A systematic analysis of common single-nucleotide polymorphisms (SNPs) in genes responsible for Mendelian diseases with kidney phenotypes has not been performed. We thus developed a comprehensive database of genes for Mendelian kidney conditions and evaluated the association between common genetic variants within these genes and kidney function in the general population. Using the Online Mendelian Inheritance in Man database, we identified 731 unique disease entries related to specific renal search terms and confirmed a kidney phenotype in 218 of these entries, corresponding to mutations in 258 genes. We interrogated common SNPs (minor allele frequency >5%) within these genes for association with the estimated GFR in 74,354 European-ancestry participants from the CKDGen Consortium. However, the top four candidate SNPs (rs6433115 at LRP2, rs1050700 at TSC1, rs249942 at PALB2, and rs9827843 at ROBO2) did not achieve significance in a stage 2 meta-analysis performed in 56,246 additional independent individuals, indicating that these common SNPs are not associated with estimated GFR. The effect of less common or rare variants in these genes on kidney function in the general population and disease-specific cohorts requires further research.

Nuclear and mitochondrial genetic variants associated with mitochondrial DNA copy number.

Mitochondrial DNA copy number (mtDNA-CN) is a biomarker for mitochondrial dysfunction associated with several diseases. Previous genome-wide association studies (GWAS) have been performed to unravel underlying mechanisms of mtDNA-CN regulation. However, the identified gene regions explain only a small fraction of mtDNA-CN variability. Most of this data has been estimated from microarrays based on various pipelines. In the present study we aimed to (1) identify genetic loci for qPCR-measured mtDNA-CN from three studies (16,130 participants) using GWAS, (2) identify potential systematic differences between our qPCR derived mtDNA-CN measurements compared to the published microarray intensity-based estimates, and (3) disentangle the nuclear from mitochondrial regulation of the mtDNA-CN phenotype. We identified two genome-wide significant autosomal loci associated with qPCR-measured mtDNA-CN: at HBS1L (rs4895440, p = 3.39 × 10-13) and GSDMA (rs56030650, p = 4.85 × 10-08) genes. Moreover, 113/115 of the previously published SNPs identified by microarray-based analyses were significantly equivalent with our findings. In our study, the mitochondrial genome itself contributed only marginally to mtDNA-CN regulation as we only detected a single rare mitochondrial variant associated with mtDNA-CN. Furthermore, we incorporated mitochondrial haplogroups into our analyses to explore their potential impact on mtDNA-CN. However, our findings indicate that they do not exert any significant influence on our results.

Replenishing IRAK-M expression in retinal pigment epithelium attenuates outer retinal degeneration.

Chronic inflammation is a constitutive component of many age-related diseases, including age-related macular degeneration (AMD). Here, we identified interleukin-1 receptor-associated kinase M (IRAK-M) as a key immunoregulator in retinal pigment epithelium (RPE) that declines during the aging process. Rare genetic variants of IRAK3, which encodes IRAK-M, were associated with an increased likelihood of developing AMD. In human samples and mouse models, IRAK-M abundance in the RPE declined with advancing age or exposure to oxidative stress and was further reduced in AMD. Irak3-knockout mice exhibited an increased incidence of outer retinal degeneration at earlier ages, which was further exacerbated by oxidative stressors. The absence of IRAK-M led to a disruption in RPE cell homeostasis, characterized by compromised mitochondrial function, cellular senescence, and aberrant cytokine production. IRAK-M overexpression protected RPE cells against oxidative or immune stressors. Subretinal delivery of adeno-associated virus (AAV)-expressing human IRAK3 rescued light-induced outer retinal degeneration in wild-type mice and attenuated age-related spontaneous retinal degeneration in Irak3-knockout mice. Our data show that replenishment of IRAK-M in the RPE may redress dysregulated pro-inflammatory processes in AMD, suggesting a potential treatment for retinal degeneration.

X-chromosome and kidney function: evidence from a multi-trait genetic analysis of 908,697 individuals reveals sex-specific and sex-differential findings in genes regulated by androgen response elements.

X-chromosomal genetic variants are understudied but can yield valuable insights into sexually dimorphic human traits and diseases. We performed a sex-stratified cross-ancestry X-chromosome-wide association meta-analysis of seven kidney-related traits (n = 908,697), identifying 23 loci genome-wide significantly associated with two of the traits: 7 for uric acid and 16 for estimated glomerular filtration rate (eGFR), including four novel eGFR loci containing the functionally plausible prioritized genes ACSL4, CLDN2, TSPAN6 and the female-specific DRP2. Further, we identified five novel sex-interactions, comprising male-specific effects at FAM9B and AR/EDA2R, and three sex-differential findings with larger genetic effect sizes in males at DCAF12L1 and MST4 and larger effect sizes in females at HPRT1. All prioritized genes in loci showing significant sex-interactions were located next to androgen response elements (ARE). Five ARE genes showed sex-differential expressions. This study contributes new insights into sex-dimorphisms of kidney traits along with new prioritized gene targets for further molecular research.

Study design of DIACORE (DIAbetes COhoRtE) - a cohort study of patients with diabetes mellitus type 2.

BACKGROUND: Diabetes mellitus type 2 (DM2) is highly associated with increased risk for chronic kidney disease (CKD), end stage renal disease (ESRD) and cardiovascular morbidity. Epidemiological and genetic studies generate hypotheses for innovative strategies in DM2 management by unravelling novel mechanisms of diabetes complications, which is essential for future intervention trials. We have thus initiated the DIAbetes COhoRtE study (DIACORE). METHODS: DIACORE is a prospective cohort study aiming to recruit 6000 patients of self-reported Caucasian ethnicity with prevalent DM2 for at least 10 years of follow-up. Study visits are performed in University-based recruiting clinics in Germany using standard operating procedures. All prevalent DM2 patients in outpatient clinics surrounding the recruiting centers are invited to participate. At baseline and at each 2-year follow-up examination, patients are subjected to a core phenotyping protocol. This includes a standardized online questionnaire and physical examination to determine incident micro- and macrovascular DM2 complications, malignancy and hospitalization, with a primary focus on renal events. Confirmatory outcome information is requested from patient records. Blood samples are obtained for a centrally analyzed standard laboratory panel and for biobanking of aliquots of serum, plasma, urine, mRNA and DNA for future scientific use. A subset of the cohort is subjected to extended phenotyping, e.g. sleep apnea screening, skin autofluorescence measurement, non-mydriatic retinal photography and non-invasive determination of arterial stiffness. DISCUSSION: DIACORE will enable the prospective evaluation of factors involved in DM2 complication pathogenesis using high-throughput technologies in biosamples and genetic epidemiological studies.