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

SUMMARY: This paper 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: GitHub: https://github.com/genepi-freiburg/metaGWASmanager. SUPPLEMENTARY INFORMATION: Supplementary materials are available at Bioinformatics online.

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.

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 Age-Related Decline of IRAK-M Expression in Retinal Pigment Epithelium Attenuates Outer Retinal Degeneration.

Unchecked, chronic inflammation is a constitutive component of age-related diseases, including age-related macular degeneration (AMD). Here we identified interleukin-1 receptor-associated kinase (IRAK)-M as a key immunoregulator in retinal pigment epithelium (RPE) that declines with age. Rare genetic variants of IRAK-M increased the likelihood of AMD. IRAK-M expression in RPE declined with age or oxidative stress and was further reduced in AMD. IRAK-M-deficient mice exhibited increased incidence of outer retinal degeneration at earlier ages, which was further exacerbated by oxidative stressors. The absence of IRAK-M disrupted RPE cell homeostasis, including compromised mitochondrial function, cellular senescence, and aberrant cytokine production. IRAK-M overexpression protected RPE cells against oxidative or immune stressors. Subretinal delivery of AAV-expressing IRAK-M rescued light-induced outer retinal degeneration in wild-type mice and attenuated age-related spontaneous retinal degeneration in IRAK-M-deficient mice. Our data support that replenishment of IRAK-M expression may redress dysregulated pro-inflammatory processes in AMD, thereby treating degeneration.

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 Risk Score Analysis Supports a Joint View of Two Classification Systems for Age-Related Macular Degeneration.

Purpose: The purpose of this study was to evaluate the utility of combining the Clinical Classification (CC) and the Three Continent age-related macular degeneration (AMD) Consortium Severity Scale (3CACSS) for classification of AMD. Methods: In two independent cross-sectional datasets of our population-based AugUR study (Altersbezogene Untersuchungen zur Gesundheit der Universität Regensburg), we graded AMD via color fundus images applying two established classification systems (CC and 3CACSS). We calculated the genetic risk score (GRS) across 50 previously identified variants for late AMD, its association via logistic regression, and area under the curve (AUC) for each AMD stage. Results: We analyzed 2188 persons aged 70 to 95 years. When comparing the two classification systems, we found a distinct pattern: CC "age-related changes" and CC "early AMD" distinguished individuals with 3CACSS "no AMD"; 3CACSS "mild/moderate/severe early AMD" stages, and distinguished CC "intermediate AMD". This suggested a 7-step scale combining the 2 systems: (i) "no AMD", (ii) "age-related changes", (iii) "very early AMD", (i.e. CC "early"), (iv) "mild early AMD", (v) "moderate early AMD", (vi) "severe early AMD", and (vii) "late AMD". GRS association and diagnostic accuracy increased stepwise by increased AMD severity in the 7-step scale and by increased restriction of controls (e.g. for CC "no AMD without age-related changes": AUC = 55.1%, 95% confidence interval [CI] = 51.6, 58.6, AUC = 62.3%, 95% CI = 59.1, 65.6, AUC = 63.8%, 95% CI = 59.3, 68.3, AUC = 78.1%, 95% CI = 73.6, 82.5, AUC = 82.2%, 95% CI = 78.4, 86.0, and AUC = 79.2%, 95% CI = 75.4, 83.0). A stepwise increase was also observed by increased drusen size and area. Conclusions: The utility of a 7-step scale is supported by our clinical and GRS data. This harmonization and full data integration provides an immediate simplification over using either CC or 3CACSS and helps to sharpen the control group.

Polygenic scores for estimated glomerular filtration rate in a population of general adults and elderly - comparative results from the KORA and AugUR study.

