Urinary metabolite profiling and risk of progression of diabetic nephropathy in 2670 individuals with type 1 diabetes.

AIMS/HYPOTHESIS: This prospective, observational study examines associations between 51 urinary metabolites and risk of progression of diabetic nephropathy in individuals with type 1 diabetes by employing an automated NMR metabolomics technique suitable for large-scale urine sample collections. METHODS: We collected 24-h urine samples for 2670 individuals with type 1 diabetes from the Finnish Diabetic Nephropathy study and measured metabolite concentrations by NMR. Individuals were followed up for 9.0 ± 5.0 years until their first sign of progression of diabetic nephropathy, end-stage kidney disease or study end. Cox regressions were performed on the entire study population (overall progression), on 1999 individuals with normoalbuminuria and 347 individuals with macroalbuminuria at baseline. RESULTS: Seven urinary metabolites were associated with overall progression after adjustment for baseline albuminuria and chronic kidney disease stage (p < 8 × 10-4): leucine (HR 1.47 [95% CI 1.30, 1.66] per 1-SD creatinine-scaled metabolite concentration), valine (1.38 [1.22, 1.56]), isoleucine (1.33 [1.18, 1.50]), pseudouridine (1.25 [1.11, 1.42]), threonine (1.27 [1.11, 1.46]) and citrate (0.84 [0.75, 0.93]). 2-Hydroxyisobutyrate was associated with overall progression (1.30 [1.16, 1.45]) and also progression from normoalbuminuria (1.56 [1.25, 1.95]). Six amino acids and pyroglutamate were associated with progression from macroalbuminuria. CONCLUSIONS/INTERPRETATION: Branched-chain amino acids and other urinary metabolites were associated with the progression of diabetic nephropathy on top of baseline albuminuria and chronic kidney disease. We found differences in associations for overall progression and progression from normo- and macroalbuminuria. These novel discoveries illustrate the utility of analysing urinary metabolites in entire population cohorts.

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

