Genome-wide association study of placental weight identifies distinct and shared genetic influences between placental and fetal growth.

A well-functioning placenta is essential for fetal and maternal health throughout pregnancy. Using placental weight as a proxy for placental growth, we report genome-wide association analyses in the fetal (n = 65,405), maternal (n = 61,228) and paternal (n = 52,392) genomes, yielding 40 independent association signals. Twenty-six signals are classified as fetal, four maternal and three fetal and maternal. A maternal parent-of-origin effect is seen near KCNQ1. Genetic correlation and colocalization analyses reveal overlap with birth weight genetics, but 12 loci are classified as predominantly or only affecting placental weight, with connections to placental development and morphology, and transport of antibodies and amino acids. Mendelian randomization analyses indicate that fetal genetically mediated higher placental weight is causally associated with preeclampsia risk and shorter gestational duration. Moreover, these analyses support the role of fetal insulin in regulating placental weight, providing a key link between fetal and placental growth.

European and multi-ancestry genome-wide association meta-analysis of atopic dermatitis highlights importance of systemic immune regulation.

Atopic dermatitis (AD) is a common inflammatory skin condition and prior genome-wide association studies (GWAS) have identified 71 associated loci. In the current study we conducted the largest AD GWAS to date (discovery N = 1,086,394, replication N = 3,604,027), combining previously reported cohorts with additional available data. We identified 81 loci (29 novel) in the European-only analysis (which all replicated in a separate European analysis) and 10 additional loci in the multi-ancestry analysis (3 novel). Eight variants from the multi-ancestry analysis replicated in at least one of the populations tested (European, Latino or African), while two may be specific to individuals of Japanese ancestry. AD loci showed enrichment for DNAse I hypersensitivity and eQTL associations in blood. At each locus we prioritised candidate genes by integrating multi-omic data. The implicated genes are predominantly in immune pathways of relevance to atopic inflammation and some offer drug repurposing opportunities.

Comprehensive genome-wide association study of different forms of hernia identifies more than 80 associated loci.

Hernias are characterized by protrusion of an organ or tissue through its surrounding cavity and often require surgical repair. In this study we identify 65,492 cases for five hernia types in the UK Biobank and perform genome-wide association study scans for these five types and two combined groups. Our results show associated variants in all scans. Inguinal hernia has the most associations and we conduct a follow-up study with 23,803 additional cases from four study groups giving 84 independently associated variants. Identified variants from all scans are collapsed into 81 independent loci. Further testing shows that 26 loci are associated with more than one hernia type, suggesting substantial overlap between the underlying genetic mechanisms. Pathway analyses identify several genes with a strong link to collagen and/or elastin (ADAMTS6, ADAMTS16, ADAMTSL3, LOX, ELN) in the vicinity of associated loci for inguinal hernia, which substantiates an essential role of connective tissue morphology.

Genome-wide association study of febrile seizures implicates fever response and neuronal excitability genes.

Febrile seizures represent the most common type of pathological brain activity in young children and are influenced by genetic, environmental and developmental factors. In a minority of cases, febrile seizures precede later development of epilepsy. We conducted a genome-wide association study of febrile seizures in 7635 cases and 83 966 controls identifying and replicating seven new loci, all with P < 5 × 10-10. Variants at two loci were functionally related to altered expression of the fever response genes PTGER3 and IL10, and four other loci harboured genes (BSN, ERC2, GABRG2, HERC1) influencing neuronal excitability by regulating neurotransmitter release and binding, vesicular transport or membrane trafficking at the synapse. Four previously reported loci (SCN1A, SCN2A, ANO3 and 12q21.33) were all confirmed. Collectively, the seven novel and four previously reported loci explained 2.8% of the variance in liability to febrile seizures, and the single nucleotide polymorphism heritability based on all common autosomal single nucleotide polymorphisms was 10.8%. GABRG2, SCN1A and SCN2A are well-established epilepsy genes and, overall, we found positive genetic correlations with epilepsies (rg = 0.39, P = 1.68 × 10-4). Further, we found that higher polygenic risk scores for febrile seizures were associated with epilepsy and with history of hospital admission for febrile seizures. Finally, we found that polygenic risk of febrile seizures was lower in febrile seizure patients with neuropsychiatric disease compared to febrile seizure patients in a general population sample. In conclusion, this largest genetic investigation of febrile seizures to date implicates central fever response genes as well as genes affecting neuronal excitability, including several known epilepsy genes. Further functional and genetic studies based on these findings will provide important insights into the complex pathophysiological processes of seizures with and without fever.