BACKGROUND: Polygenic scores (PGSs) combining genetic variants found to be associated with creatinine-based estimated glomerular filtration rate (eGFRcrea) have been applied in various study populations with different age ranges. This has shown that PGS explain less eGFRcrea variance in the elderly. Our aim was to understand how differences in eGFR variance and the percentage explained by PGS varies between population of general adults and elderly. RESULTS: We derived a PGS for cystatin-based eGFR (eGFRcys) from published genome-wide association studies. We used the 634 variants known for eGFRcrea and the 204 variants identified for eGFRcys to calculate the PGS in two comparable studies capturing a general adult and an elderly population, KORA S4 (n = 2,900; age 24-69 years) and AugUR (n = 2,272, age ≥ 70 years). To identify potential factors determining age-dependent differences on the PGS-explained variance, we evaluated the PGS variance, the eGFR variance, and the beta estimates of PGS association on eGFR. Specifically, we compared frequencies of eGFR-lowering alleles between general adult and elderly individuals and analyzed the influence of comorbidities and medication intake. The PGS for eGFRcrea explained almost twice as much (R2 = 9.6%) of age-/sex adjusted eGFR variance in the general adults compared to the elderly (4.6%). This difference was less pronounced for the PGS for eGFRcys (4.7% or 3.6%, respectively). The beta-estimate of the PGS on eGFRcrea was higher in the general adults compared to the elderly, but similar for the PGS on eGFRcys. The eGFR variance in the elderly was reduced by accounting for comorbidities and medication intake, but this did not explain the difference in R2-values. Allele frequencies between general adult and elderly individuals showed no significant differences except for one variant near APOE (rs429358). We found no enrichment of eGFR-protective alleles in the elderly compared to general adults. CONCLUSIONS: We concluded that the difference in explained variance by PGS was due to the higher age- and sex-adjusted eGFR variance in the elderly and, for eGFRcrea, also by a lower PGS association beta-estimate. Our results provide little evidence for survival or selection bias.

Differential and shared genetic effects on kidney function between diabetic and non-diabetic individuals.

Reduced glomerular filtration rate (GFR) can progress to kidney failure. Risk factors include genetics and diabetes mellitus (DM), but little is known about their interaction. We conducted genome-wide association meta-analyses for estimated GFR based on serum creatinine (eGFR), separately for individuals with or without DM (nDM = 178,691, nnoDM = 1,296,113). Our genome-wide searches identified (i) seven eGFR loci with significant DM/noDM-difference, (ii) four additional novel loci with suggestive difference and (iii) 28 further novel loci (including CUBN) by allowing for potential difference. GWAS on eGFR among DM individuals identified 2 known and 27 potentially responsible loci for diabetic kidney disease. Gene prioritization highlighted 18 genes that may inform reno-protective drug development. We highlight the existence of DM-only and noDM-only effects, which can inform about the target group, if respective genes are advanced as drug targets. Largely shared effects suggest that most drug interventions to alter eGFR should be effective in DM and noDM.

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.

Discovery and prioritization of variants and genes for kidney function in >1.2 million individuals.

Genes underneath signals from genome-wide association studies (GWAS) for kidney function are promising targets for functional studies, but prioritizing variants and genes is challenging. By GWAS meta-analysis for creatinine-based estimated glomerular filtration rate (eGFR) from the Chronic Kidney Disease Genetics Consortium and UK Biobank (n = 1,201,909), we expand the number of eGFRcrea loci (424 loci, 201 novel; 9.8% eGFRcrea variance explained by 634 independent signal variants). Our increased sample size in fine-mapping (n = 1,004,040, European) more than doubles the number of signals with resolved fine-mapping (99% credible sets down to 1 variant for 44 signals, ≤5 variants for 138 signals). Cystatin-based eGFR and/or blood urea nitrogen association support 348 loci (n = 460,826 and 852,678, respectively). Our customizable tool for Gene PrioritiSation reveals 23 compelling genes including mechanistic insights and enables navigation through genes and variants likely relevant for kidney function in human to help select targets for experimental follow-up.

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.

Target genes, variants, tissues and transcriptional pathways influencing human serum urate levels.

Elevated serum urate levels cause gout and correlate with cardiometabolic diseases via poorly understood mechanisms. We performed a trans-ancestry genome-wide association study of serum urate in 457,690 individuals, identifying 183 loci (147 previously unknown) that improve the prediction of gout in an independent cohort of 334,880 individuals. Serum urate showed significant genetic correlations with many cardiometabolic traits, with genetic causality analyses supporting a substantial role for pleiotropy. Enrichment analysis, fine-mapping of urate-associated loci and colocalization with gene expression in 47 tissues implicated the kidney and liver as the main target organs and prioritized potentially causal genes and variants, including the transcriptional master regulators in the liver and kidney, HNF1A and HNF4A. Experimental validation showed that HNF4A transactivated the promoter of ABCG2, encoding a major urate transporter, in kidney cells, and that HNF4A p.Thr139Ile is a functional variant. Transcriptional coregulation within and across organs may be a general mechanism underlying the observed pleiotropy between urate and cardiometabolic traits.

Genome-wide association meta-analyses and fine-mapping elucidate pathways influencing albuminuria.