AIMS/HYPOTHESIS: Diabetic kidney disease (DKD) is the leading cause of kidney failure and has a substantial genetic component. Our aim was to identify novel genetic factors and genes contributing to DKD by performing meta-analysis of previous genome-wide association studies (GWAS) on DKD and by integrating the results with renal transcriptomics datasets. METHODS: We performed GWAS meta-analyses using ten phenotypic definitions of DKD, including nearly 27,000 individuals with diabetes. Meta-analysis results were integrated with estimated quantitative trait locus data from human glomerular (N=119) and tubular (N=121) samples to perform transcriptome-wide association study. We also performed gene aggregate tests to jointly test all available common genetic markers within a gene, and combined the results with various kidney omics datasets. RESULTS: The meta-analysis identified a novel intronic variant (rs72831309) in the TENM2 gene associated with a lower risk of the combined chronic kidney disease (eGFR<60 ml/min per 1.73 m2) and DKD (microalbuminuria or worse) phenotype (p=9.8×10-9; although not withstanding correction for multiple testing, p>9.3×10-9). Gene-level analysis identified ten genes associated with DKD (COL20A1, DCLK1, EIF4E, PTPRN-RESP18, GPR158, INIP-SNX30, LSM14A and MFF; p<2.7×10-6). Integration of GWAS with human glomerular and tubular expression data demonstrated higher tubular AKIRIN2 gene expression in individuals with vs without DKD (p=1.1×10-6). The lead SNPs within six loci significantly altered DNA methylation of a nearby CpG site in kidneys (p<1.5×10-11). Expression of lead genes in kidney tubules or glomeruli correlated with relevant pathological phenotypes (e.g. TENM2 expression correlated positively with eGFR [p=1.6×10-8] and negatively with tubulointerstitial fibrosis [p=2.0×10-9], tubular DCLK1 expression correlated positively with fibrosis [p=7.4×10-16], and SNX30 expression correlated positively with eGFR [p=5.8×10-14] and negatively with fibrosis [p<2.0×10-16]). CONCLUSIONS/INTERPRETATION: Altogether, the results point to novel genes contributing to the pathogenesis of DKD. DATA AVAILABILITY: The GWAS meta-analysis results can be accessed via the type 1 and type 2 diabetes (T1D and T2D, respectively) and Common Metabolic Diseases (CMD) Knowledge Portals, and downloaded on their respective download pages ( https://t1d.hugeamp.org/downloads.html ; https://t2d.hugeamp.org/downloads.html ; https://hugeamp.org/downloads.html ).

Association of Coding Variants in Hydroxysteroid 17-beta Dehydrogenase 14 (HSD17B14) with Reduced Progression to End Stage Kidney Disease in Type 1 Diabetes.

BACKGROUND: Rare variants in gene coding regions likely have a greater impact on disease-related phenotypes than common variants through disruption of their encoded protein. We searched for rare variants associated with onset of ESKD in individuals with type 1 diabetes at advanced kidney disease stage. METHODS: Gene-based exome array analyses of 15,449 genes in five large incidence cohorts of individuals with type 1 diabetes and proteinuria were analyzed for survival time to ESKD, testing the top gene in a sixth cohort (n=2372/1115 events all cohorts) and replicating in two retrospective case-control studies (n=1072 cases, 752 controls). Deep resequencing of the top associated gene in five cohorts confirmed the findings. We performed immunohistochemistry and gene expression experiments in human control and diseased cells, and in mouse ischemia reperfusion and aristolochic acid nephropathy models. RESULTS: Protein coding variants in the hydroxysteroid 17-ß dehydrogenase 14 gene (HSD17B14), predicted to affect protein structure, had a net protective effect against development of ESKD at exome-wide significance (n=4196; P value=3.3 × 10-7). The HSD17B14 gene and encoded enzyme were robustly expressed in healthy human kidney, maximally in proximal tubular cells. Paradoxically, gene and protein expression were attenuated in human diabetic proximal tubules and in mouse kidney injury models. Expressed HSD17B14 gene and protein levels remained low without recovery after 21 days in a murine ischemic reperfusion injury model. Decreased gene expression was found in other CKD-associated renal pathologies. CONCLUSIONS: HSD17B14 gene is mechanistically involved in diabetic kidney disease. The encoded sex steroid enzyme is a druggable target, potentially opening a new avenue for therapeutic development.

Whole-Genome Sequencing of Finnish Type 1 Diabetic Siblings Discordant for Kidney Disease Reveals DNA Variants associated with Diabetic Nephropathy.

BACKGROUND: Several genetic susceptibility loci associated with diabetic nephropathy have been documented, but no causative variants implying novel pathogenetic mechanisms have been elucidated. METHODS: We carried out whole-genome sequencing of a discovery cohort of Finnish siblings with type 1 diabetes who were discordant for the presence (case) or absence (control) of diabetic nephropathy. Controls had diabetes without complications for 15-37 years. We analyzed and annotated variants at genome, gene, and single-nucleotide variant levels. We then replicated the associated variants, genes, and regions in a replication cohort from the Finnish Diabetic Nephropathy study that included 3531 unrelated Finns with type 1 diabetes. RESULTS: We observed protein-altering variants and an enrichment of variants in regions associated with the presence or absence of diabetic nephropathy. The replication cohort confirmed variants in both regulatory and protein-coding regions. We also observed that diabetic nephropathy-associated variants, when clustered at the gene level, are enriched in a core protein-interaction network representing proteins essential for podocyte function. These genes include protein kinases (protein kinase C isoforms ε and ι) and protein tyrosine kinase 2. CONCLUSIONS: Our comprehensive analysis of a diabetic nephropathy cohort of siblings with type 1 diabetes who were discordant for kidney disease points to variants and genes that are potentially causative or protective for diabetic nephropathy. This includes variants in two isoforms of the protein kinase C family not previously linked to diabetic nephropathy, adding support to previous hypotheses that the protein kinase C family members play a role in diabetic nephropathy and might be attractive therapeutic targets.

A Targeted Multiomics Approach to Identify Biomarkers Associated with Rapid eGFR Decline in Type 1 Diabetes.

BACKGROUND: Individuals with type 1 diabetes (T1D) demonstrate varied trajectories of estimated glomerular filtration rate (eGFR) decline. The molecular pathways underlying rapid eGFR decline in T1D are poorly understood, and individual-level risk of rapid eGFR decline is difficult to predict. METHODS: We designed a case-control study with multiple exposure measurements nested within 4 well-characterized T1D cohorts (FinnDiane, Steno, EDC, and CACTI) to identify biomarkers associated with rapid eGFR decline. Here, we report the rationale for and design of these studies as well as results of models testing associations of clinical characteristics with rapid eGFR decline in the study population, upon which "omics" studies will be built. Cases (n = 535) and controls (n = 895) were defined as having an annual eGFR decline of ≥3 and <1 mL/min/1.73 m2, respectively. Associations of demographic and clinical variables with rapid eGFR decline were tested using logistic regression, and prediction was evaluated using area under the curve (AUC) statistics. Targeted metabolomics, lipidomics, and proteomics are being performed using high-resolution mass-spectrometry techniques. RESULTS: At baseline, the mean age was 43 years, diabetes duration was 27 years, eGFR was 94 mL/min/1.73 m2, and 62% of participants were normoalbuminuric. Over 7.6-year median follow-up, the mean annual change in eGFR in cases and controls was -5.7 and 0.6 mL/min/1.73 m2, respectively. Younger age, longer diabetes duration, and higher baseline HbA1c, urine albumin-creatinine ratio, and eGFR were significantly associated with rapid eGFR decline. The cross-validated AUC for the predictive model incorporating these variables plus sex and mean arterial blood pressure was 0.74 (95% CI: 0.68-0.79; p < 0.001). CONCLUSION: Known risk factors provide moderate discrimination of rapid eGFR decline. Identification of blood and urine biomarkers associated with rapid eGFR decline in T1D using targeted omics strategies may provide insight into disease mechanisms and improve upon clinical predictive models using traditional risk factors.

ErbB2-driven breast cancer cell invasion depends on a complex signaling network activating myeloid zinc finger-1-dependent cathepsin B expression.

Aberrant ErbB2 receptor tyrosine kinase activation in breast cancer is strongly linked to an invasive disease. The molecular basis of ErbB2-driven invasion is largely unknown. We show that cysteine cathepsins B and L are elevated in ErbB2 positive primary human breast cancer and function as effectors of ErbB2-induced invasion in vitro. We identify Cdc42-binding protein kinase beta, extracellular regulated kinase 2, p21-activated protein kinase 4, and protein kinase C alpha as essential mediators of ErbB2-induced cysteine cathepsin expression and breast cancer cell invasiveness. The identified signaling network activates the transcription of cathepsin B gene (CTSB) via myeloid zinc finger-1 transcription factor that binds to an ErbB2-responsive enhancer element in the first intron of CTSB. This work provides a model system for ErbB2-induced breast cancer cell invasiveness, reveals a signaling network that is crucial for invasion in vitro, and defines a specific role and targets for the identified serine-threonine kinases.

Genome-Wide Association Study of Diabetic Kidney Disease Highlights Biology Involved in Glomerular Basement Membrane Collagen.

BACKGROUND: Although diabetic kidney disease demonstrates both familial clustering and single nucleotide polymorphism heritability, the specific genetic factors influencing risk remain largely unknown. METHODS: To identify genetic variants predisposing to diabetic kidney disease, we performed genome-wide association study (GWAS) analyses. Through collaboration with the Diabetes Nephropathy Collaborative Research Initiative, we assembled a large collection of type 1 diabetes cohorts with harmonized diabetic kidney disease phenotypes. We used a spectrum of ten diabetic kidney disease definitions based on albuminuria and renal function. RESULTS: Our GWAS meta-analysis included association results for up to 19,406 individuals of European descent with type 1 diabetes. We identified 16 genome-wide significant risk loci. The variant with the strongest association (rs55703767) is a common missense mutation in the collagen type IV alpha 3 chain (COL4A3) gene, which encodes a major structural component of the glomerular basement membrane (GBM). Mutations in COL4A3 are implicated in heritable nephropathies, including the progressive inherited nephropathy Alport syndrome. The rs55703767 minor allele (Asp326Tyr) is protective against several definitions of diabetic kidney disease, including albuminuria and ESKD, and demonstrated a significant association with GBM width; protective allele carriers had thinner GBM before any signs of kidney disease, and its effect was dependent on glycemia. Three other loci are in or near genes with known or suggestive involvement in this condition (BMP7) or renal biology (COLEC11 and DDR1). CONCLUSIONS: The 16 diabetic kidney disease-associated loci may provide novel insights into the pathogenesis of this condition and help identify potential biologic targets for prevention and treatment.

Biomarker panels associated with progression of renal disease in type 1 diabetes.

AIMS/HYPOTHESIS: We aimed to identify a sparse panel of biomarkers for improving the prediction of renal disease progression in type 1 diabetes. METHODS: We considered 859 individuals recruited from the Scottish Diabetes Research Network Type 1 Bioresource (SDRNT1BIO) and 315 individuals from the Finnish Diabetic Nephropathy (FinnDiane) study. All had an entry eGFR between 30 and 75 ml min-1[1.73 m]-2, with those from FinnDiane being oversampled for albuminuria. A total of 297 circulating biomarkers (30 proteins, 121 metabolites, 146 tryptic peptides) were measured in non-fasting serum samples using the Luminex platform and LC electrospray tandem MS (LC-MS/MS). We investigated associations with final eGFR adjusted for baseline eGFR and with rapid progression (a loss of more than 3 ml min-1[1.73 m]-2 year-1) using linear and logistic regression models. Panels of biomarkers were identified using a penalised Bayesian approach, and their performance was evaluated through 10-fold cross-validation and compared with using clinical record data alone. RESULTS: For final eGFR, 16 proteins and 30 metabolites or tryptic peptides showed significant association in SDRNT1BIO, and nine proteins and five metabolites or tryptic peptides in FinnDiane, beyond age, sex, diabetes duration, study day eGFR and length of follow-up (all at p < 10-4). The strongest associations were with CD27 antigen (CD27), kidney injury molecule 1 (KIM-1) and α1-microglobulin. Including the Luminex biomarkers on top of baseline covariates increased the r2 for prediction of final eGFR from 0.47 to 0.58 in SDRNT1BIO and from 0.33 to 0.48 in FinnDiane. At least 75% of the increment in r2 was attributable to CD27 and KIM-1. However, using the weighted average of historical eGFR gave similar performance to biomarkers. The LC-MS/MS platform performed less well. CONCLUSIONS/INTERPRETATION: Among a large set of associated biomarkers, a sparse panel of just CD27 and KIM-1 contains most of the predictive information for eGFR progression. The increment in prediction beyond clinical data was modest but potentially useful for oversampling individuals with rapid disease progression into clinical trials, especially where there is little information on prior eGFR trajectories.

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.

Early Trichinella spiralis and Trichinella nativa infections induce similar gene expression profiles in rat jejunal mucosa.

Trichinella spiralis causes a significantly higher parasite burden in rat muscle than Trichinella nativa. To assess whether the difference in infectivity is due to the early intestinal response, we analyzed gene expression changes in the rat jejunum during Trichinella infection with a whole-genome microarray. The rats were euthanized on day five of infection, and their jejunal mucosa was sampled for microarray analysis. In addition, intestinal histology and hematology were examined. Against our expectations, the gene expression changes were similar in both T.nativa- and T. spiralis-infected groups. The two groups were hence pooled, and in the combined Trichinella-infected group, 551 genes were overexpressed and 427 underexpressed when compared to controls (false discovery rate ≤ 0.001 and fold change at least 2 in either direction). Pathway analysis identified seven pathways significantly associated with Trichinella infection (p < 0.05). The microarray data suggested nonspecific damage and an inflammatory response in the jejunal mucosa. Histological findings, including hyperemia, hemorrhage and a marked infiltration of inflammatory cells, supported the microarray data. Trichinella infection caused complex gene expression changes that indicate a host response to tissue damage in the mucosa of the jejunum, but the changes were not notably dependent on the studied species of Trichinella.

Rate of Kidney Function Decline is Associated With Kidney and Heart Failure in Individuals With Type 1 Diabetes.

Introduction: Diabetes is the most common cause of chronic kidney disease (CKD). Urinary albumin excretion rate (AER) and estimated glomerular filtration rate (eGFR) are commonly used to monitor the onset and progression of diabetic kidney disease (DKD). We studied if the preceding rate of kidney function decline, that is, the eGFR slope, is independently associated with incident clinical cardiorenal events. Methods: This study included longitudinal data for 2498 Finnish individuals with type 1 diabetes (T1D). The eGFR slope was calculated from 5 years preceding the study visit. Data on kidney failure, coronary heart disease (CHD), stroke, 3-point major adverse cardiovascular events (MACE), heart failure, and death were obtained from national registries. The associations between the eGFR slope and incident events were assessed with multivariable competing risk models during the average follow-up of 9.2 years. Results: The eGFR slopes were associated (P ≤ 0.001) with all outcomes when adjusted for age, sex, and HbA1c. However, eGFR slope remained associated only with the composite outcome of kidney failure or death when the albuminuria group and eGFR at the study visit were included in the model (P = 0.041). In addition, eGFR slope was independently associated with kidney failure in individuals without CKD (eGFR > 60 ml/min per 1.73 m2; P = 0.044), and with heart failure in those with CKD (P = 0.033). However, eGFR slope did not markedly improve the model C-index. Conclusion: The eGFR slope was independently associated with kidney failure in those without CKD, and with heart failure in those with CKD. However, it is unlikely to have major relevance for clinical practice when the current eGFR and albuminuria status are known.

Medication profiling in women with type 1 diabetes highlights the importance of adequate, guideline-based treatment in low-risk groups.

Effective treatment may prevent kidney complications, but women might be underprescribed. Novel, data-driven insights into prescriptions and their relationship with kidney health in women with type 1 diabetes may help to optimize treatment. We identified six medication profiles in 1164 women from the FinnDiane Study with normal albumin excretion rate based on clusters of their baseline prescription data using a self-organizing map. Future rapid kidney function decline was defined as an annual estimated glomerular filtration rate (eGFR) loss > 3 ml/min/1.73 m2 after baseline. Two profiles were associated with future decline: Profile ARB with the highest proportion of angiotensin receptor blockers (odds ratio [OR] 2.75, P = 0.02) and highly medicated women in profile HighMed (OR 2.55, P = 0.03). Compared with profile LowMed (low purchases of all), profile HighMed had worse clinical characteristics, whereas in profile ARB only systolic blood pressure was elevated. Importantly, the younger women in profile ARB with fewer kidney protective treatments developed a rapid decline despite otherwise similar baseline characteristics to profile ACE & Lipids (the highest proportions of ACE inhibitors and lipid-modifying agents) without a future rapid decline. In conclusion, medication profiles identified different future eGFR trajectories in women with type 1 diabetes revealing potential treatment gaps for younger women.

Genome-wide mRNA profiling in urinary extracellular vesicles reveals stress gene signature for diabetic kidney disease.

Urinary extracellular vesicles (uEV) are a largely unexplored source of kidney-derived mRNAs with potential to serve as a liquid kidney biopsy. We assessed ∼200 uEV mRNA samples from clinical studies by genome-wide sequencing to discover mechanisms and candidate biomarkers of diabetic kidney disease (DKD) in Type 1 diabetes (T1D) with replication in Type 1 and 2 diabetes. Sequencing reproducibly showed >10,000 mRNAs with similarity to kidney transcriptome. T1D DKD groups showed 13 upregulated genes prevalently expressed in proximal tubules, correlated with hyperglycemia and involved in cellular/oxidative stress homeostasis. We used six of them (GPX3, NOX4, MSRB, MSRA, HRSP12, and CRYAB) to construct a transcriptional "stress score" that reflected long-term decline of kidney function and could even identify normoalbuminuric individuals showing early decline. We thus provide workflow and web resource for studying uEV transcriptomes in clinical urine samples and stress-linked DKD markers as potential early non-invasive biomarkers or drug targets.

Genetic Risk Score Enhances Coronary Artery Disease Risk Prediction in Individuals With Type 1 Diabetes.

OBJECTIVE: Individuals with type 1 diabetes are at a high lifetime risk of coronary artery disease (CAD), calling for early interventions. This study explores the use of a genetic risk score (GRS) for CAD risk prediction, compares it to established clinical markers, and investigates its performance according to the age and pharmacological treatment. RESEARCH DESIGN AND METHODS: This study in 3,295 individuals with type 1 diabetes from the Finnish Diabetic Nephropathy Study (467 incident CAD, 14.8 years follow-up) used three risk scores: a GRS, a validated clinical score, and their combined score. Hazard ratios (HR) were calculated with Cox regression, and model performances were compared with the Harrell C-index (C-index). RESULTS: A HR of 6.7 for CAD was observed between the highest and the lowest 5th percentile of the GRS (P = 1.8 × 10-6). The performance of GRS (C-index = 0.562) was similar to HbA1c (C-index = 0.563, P = 0.96 for difference), HDL (C-index = 0.571, P = 0.6), and total cholesterol (C-index = 0.594, P = 0.1). The GRS was not correlated with the clinical score (r = -0.013, P = 0.5). The combined score outperformed the clinical score (C-index = 0.813 vs. C-index = 0.820, P = 0.003). The GRS performed better in individuals below the median age (38.6 years) compared with those above (C-index = 0.637 vs. C-index = 0.546). CONCLUSIONS: A GRS identified individuals at high risk of CAD and worked better in younger individuals. GRS was also an independent risk factor for CAD, with a predictive power comparable to that of HbA1c and HDL and total cholesterol, and when incorporated into a clinical model, modestly improved the predictions. The GRS promises early risk stratification in clinical practice by enhancing the prediction of CAD.

Effect of serum sample storage temperature on metabolomic and proteomic biomarkers.

Prospective biomarker studies can be used to identify biomarkers predictive of disease onset. However, if serum biomarkers are measured years after their collection, the storage conditions might affect analyte concentrations. Few data exists concerning which metabolites and proteins are affected by storage at - 20 °C vs - 80 °C. Our objectives were to document analytes affected by storage of serum samples at - 20 °C vs - 80 °C, and to identify those indicative of the storage temperature. We utilized liquid chromatography tandem mass spectrometry and Luminex to quantify 300 analytes from serum samples of 16 Finnish individuals with type 1 diabetes, with split-aliquot samples stored at - 80 °C and - 20 °C for a median of 4.2 years. Results were validated in 315 Finnish and 916 Scottish individuals with type 1 diabetes, stored at - 20 °C and at - 80 °C, respectively. After quality control, we analysed 193 metabolites and proteins of which 120 were apparently unaffected and 15 clearly susceptible to storage at - 20 °C vs - 80 °C. Further, we identified serum glutamate/glutamine ratio greater than 0.20 as a biomarker of storage at - 20 °C vs - 80 °C. The results provide a catalogue of analytes unaffected and affected by storage at - 20 °C vs - 80 °C and biomarkers indicative of sub-optimal storage.

Whole-exome sequencing identifies novel protein-altering variants associated with serum apolipoprotein and lipid concentrations.

BACKGROUND: Dyslipidemia is a major risk factor for cardiovascular disease, and diabetes impacts the lipid metabolism through multiple pathways. In addition to the standard lipid measurements, apolipoprotein concentrations provide added awareness of the burden of circulating lipoproteins. While common genetic variants modestly affect the serum lipid concentrations, rare genetic mutations can cause monogenic forms of hypercholesterolemia and other genetic disorders of lipid metabolism. We aimed to identify low-frequency protein-altering variants (PAVs) affecting lipoprotein and lipid traits. METHODS: We analyzed whole-exome (WES) and whole-genome sequencing (WGS) data of 481 and 474 individuals with type 1 diabetes, respectively. The phenotypic data consisted of 79 serum lipid and apolipoprotein phenotypes obtained with clinical laboratory measurements and nuclear magnetic resonance spectroscopy. RESULTS: The single-variant analysis identified an association between the LIPC p.Thr405Met (rs113298164) and serum apolipoprotein A1 concentrations (p=7.8×10-8). The burden of PAVs was significantly associated with lipid phenotypes in LIPC, RBM47, TRMT5, GTF3C5, MARCHF10, and RYR3 (p<2.9×10-6). The RBM47 gene is required for apolipoprotein B post-translational modifications, and in our data, the association between RBM47 and apolipoprotein C-III concentrations was due to a rare 21 base pair p.Ala496-Ala502 deletion; in replication, the burden of rare deleterious variants in RBM47 was associated with lower triglyceride concentrations in WES of >170,000 individuals from multiple ancestries (p=0.0013). Two PAVs in GTF3C5 were highly enriched in the Finnish population and associated with cardiovascular phenotypes in the general population. In the previously known APOB gene, we identified novel associations at two protein-truncating variants resulting in lower serum non-HDL cholesterol (p=4.8×10-4), apolipoprotein B (p=5.6×10-4), and LDL cholesterol (p=9.5×10-4) concentrations. CONCLUSIONS: We identified lipid and apolipoprotein-associated variants in the previously known LIPC and APOB genes, as well as PAVs in GTF3C5 associated with LDLC, and in RBM47 associated with apolipoprotein C-III concentrations, implicated as an independent CVD risk factor. Identification of rare loss-of-function variants has previously revealed genes that can be targeted to prevent CVD, such as the LDL cholesterol-lowering loss-of-function variants in the PCSK9 gene. Thus, this study suggests novel putative therapeutic targets for the prevention of CVD.

Urinary Proteomics Identifies Cathepsin D as a Biomarker of Rapid eGFR Decline in Type 1 Diabetes.

OBJECTIVE: Understanding mechanisms underlying rapid estimated glomerular filtration rate (eGFR) decline is important to predict and treat kidney disease in type 1 diabetes (T1D). RESEARCH DESIGN AND METHODS: We performed a case-control study nested within four T1D cohorts to identify urinary proteins associated with rapid eGFR decline. Case and control subjects were categorized based on eGFR decline ≥3 and <1 mL/min/1.73 m2/year, respectively. We used targeted liquid chromatography-tandem mass spectrometry to measure 38 peptides from 20 proteins implicated in diabetic kidney disease. Significant proteins were investigated in complementary human cohorts and in mouse proximal tubular epithelial cell cultures. RESULTS: The cohort study included 1,270 participants followed a median 8 years. In the discovery set, only cathepsin D peptide and protein were significant on full adjustment for clinical and laboratory variables. In the validation set, associations of cathepsin D with eGFR decline were replicated in minimally adjusted models but lost significance with adjustment for albuminuria. In a meta-analysis with combination of discovery and validation sets, the odds ratio for the association of cathepsin D with rapid eGFR decline was 1.29 per SD (95% CI 1.07-1.55). In complementary human cohorts, urine cathepsin D was associated with tubulointerstitial injury and tubulointerstitial cathepsin D expression was associated with increased cortical interstitial fractional volume. In mouse proximal tubular epithelial cell cultures, advanced glycation end product-BSA increased cathepsin D activity and inflammatory and tubular injury markers, which were further increased with cathepsin D siRNA. CONCLUSIONS: Urine cathepsin D is associated with rapid eGFR decline in T1D and reflects kidney tubulointerstitial injury.

Posttranscriptional regulation of angiotensin II type 1 receptor expression by glyceraldehyde 3-phosphate dehydrogenase.

Regulation of angiotensin II type 1 receptor (AT1R) has a pathophysiological role in hypertension, atherosclerosis and heart failure. We started from an observation that the 3'-untranslated region (3'-UTR) of AT1R mRNA suppressed AT1R translation. Using affinity purification for the separation of 3'-UTR-binding proteins and mass spectrometry for their identification, we describe glyceraldehyde 3-phosphate dehydrogenase (GAPDH) as an AT1R 3'-UTR-binding protein. RNA electrophoretic mobility shift analysis with purified GAPDH further demonstrated a direct interaction with the 3'-UTR while GAPDH immunoprecipitation confirmed this interaction with endogenous AT1R mRNA. GAPDH-binding site was mapped to 1-100 of 3'-UTR. GAPDH-bound target mRNAs were identified by expression array hybridization. Analysis of secondary structures shared among GAPDH targets led to the identification of a RNA motif rich in adenines and uracils. Silencing of GAPDH increased the expression of both endogenous and transfected AT1R. Similarly, a decrease in GAPDH expression by H(2)O(2) led to an increased level of AT1R expression. Consistent with GAPDH having a central role in H(2)O(2)-mediated AT1R regulation, both the deletion of GAPDH-binding site and GAPDH overexpression attenuated the effect of H(2)O(2) on AT1R mRNA. Taken together, GAPDH is a translational suppressor of AT1R and mediates the effect of H(2)O(2) on AT1R mRNA.