Genetic regulation of spermine oxidase activity and cancer risk: a Mendelian randomization study.

Spermine oxidase (SMOX) catalyzes the oxidation of spermine to spermidine. Observational studies have reported SMOX as a source of reactive oxygen species associated with cancer, implying that inhibition of SMOX could be a target for chemoprevention. Here we test causality of SMOX levels with cancer risk using a Mendelian randomization analysis. We performed a GWAS of spermidine/spermine ratio to identify genetic variants associated with regulation of SMOX activity. Replication analysis was performed in two datasets of SMOX gene expression. We then did a Mendelian randomization analysis by testing the association between the SMOX genetic instrument and neuroblastoma, gastric, lung, breast, prostate, and colorectal cancers using GWAS summary statistics. GWAS of spermidine/spermine ratio identified SMOX locus (P = 1.34 × 10-49) explaining 32% of the variance. The lead SNP rs1741315 was also associated with SMOX gene expression in newborns (P = 8.48 × 10-28) and adults (P = 2.748 × 10-8) explaining 37% and 6% of the variance, respectively. Genetically determined SMOX activity was not associated with neuroblastoma, gastric, lung, breast, prostate nor colorectal cancer (P > 0.05). A PheWAS of rs1741315 did not reveal any relevant associations. Common genetic variation in the SMOX gene was strongly associated with SMOX activity in newborns, and less strongly in adults. Genetic down-regulation of SMOX was not significantly associated with lower odds of neuroblastoma, gastric, lung, breast, prostate and colorectal cancer. These results may inform studies of SMOX inhibition as a target for chemoprevention.

Shared genetic etiology between idiopathic pulmonary fibrosis and COVID-19 severity.

BACKGROUND: Idiopathic pulmonary fibrosis (IPF) is a complex lung disease, characterized by progressive lung scarring. Severe COVID-19 is associated with substantial pneumonitis and has a number of shared major risk factors with IPF. This study aimed to determine the genetic correlation between IPF and severe COVID-19 and assess a potential causal role of genetically increased risk of IPF on COVID-19 severity. METHODS: The genetic correlation between IPF and COVID-19 severity was estimated with linkage disequilibrium (LD) score regression. We performed a Mendelian randomization (MR) study for IPF causality in COVID-19. Genetic variants associated with IPF susceptibility (P<5 × 10-8) in previous genome-wide association studies (GWAS) were used as instrumental variables (IVs). Effect estimates of those IVs on COVID-19 severity were gathered from the GWAS meta-analysis by the COVID-19 Host Genetics Initiative (4,336 cases & 623,902 controls). FINDINGS: We detected a positive genetic correlation of IPF with COVID-19 severity (rg=0·31 [95% CI 0·04-0·57], P = 0·023). The MR estimates for severe COVID-19 did not reveal any genetic association (OR 1·05, [95% CI 0·92-1·20], P = 0·43). However, outlier analysis revealed that the IPF risk allele rs35705950 at MUC5B had a different effect compared with the other variants. When rs35705950 was excluded, MR results provided evidence that genetically increased risk of IPF has a causal effect on COVID-19 severity (OR 1·21, [95% CI 1·06-1·38], P = 4·24 × 10-3). Furthermore, the IPF risk-allele at MUC5B showed an apparent protective effect against COVID-19 hospitalization only in older adults (OR 0·86, [95% CI 0·73-1·00], P = 2·99 × 10-2) . INTERPRETATION: The strongest genetic determinant of IPF, rs35705950 at MUC5B, seems to confer protection against COVID-19, whereas the combined effect of all other IPF risk loci seem to confer risk of COVID-19 severity. The observed effect of rs35705950 could either be due to protective effects of mucin over-production on the airways or a consequence of selection bias due to (1) a patient group that is heavily enriched for the rs35705950 T undertaking strict self-isolation and/or (2) due to survival bias of the rs35705950 non-IPF risk allele carriers. Due to the diverse impact of IPF causal variants on SARS-CoV-2 infection, with a possible selection bias as an explanation, further investigation is needed to address this apparent paradox between variance at MUC5B and other IPF genetic risk factors. FUNDING: Novo Nordisk Foundation and Oak Foundation.

Integrating genetics with newborn metabolomics in infantile hypertrophic pyloric stenosis.