Increased levels of the urinary albumin-to-creatinine ratio (UACR) are associated with higher risk of kidney disease progression and cardiovascular events, but underlying mechanisms are incompletely understood. Here, we conduct trans-ethnic (n = 564,257) and European-ancestry specific meta-analyses of genome-wide association studies of UACR, including ancestry- and diabetes-specific analyses, and identify 68 UACR-associated loci. Genetic correlation analyses and risk score associations in an independent electronic medical records database (n = 192,868) reveal connections with proteinuria, hyperlipidemia, gout, and hypertension. Fine-mapping and trans-Omics analyses with gene expression in 47 tissues and plasma protein levels implicate genes potentially operating through differential expression in kidney (including TGFB1, MUC1, PRKCI, and OAF), and allow coupling of UACR associations to altered plasma OAF concentrations. Knockdown of OAF and PRKCI orthologs in Drosophila nephrocytes reduces albumin endocytosis. Silencing fly PRKCI further impairs slit diaphragm formation. These results generate a priority list of genes and pathways for translational research to reduce albuminuria.

Mapping eGFR loci to the renal transcriptome and phenome in the VA Million Veteran Program.

Chronic kidney disease (CKD), defined by low estimated glomerular filtration rate (eGFR), contributes to global morbidity and mortality. Here we conduct a transethnic Genome-Wide Association Study of eGFR in 280,722 participants of the Million Veteran Program (MVP), with replication in 765,289 participants from the Chronic Kidney Disease Genetics (CKDGen) Consortium. We identify 82 previously unreported variants, confirm 54 loci, and report interesting findings including association of the sickle cell allele of betaglobin among non-Hispanic blacks. Our transcriptome-wide association study of kidney function in healthy kidney tissue identifies 36 previously unreported and nine known genes, and maps gene expression to renal cell types. In a Phenome-Wide Association Study in 192,868 MVP participants using a weighted genetic score we detect associations with CKD stages and complications and kidney stones. This investigation reinterprets the genetic architecture of kidney function to identify the gene, tissue, and anatomical context of renal homeostasis and the clinical consequences of dysregulation.

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

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

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.

Poor risk factor control in outpatients with diabetes mellitus type 2 in Germany: The DIAbetes COhoRtE (DIACORE) study.

INTRODUCTION: Patients with diabetes mellitus type 2 (DM2) are at high risk for micro- and macrovascular disease. Here, we explore the degree of traditional risk factor control in the baseline visit of a cohort of DM2 outpatients. METHODS: DIACORE (DIAbetes COhoRtE) is a prospective cohort study of 3000 adult DM2 outpatients. Here, we present results from the baseline visit. Sociodemographic and anthropometric variables, cardiovascular risk factors, comorbidities and medication were assessed by interview and medical exams. Serum-creatinine based estimated glomerular filtration rate (eGFRcrea) and urinary albumin-creatinine ratio (UACR) were determined for classification of chronic kidney disease (CKD). The proportion of patients with adequate control of traditional risk factors (blood pressure<140/90mmHg, HbA1c<7.5%, LDL<100mg/dl) was calculated in 2892 patients with non-missing data in 9 relevant variables within each KDIGO 2012 CKD class. RESULTS: In the analyzed baseline data (n = 2892, 60.2% men), mean (standard deviation) values for age, DM2 duration and HbA1c were 65.3 (9.3) years, 10.3 (8.4) years and 6.9% (1.1) respectively. Of these 2892 patients, 18.7% had CKD stage 3 or higher, 25.7% had UACR≥30mg/g. Adequate blood pressure, HbA1c and LDL control was achieved in 55.7%, 78.5% and 34.4%, respectively. In 16.4% of patients (473), all three risk factors were below recommended targets. The proportion of adequate risk factor control was similar across KDIGO eGFRcrea classes. Adequate blood pressure and HbA1c control were significantly associated with lower UACR category without and with controlling for other risk factors (p<0.0001, p = 0.0002, respectively). CONCLUSION: In our study of patients with diabetes mellitus type 2, we observed a low level of risk factor control indicating potential for risk reduction.

Daytime napping and diabetes-associated kidney disease.

BACKGROUND: Diabetes-associated Kidney Disease (DKD) is a common comorbidity in patients with type 2 diabetes. The present study investigates whether daytime sleeping duration in patients, ill with type 2 diabetes, is associated with DKD. METHODS: A total of 733 outpatients of the cross-sectional baseline survey of the DIACORE study were analyzed with respect to their self-reported daytime sleeping duration, assessed by a standardized questionnaire. DKD was defined as eGFR <60 ml/min/1.73 m2 and/or urinary albumin-to-creatinine-ratio (UACR) > 30 mg/g. RESULTS: Mean daytime sleeping duration was 17 ± 27 min. With increasing daytime sleeping duration a statistically significant decrease in eGFR (p = 0.002) and increase in UACR (p < 0.001) were found, respectively. Prevalence of DKD was significantly higher in patients with longer daytime sleeping duration (31% in patients not napping, 40% in patients napping less than 30 min, 47% in patients napping 30-60 min, 56% in patients napping 60 min or more; p = 0.001). After accounting for known modulators (Age, sex, BMI, waist-hip-ratio, systolic and diastolic blood pressure, physical activity, diabetes duration, HbA1c, homeostasis model assessment (HOMA-Index), nighttime sleeping duration, apnea-hypopnea-index (AHI), Epworth Sleepiness Scale (ESS)), longer daytime sleeping duration was significantly associated with impaired eGFR [B (95% CI) = -0.05 (-0.09; 0.00), p = 0.044] and increased UACR [B (95% CI) = 0.01 (0.01; 0.02), p < 0.001], respectively. CONCLUSION: Increased daytime sleeping duration is significantly associated with reduced eGFR and higher UACR, independent of known modulators of DKD. The direction of this relationship remains unclear.