Novel Linkage Peaks Discovered for Diabetic Nephropathy in Individuals With Type 1 Diabetes.

Genome-wide association studies (GWAS) and linkage studies have had limited success in identifying genome-wide significantly linked regions or risk loci for diabetic nephropathy (DN) in individuals with type 1 diabetes (T1D). As GWAS cohorts have grown, they have also included more documented and undocumented familial relationships. Here we computationally inferred and manually curated pedigrees in a study cohort of >6,000 individuals with T1D and their relatives without diabetes. We performed a linkage study for 177 pedigrees consisting of 452 individuals with T1D and their relatives using a genome-wide genotyping array with >300,000 single nucleotide polymorphisms and PSEUDOMARKER software. Analysis resulted in genome-wide significant linkage peaks on eight chromosomal regions from five chromosomes (logarithm of odds score >3.3). The highest peak was localized at the HLA region on chromosome 6p, but whether the peak originated from T1D or DN remained ambiguous. Of other significant peaks, the chromosome 4p22 region was localized on top of ARHGAP24, a gene associated with focal segmental glomerulosclerosis, suggesting this gene may play a role in DN as well. Furthermore, rare variants have been associated with DN and chronic kidney disease near the 4q25 peak, localized on top of CCSER1.

Genome-wide association study on coronary artery disease in type 1 diabetes suggests beta-defensin 127 as a risk locus.