INTRODUCTION: Infantile hypertrophic pyloric stenosis (IHPS) is caused by hypertrophy of the pyloric sphincter muscle. OBJECTIVES: Since previous reports have implicated lipid metabolism, we aimed to (1) investigate associations between IHPS and a wide array of lipid-related metabolites in newborns, and (2) address whether detected differences in metabolite levels were likely to be driven by genetic differences between IHPS cases and controls or by differences in early life feeding patterns. METHODS: We used population-based random selection of IHPS cases and controls born in Denmark between 1997 and 2014. We randomly took dried blood spots of newborns from 267 pairs of IHPS cases and controls matched by sex and day of birth. We used a mixed-effects linear regression model to evaluate associations between 148 metabolites and IHPS in a matched case-control design. RESULTS: The phosphatidylcholine PC(38:4) showed significantly lower levels in IHPS cases (P = 4.68 × 10-8) as did six other correlated metabolites (four phosphatidylcholines, acylcarnitine AC(2:0), and histidine). Associations were driven by 98 case-control pairs born before 2009, when median age at sampling was 6 days. No association was seen in 169 pairs born in 2009 or later, when median age at sampling was 2 days. More IHPS cases than controls had a diagnosis for neonatal difficulty in feeding at breast (P = 6.15 × 10-3). Genetic variants known to be associated with PC(38:4) levels did not associate with IHPS. CONCLUSIONS: We detected lower levels of certain metabolites in IHPS, possibly reflecting different feeding patterns in the first days of life.

Sex-dimorphic genetic effects and novel loci for fasting glucose and insulin variability.

Differences between sexes contribute to variation in the levels of fasting glucose and insulin. Epidemiological studies established a higher prevalence of impaired fasting glucose in men and impaired glucose tolerance in women, however, the genetic component underlying this phenomenon is not established. We assess sex-dimorphic (73,089/50,404 women and 67,506/47,806 men) and sex-combined (151,188/105,056 individuals) fasting glucose/fasting insulin genetic effects via genome-wide association study meta-analyses in individuals of European descent without diabetes. Here we report sex dimorphism in allelic effects on fasting insulin at IRS1 and ZNF12 loci, the latter showing higher RNA expression in whole blood in women compared to men. We also observe sex-homogeneous effects on fasting glucose at seven novel loci. Fasting insulin in women shows stronger genetic correlations than in men with waist-to-hip ratio and anorexia nervosa. Furthermore, waist-to-hip ratio is causally related to insulin resistance in women, but not in men. These results position dissection of metabolic and glycemic health sex dimorphism as a steppingstone for understanding differences in genetic effects between women and men in related phenotypes.

Genetic predisposition to hypertension is associated with preeclampsia in European and Central Asian women.

Preeclampsia is a serious complication of pregnancy, affecting both maternal and fetal health. In genome-wide association meta-analysis of European and Central Asian mothers, we identify sequence variants that associate with preeclampsia in the maternal genome at ZNF831/20q13 and FTO/16q12. These are previously established variants for blood pressure (BP) and the FTO variant has also been associated with body mass index (BMI). Further analysis of BP variants establishes that variants at MECOM/3q26, FGF5/4q21 and SH2B3/12q24 also associate with preeclampsia through the maternal genome. We further show that a polygenic risk score for hypertension associates with preeclampsia. However, comparison with gestational hypertension indicates that additional factors modify the risk of preeclampsia.

Genetic variant effects on gene expression in human pancreatic islets and their implications for T2D.

Most signals detected by genome-wide association studies map to non-coding sequence and their tissue-specific effects influence transcriptional regulation. However, key tissues and cell-types required for functional inference are absent from large-scale resources. Here we explore the relationship between genetic variants influencing predisposition to type 2 diabetes (T2D) and related glycemic traits, and human pancreatic islet transcription using data from 420 donors. We find: (a) 7741 cis-eQTLs in islets with a replication rate across 44 GTEx tissues between 40% and 73%; (b) marked overlap between islet cis-eQTL signals and active regulatory sequences in islets, with reduced eQTL effect size observed in the stretch enhancers most strongly implicated in GWAS signal location; (c) enrichment of islet cis-eQTL signals with T2D risk variants identified in genome-wide association studies; and (d) colocalization between 47 islet cis-eQTLs and variants influencing T2D or glycemic traits, including DGKB and TCF7L2. Our findings illustrate the advantages of performing functional and regulatory studies in disease relevant tissues.

Genome-wide Study Identifies Association between HLA-B∗55:01 and Self-Reported Penicillin Allergy.