Investigating the modulation of genetic effects on late AMD by age and sex: Lessons learned and two additional loci.

Late-stage age-related macular degeneration (AMD) is the leading cause of visual impairment in the elderly with a complex etiology. The most important non-modifiable risk factors for onset and progression of late AMD are age and genetic risk factors, however, little is known about the interplay between genetics and age or sex. Here, we conducted a large-scale age- and sex-stratified genome-wide association study (GWAS) using 1000 Genomes imputed genome-wide and ExomeChip data (>12 million variants). The data were established by the International Age-related Macular Degeneration Genomics Consortium (IAMDGC) from 16,144 late AMD cases and 17,832 controls. Our systematic search for interaction effects yielded significantly stronger effects among younger individuals at two known AMD loci (near CFH and ARMS2/HTRA1). Accounting for age and gene-age interaction using a joint test identified two additional AMD loci compared to the previous main effect scan. One of these two is a novel AMD GWAS locus, near the retinal clusterin-like protein (CLUL1) gene, and the other, near the retinaldehyde binding protein 1 (RLBP1), was recently identified in a joint analysis of nuclear and mitochondrial variants. Despite considerable power in our data, neither sex-dependent effects nor effects with opposite directions between younger and older individuals were observed. This is the first genome-wide interaction study to incorporate age, sex and their interaction with genetic effects for late AMD. Results diminish the potential for a role of sex in the etiology of late AMD yet highlight the importance and existence of age-dependent genetic effects.

HDL-cholesterol levels and risk of age-related macular degeneration: a multiethnic genetic study using Mendelian randomization.

Background: Dyslipidemia, particularly high-density lipoprotein cholesterol (HDL-C), has recently been implicated in the pathogenesis of age-related macular degeneration (AMD), the leading cause of vision loss. However, epidemiological studies have yielded conflicting results. Methods: We investigated the causal role of plasma lipid levels in AMD in multiethnic populations comprising 16 144 advanced AMD cases and 17 832 controls of European descent, together with 2219 cases and 5275 controls of Asian descent, using Mendelian randomization in three models. Model 1 is a conventional meta-analysis which does not account for pleiotropy of instrumental variable (IV) effects. Model 2 is a univariate, inverse variance weighted regression analysis that accounts for potential unbalanced pleiotropy using MR-Egger method. Finally, Model 3 is a multivariate regression analysis that addresses pleiotropy by MR-Egger method and by adjusting for effects on other lipid traits. Results: A 1 standard deviation (SD) higher HDL-cholesterol level was associated with an odds ratio (OR) for AMD of 1.17 (95% confidence interval: 1.07-1.29) in Europeans (P = 6.88 × 10-4) and of 1.58 (1.24-2.00) in Asians (P = 2.92 × 10-4) in Model 3. The corresponding OR estimates were 1.30 (1.09-1.55) in Europeans (P = 3.18 × 10-3) and 1.42 (1.11-1.80) in Asians (P = 4.42 × 10-3) in Model 1, and 1.21 (1.11-1.31) in Europeans (P = 3.12 × 10-5) and 1.51 (1.20-1.91) in Asians (P = 7.61 × 10-4) in Model 2. Conversely, neither LDL-C (Europeans: OR = 0.96, P = 0.272; Asians: OR = 1.02, P = 0.874; Model 3) nor triglyceride levels (Europeans: OR = 0.91, P = 0.102; Asians: OR = 1.06, P = 0.613) were associated with AMD. We also assessed the association between lipid levels and polypoidal choroidal vasculopathy (PCV) in Asians, a subtype of AMD, and found a similar trend for association of PCV with HDL-C levels. Conclusions: Our study shows that high levels of plasma HDL-C are causally associated with an increased risk for advanced AMD in European and Asian populations, implying that strategies reducing HDL-C levels may be useful to prevent and treat AMD.