AIMS: Diabetes is a known risk factor for coronary artery disease (CAD). There is accumulating evidence that CAD pathogenesis differs for individuals with type 1 diabetes (T1D). However, the genetic background has not been extensively studied. We aimed to discover genetic loci increasing CAD susceptibility, especially in T1D, to examine the function of these discoveries and to study the role of the known risk loci in T1D. METHODS AND RESULTS: We performed the largest genome-wide association study to date for CAD in T1D, comprising 4869 individuals with T1D (cases/controls: 941/3928). Two loci reached genome-wide significance, rs1970112 in CDKN2B-AS1 [odds ratio (OR) = 1.32, P = 1.50 × 10-8], and rs6055069 on DEFB127 promoter (OR = 4.17, P = 2.35 × 10-9), with consistent results in survival analysis. The CDKN2B-AS1 variant replicated (P = 0.04) when adjusted for diabetic kidney disease in three additional T1D cohorts (cases/controls: 434/3123). Furthermore, we explored the function of the lead discoveries with a cardio-phenome-wide analysis. Among the eight suggestive loci (P < 1 × 10-6), rs70962766 near B3GNT2 associated with central blood pressure, rs1344228 near CNTNAP5 with intima media thickness, and rs2112481 on GRAMD2B promoter with serum leucocyte concentration. Finally, we calculated genetic risk scores for individuals with T1D with the known susceptibility loci. General population risk variants were modestly but significantly associated with CAD also in T1D (P = 4.21 × 10-7). CONCLUSION: While general population CAD risk loci had limited effect on the risk in T1D, for the first time, variants at the CDKN2B-AS1 locus were robustly associated with CAD in individuals with T1D. The novel finding on ß-defensin DEFB127 promoter provides a link between diabetes, infection susceptibility, and CAD, although pending on future confirmation.