Hypersensitivity reactions to drugs are often unpredictable and can be life threatening, underscoring a need for understanding their underlying mechanisms and risk factors. The extent to which germline genetic variation influences the risk of commonly reported drug allergies such as penicillin allergy remains largely unknown. We extracted data from the electronic health records of more than 600,000 participants from the UK, Estonian, and Vanderbilt University Medical Center's BioVU biobanks to study the role of genetic variation in the occurrence of self-reported penicillin hypersensitivity reactions. We used imputed SNP to HLA typing data from these cohorts to further fine map the human leukocyte antigen (HLA) association and replicated our results in 23andMe's research cohort involving a total of 1.12 million individuals. Genome-wide meta-analysis of penicillin allergy revealed two loci, including one located in the HLA region on chromosome 6. This signal was further fine-mapped to the HLA-B∗55:01 allele (OR 1.41 95% CI 1.33-1.49, p value 2.04 × 10-31) and confirmed by independent replication in 23andMe's research cohort (OR 1.30 95% CI 1.25-1.34, p value 1.00 × 10-47). The lead SNP was also associated with lower lymphocyte counts and in silico follow-up suggests a potential effect on T-lymphocytes at HLA-B∗55:01. We also observed a significant hit in PTPN22 and the GWAS results correlated with the genetics of rheumatoid arthritis and psoriasis. We present robust evidence for the role of an allele of the major histocompatibility complex (MHC) I gene HLA-B in the occurrence of penicillin allergy.

Glucocorticoid induces human beta cell dysfunction by involving riborepressor GAS5 LincRNA.

OBJECTIVE: A widely recognized metabolic side effect of glucocorticoid (GC) therapy is steroid-induced diabetes mellitus (DM). However, studies on the molecular basis of GC-induced pancreatic beta cell dysfunction in human beta cells are lacking. The significance of non-coding RNAs in various cellular processes is emerging. In this study, we aimed to show the direct negative impact of GC on beta cell function and elucidate the role of riborepressor GAS5 lincRNA in the GC signaling pathway in human pancreatic beta cells. METHODS: Patients undergoing two weeks of high-dose prednisolone therapy were monitored for C-peptide levels. Human pancreatic islets and the human beta cell line EndoC-ßH1 were incubated in pharmacological concentrations of dexamethasone. The GAS5 level was modulated using anti-sense LNA gapmeR or short oligonucleotides with GAS5 HREM (hormone response element motif). Immunoblotting and/or real-time PCR were used to assess changes in protein and RNA expression, respectively. Functional characterization included glucose-stimulated insulin secretion and apoptosis assays. Correlation analysis was performed on RNAseq data of human pancreatic islets. RESULTS: We found reduced C-peptide levels in patients undergoing high-dose GC therapy. Human islets and the human beta cell line EndoC-ßH1 exposed to GC exhibited reduced insulin secretion and increased apoptosis. Concomitantly, reduced expression of important beta cell transcription factors, PDX1 and NKX6-1, as well as exocytotic protein SYT13 were observed. The expression of the glucocorticoid receptor was decreased, while that of serum and glucocorticoid-regulated kinase 1 (SGK1) was elevated. The expression of these genes was found to significantly correlate with GAS5 in human islet transcriptomics data. Increasing GAS5 levels using GAS5 HREM alleviated the inhibitory effects of dexamethasone on insulin secretion. CONCLUSIONS: The direct adverse effect of glucocorticoid in human beta cell function is mediated via important beta cell proteins and components of the GC signaling pathway in an intricate interplay with GAS5 lincRNA, a potentially novel therapeutic target to counter GC-mediated beta cell dysfunction.

Variants in the fetal genome near pro-inflammatory cytokine genes on 2q13 associate with gestational duration.

The duration of pregnancy is influenced by fetal and maternal genetic and non-genetic factors. Here we report a fetal genome-wide association meta-analysis of gestational duration, and early preterm, preterm, and postterm birth in 84,689 infants. One locus on chromosome 2q13 is associated with gestational duration; the association is replicated in 9,291 additional infants (combined P = 3.96 × 10-14). Analysis of 15,588 mother-child pairs shows that the association is driven by fetal rather than maternal genotype. Functional experiments show that the lead SNP, rs7594852, alters the binding of the HIC1 transcriptional repressor. Genes at the locus include several interleukin 1 family members with roles in pro-inflammatory pathways that are central to the process of parturition. Further understanding of the underlying mechanisms will be of great public health importance, since giving birth either before or after the window of term gestation is associated with increased morbidity and mortality.

GNAS gene is an important regulator of insulin secretory capacity in pancreatic ß-cells.

BACKGROUND: Type 2 diabetes (T2D) is a complex polygenic disease with unclear mechanism. In an attempt to identify novel genes involved in ß-cell function, we harness a bioinformatics method called Loss-of-function tool (LoFtool) gene score. METHODS: RNA-sequencing data from human islets were used to cross-reference genes within the 1st quartile of most intolerant LoFtool score with the 100th most expressed genes in human islets. Out of these genes, GNAS and EEF1A1 genes were selected for further investigation in diabetic islets, metabolic tissues along with their correlation with diabetic phenotypes. The influence of GNAS and EEF1A1 on insulin secretion and ß-cell function were validated in INS-1 cells. RESULTS: A comparatively higher expression level of GNAS and EEF1A1 was observed in human islets than fat, liver and muscle tissues. Furthermore, diabetic islets displayed a reduced expression of GNAS, but not of EEF1A, compared to non-diabetic islets. The expression of GNAS was positively correlated with insulin secretory index, GLP1R, GIPR and inversely correlated with HbA1c. Diabetic human islets displayed a reduced cAMP generation and insulin secretory capacity in response to glucose. Moreover, siRNA silencing of GNAS in INS-1 cells reduced insulin secretion, insulin content, and cAMP production. In addition, the expression of Insulin, PDX1, and MAFA was significantly down-regulated in GNAS-silenced cells. However, cell viability and apoptosis rate were unaffected. CONCLUSION: LoFtool is a powerful tool to identify genes associated with pancreatic islets dysfunction. GNAS is a crucial gene for the ß-cell insulin secretory capacity.

Genome-wide meta-analysis identifies BARX1 and EML4-MTA3 as new loci associated with infantile hypertrophic pyloric stenosis.

Infantile hypertrophic pyloric stenosis (IHPS) is a disorder of young infants with a population incidence of ∼2/1000 live births, caused by hypertrophy of the pyloric sphincter smooth muscle. Reported genetic loci associated with IHPS explain only a minor proportion of IHPS risk. To identify new risk loci, we carried out a genome-wide meta-analysis on 1395 surgery-confirmed cases and 4438 controls, with replication in a set of 2427 cases and 2524 controls. We identified and replicated six independent genomic loci associated with IHPS risk at genome wide significance (P < 5 × 10-8), including novel associations with two single nucleotide polymorphisms (SNPs). One of these SNPs, rs6736913 [odds ratio (OR) = 2.32; P = 3.0 × 10-15], is a low frequency missense variant in EML4 at 2p21. The second SNP, rs1933683 (OR = 1.34; P = 3.1 × 10-9) is 1 kb downstream of BARX1 at 9q22.32, an essential gene for stomach formation in embryogenesis. Using the genome-wide complex trait analysis method, we estimated the IHPS SNP heritability to be 30%, and using the linkage disequilibrium score regression method, we found support for a previously reported genetic correlation of IHPS with lipid metabolism. By combining the largest collection of IHPS cases to date (3822 cases), with results generalized across populations of different ancestry, we elucidate novel mechanistic avenues of IHPS disease architecture.

Controllability in an islet specific regulatory network identifies the transcriptional factor NFATC4, which regulates Type 2 Diabetes associated genes.

Probing the dynamic control features of biological networks represents a new frontier in capturing the dysregulated pathways in complex diseases. Here, using patient samples obtained from a pancreatic islet transplantation program, we constructed a tissue-specific gene regulatory network and used the control centrality (Cc) concept to identify the high control centrality (HiCc) pathways, which might serve as key pathobiological pathways for Type 2 Diabetes (T2D). We found that HiCc pathway genes were significantly enriched with modest GWAS p-values in the DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) study. We identified variants regulating gene expression (expression quantitative loci, eQTL) of HiCc pathway genes in islet samples. These eQTL genes showed higher levels of differential expression compared to non-eQTL genes in low, medium, and high glucose concentrations in rat islets. Among genes with highly significant eQTL evidence, NFATC4 belonged to four HiCc pathways. We asked if the expressions of T2D-associated candidate genes from GWAS and literature are regulated by Nfatc4 in rat islets. Extensive in vitro silencing of Nfatc4 in rat islet cells displayed reduced expression of 16, and increased expression of four putative downstream T2D genes. Overall, our approach uncovers the mechanistic connection of NFATC4 with downstream targets including a previously unknown one, TCF7L2, and establishes the HiCc pathways' relationship to T2D.

Genome-wide association study of Hirschsprung disease detects a novel low-frequency variant at the RET locus.

Hirschsprung disease (HSCR) is a congenital disorder with a population incidence of ~1/5000 live births, defined by an absence of enteric ganglia along variable lengths of the colon. HSCR genome-wide association studies (GWAS) have found common associated variants at RET, SEMA3, and NRG1, but they still fail to explain all of its heritability. To enhance gene discovery, we performed a GWAS of 170 cases identified from the Danish nationwide pathology registry with 4717 controls, based on 6.2 million variants imputed from the haplotype reference consortium panel. We found a novel low-frequency variant (rs144432435), which, when conditioning on the lead RET single-nucleotide polymorphism (SNP), was of genome-wide significance in the discovery analysis. This conditional association signal was replicated in a Swedish HSCR cohort with discovery plus replication meta-analysis conditional odds ratio of 6.6 (P = 7.7 × 10-10; 322 cases and 4893 controls). The conditional signal was, however, not replicated in two HSCR cohorts from USA and Finland, leading to the hypothesis that rs144432435 tags a rare haplotype present in Denmark and Sweden. Using the genome-wide complex trait analysis method, we estimated the SNP heritability of HSCR to be 88%, close to estimates based on classical family studies. Moreover, by using Lasso (least absolute shrinkage and selection operator) regression we were able to construct a genetic HSCR predictor with a area under the receiver operator characteristics curve of 76% in an independent validation set. In conclusion, we combined the largest collection of sporadic Hirschsprung cases to date (586 cases) to further elucidate HSCR's genetic architecture.

Glucose-Induced Changes in Gene Expression in Human Pancreatic Islets: Causes or Consequences of Chronic Hyperglycemia.

Dysregulation of gene expression in islets from patients with type 2 diabetes (T2D) might be causally involved in the development of hyperglycemia, or it could develop as a consequence of hyperglycemia (i.e., glucotoxicity). To separate the genes that could be causally involved in pathogenesis from those likely to be secondary to hyperglycemia, we exposed islets from human donors to normal or high glucose concentrations for 24 h and analyzed gene expression. We compared these findings with gene expression in islets from donors with normal glucose tolerance and hyperglycemia (including T2D). The genes whose expression changed in the same direction after short-term glucose exposure, as in T2D, were considered most likely to be a consequence of hyperglycemia. Genes whose expression changed in hyperglycemia but not after short-term glucose exposure, particularly those that also correlated with insulin secretion, were considered the strongest candidates for causal involvement in T2D. For example, ERO1LB, DOCK10, IGSF11, and PRR14L were downregulated in donors with hyperglycemia and correlated positively with insulin secretion, suggesting a protective role, whereas TMEM132C was upregulated in hyperglycemia and correlated negatively with insulin secretion, suggesting a potential pathogenic role. This study provides a catalog of gene expression changes in human pancreatic islets after exposure to glucose.

The Genetic Landscape of Renal Complications in Type 1 Diabetes.

Diabetes is the leading cause of ESRD. Despite evidence for a substantial heritability of diabetic kidney disease, efforts to identify genetic susceptibility variants have had limited success. We extended previous efforts in three dimensions, examining a more comprehensive set of genetic variants in larger numbers of subjects with type 1 diabetes characterized for a wider range of cross-sectional diabetic kidney disease phenotypes. In 2843 subjects, we estimated that the heritability of diabetic kidney disease was 35% (P=6.4×10-3). Genome-wide association analysis and replication in 12,540 individuals identified no single variants reaching stringent levels of significance and, despite excellent power, provided little independent confirmation of previously published associated variants. Whole-exome sequencing in 997 subjects failed to identify any large-effect coding alleles of lower frequency influencing the risk of diabetic kidney disease. However, sets of alleles increasing body mass index (P=2.2×10-5) and the risk of type 2 diabetes (P=6.1×10-4) associated with the risk of diabetic kidney disease. We also found genome-wide genetic correlation between diabetic kidney disease and failure at smoking cessation (P=1.1×10-4). Pathway analysis implicated ascorbate and aldarate metabolism (P=9.0×10-6), and pentose and glucuronate interconversions (P=3.0×10-6) in pathogenesis of diabetic kidney disease. These data provide further evidence for the role of genetic factors influencing diabetic kidney disease in those with type 1 diabetes and highlight some key pathways that may be responsible. Altogether these results reveal important biology behind the major cause of kidney disease.